data <- read.csv("data.csv")
result <- data.frame(estimate=0,se=0,p=0)
O3_lag0 <- data[,1]
for (i in 1:5) {
linear <- lm(O3_lag0~data[,i+4])
result[i,1]<-summary(linear)$coefficients[2,1]
result[i,2]<-summary(linear)$coefficients[2,3]
result[i,3]<-summary(linear)$coefficients[2,4]
}
result
我的问题,我想根据这个代码作进一步的改动,就是多个因变量和多个自变量分别做线性回归,目前这个我通过2个for循环实现了,可是问题是关于提取系数部分的结果不会写w(゚Д゚)w,不知道怎么才能依次把得到的结果提取出来,下面是我写的渣渣代码,就是标红部分有问题,不会写o(╥﹏╥)o:
data <- read.csv("data.csv")
result <- data.frame(estimate=0,se=0,p=0)
for (j in 1:4) {
for (i in 1:5) {
linear <- lm(data[,j]~data[,i+4])
result[i,1]<-summary(linear)$coefficients[2,1]
result[i,2]<-summary(linear)$coefficients[2,3]
result[i,3]<-summary(linear)$coefficients[2,4]
}
}
result
附上我的数据~
- data.csv