R/include.2l.predictors_v1.R

In miceadds: Some Additional Multiple Imputation Functions, Especially for 'mice'

## File Name: include.2l.predictors_v1.R
## File Version: 0.12


# function for inclusion of group level predictors
include.2l.predictors_v1 <- function( y, x, ry, type, vname, newstate, ... )
{
    X <- as.matrix( x[, type %in% c(1,2)] )
    X <- cbind( 1, X )
    # group level predictors
    if ( sum( type==-2 ) > 0 ){
        cluster <- x[, type==-2 ]
        if ( sum( type==2 ) > 0 ){
            x1a <- as.matrix( cbind( x[, type==- 2 ], x[, type==2 ]  ) )
            colnames(x1a) <- c( colnames(x)[ type==-2 ], colnames(x)[ type==2 ] )
            gm0 <- mice.impute.2l.groupmean(y=y, ry=FALSE * ry, x=x1a,
                            type=c( -2, rep(1, ncol(x1a)-1 ) ),
                            grmeanwarning=FALSE, ... )
            gm0 <- as.matrix(gm0)
            colnames(gm0) <- paste0("M.", colnames(x1a)[-1] )
            X <- as.matrix( cbind( X, gm0 ))
        }
    } else {
            cluster <- NULL
    }
    res <- mice_imputation_get_states()
    vname <- res$vname
    res$X <- X
    res$cluster <- cluster
    return(res)
}