Methods are proposed for the analysis of diseases with variable age at onset. The Cox proportional hazards model, widely used for epidemiologic analysis, is modified to include both measured (environmental) covariates and latent (genetic) variables, as well as their interactions. A Monte Carlo technique known as Gibbs sampling is utilized to generate observations from the posterior distributions of all model parameters. A correction to account for single ascertainment of pedigrees is also described. Simulation studies show that parameter estimation is nearly unbiased for a wide variety of models, and that moderate gene-environment interaction effects can be detected.