R/checkModelQuality.r

#' Function checkModelQuality
#' 
#' Returns information on model fit and psychomeric quality based on data.
#' @export


checkModelQuality <- function (data) 
{
    library(psych)
    library(TAM)
    a = suppressWarnings(alpha(data))
    print("Model Fit")
    ta = tam.mml(data, control = list(progress = F))
    mf = tam.modelfit(ta, progress = F)
    iff = tam.fit(ta, progress = F)
    l = c(length(data[1, ]), ta$EAP, tam.wle(ta, progress = F)$WLE.rel[1], 
        a$total$std.alpha, min(a$item.stats$r.cor), mf$statlist$pmaxX2, 
        mf$stat.MADaQ3$p, logLik(ta), -AIC(ta), -BIC(ta), (-mf$fitstat[[1]]), 
        (-mf$fitstat[[2]]), (-mf$fitstat[[3]]))
    names(l) = c("items", "EAPrel", "WLErel", "alpha", "minCPWC", 
        "pmaxX2", "MADaQ3$p", "logLik", "AIC", "BIC", "-MADCOV", 
        "-SRMR", "-SRMSR")
    return(cbind(`Goodness of fit` = l))
}
AndreasFischer1985/TAMirt documentation built on May 26, 2019, 7:01 p.m.