R/summary.mlnormal.R

In alexanderrobitzsch/LAM: Some Latent Variable Models

## File Name: summary.mlnormal.R
## File Version: 0.30

#*******************************************************
# Summary for mlnormal object
summary.mlnormal <- function( object, digits=4, file=NULL, ...)
{
    # open sink
    CDM::osink( file=file, suffix=paste0( "__SUMMARY.Rout") )

    cat("-----------------------------------------------------------------\n")
    #- package and R session
    sirt::sirt_summary_print_package_rsession(pack="LAM")

    cat( "Date of Analysis:", paste( object$s2 ), "\n" )
    cat("Computation Time:", print(object$s2 - object$s1), "\n\n")

    #- print call
    sirt::sirt_summary_print_call(CALL=object$CALL)

    cat( object$descriptions["des_method"], "\n\n")

    if ( object$REML & ! ( object$vcov_REML ) ){
        cat(" * Standard errors for random effects are not properly calculated! \n")
        cat(" * Use vcov=TRUE as an argument in 'mlnormal' if valid standard errors\n")
        cat("   for random effects parameters are requested.\n\n")
    }

    cat("-----------------------------------------------------------------\n")
    cat( "Number of observations=", object$ic$N, "\n" )
    cat( "Number of clusters=", object$ic$G, "\n\n" )

    cat( "Number of iterations=", object$iter, "\n\n" )
    cat( "Deviance=", round( object$deviance, 2 ), "\n" )
    cat( paste0( object$descriptions["log_like_verbose"]
                    , "=" ), round( object$ic$loglike, 2 ), "\n" )

    #---- print posterior
    if ( object$posterior_obj$display_posterior ){
        cat( "Log prior=",
                round( object$posterior_obj$log_prior, 2 ), "\n" )
        cat( "Log posterior=",
                round( object$posterior_obj$log_posterior, 2 ), "\n" )
    }

    cat("\n")
    cat( "Number of estimated parameters=", object$ic$np, "\n" )
    cat( "  Number of estimated beta parameters=", object$ic$np.beta,
                "\n" )
    cat( "  Number of estimated theta parameters=", object$ic$np.theta,
                "\n\n" )

    if ( ! object$REML ){
        cat( "AIC=", round( object$ic$AIC, 1 ), " | penalty=",
                round( object$ic$AIC - object$ic$deviance,2 ),
                "   | AIC=-2*LL + 2*p  \n" )
    }

#    cat( "AICc=", round( object$ic$AICc, 0 )," | penalty=", round( object$ic$AICc - object$ic$deviance,2 ) )
#        cat("    | AICc=-2*LL + 2*p + 2*p*(p+1)/(n-p-1)  (bias corrected AIC)\n" )
#    cat( "BIC=", round( object$ic$BIC, 0 ), " | penalty=", round( object$ic$BIC - object$ic$deviance,2 ),
#            "   | BIC=-2*LL + log(n)*p  \n" )
#    cat( "CAIC=", round( object$ic$CAIC, 0 )," | penalty=", round( object$ic$CAIC - object$ic$deviance,2 ) )
#        cat("   | CAIC=-2*LL + [log(n)+1]*p  (consistent AIC)\n\n" )

    cat("-----------------------------------------------------------------\n")
    cat("Beta Parameters\n")
    excl <- c("parm","prior")
    obji <- object$beta_summary
    mlnormal_summary_round_helper(obji, digits=digits, exclude=excl, print=TRUE)

    cat("-----------------------------------------------------------------\n")
    cat("Theta Parameters\n")
    obji <- object$theta_summary
    mlnormal_summary_round_helper(obji, digits=digits, exclude=excl, print=TRUE)

    # close sink
    CDM::csink( file=file )
}
#*******************************************************