Nothing
print.penmodel_c <- function (x, digits = max(3, getOption("digits") - 3), ...)
{
savedig <- options(digits = digits)
on.exit(options(savedig))
frailty.dist <- attr(x, "frailty.dist")
base.dist <- attr(x,"base.dist")
if(base.dist[1]==base.dist[2]) base.dist <- base.dist[1]
cat("Call: \n")
if(frailty.dist=="none" | is.null(frailty.dist)){
cat("Penetrance model for",attr(x,"design"),"design using",
base.dist, "baseline distributions. \n")
}
else{
f <- match(frailty.dist, c("gamma", "lognormal", "cgamma", "clognormal"))
fname <- c("Gamma", "Log-normal", "Correlated gamma", "Correlated log-normal")[f]
cat("Penetrance model for",attr(x,"design"),"design using \n", fname, "frailty distribution with",
base.dist, "baseline distributions. \n")
}
cat("Minimum age at onset used: ", attr(x, "agemin"))
cat("\n")
cat("\nEstimates: \n")
print(x$estimates)
if(attr(x, "robust")) cat("Robust standard errors was obtained. \n")
invisible(x)
}
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.