我有一个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),所以我知道计算应该全部解决。
编辑:我尝试了评论中的建议
listslope< -lapply(score1,function(x)coef(lm(x~obs))[2])
我正在
lm.fit中的错误(x,y,offset = offset,singular.ok = singular.ok,...): 0(非NA)案件
我将尝试弄清楚是什么导致了这个问题,但我想更新我的帖子并说声谢谢。
答案 0 :(得分:1)
lapply(score1, function(x) coef(lm(x~obs))[2])
应该在列表中给你。
如果您更喜欢矢量,请使用sapply
。