Modeling hazard functions in families

A genetic frailty model is presented for censored age of onset data in nuclear families where individuals carrying a genetic susceptibility gene have an increased risk of becoming affected. We use maximum likelihood via the EM algorithm to estimate the genetic relative risk and the allele frequency under a dominant susceptibility type and a proportional hazards model. When sampling is from a disease registry, likelihood corrections are necessary for reducing bias in the parameter estimates. In these biased samples, the full conditional likelihood is approximated by a likelihood conditional on the proband's age of onset. For unbiased samples, simulations show the distributions of the estimates are similar under both a semiparametric and the correctly specified parametric likelihoods. For biased samples, simulations under the approximate conditional likelihood show the median estimates of the allele frequency and genetic relative risk tend to under- and overestimate, respectively, the true values; however, the approximation is better for rarer allele frequencies (0.0033 vs. 0.01). In practice, large samples or more complex ascertainment corrections are recommended. Using the approximate conditional likelihood on familial breast cancer onset data collected as part of a case-control study at the Fred Hutchinson Cancer Research Center in Seattle, Washington, we estimate an allele frequency of 0.0009 (approximate 95% CI 0.0006-0.002) and a genetic relative risk of 104 (approximate 95% CI 55-181).