search for: vorvb

...l("http://www.anicca-vijja.de/lg/dframe.Rdata")) # load object "dframe" dim(dframe) apply(dframe,2,function(x) sum(is.na(x))) # check how many NAs in the dataset #dframe <- subset(dframe,select=-c(alter,grpzugeh,is1,is4,is6,klassenstufe,mmit,vorai,vorap,voras,vorkf,vorsg,vorvb)) s1 <- prelim.norm(dframe) s1$nmis # re-check of NAs should be identical to above s2 <- prelim.norm(dframe[,1:32])# see below -> still NAs are available - _not_ imputed thetahat1 <- em.norm(s1) theta1 <- da.norm(s1,thetahat1,steps=20,showits=TRUE) # error:...