R/summary.immer_cml.R

In immer: Item Response Models for Multiple Ratings

## File Name: summary.immer_cml.R
## File Version: 0.192


#*******************************************************
# Summary for immer object
summary.immer_cml <- function( object, digits=3, file=NULL, ... )
{
    #-- open sink
    immer_osink( file=file )

    #-- package description
    cat("-----------------------------------------------------------------\n")
    immer_summary_print_package_rsession(pack="immer")

    cat( object$description, "\n\n")

    #-- call
    immer_summary_print_call(CALL=object$CALL)

    #-- computation time
    immer_summary_print_computation_time( object=object )

    cat("-----------------------------------------------------------------\n")

    cat( "Deviance=", round( object$dev, 2 ), " | " )
    cat( "Log Likelihood=", round( object$loglike, 2 ), "\n" )

    cat( "Number of persons=", object$N, "\n" )
    cat( "Number of missing data patterns=", object$NP, "\n" )
    cat( "Number of items=", object$I, "\n" )
    cat( "Number of estimated parameters=", object$npars, "\n" )

    if (FALSE){
        cat( "AIC=", round( object$ic$AIC, 2 ), " | penalty=",
                round( object$ic$AIC - object$ic$dev,2 ),
                "   | AIC=-2*LL + 2*p  \n" )
        cat( "AIC3=", round( object$ic$AIC3, 2 ), " | penalty=",
                    round( object$ic$AIC3 - object$ic$dev,2 ),
                    "   | AIC3=-2*LL + 3*p  \n" )
        cat( "AICc=", round( object$ic$AICc, 2 )," | penalty=",
                round( object$ic$AICc - object$ic$dev,2 ) )
            cat("    | AICc=-2*LL + 2*p + 2*p*(p+1)/(n-p-1)  (bias corrected AIC)\n" )
        cat( "BIC=", round( object$ic$BIC, 2 ), " | penalty=",
                    round( object$ic$BIC - object$ic$dev,2 ),
                "   | BIC=-2*LL + log(n)*p  \n" )
        cat( "aBIC=", round( object$ic$aBIC, 2 ), " | penalty=",
                    round( object$ic$aBIC - object$ic$dev,2 ),
                "   | aBIC=-2*LL + log((n-2)/24)*p  (adjusted BIC) \n" )
        cat( "CAIC=", round( object$ic$CAIC, 2 )," | penalty=",
                        round( object$ic$CAIC - object$ic$dev,2 ) )
            cat("   | CAIC=-2*LL + [log(n)+1]*p  (consistent AIC)\n\n" )
    }

    cat("-----------------------------------------------------------------\n")
    cat("Item-Category Parameters \n")
    immer_summary_print_objects(obji=object$item, from=2, digits=digits,
                rownames_null=TRUE)

    cat("-----------------------------------------------------------------\n")
    cat("Estimated Basis Parameters \n")
    immer_summary_print_objects(obji=object$par_summary, from=2, digits=digits,
                rownames_null=TRUE)

    #-- close sink
    immer_csink( file=file )
}
#*******************************************************