Displaying 1 result from an estimated 1 matches for "prtyid".
Did you mean:
partyid
2004 Dec 16
1
help with multiple imputation using imp.mix
...vard.edu/~gov2001/Data/immigrati
on.dat',header=TRUE)
library(mix)
rngseed(12345678)
# Preare data for imputation
gender1<-c()
gender1<-as.integer(data$gender)
gender1[gender1==1]<-2
gender1[gender1==0]<-1
data$gender<-gender1
x<-cbind(data$gender,data$ipip,data$ideol,data$prtyid, data$wage1992)
colnames(x)<-c("gender","ipip", "ideol", "prtyid","wage")
# start imputation
s <- prelim.mix(x,4)
thetahat <- em.mix(s)
And here comes the error message:
> newtheta <- da.mix(s,thetahat, steps=100,showits=TRUE)
Steps...