Semiparametric estimation of major gene effects for age of onset

Analysis of age of onset is a key factor in the segregation and linkage analysis of some complex genetic traits. Previous work in the genetics literature has used parametric distributional assumptions on age of onset. In this paper, a Cox model with latent major gene effects is used: a semiparametric model with unspecified baseline hazard. A Monte Carlo EM procedure is used to obtain maximum likelihood estimates. Markov chain Monte Carlo is used to realize genotypic configurations from the posterior distribution given the current model and the observed data, and these genotypic configurations are used to estimate the expectations in the EM algorithm. Simulated data sets indicate that the parameters can be estimated well, and one real data set shows the practical applicability of the proposed method.