Amelia を使用して欠損値を代入しています。
Zelig と Amelia を使用して計算を行うことはできますが...
これらのパッケージを使用して、新しく帰属されたデータのプールされた平均と標準偏差を見つけるにはどうすればよいですか?
library(Amelia)
library(Zelig)
n= 100
x1= rnorm(n,0,1) #random normal distribution
x2= .4*x1+rnorm(n,0,sqrt(1-.4)^2) #x2 is correlated with x1, r=.4
x1= ifelse(rbinom(n,1,.2)==1,NA,x1) #randomly creating missing values
d= data.frame(cbind(x1,x2))
m=5 #set 5 imputed data frames
d.imp=amelia(d,m=m) #imputed data
summary(d.imp) #provides summary of imputation process