在R: Gdd计算中创建自定义函数

# 定义计算GDD的函数
calculate_GDD <- function(temperatures, base_temp = 5) {
  # 计算每日GDD
  daily_gdd <- pmax(temperatures - base_temp, 0)
  # 返回累积GDD
  return(sum(daily_gdd))
}

# 示例数据
dates <- seq(as.Date("2023-01-01"), as.Date("2023-01-31"), by="days")
temperatures <- runif(length(dates), min = -5, max = 15)  # 随机生成温度数据

# 计算GDD
gdd_result <- calculate_GDD(temperatures)
print(paste("累积GDD:", gdd_result))

temperatures[is.na(temperatures)] <- mean(temperatures, na.rm = TRUE)  # 填充缺失值

# 使用IQR方法识别异常值
Q1 <- quantile(temperatures, 0.25)
Q3 <- quantile(temperatures, 0.75)
IQR <- Q3 - Q1
lower_bound <- Q1 - 1.5 * IQR
upper_bound <- Q3 + 1.5 * IQR
temperatures[temperatures < lower_bound] <- lower_bound
temperatures[temperatures > upper_bound] <- upper_bound