在循环中迭代计算线性回归

Question

我在R工作，我有2个向量，LE和NEE（长度为5265），我想用它计算两者之间的线性回归，但是一次得到5个点，并将每个5点回归系数值放入另一个名为WUE_5的向量。 这是我计算所有点的回归系数的函数：

WUE_function <- function(NEE, LE) {
  return(lm(NEE ~ LE)$coefficients[2])
}

我无法弄清楚如何编写一个基本上由以下组成的循环

WUE_5 <- c(lm(NEE[1:5] ~ LE[1:5])$coefficients[2],
           lm(NEE[2:6] ~ LE[2:6])$coefficients[2],
           ..., etc)

我尝试了一些任意向量x和y，但这是结果：

x <- c(1:10)
y <- c(2, 5, 3, 6, 7, 8, 13, 6, 3, 8)
n <- 7
i <- 1
z <- NULL


while(i < n){
  z[i] <- lm(x[i:i+4] ~ y[i:i+4])$coefficients[2]
  i <- i+1 
} 

z
# [1] NA NA NA NA NA NA 

我不明白我做错了什么让z充满NA的感谢你提前感谢你的帮助！

Answer 1

序列运算符:具有高优先级，然后是加法运算符+ 。 这意味着i:i+4转换为i+4 （这是单个元素）。

如果你这样做了

z[i]<- lm(x[i:(i+4)]~y[i:(i+4)])$coefficients[2]

相反，我相信你会得到你所追求的结果