Plot estimated penetrance curves from a model fit

Share:

Description

Plot method for penmodel. Plots penetrance curves estimated from the model fit and overlays non-parametric penetrance curves estimated from the data without proabands.

Usage

1
2
## S3 method for class 'penmodel'
plot(x, ...)

Arguments

x

An object class of penmodel, a fitted model by penmodel function.

...

additional arguments

Value

It displays the following summary values:

parameters

Parameter estimates of baseline parameters (λ, ρ) and regression coefficients for gender and mutation status (β_s, β_g) including their standard errors and also robust standard errors.

pen70

Penetrance estimates by age 70 specific to gender and mutation-status subgroups including their standard errors and 95% confidence intervals.

Author(s)

Yun-Hee Choi

See Also

penmodel, summary.penmodel, simfam

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
# Simulated family data  

fam <- simfam(N.fam=300, design="pop+", base.dist="Weibull", variation="none",
       base.parms=c(0.01,3), vbeta=c(-1.13, 2.35), allelefreq=0.02, agemin=20)
 
# Fit family data

fit <- penmodel(parms=c(0.01, 3), vbeta=c(-1.13, 2.35), data=fam, 
       design="pop+", base.dist="Weibull")
       
# Plot penetrance function curves

plot(fit)