请参见之前回答... 使用R语言实现。
https://bbs.pinggu.org/forum.php?mod=viewthread&tid=6783316&page=1#pid55118498
#### change the following address by yourslf ####
data <- read.csv("D:/path/data.csv",)
result <- data.frame(estimate=0,se=0,p=0)
O3_lag0 <- data[,1]
for (j in 1:4) {
for (i in 1:5) {
linear <- lm(data[,i]~data[,i+4])
for (k in 1:3) {
result[(5*j)+(i)-5,k] <- summary(linear)$coefficients[2,k]
}
}
}
result