【发布时间】:2018-10-22 11:56:17
【问题描述】:
大家好,我有一些数据,想根据这段代码自动运行一些回归:
library(dplyr)
library(plyr)
set.seed(123)
v1<-c(rnorm(340,0,1))
v2<-c(rnorm(340,1,1))
v3<-rnorm(340,2,1)
v4<-rnorm(340,4,1)
v5<-rnorm(340,8,1)
v6<-rnorm(340,0,3)
v7<-rnorm(340,2,6)
v8<-rnorm(340,3,2)
v9<-rnorm(340,9,1)
v10<-rnorm(340,0,2)
date<-seq(as.Date("1990-01-01"), as.Date("2018-04-01"),by="months")
data<-data.frame(date,v1,v2,v3,v4,v5,v6,v7,v8,v9,v10)
rownames(data)<-data$date
data<- data[-1]
rm(v1,v2,v3,v4,v5,v6,v7,v8,v9,v10)
all.combs<- function(dataframe,y){
df<-dataframe
df[[substitute(y)]]<- NULL
m1<- data.frame(t(combn(colnames(df),1)))
colnames(m1)<-"x1"
m2<- data.frame(t(combn(colnames(df),2)))
colnames(m2)<-c("x1","x2")
m3<- data.frame(t(combn(colnames(df),3)))
colnames(m3)<-c("x1","x2","x3")
m4<- data.frame(t(combn(colnames(df),4)))
colnames(m4)<-c("x1","x2","x3","x4")
liste<- rbind.fill(m1,m2,m3,m4)
return(liste)
}
test<-all.combs(data,v3)
其中test 是一个数据框,其中每一行都是自变量列表,v3 是我的因变量。现在我基本上想用所有测试行运行回归。即类似的东西
for (i in 1:nrow(test){
reg<-lm(data$v3~data$test[i]
regs$"i"<-predict(reg)
}
所以我有一个数据框,其中包含我在 v3 上使用我的函数生成的所有组合的所有回归的所有预测值。
提前致谢, 彼得
【问题讨论】:
-
你到底在做什么,为什么?
-
@user2974951 最后,我想运行向量自回归来检查哪些变量对预测通货膨胀有用以及它们如何影响通货膨胀,基于以下论文:snb.ch/de/mmr/studies/id/economic_studies_2006_02。我已经有一个用 Gretl 编写的代码,我已经设置好了所有东西,除了我必须手写所有模型。
标签: r