R/polr.mi.R

In mi: Missing Data Imputation and Model Checking

polr.mi <- function (formula, mi.object, ... ) 
{
    call   <- match.call( )
    m      <- m(mi.object)
    result <- vector( "list", m )
    names( result ) <- as.character(paste( "Chain", seq( m ), sep = "" ))
    mi.data <- mi.completed(mi.object)
    
    for ( i in 1:m ) {
      result[[i]] <- update(polr( formula, data = mi.data[[i]], ... ), Hess=TRUE)
    }
    coef   <- vector( "list", m )
    se     <- vector( "list", m )
    for( j in 1:m ) {
      coef[[j]]<- lapply(result, summary)[[j]]$coefficients[ ,1]
      se[[j]]  <- lapply(result, summary)[[j]]$coefficients[ ,2]
    }
    
    pooled <- mi.pooled(coef, se)
    mi.pooled.object <- new("mi.pooled",
                          call = call, 
                          mi.pooled = pooled,
                          mi.fit = result)
    return( mi.pooled.object )
}