ループ内の回帰ループ

Question

このコードをループに書き直す最良の方法は何ですか？

for (i in 1:10) { 
a.data[i] <- read.csv('outdata[i].csv') 
growth.sub.QOG[i] <- merge(QOG, a.data[i], by = c('year', 'country'), all = F) 
growth.re[i] <- plm(NY.GDP.PCAP.KD.ZG ~ log(Enrolment.in.all.programmes..Tertiary..Total) + law + 
        engineering + log(SP.POP.TOTL) + lp.legor 
        ,data=growth.sub.QOG[i], model="random") 
summary(growth.re[i]) 
eststo(growth.re[i]) 
} 

が、それは動作しませんでした、それは私が間違ってやっているということです：

a.data1 <- read.csv('outdata1.csv') 
growth.sub.QOG1 <- merge(QOG, a.data1, by = c('year', 'country'), all = F) 
growth.re1 <- plm(NY.GDP.PCAP.KD.ZG ~ log(Enrolment.in.all.programmes..Tertiary..Total) + law + engineering + log(SP.POP.TOTL) + lp.legor 
,data=growth.sub.QOG1, model="random") 
summary(growth.re1) 
eststo(growth.re1) 


a.data2 <- read.csv('outdata2.csv') 
growth.sub.QOG2 <- merge(QOG, a.data2, by = c('year', 'country'), all = F) 
growth.re2 <- plm(NY.GDP.PCAP.KD.ZG ~ log(Enrolment.in.all.programmes..Tertiary..Total) + law + 
        engineering + log(SP.POP.TOTL) + lp.legor 
        ,data=growth.sub.QOG2, model="random") 
summary(growth.re2) 
eststo(growth.re2) 

a.data3 <- read.csv('outdata3.csv') 
growth.sub.QOG3 <- merge(QOG, a.data3, by = c('year', 'country'), all = F) 
growth.re3 <- plm(NY.GDP.PCAP.KD.ZG ~ log(Enrolment.in.all.programmes..Tertiary..Total) + law + 
        engineering + log(SP.POP.TOTL) + lp.legor 
        ,data=growth.sub.QOG3, model="random") 
summary(growth.re3) 
eststo(growth.re3) 

私はこのような何かをしようとしましたか？

Answer 1

があなたのファイル名を構築お知らせください。

残りのコードをどのように構造化するかは、それらの中間変数を気にするかどうかによって異なります。私はあなたがそうすると仮定していたので、たくさんの別々のループがあります。気にしない場合は、lapply文をマージしてください。

データを読み込みます。

all_data <- lapply(file, read.csv) 

マージ。

merged <- lapply(all_data, function(data) 
{ 
    merge(QOG, data, by = c('year', 'country'), all = FALSE) 
}) 

モデル。

models <- lapply(merged, function(data) 
{ 
    plm(
    NY.GDP.PCAP.KD.ZG ~ log(Enrolment.in.all.programmes..Tertiary..Total) + law + engineering + log(SP.POP.TOTL) + lp.legor, 
    data, 
    model = "random" 
) 
}) 

一部の出力を表示します。

(summaries <- lapply(models, summary)) 
(eststos <- lapply(models, eststo)) 

Answer 2

いくつかのサンプルデータは素晴らしかったですが、自発的に私はあなたがそのようなファイルで読むことができないというエラーを見ています。試してみてください。

file.name <- paste('outdata', i, '.csv', sep='') 
    variable <- paste('a.data', i, sep='') 
    data.in <- read.csv(file.name) 

を使用すると、動的に作成変数に格納したい場合は、これは次のように動作します。

assign(variable, data.in) 

これは最初の部分を修正する必要があります！

Answer 3

私は、これは

#instance of your directory 
datadir <-"D:/Regression" 
# set working directory, i.e. R knows where to get the data files 
setwd(datadir) 

csvfiles <- list.files(datadir,".csv$") 

#read data from datadir 
for(x in csvfiles) 
{ 
    assign(gsub(" ","",sub(".csv","",x)),read.csv(x,header=TRUE,stringsAsFactors=F,sep=";")) 
} 

data<-c("outdata1,outdata2,outdata3,...") 

i<-1 
for(x in data) 
{ 
    tmp <- eval(parse(text=x)) 
    growth.sub.QOG[i]<- merge(QOG,tmp, by = c('year', 'country'), all = F) 
    growth.re[i] <- plm(NY.GDP.PCAP.KD.ZG ~ log(Enrolment.in.all.programmes..Tertiary..Total) 
        + law + engineering + log(SP.POP.TOTL) + lp.legor, 
        data=tmp, model="random") 
    Summary[i]<-summary(growth.re[i]) 
    Est[i]<-eststo(growth.re[i]) 
    rm(tmp) 
    i<-i+1 
} 

幸運をうまくいくと思うし、私はあなたには、いくつかのエラーが発生した場合...