print.glh.rdc: Print a 'glh' object tested in RDC

In gmonette/spida15: Collection of miscellaneous functions for mixed models etc. prepared for SPIDA 2009+ (development version)

Usage

print.glh.rdc(x, round = 6, pround = round - 1, L = T, cov = T, ...)

Arguments

x
round
pround
L
cov
...

Examples

##---- Should be DIRECTLY executable !! ----
##-- ==>  Define data, use random,
##--	or do  help(data=index)  for the standard data sets.

## The function is currently defined as
function (x, round = 6, pround = round - 1, L = T, cov = T, ...) 
{
    rnd <- function(x, digits) {
        if (is.numeric(x)) 
            x <- round(x, digits = digits)
        format(x)
    }
    for (ii in 1:length(x)) {
        nn <- names(x)[[ii]]
        tt <- x[[ii]]
        ta <- tt$anova
        tap <- array("", dim = dim(ta), dimnames = dimnames(ta))
        cat("\n", nn, "\n", sep = "")
        print(formatCoefmat(ta, digits = round, pdigits = pround), 
            quote = F, right = T)
        cat("\n")
        te <- tt$estimate
        if (!is.null(zhead <- attr(tt, "heading"))) 
            cat(zhead, "\n")
        print(formatCoefmat(te, digits = round, pdigits = pround), 
            quote = F, right = T)
        if (L == T) {
            cat("\nL:\n")
            print(tt$L)
            if (dim(tt$L.full)[1] < dim(tt$L)[1]) {
                cat("\nL (full rank):\n")
                print(tt$L.full)
            }
        }
        if (cov == T) {
            cat("\nVar-Cov of estimates:\n")
            print(tt$vcov)
            cat("\nCorrelations:\n")
            print(tt$vcor)
        }
    }
    invisible(x)
  }

gmonette/spida15 documentation built on May 17, 2019, 7:26 a.m.