Nothing
## File Name: summary.rasch.copula2.R
## File Version: 0.205
#** Summary for rasch.copula object
summary.rasch.copula2 <- function( object, file=NULL, digits=3, ... )
{
is_copula2 <- inherits(object,"rasch.copula2")
is_copula3 <- inherits(object,"rasch.copula3")
# open sink for a file
sirt_osink(file=file)
display_string <- sirt_summary_print_display(symbol="-", len=65)
cat(display_string)
#- package and R session
sirt_summary_print_package_rsession(pack="sirt")
#- print call
sirt_summary_print_call(CALL=object$CALL)
#-- print computation time
sirt_summary_print_computation_time_s1(object=object)
cat("Marginal Maximum Likelihood Estimation \n")
cat(paste( "Raschtype Copula Model with generalized logistic link function: Estimation of alpha1 and alpha2 \n") )
cat( paste0("Function '", class(object), "'\n"))
cat("alpha1","=", round(object$alpha1, digits)," alpha2", "=",
round(object$alpha2, digits), "\n")
cat(display_string)
cat( "Deviance", "=", round( object$deviance, 2 ), "\n" )
cat( "Number of persons", "=", object$ic$n, " (", nrow(object$pattern), " Response Patterns)\n" )
cat( "Number of estimated parameters", "=", object$ic$np, "\n" )
cat( "Number of iterations", "=", object$iter, "\n" )
cat( "AIC", "=", round( object$ic$AIC, 2 ), " | penalty", "=",
round( object$ic$AIC - object$ic$deviance,2 ), "\n" )
cat( "AICc", "=", round( object$ic$AICc, 2 ), " | penalty", "=", round( object$ic$AICc - object$ic$deviance,2 ) )
cat(" (bias corrected AIC)\n" )
cat( "BIC", "=", round( object$ic$BIC, 2 ), " | penalty", "=", round( object$ic$BIC - object$ic$deviance,2 ), "\n" )
cat( "CAIC", "=", round( object$ic$CAIC, 2 )," | penalty", "=", round( object$ic$CAIC - object$ic$deviance,2 ) )
cat( " (consistent AIC) \n\n" )
#--- trait distribution
if (is_copula2){
cat( "Trait Distribution (", length(object$theta.k), " Knots )\n",
"Mean", "=", 0, " SD", "=", 1, "\n")
cat(paste("\nEAP Reliability:", round( object$EAP.Rel, digits)),"\n\n")
}
if (is_copula3){
cat( "Trait Distribution (", length(object$theta.k), " Knots )\n" )
cat("\nMean Vector\n")
print( round( object$mu, digits))
cat("\nCorrelation Matrix\n")
print( round( object$sigma, digits))
}
# item paramters
cat(display_string)
cat("Item Parameter \n")
obji <- object$item
sirt_summary_print_objects(obji=obji, digits=digits, from=2, rownames_null=TRUE)
# dependency parameters
cat("\nDependency parameters\n")
obji <- object$summary.delta
sirt_summary_print_objects(obji=obji, digits=digits, from=1, rownames_null=TRUE)
sirt_csink(file=file)
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.