写一段R代码用来计算数据文件中多年冬小麦生长季期间温度的平均值降雨的和并提取成熟时的产量冬小麦生长季定义为第一年播种到次年成熟时

以下是一段R代码，用于计算数据文件中多年冬小麦生长季期间温度的平均值，降雨的和，并提取成熟时的产量。冬小麦生长季定义为第一年播种到次年成熟时。

# 读取数据文件
data <- read.csv("data.csv")

# 计算冬小麦生长季的开始年份和结束年份
start_year <- min(data$Year)
end_year <- max(data$Year) + 1

# 创建空向量存储每年的温度和降雨数据
temperature <- numeric(length = end_year - start_year)
rainfall <- numeric(length = end_year - start_year)

# 计算每年的温度和降雨数据
for (i in start_year:end_year) {
  year_data <- subset(data, Year == i)
  temperature[i - start_year + 1] <- mean(year_data$Temperature)
  rainfall[i - start_year + 1] <- sum(year_data$Rainfall)
}

# 提取成熟时的产量
harvest_yield <- subset(data, Event == "Harvest")$Yield

# 输出结果
mean_temperature <- mean(temperature)
total_rainfall <- sum(rainfall)
mean_yield <- mean(harvest_yield)

cat("冬小麦生长季期间的平均温度:", mean_temperature, "\n")
cat("冬小麦生长季期间的总降雨量:", total_rainfall, "\n")
cat("冬小麦生长季期间成熟时的平均产量:", mean_yield, "\n")

请注意，上述代码假设数据文件名为"data.csv"，且该文件应包含以下列：Year（年份）、Temperature（温度）、Rainfall（降雨）和Yield（产量）。你需要将代码中的"data.csv"替换为实际数据文件的名称，并确保数据文件中包含所需的列