Combining innovations in hazard modeling with those in multilevel modeling, we develop a method to estimate discrete-time multilevel hazard models. We derive the likelihood of and formulate assumptions for a discrete-time multilevel hazard model with time-varying covariates at two levels. We pay special attention to assumptions justifying the estimation method. Next, we demonstrate file construction and estimation of the models using two common software packages, HLM and MLN. We also illustrate the use of both packages by estimating a model of the hazard of contraceptive use in rural Nepal using time-varying covariates at both individual and neighborhood levels.