如何通过循环将linearmodel2摘要结果添加到测量的时间序列中？

Question

我有一些数据，其中包含两个设备在不同月份的二氧化碳通量测量值。我已经使用此循环创建了时间序列图：

k

我已经从模型（MA）-[行] [列]中提取了所需的值

library(lmodel2)
layout(matrix(1:12, 4,3, byrow = TRUE))
par(mar=c(4,4,0,0)+0.5)

for (i in 1:12){
  plot(data_Hainich_2017$co2_flux_6262_R3[data_Hainich_2017$Month==i], 
       data_Hainich_2017$co2_flux_7200_HS[data_Hainich_2017$Month==i],
       xlab="CO2_flux_6262_R3", ylab = "CO2_flux_7200_HS")

  lm_co2_flux_6262_R3_7200_HS <- lmodel2(co2_flux_7200_HS[Month==i] ~ co2_flux_6262_R3[Month==i], 
                                         data=data_Hainich_2017, nperm=99)
}

现在，我想在时间序列中添加截距，斜率，r²和斜率，但是我不知道如何做。因此，应该有12个截距，12个坡度和12个R²（每个地块）。这是我的lmodel2的输出，它为我提供了3种不同类型的模型，这就是[row] [column]的原因。

intercept <- lm_co2_flux_6262_R3_7200_HS$regression.results[[2]][2]
slope <- lm_co2_flux_6262_R3_7200_HS$regression.results[[3]][2]
r.square <- lm_co2_flux_6262_R3_7200_HS$rsquare

非常感谢您的回答！