Our estimation strategy uses sequences of conditional probability functions, similar to those used in discrete time hazard rate analyses, to construct a discrete approximation to the density function of an outcome of interest conditional on exogenous explanatory variables. Once the conditional density function has been constructed, we can examine expectations of arbitrary functions of the outcome of interest and evaluate how these expectations vary with observed exogenous covariates. We demonstrate the features and precision of the conditional density estimation method (and compare it to other commonly used methods) through Monte Carlo experiments and an application to health expenditures using the RAND Health Insurance Experiment data. Overall, we find that the approximate conditional density estimator provides accurate and precise estimates of derivatives of expected outcomes for a wide range of types of explanatory variables.