Value-added models (VAMs) are used by many states to assess contributions of individual teachers and schools to students' academic growth. The generalized persistence VAM, one of the most flexible in the literature, estimates the "value added" by individual teachers to their students' current and future test scores by employing a mixed model with a longitudinal database of test scores. There is concern, however, that missing values that are common in the longitudinal student scores can bias value-added assessments, especially when the models serve as a basis for personnel decisions -such as promoting or dismissing teachers -as they are being used in some states. Certain types of missing data require that the VAM be modeled jointly with the missingness process in order to obtain unbiased parameter estimates.This dissertation studies two problems. First, the flexibility and multimembership random effects structure of the generalized persistence model lead to computational challenges that have limited the model's availability. To this point, no methods have been developed for scalable maximum likelihood estimation of the model. An EM algorithm to compute maximum likelihood estimates efficiently is developed, making use of the sparse structure of the random effects and error covariance matrices. The algorithm is implemented in the package GPvam in R statistical software. Illustrations of the gains in computational efficiency achieved by the estimation procedure are given.Furthermore, to address the presence of potentially nonignorable missing data, a flexible correlated random effects model is developed that extends the generalized persistence model to jointly model the test scores and the missingness process, allowing the process to depend on both students and teachers. The joint model gives the ability to test the sensitivity of the VAM to the presence of noni ignorable missing data. Estimation of the model is challenging due to the nonhierarchical dependence structure and the resulting intractable high-dimensional integrals. Maximum likelihood estimation of the model is performed using an EM algorithm with fully exponential Laplace approximations for the E step. The methods are illustrated with data from university calculus classes and with data from standardized test scores from an urban school district.ii To Laura, Oliver, and the rest of my family.iii ACKNOWLEDGEMENTS I owe a heartfelt thanks to my advisors, Dr. Lohr and Dr. Yang, who have dedicated so much of their time and effort to helping me develop this dissertation.They were selfless in their willingness to respond quickly to all of my emails and to meet with me whenever I needed their assistance. I am indebted to them for significant portions of my graduate education that they went out of their way to provide me with, including encouragement and funding to travel to JSM in Vancouver and Miami Beach, help with job searching, guiding me through the publication and peer-review process, and in general teaching me how to work in an academic environme...