## File Name: summary.btm.R
## File Version: 0.268
#--- summary.btm
summary.btm <- function( object, file=NULL, digits=4,... )
{
# open sink
sirt_osink( file=file )
res <- object
cat("------------------------------------------------------------\n")
#- package and R session
sirt_summary_print_package_rsession(pack="sirt")
#-- print computation time
sirt_summary_print_computation_time_s1(object=object)
cat("Computation Time Algorithm", print(object$time_alg), "\n\n")
#- print call
sirt_summary_print_call(CALL=object$CALL)
cat("Bradley-Terry Model with Ties and Home Advantage Parameters\n")
cat("------------------------------------------------------------\n")
cat( "Log-likelihood value", "=", round(object$ll,2), "\n" )
cat( "Number of iterations", "=", object$iter, "\n" )
# cat( "Deviance=", round( object$deviance, 2 ), " | " )
# cat( "Log Likelihood=", round( -object$deviance/2, 2 ), "\n" )
cat( "Number of individuals", "=", object$ic$n, "\n" )
cat( "Number of pairwise comparisons", "=", object$ic$D, "\n" )
cat( "Epsilon value", "=", object$eps, "\n" )
cat( "ignore.ties", "=", object$ignore.ties, "\n" )
cat( "wgt.ties", "=", round(object$wgt.ties,3), "\n" )
# cat( "Number of estimated parameters=", object$ic$np, "\n" )
cat("------------------------------------------------------------\n")
cat("Ties and Home advantage parameters\n")
obji <- res$pars
sirt_summary_print_objects(obji=obji, digits=digits, from=3)
cat("------------------------------------------------------------\n")
cat("Summary of individual effects parameters\n")
obji <- res$summary.effects
sirt_summary_print_objects(obji=obji, digits=digits, from=1, rownames_null=TRUE)
cat("------------------------------------------------------------\n")
cat("MLE reliability (separation reliability)\n")
cat(paste0("MLE Rel", "=", round( res$mle.rel, digits ), "\n") )
cat(paste0("Separation index", "=", round( res$sepG, digits ), "\n") )
cat("------------------------------------------------------------\n")
cat("Individual effects parameters\n")
obji <- res$effects
sirt_summary_print_objects(obji=obji, digits=digits, from=3)
if (res$multiple_judges){
cat("------------------------------------------------------------\n")
cat("Fit statistics judges\n")
obji <- res$fit_judges
sirt_summary_print_objects(obji=obji, digits=digits, from=2)
}
# close sink
sirt_csink( file=file )
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.