是否有一种简单的方法可以基于多个点计算线性回归斜率？

Question

data <- data.frame(Year1 = c(4, 2), Year2 = c(0, 10), Year3 = c(5, 6), Year4 = c(18, 4))

给定上面包含2年/列的2行/观测值。 我想让R基于每个观察的所有4个点计算线性回归线的斜率。

斜率值应出现在数据框右侧的新列中。

slope_data <- data.frame(Year1 = c(4, 2), Year2 = c(0, 10), Year3 = c(5, 6), Year4 = c(18, 4), Slope = c(5, 5))

我想知道是否有一个简单的函数或程序包可以实现这一目标？

Answer 1

根据评论，这将需要两行以上。 您可以使用lm获取一行的所需斜率，并使用sapply获取所有行的斜率。

Slope = function(x) {
    TempDF = data.frame(x, year=1:ncol(data))
    lm(x ~ year, data=TempDF)$coefficients[2]
}

TData = as.data.frame(t(data))
data$slope = sapply(TData, Slope)
data
  Year1 Year2 Year3 Year4 slope
1     4     0     5    18   4.7
2     2    10     6     4   0.2