R/print.aareg.R

In survival: Survival Analysis

print.aareg <- function(x, maxtime, test=c('aalen', 'nrisk'), scale=1, ...) {
    if (!inherits(x, 'aareg')) stop ("Must be an addreg object")
    if (!is.null(cl<- x$call)) {
	cat("Call:\n")
	dput(cl)
	cat("\n")
	}

    if (missing(test)) test <- x$test
    else test <- match.arg(test)

    if (missing(maxtime)) summ <- summary(x, test=test, scale=scale)
    else                  summ <- summary(x, maxtime=maxtime, test=test,
                                          scale=scale)

    omit <- x$na.action
    if (length(omit))
	cat("  n=", x$n[1], " (", naprint(omit), ")\n", sep="")
    else cat("  n=", x$n[1], "\n")
    cat("   ", summ$n[2], "out of", x$n[3], "unique event times used\n\n")
    print(signif(summ$table,3))
    chi <- summ$chisq
    df <- nrow(summ$table) -1
    pdig <- max(1, getOption("digits")-4)  # default it too high IMO
    cat("\nChisq=", format(round(chi,2)), " on ", df, " df, p=",
	            format.pval(pchisq(chi,df, lower.tail=FALSE), digits=pdig), 
	            "; test weights=", x$test, "\n", sep="")
    invisible(x)
    }