This function fits the penetrance model to the family data based on a prospective likelihood with ascertainment correction and provides parameter estimates as well as the gender and mutationspecific penetrance estimates.
1 
parms 
Vector of initial values for baseline parameters.

vbeta 
Vector of initial values for regression coefficients for gender and majorgene.

data 
Family data structure should follow the format of the data generated from 
design 
Study design of the family data. Possible choices are: 
base.dist 
Choice of baseline hazard distribution to fit. Possible choices are: 
The penetrance model is fitted to family data with a specified baseline hazard distribution,
h(tx_s, x_g) = h_0(t) \exp(β_s x_s+β_g x_g)
where h_0(t) is the baseline hazards function specified by base.dist
, which depends on the shape and scale parameters, λ and ρ; x_s indicates male (1) and female (0) and x_g indicates carrier (1) or noncarrier (0) of a gene of interest (major gene).
For family data arising from population or clinicbased study designs (design="pop", "pop+"
, "cli"
, or "cli+"
), the parameters of the penetrance model are estimated from the ascertainmentcorrected prospective likelihood approach (Choi, Kopciuk and Briollais, 2008).
For family data arising from a twostage study design (design="twostage"
), model parameters are estimated based on the composite likelihood approach (Choi and Briollais, 2011)
An object of class penmodel
, a list including elements
parms.est 
Parameter estimates of baseline parameters (λ, ρ) and regression coefficients for gender and mutation status (β_s, β_g) including their standard errors and also robust standard errors. 
parms.cov 
Covariance matrix of parameter estimates. 
parms.se 
Standard errors of parameter estimates. 
parms.rcov 
Robust (sandwich) covariance matrix of parameter estimates. 
parms.rse 
Robust standard errors of parameter estimates. 
pen70.est 
Penetrance estimates by age 70 specific to gender and mutationstatus subgroups. 
pen70.se 
Standard errors of penetrance estimates by age 70 specific to gender and mutationstatus subgroups. 
pen70.ci 
95% confidence interval for penetrance estimates by age 70 specific to gender and mutationstatus groups. 
ageonset 
Vector of ages of onset ranging from 
pen.maleCarr 
Vector of penetrance estimates for male carriers from 
pen.femaleCarr 
Vector of penetrance estimates for female carriers from 
pen.maleNonCarr 
Vector of penetrance estimates for male noncarriers from 
pen.femaleNonCarr 
Vector of penetrance estimates for female noncarriers from 
YunHee Choi
Choi, Y.H., Kopciuk, K. and Briollais, L. (2008) Estimating Disease Risk Associated Mutated Genes in FamilyBased Designs, Human Heredity 66, 238251
Choi, Y.H. and Briollais (2011) An EM Composite Likelihood Approach for Multistage Sampling of Family Data with Missing Genetic Covariates, Statistica Sinica 21, 231253
simfam, penplot, summary.penmodel, plot.penmodel
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18  # Family data simulated from populationbased design using a Weibull baseline hazard
fam < simfam(N.fam=300, design="pop+", variation="none", base.dist="Weibull",
base.parms=c(0.01,3), vbeta=c(1.13, 2.35), agemin=20, allelefreq=0.02)
# Penetrance model fit for simulated family data
fit < penmodel(parms=c(0.01, 3), vbeta=c(1.13, 2.35), data=fam,
design="pop+", base.dist="Weibull")
# Summary of the model parameter and penetrance estimates from model fit
summary(fit)
# Generate the lifetime penetrance curves from model fit for specific gender and
# mutation status groups along with their nonparametric penetrance curves
plot(fit)

Questions? Problems? Suggestions? Tweet to @rdrrHQ or email at ian@mutexlabs.com.
All documentation is copyright its authors; we didn't write any of that.