R/summary.amh.R

In LAM: Some Latent Variable Models

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


##*** 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[1L:3,]
    obji[,-1] <- round( obji[,-1], 2 )
    rownames(obji) <- NULL
    print(obji)

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

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