R/summaryMIwaldtest_linear_hyp.R

Defines functions summaryMIwaldtest_linear_hyp

## File Name: summaryMIwaldtest_linear_hyp.R
## File Version: 0.09



summaryMIwaldtest_linear_hyp <- function(object, digits)
{
    require_namespace(pkg="CDM")
    x <- object
    table <- array( x$qbar, dim=c(length(x$qbar), 10) )
    dimnames(table) <- list(labels(x$qbar),
            c("est", "se", "t", "df", "Pr(>|t|)", "lo 95", "hi 95",
                    "fmi", "fmi_Betw", "fmi_Within"))
    table[, 2] <- sqrt( diag(x$Tm) )
    table[, 3] <- table[, 1]/table[, 2]
    table[, 4] <- x$df
    table[, 5] <- if (all(x$df > 0))
                    2 * (1 - stats::pt(abs(table[, 3]), x$df)) else NA
    table[, 6] <- table[, 1] - stats::qt(0.975, x$df) * table[, 2]
    table[, 7] <- table[, 1] + stats::qt(0.975, x$df) * table[, 2]
    table[, "fmi"] <- x$lambda
    table[, "fmi_Betw"] <- x$lambda_Between
    table[, "fmi_Within"] <- x$lambda_Within
    table <- as.data.frame(table)
    if ( is.na(table$se)[1] ){
        table$df <- NA
    }
    table$fmi_Betw <- NULL
    table$fmi_Within <- NULL
    CDM::cdm_print_summary_data_frame(table, from=1, digits=digits)
    invisible(table)
}

Try the miceadds package in your browser

Any scripts or data that you put into this service are public.

miceadds documentation built on May 29, 2024, 11:05 a.m.