R/summary.amh.R

In alexanderrobitzsch/LAM: Some Latent Variable Models

## File Name: summary.amh.R
## File Version: 0.356


##*** summary amh
summary.amh <- function( object, digits=3, file=NULL,... )
{
    # open sink
    CDM::osink( file=file, suffix="__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( "Number of iterations", "=", object$n.iter, "\n" )
    cat( "Number of burnin iterations", "=", object$n.burnin, "\n" )
    cat( "Number of saved iterations", "=", object$n.saved, "\n\n" )

    cat("-----------------------------------------------------------------\n")
    cat("Marginal MAP Estimation\n")

    cat( "Deviance", "=", round( object$deviance, 2 ), "\n" )
    cat( "Log Likelihood", "=", round( -object$deviance/2, 2 ), "\n" )
    cat( "Log Prior", "=", round( object$ic$logprior, 2 ), "\n" )
    cat( "Log Posterior", "=", round( object$ic$logpost, 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", "=",
            lam_print_summary_ic_penalty(ic=object$ic, crit="AIC", digits=2),
            "   | AIC=-2*LL + 2*p  \n" )
    cat( "AICc", "=", round( object$ic$AICc, 0 )," | penalty", "=",
            lam_print_summary_ic_penalty(ic=object$ic, crit="AICc", digits=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", "=",
            lam_print_summary_ic_penalty(ic=object$ic, crit="BIC", digits=2),
            "   | BIC=-2*LL + log(n)*p  \n" )
    cat( "DIC", "=", round( object$ic$DIC, 0 ), " | penalty", "=",
                round( 2*object$ic$pD,2 ),
            "   | DIC=-2*LL + 2*pD  \n" )
    cat( "CAIC", "=", round( object$ic$CAIC, 0 )," | penalty", "=",
        lam_print_summary_ic_penalty(ic=object$ic, crit="CAIC", digits=2) )
        cat("   | CAIC=-2*LL + [log(n)+1]*p  (consistent AIC)\n\n" )

    cat("-----------------------------------------------------------------\n")
    cat("Prior Summary \n")

    print(object$prior_summary)

    cat("-----------------------------------------------------------------\n")
    cat("Parameter Summary (Marginal MAP estimation) \n")

    obji <- object$amh_summary
    vars <- c("parameter","MAP","SD", "Q2.5", "Q97.5", "Rhat","SERatio",
                    "effSize", "accrate")
    obji <- obji[-1,vars]
    NO <- ncol(obji) - 1
    for (vv in c( 2:(NO-1), NO+1) ){
        obji[,vv] <- round( obji[,vv], digits )
        }
    obji[,NO] <- round( obji[,NO] )
    rownames(obji) <- NULL
    print(obji)

    cat("-----------------------------------------------------------------\n")
    cat("Comparison of Different Estimators\n\n")

    cat("MAP: Univariate marginal MAP estimation\n")
    cat("mMAP: Multivariate MAP estimation (penalized likelihood estimate)\n")
    cat("Mean: Mean of posterior distributions\n\n")

    cat("Comparison Posterior:\n")
    obji <- object$comp_estimators[1:3,]
    obji[,-1] <- round( obji[,-1], 2 )
    rownames(obji) <- NULL
    print(obji)

    cat("\nParameter Summary:\n")
    obji <- object$comp_estimators[-c(1:3),]
    obji[,-1] <- round( obji[,-1], digits )
    rownames(obji) <- NULL
    print(obji)

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