s <- prelim.mix(hep3,3) # do preliminary manipulations
margins <- c(1,2,3) # saturated contingency table model
design <- diag(rep(1,12)) # identity matrix D=no of cells
thetahat <- ecm.mix(s,margins,design) # find ML estimate
rngseed(1234567) # random generator seed
newtheta <- dabipf.mix(s,margins,design,thetahat,steps=200)
ximp <- imp.mix(s,newtheta,hep3) # impute under newtheta
语句中得 margin 代表什么？
算到ecm那总是出现
Steps of ECM:
1...2...错误: 外接函数调用时不能有NA/NaN/Inf(arg5)