R/immer_ic_hrm.R

In immer: Item Response Models for Multiple Ratings

## File Name: immer_ic_hrm.R
## File Version: 0.102


#-- information criteria
immer_ic_hrm <- function( ic, summary.mcmcobj )
{
    ic$n <- ic$N
    pars <- paste(summary.mcmcobj$parameter)
    vars <- c('mu', 'sigma', 'a', 'b',  'phi', 'psi')
    VV <- length(vars)
    Npars <- NULL
    for (vv in 1L:VV){
        ind <- which( substring( pars, 1, nchar( vars[vv] ) )==vars[vv] )
        Npars[ vars[vv] ] <- length(ind)
    }
    ic$Npars <- Npars
    ic$np <- sum(Npars)
    ic <- immer_IC_calc(ic)
    return(ic)
}