0

我有一个 data.frame 和一个观察列表 http://imgur.com/7hL2bDC

以及观察时间列表 http://imgur.com/w1P0qsg

我想创建一个列表,其中包含该数据框每一列的线性回归线斜率和观察时间。

如果我是手动完成这一切,我会这样做:

raw_slope <- lm(score1$AKR ~ obs)

坡度 <- raw_slope$coefficients[[2]]

然后“斜率”是我想要的值。我试图了解如何运行某些东西并为我的 data.frame 中的 19 个公司符号中的每一个吐出斜率。

我是 R 的新手,但我正在努力学习它,我感谢任何帮助指出我正确方向的人。

EDIT2:我为每一行手动运行了 lm(score1$... ~ obs) ,所以我知道计算应该都可以解决。

编辑:我尝试了评论中的建议

listlope<-lapply(score1, function(x) coef(lm(x~obs))[2])

我越来越

lm.fit(x,y,offset = offset,singular.ok =singular.ok,...)中的错误:0(非 NA)案例

我将尝试找出导致问题的原因,但我想更新我的帖子并表示感谢。

4

1 回答 1

1
lapply(score1, function(x) coef(lm(x~obs))[2])

应该在列表中给你。

sapply如果您更喜欢矢量,请使用。

于 2013-11-06T01:04:13.830 回答