R/summary.pmle.R

In alexanderrobitzsch/LAM: Some Latent Variable Models

## File Name: summary.pmle.R
## File Version: 0.36

#*** summary pmle
summary.pmle <- function( object, digits=3, 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( object$description, "\n\n")

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

    cat( "Date of Analysis:", "\n" )
    cat( "   Start:", paste( object$time$start ), "\n" )
    cat( "   End  :", paste( object$time$end ), "\n" )
    cat("Computation time:", print(object$time$end - object$time$start), "\n\n")

    cat( "Optimization function", "=", object$optim_fct, "\n" )
    cat( "Gradient provided", "=", object$use_grad, "\n" )
    cat( "Convergence Code", "=", object$results_optim$convergence, "\n" )
    cat( "CONVERGED", "=", object$converged, "\n" )

    cat("-----------------------------------------------------------------\n")
    cat( "Deviance", "=", round( object$deviance, 2 ), "\n" )
    cat( "Log Likelihood", "=", round( object$ic$loglike, 2 ), "\n" )
    cat( "Log Prior", "=", round( object$ic$prior, 2 ), "\n" )
    cat( "Log Posterior", "=", round( object$ic$post, 2 ), "\n\n" )

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

    cat( "AIC", "=", round( object$ic$AIC, 0 ), " | 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("Prior Summary \n")
    obji <- object$prior_summary
    print(obji)

    cat("-----------------------------------------------------------------\n")
    cat("Parameter Summary \n")
    obji <- object$pmle_summary
    mlnormal_summary_round_helper(obji, digits=digits, print=TRUE, start_index=2)

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