从OLS中提取系数

时间:2019-09-16 07:12:53

标签: r for-loop

我正在使用数据集进行建模,并使用循环执行线性回归。您可以在下面看到代码:

#CODE

DATA<-data.frame(
  Time=c(1,2,3,4,5,6,7,8,9,10),
  NACE1=c(1,2,3,4,5,6,7,8,9,10),
  NACE2=c(2,4,6,8,10,12,14,18,20,22),
  NACE3=c(8,4,10,8,15,17,12,12,24,25)
)

MY_DATA<-data.frame(DATA)
MY_OLS=apply(MY_DATA,2,function(x){
  return(
    list(
      #Constant 
      summary(lm(log(x)~DATA$Time))
    )
  )
})

因此,下一步应该是使用此代码将系数提取到一个表中。

FINAL_TABLE<-lapply(MY_OLS, function(x) t(sapply(x, function(y) 
  c(testset = y[[1]]$coefficients))))

FINAL_TABLE<-data.frame(FINAL_TABLE)

但是此代码不起作用。有人可以帮助解决此问题并将系数放入一张表吗?

0 个答案:

没有答案