A new semiparametric estimator for estimating conditional expectation functions from incomplete data is proposed, which integrates parametric regression with nonparametric matching estimators. Besides its applicability to missing data situations due to non-response or attrition, the estimator can also be used for analyzing treatment effect heterogeneity and statistical treatment rules, where data on potential outcomes is missing by definition. By combining moments from a parametric specification with nonparametric estimates of mean outcomes in the non-responding population within a GMM framework, the estimator seeks to balance a good fit in the responding population with low bias in the non-responding population. The estimator is applied to analyzing treatment effect heterogeneity among Swedish rehabilitation programmes.