Partially linear models are commonly used in observational studies of the causal effect of treatment and/or exposure when there are observed confounding variables. The models are robust and asymptotically distribution-free for testing the causal null hypothesis. In this research, we investigate methods for estimating the partially linear models with data missing at random in all the variables, including the response, the treatment, and the confounding variables. We develop a general estimation method for inference in partially linear models with nonmonotone missing at random data. It proposes using partially linear working models to improve the estimation efficiency of the standard complete case method. It can be shown that the new estimator is consistent, which does not depend on the correctness of the working models. In addition, we recommend bootstrap estimates for the asymptotic variances and semiparametric models for the missing data probabilities. It is computationally simple and can be directly implemented in standard software. Simulation studies are provided to examine its performance. A real data example with sparsely observed missingness patterns is used to illustrate the method.