请大家帮我看一下这个程序，写的有点发晕了，希望可以得到大家的点拨！谢谢

首先R不会智能的把数据累加或者合并起来，其次某个数据集中的某个变量可以用$符号调用，如Rdata$d。当然还有其他问题。。。

这一句，

for (x1 in a1) {<br />
    Rdata1 <- read.csv(x1)<br />
    print(Rdata1)<br />
}

每一次循环都把Rdata1覆盖掉了，合理的做法是将几个data.frame存到一个list中。

Rdata <- list()<br />
for (i in seq_along(a1)) {<br />
    Rdata[[i]] <- read.csv(a1[i])<br />
    print(Rdata[[i]])<br />
}

在后面的循环中，你应该把数据拆成d = 0和d = 1的两部分，然后对d = 1的部分进行插补。

比如

Rdata.out <- list()<br />
for (i in seq_along(a1)) {<br />
    d1 <- subset(Rdata[[i]], d == 1)<br />
    d0 <- subset(Rdata[[i]], d == 0)<br />
    Rdata.out[[i]] <- amelia(x = d1, ...)<br />
    ...<br />
}

</p>

回复 第6楼 的 zjk0507：看起来好复杂，有木有试过rbind直接堆起来呢？

话说4楼回答的很清楚了，最后如果一起分析，7楼的方法也应该行得通。

回复 第10楼 的 zjk0507：一种方法是你在你可以在subset之后直接rbind成列表：

<br />
Rdata.out <- list()<br />
d1=d0=NULL<br />
for (i in seq_along(a1)) {<br />
    d1 <- rbind(d1,subset(Rdata[[i]], d == 1))<br />
    d0 <- rbind(d0,subset(Rdata[[i]], d == 0))<br />
    Rdata.out <- amelia(x = d1, ...)<br />
    ...<br />
}<br />

或者最后一起来：

<br />
Rdata.out <- list()<br />
Rdata.rest <- list()<br />
for (i in seq_along(a1)) {<br />
    d1 <- subset(Rdata[[i]], d == 1)<br />
    Rdata.rest[[i]] <- subset(Rdata[[i]], d == 0)<br />
    Rdata.out[[i]] <- amelia(x = d1, ...)<br />
    ...<br />
}<br />
Rdata.imp <- rbind(do.call(rbind,Rdata.out),do.call(rbind,Rdata.rest))<br />

话说你的这句

<br />
Rdata.imp[[i]] <- rbind(Rdata.imp[[i]],d0)<br />

应该是

<br />
Rdata.imp[[i]] <- rbind(Rdata.out[[i]],d0)<br />

且应该在for循环中。

如果还是有同样错误应该是数据问题了。
</p>

回复 第15楼 的 zjk0507：话说你可以先传数据的一小部分上来，你的要求我也大致了解了。还是实战演练的好~

回复 第18楼 的 zjk0507：有数据就好搞了：(PS：水平有限可能烦了点)

<br />
trans <- function(filepath="E:/SAS data"){<br />
  setwd(filepath)<br />
  a1 <- list.files(pattern="*.csv")<br />
  Rdata1=Rdata2=list()<br />
  for (i in a1) {<br />
    Rdata1[[i]] <- subset(read.csv(i),d==1)<br />
    Rdata2[[i]] <- subset(read.csv(i),d==0)<br />
    #print(Rdata1)<br />
  }</p>
<p>  for(i in seq_along(a1)){<br />
    Rdata.out <- amelia(x=Rdata1[[i]],m=5,p2s=0,startvals=0,logs="y",<br />
                        idvars = c("rho", "nu","Rnum","y0","d0","c","d","y1","rind"))<br />
    for(j in 1:5){<br />
      write.csv(rbind(Rdata2[[i]],Rdata.out[[1]][[j]]),file=paste(c("Trans-final",i,"-               ",j,".csv"),collapse=''),row.names=F)<br />
    }<br />
  }<br />
}<br />

应该可以了。

应该有更简单的方法，代码有问题可以问我。[s:11]

主要是你Amelia后数据变成了列表，我不太擅长哈~
</p>