metadata: authors: 0: author: Andrabi, Tahir
1: author: Das, Jishnu
2: author: Khwaja, Asim Ijaz
3: author: Zajonc, Tristan
count: Pakistan
docna: 0: docna: Do value-added estimates add value ?
            accounting for learning dynamics
docty: Policy Research Working Paper
theme: Education for all
majtheme: Human development
subtopic: Health Service Management and Delivery,Economics of Education,Educational Sciences,Educational Institutions & Facilities,Education Finance,Effective Schools and Teachers
teratopic: Education
abstracts: cdata!: Evaluations of educational programs
            commonly assume that what children learn persists over time.
            The authors compare learning in Pakistani public and private
            schools using dynamic panel methods that account for three
            key empirical challenges to widely used value-added models:
            imperfect persistence, unobserved student heterogeneity, and
            measurement error. Their estimates suggest that only a fifth
            to a half of learning persists between grades and that
            private schools increase average achievement by 0.25
            standard deviations each year. In contrast, estimates from
            commonly used value-added models significantly understate
            the impact of private schools on student achievement and/or
            overstate persistence. These results have implications for
            program evaluation and value-added accountability system design.
geo_reg_key: 80475,80451,80470
repnme: repnme: Do value-added estimates add value ?
            accounting for learning dynamics
display_title: Do value-added estimates add value
            accounting for learning dynamics
listing_relative_url: /research/2009/09/11136498/value-added-estimates-add-value-accounting-learning-dynamics
topicv3: Teachers Management m1326119 1523,Education m1327048 470,Private Education m1327271 1490,Psychiatry, Psychology m1504686 3613,Diagnostics m1504698 3625
publishtoextweb_dt: 2010-07-01T00:00:00Z
historic_topic: Education
no_of_pages: 44
model_extractions: WPS5066 Policy Research Working Paper 5066

# Do Value-Added Estimates Add Value? Accounting for Learning Dynamics

_Tahir Andrabi Jishnu Das Asim Ijaz Khwaja Tristan Zajonc_ The World Bank Development Research Group Human Development and Public Services Team September 2009
Policy Research Working Paper 5066

# **Abstract**

Evaluations of educational programs commonly assume that what children learn persists over time. The authors compare learning in Pakistani public and private schools using dynamic panel methods that account for three key empirical challenges to widely used value-added models: imperfect persistence, unobserved student heterogeneity, and measurement error. Their estimates suggest that only a fifth to a half of learning persists between grades and that private schools increase average achievement by 0.25 standard deviations each year. In contrast, estimates from commonly used value-added models significantly understate the impact of private schools’ on student achievement and/or overstate persistence. These results have implications for program evaluation and valueadded accountability system design.

This paper—a product of the Human Development and Public Services Team, Development Research Group—is part of a larger effort in the department to expand our knowledge of child learning and test scores as a broad measure of educational outcomes. Policy Research Working Papers are also posted on the Web at http://econ.worldbank.org. The author may be contacted at jdas1@worldbank.org.

_The Policy Research Working Paper Series disseminates the findings of work in progress to encourage the exchange of ideas about development issues. An objective of the series is to get the findings out quickly, even if the presentations are less than fully polished. The papers carry the names of the authors and should be cited accordingly. The findings, interpretations, and conclusions expressed in this paper are entirely those of the authors. They do not necessarily represent the views of the International Bank for Reconstruction and Development/World Bank and its affiliated organizations, or those of the Executive Directors of the World Bank or the governments they represent._ Produced by the Research Support Team
# Do Value-Added Estimates Add Value? Accounting for Learning Dynamics<sup>_∗_</sup>

Tahir Andrabi Jishnu Das Asim Ijaz Khwaja Tristan Zajonc<sup>_†_</sup>

## 1 Introduction

Models of learning often assume that children’s achievement persists between grades—what a child learns today largely stays with her tomorrow. Yet recent research highlights that treatment effects measured by test scores fade rapidly, both in randomized interventions and observational studies. Jacob, Lefgren and Sims (2008), Kane and Staiger (2008), and Rothstein (2008) find that teacher effects dissipate by between 50 and 80 percent over one year. The same pattern holds in several studies of supplemental education programs in developed and developing countries. Currie and Thomas (1995) document the rapid fade out of Head Start’s impact in the United States, and Glewwe, Ilias and Kremer (2003) and Banerjee et al. (2007) report on education experiments in Kenya and India where over 70 percent of the one-year treatment effect is lost after an additional year. Low persistence may in fact be the norm rather than the exception, and a central feature of learning.

Low persistence has critical implications for commonly used program evaluation strategies that rest heavily on assumptions about or estimation of persistence. Using primary data on public and private schools in Pakistan, this paper addresses the challenges to value-added evaluation strategies posed by 1) imperfect persistence of achievement, 2) heterogeneity in

> _∗_ An earlier version of this paper also circulated under the title “Here Today, Gone Tomorrow? Examining the Extent and Implications of Low Persistence in Child Learning”.

> _†tandrabi@pomona.edu,_ Pomona College. _jdas1@worldbank.org,_ World Bank, Washington DC and Center for Policy Research, New Delhi; _akhwaja@ksg.harvard.edu,_ Kennedy School of Government, Harvard University, BREAD, NBER; _tzajonc@fas.harvard.edu,_ Kennedy School of Government, Harvard University. We are grateful to Alberto Abadie, Chris Avery, David Deming, Pascaline Dupas, Brian Jacob, Dale Jorgenson, Elizabeth King, Karthik Muralidharan, David McKenzie, Rohini Pande, Lant Pritchett, Jesse Rothstein, Douglas Staiger, Tara Vishwanath, and seminar participants at Harvard, NEUDC and BREAD for helpful comments on drafts of this paper. This research was funded by grants from the Poverty and Social Impact Analysis and Knowledge for Change Program Trust Funds and the South Asia region of the World Bank. The findings, interpretations, and conclusions expressed here are those of the authors and do not necessarily represent the views of the World Bank, its Executive Directors, or the governments they represent.

1
learning, and 3) measurement error in test scores. We find that ignoring any of these learning dynamics biases estimates of persistence and can dramatically affect estimates of the valueadded of private schools.

To fix concepts, consider a simple model of learning evolution, _yit_ = _αTit_ + _βyi,t−_ 1 + _ηi_ + _υit,_ where _yit_ is child achievement measured by test scores in period _t_ , _Tit_ is the treatment or program effect in period _t_ , and _ηi_ is unobserved student ability that speeds learning each period. We refer to _β_ , the parameter that links test scores across periods, as persistence.<sup>1</sup> The canonical restricted value-added or gain-score model assumes that _β_ = 1 (for examples, see Hanushek (2003)). When _β <_ 1 _,_ test scores exhibit mean reversion. Estimates of the treatment or program effect, _α_ , that assume _β_ = 1 will be biased if the baseline achievement of the treatment and control groups differs and persistence is imperfect. This has led many researchers to advocate leaving lagged achievement on the right-hand side. However doing so is not entirely straightforward: if estimated by OLS, omitted heterogeneity that speeds learning, _ηi_ , will generally bias _β_ upward and any measurement in test scores _yi,t−_ 1 will bias _β_ downward. Both the estimate of persistence _β_ and the treatment effect _α_ may remain biased when estimated by standard methods.

To address these concerns, we introduce techniques from the dynamic panel literature (Arellano and Honore, 2001; Arellano, 2003) that require three years of data. There are several findings. First, we find that learning persistence is low: only a fifth to a half of achievement persists between grades. That is, _β_ is between 0.2 and 0.5 rather than closer to 1. These estimates are remarkably similar to those obtained in the United States (Jacob, Lefgren and Sims, 2008; Kane and Staiger, 2008; Rothstein, 2008). Second, OLS estimates of _β_ are contaminated both by measurement error in test scores _and_ unobserved student-level heterogeneity in learning. Ignoring both biases leads to higher persistence estimates between 0.5 and 0.6; correcting only for measurement error results in estimates between 0.7 and 0.8. For persistence, the upward bias from omitted heterogeneity outweighs measurement error attenuation.

Third, our estimate of the private schooling effect is highly sensitive to the persistence parameter. Since private schooling is a school input that that is continually applied and leads to a large baseline gap in achievement, this is expected. Indeed, we find that incorrectly assuming _β_ = 1 significantly understates and occasionally yields the wrong sign for private schools’ impact on achievement—providing a compelling example of Lord’s paradox (Lord, 1967). The restricted value-added model suggests that private schools contribute no more than public schools; in contrast, our dynamic panel estimates suggest large and significant

> 1There are several different uses of the term persistence in the education literature. We refer to persistence as the fraction of knowledge that persist from one period to the next. The education literature, however, also uses the term “persistence” to indicate the probability of continuation from grade to grade (as opposed to dropping out), or to indicate a child’s motivation or propensity to complete tasks in the face of adversity.

2
contributions ranging from 0.19 to 0.32 standard-deviations a year.<sup>2</sup> Notably, the lagged valueadded model estimated by OLS gives similar results for the private school effect as our more data intensive dynamic panel methods. This is due to the countervailing heterogeneity and measurement error biases on _β_ and because lagged achievement can act as a partial proxy for omitted heterogeneity in learning.<sup>3</sup> Finally, towards an economic interpretation of low persistence, we use question-level exam responses as well as household expenditure and time-use data to explore whether psychometric testing issues, behavioral responses, or forgetting contribute to low persistence—causes that have different welfare implications. This investigation suggests that measurement error, mechanical psychometric testing, and behavioral response based explanations are insufficient. Understanding the behavioral or technological reasons for low persistence remains a critical issue in the literature.

The value-added of our contribution is several fold. To begin with, we show that the restricted value-added estimates based on longitudinal data may be _worse_ than the naive crosssectional OLS estimates. Second, we demonstrate how unobserved heterogeneity in learning and measurement error in test scores can bias estimates of persistence. The low persistence we find implies that long-run extrapolations from short-run impacts are fraught with danger. In the model above, the long-run impact of continued treatment is _α/_ (1 _− β_ ); with estimates of _β_ around 0.2 to 0.5, these gains may be much smaller than those obtained by assuming that _β_ is close to 1.<sup>4</sup> Third, we find that students in Pakistan’s private schools learn significantly more each year than their public school counterparts but that the popular gain-score or restricted value-added model would have detected no effect of private schools on learning.<sup>5</sup> From a public 2Harris and Sass (2006) find the that the persistence parameter makes little difference to estimates of teacher effects, while we find it starkly affects the estimates of school type. This can be explained by the relative gaps in baseline achievement. It is likely that a child does not continue with the same teacher, or an equally good teacher, over time. Hence, even if we don’t observe children’s educational history, two children who currently have different teachers may have been exposed to a similar quality teachers in the past. As such, children with different teachers often do not differ substantially in their baseline learning levels. In contrast, given that there is little switching across school types, children currently in different schools differ substantially in baseline learning levels.

> 3This results suggests that correcting for measurement error alone may do more harm than good. For example, Ladd and Walsh (2002) correct for measurement error in the lagged value-added model of school effects by instrumenting using double-lagged test scores but don’t address potential omitted heterogeneity. They show this correction significantly changes school rankings and benefits poorly performing districts. Given that we find unobserved heterogeneity in learning rates, rankings that correct for measurement error may be poorer than those that do not.

> 4For example, Krueger and Whitmore (2001), Angrist et al. (2002), Krueger (2003), and Gordon, Kane and Staiger (2006) calculate the economic return of various educational interventions by citing research linking test scores to earnings of young adults (e.g. Murnane, Willett and Levy, 1995; Neal and Johnson, 1996). Although effects on learning as measured by test-scores may fade, non-cognitive skills that are rewarded in the labor market could persist. For instance, Currie and Thomas (1995), Deming (2008) and Schweinhart et al. (2005) provide evidence of long run effects of Head Start and the Perry Preschool Project, even though cognitive gains largely fade after children enroll in regular classes.

> 5An alternative identification strategy to value-added models that we do not pursue here is instrumental

3
finance point of view, these different estimates matter particularly since per pupil expenditures are lower in private relative to public schools.<sup>6</sup> Our results are consistent with growing evidence that relatively inexpensive, mainstream, private schools hold potential in the developing country context (Jimenez, Lockheed and Paqueo, 1991; Alderman, Orazem and Paterno, 2001; Angrist et al., 2002; Alderman, Kim and Orazem, 2003; Tooley and Dixon, 2003; Muralidharan and Kremer, forthcoming; Andrabi, Das and Khwaja, 2008). Overall, our general results support a movement towards long-run, experimental evaluations of educational interventions.

A final contribution of our work is that it applies a wider set of econometric tools from the dynamic panel data literature than have been typically used in the education literature. In the use of dynamic panel methods, our estimators bear greatest resemblance to those discussed by Schwartz and Zabel (2005) and Sass (2006). Both use simple dynamic panel estimators, in the first case using school-level data and in the second using the Arellano and Bond (1991) differences GMM approach. Santibanez (2006) also uses the Arellano and Bond (1991) estimator to analyze the impact of teacher quality. Our efforts extend to include system GMM approaches (Arellano and Bover, 1995) and to address measurement error in test scores and alternative assumptions regarding omitted heterogeneity. We find both are important.

The rest of the paper is organized as follows: Section 2 presents the basic education production function analogy and discusses the specification and estimation of the value-added approximations to it. Section 3 summarizes our data. Section 4 reports our main results, several robustness checks, and a preliminary exploration of the economic interpretation of persistence. Section 5 concludes by discussing implications for experimental and non-experimental program evaluation.

# 2 Empirical Learning Framework

The “education production function” approach to learning relates current achievement to all previous inputs. Boardman and Murnane (1979) and Todd and Wolpin (2003) provide two accounts of this approach and the assumptions it requires; the following is a brief summary.<sup>7</sup> Using notation consistent with the dynamic panel literature, we aggregate all inputs into a single vector x _it_ and exclude interactions between past and present inputs. Achievement for variables. We are exploring such strategies using plausible exogenous geographical variation in a companion paper focused on the returns to private school education in Pakistan and competition between public and private schools. Our emphasis here is on the challenges faced by popular value-added strategies.

> 6For details on the costs of private schooling in Pakistan see Andrabi, Das and Khwaja (2008).

> 7Researchers generally assume that the model is additively separable across time and that input interactions can be captured by separable linear interactions. Cunha, Heckman and Schennach (2006) and Cunha and Heckman (2007) are two exceptions to this pattern, where dynamic complementarity between early and late investments and between cognitive and non-cognitive skills are permitted.

4
child _i_ at time (grade) _t_ is therefore where _yit_<sup>_∗_istrueachievement,measuredwithouterror,andthesummed</sup><sup>_μis_arecumulative</sup> productivity shocks.<sup>8</sup> Estimating (1) is generally impossible because researchers do not observe the full set of inputs, past and present. The value-added strategy makes estimation feasible by rewriting (1) to avoid the need for past inputs. Adding and subtracting _βyit_<sup>_∗_,normalizing</sup><sup>_θ_1to</sup> unity, and assuming that coefficients decline geometrically ( _αj_ = _βαj−_ 1 and _θj_ = _βθj−_ 1 for all _j_ ) yields the _lagged value-added model_ The basic idea behind this specification is that lagged achievement will capture the contribution of all previous inputs and any past unobservable endowments or shocks. As before, we refer to _α_ as the _input coefficient_ and _β_ as the _persistence coefficient_ . Finally, imposing the restriction that _β_ = 1 yields the gain-score or _restricted value-added model_ that is often used in the education literature: This model asserts that past achievement contains no information about future gains, or equivalently, that an input’s effect on any subsequent level of achievement does not depend on how long ago it was applied. As we will see from our results, the assumption that _β_ = 1 is clearly violated in the data, and increasingly it appears, in the literature, as well. As a result, we will focus primarily on estimating (2).

There are two potential problems with estimating (2). First, the error term _μit_ could include individual (child-level) heterogeneity in _learning_ (e.g., _μit ≡ ηi_ + _υit_ ). Lagged achievement only captures individual heterogeneity if it enters through a one-time process or endowment, but talented children may also _learn_ faster. Since this unobserved heterogeneity enters in each period, Cov( _yi,t_<sup>_∗_</sup> _−_ 1<sup>_, μit_)</sup><sup>_>_0and</sup><sup>_β_willbebiasedupwards.</sup> The second likely problem is that test scores are inherently a noisy measure of latent achievement. Letting _yit_ = _yit_<sup>_∗_+</sup><sup>_εit_denoteobservedachievement,wecanrewritethelatentlagged</sup> value-added model (2) in terms of observables. The full error term now includes measurement

> 8This starting point is more restrictive than the more general starting framework presented by Todd and Wolpin (2003). In particular, it assumes an input applied in first grade has the same effect on first grade scores as an input applied in second grade has on second grade scores.

5
error, _μit_ + _εit − βεi,t−_ 1 .

Dropping all the inputs to focus solely on the persistence coefficient, the expected bias due to both of these sources is The coefficient is biased upward by learning heterogeneity and downward by measurement error. These effects only cancel exactly when Cov( _ηi, yi,t_<sup>_∗_</sup> _−_ 1<sup>) =</sup><sup>_σ_</sup> _ε_<sup>2</sup><sup>_β_(Arellano,2003).</sup> Furthermore, bias in the persistence coefficient leads to bias in the input coefficients, _α_ . To see this, consider imposing a biased _β_<sup>ˆ</sup> and estimating the resulting model The error term now includes ( _β − β_<sup>ˆ</sup> ) _yi,t−_ 1. Since inputs and lagged achievement are generally positively correlated, the input coefficient will, in general, by biased downward if _β_<sup>ˆ</sup> _> β_ . The precise bias, however, depends on the degree of serial correlation of inputs and on the potential correlation between inputs and learning heterogeneity that remains in _μit_ .

This is more clearly illustrated in the case of the restricted value-added model (assuming that _β_ = 1) where: Therefore, if indeed there is perfect persistence as assumed and if inputs are uncorrelated with _ηi_ , OLS yields consistent estimates of the parameters _α_ . However, if _β <_ 1, OLS estimation of _α_ now results in two competing biases. By assuming an incorrect persistence coefficient we leave a portion of past achievement in the error term. This misspecification biases the input coefficient downward by the first term in (4). The second term captures possible correlation between current inputs and omitted learning heterogeneity. If there is none, then the second term is zero, and the bias will be unambiguously negative.

6
# 2.1 Addressing Child-Level Heterogeneity: Dynamic Panel Approaches to the Education Production Function

Dynamic panel approaches can address omitted child-level heterogeneity in value-added approximations of the education production function. We interpret the value-added model (2) as an autoregressive dynamic panel model with unobserved student-level effects: Identification of _β_ and _α_ is achieved by imposing appropriate moment conditions. Following Arellano and Bond (1991) and Arellano and Bover (1995), we focus on linear moment conditions and split our analysis into three groups: “differences” GMM, “differences and levels” GMM, and “levels only” GMM, which respectively refer to whether the estimates are based on the undifferenced “levels” equation (5), a differenced equation (see equation (7) below), or both. The section below provides a brief overview of the estimators we explore. Table 1 summarizes the estimators, including the standard static value-added estimators (M1-M4) and dynamic panel estimators (M5-M10). For more complete descriptions, Arellano and Honore (2001) and Arellano (2003) provide excellent reviews of these and other panel models.

## 2.1.1 GMM: Switching estimators

As noted previously, the value-added model differences out omitted endowments that might be correlated with the inputs. It does not, however, difference out heterogeneity that speeds learning. To accomplish this, the basic intuition behind the Arellano and Bond (1991) difference GMM estimator is to difference again. Differencing the dynamic panel specification of the lagged value-added model (5) yields Here, the differenced model eliminates the unobserved fixed effect _ηi_ . However, (7) cannot be estimated by OLS because _yi,t_<sup>_∗_</sup> _−_ 1<sup>iscorrelatedbyconstructionwith</sup><sup>_υi,t−_1intheerrorterm.</sup> Arellano and Bond (1991) propose instrumenting for _yi,t_<sup>_∗_</sup> _−_ 1<sup>_−y_</sup> _i,t_<sup>_∗_</sup> _−_ 2<sup>usinglagstwoperiodsand</sup> beyond, such as _yi,t_<sup>_∗_</sup> _−_ 2<sup>,orcertaininputs,dependingontheexogeneityconditions.Theselags</sup> are uncorrelated with the error term but are correlated with the change in lagged achievement, provided _β <_ 1. The input coefficient, in our case the added contribution of private schools, is primarily identified from the set of children who switch schools in the observation period.

The implementation of the GMM approach depends on the precise assumptions 7
about inputs. We consider two candidate assumptions: strictly exogenous inputs (M5) and predetermined inputs (M6). Strict exogeneity assumes past disturbances do not affect current and future inputs, ruling out feedback effects. In the educational context, this is a strong assumption. A child who experiences a positive or negative shock may adjust inputs in response. In our case, a shock may cause a child to switch schools.

To account for this possibility, we also consider the weaker case where inputs are predetermined but not strictly exogenous. Specifically, the predetermined inputs case assumes that inputs are uncorrelated with present and future disturbances but are potentially correlated with past disturbances. This case also assumes lagged achievement is uncorrelated with present and future disturbances. Compared to strict exogeneity, this approach uses only lagged inputs as instruments. Switching schools is instrumented by the original school type, allowing switches to depend on previous shocks. This estimator remains consistent if a child switches school at the same time as an achievement shock but still rules out parents anticipating and adjusting to future expected shocks.

# 2.1.2 Levels and GMM: Uncorrelated or constantly correlated

One difficulty with the differences GMM approach (M5 and M6) is that time-invariant inputs drop out of the estimated equation and their effects are therefore not identified. In our case, this means that the identification of the private school effect is based on the five percent of children who switch between public and private schools. We address the limited time-series variation using the levels and differences GMM framework proposed by Arellano and Bover (1995) and extended by Blundell and Bond (1998). Levels and differences GMM estimates a system of equations, one for the undifferenced levels equation (5) and another for the differenced equation (7). Further assumptions regarding the correlation between inputs and heterogeneity (though not necessarily between heterogeneity and lagged achievement) yield additional instruments.

We first consider predetermined inputs that have a constant correlation with the individual effects (M7). While inputs may be correlated with the omitted effects, constant correlation implies switching is not. The constant correlation assumption implies that _△_ x _it_ are available as instruments in the levels equation (Arellano and Bover, 1995). In the context of estimating school type, this estimator can be viewed as a levels and differences switching estimator since it relies on children switching school types in both the levels and differences equations. In practice, we often must assume that any time-invariant inputs are uncorrelated with the fixed effect or the levels equation, which includes the time-invariant inputs, is not fully identified.

A second possibility is that inputs are predetermined but are also uncorrelated with the omitted effects (M8). This allows using inputs x<sup>_t_</sup> _i_<sup>asinstrumentsinthelevelsmodel(5).The</sup> required assumption is fairly strong; it is natural to believe that inputs are correlated with 8
the omitted effect. Certainly, the decision to attend private school may be correlated with the child’s ability to learn. At the same time, the assumption is weaker than OLS estimation of lagged value-added model since the model (M8) allows for the omitted fixed effect to be correlated with lagged achievement.

# 2.1.3 Levels GMM: Conditional mean stationarity

In some instances, it may be reasonable to assume that, while learning heterogeneity exists, it does not affect achievement gains. A talented child may be so far ahead that imperfect persistence cancels the benefit of faster learning. That is, individual heterogeneity may be uncorrelated with gains, _yit_<sup>_∗−y_</sup> _it_<sup>_∗_</sup> _−_ 1<sup>, but not necessarily with</sup><sup>_learning, y_</sup> _it_<sup>_∗−βy_</sup> _it_<sup>_∗_</sup> _−_ 1<sup>.This situation</sup> arises when the initial conditions have reached a convergent level with respect to the fixed effect such that where _t_ = 1 is the first observed period and not the first period in the learning life-cycle. Blundell and Bond (1998) discuss this type of conditional mean stationarity restriction in considerable depth. As they point out, the key assumption is that initial deviations, _di_ , are uncorrelated with the level of _ηi/_ (1 _− β_ ). It does not imply that the achievement path, _{yi_<sup>_∗_</sup> 1<sup>_, y_</sup> _i_<sup>_∗_</sup> 2<sup>_, . . . , y_</sup> _iT_<sup>_∗}_,is</sup> stationary; inputs, including time dummies, continue to spur achievement and can be nonstationary. The assumption only requires that, conditional on the full set of controls and common time dummies, the individual effect does not influence achievement gains.

While this assumption seems too strong in the context of education, we discuss it because the dynamic panel literature has documented large downward biases of other estimators when the instruments are weak (e.g. Blundell and Bond, 1998). This occurs when persistence is perfect ( _β_ = 1) since the lagged value-added model then exhibits a unit root and lagged tests scores become weak instruments in the differenced model. The conditional mean stationarity assumption provides an additional _T −_ 2 non-redundant moment conditions that can augment the system GMM estimators. While a fully efficient approach uses these additional moments along with typical moments in the differenced equation, the conditional mean stationarity assumption ensures strong instruments in the levels equation to identify _β_ . Thus, if we prefer simplicity over efficiency, we can estimate the model using levels GMM or 2SLS and avoid the need to use a system estimator. In this simpler approach, we instrument the undifferenced value-added model (5) using lagged changes in achievement, _△yi_<sup>_∗t−_1</sup> , and either changes in inputs, _△_ x<sup>_t_</sup> _i_<sup>, or inputs directly, x</sup><sup>_t_</sup> _i_<sup>, depending on whether inputs are constantly correlated (M9)</sup> or are uncorrelated with the individual effect (M10).

9
## 2.2 Addressing Measurement Error in Test Scores

Measurement error in test scores is a central feature of educational program evaluation. Ladd and Walsh (2002), Kane and Staiger (2002), and Chay, McEwan and Urquiola (2005) all document how test-score measurement error can pose difficulties for program evaluation and valueadded accountability systems. In the context of value-added estimation, measurement error attenuates the coefficient on lagged achievement and can bias the input coefficient in the process. Dynamic panel estimators do not address measurement error on their own. For instance, if we replace true achievement with observed achievement in the standard Arellano and Bond (1991) setup, (7) becomes The standard potential instrument, _yi,t−_ 2, is uncorrelated with _△υit_ but is correlated with _△εi,t−_ 1 = _εi,t−_ 1 _− εi,t−_ 2 by construction.

The easiest solution is to use either three-period lagged test scores or alternate subjects as instruments. In the dynamic panel models discussed above, correcting for measurement error using additional lags requires four years of data for each child—a difficult requirement in most longitudinal datasets, including ours. We therefore use alternate subjects, although doing so does not address the possibility of correlated measurement error across subjects.

An alternative to instrumental variables strategies is to correct for measurement error analytically using the standard error of each test score, available from Item Response Theory.<sup>9</sup> Because the standard error is heteroscedastic—tests discriminate poorly between children at the tails of the ability distribution—one can gain efficiency by using the heteroscedastic errorsin-variables (HEIV) procedure outlined in Sullivan (2001) and followed by Jacob and Lefgren (2005), among others. Appendix A provides a detailed explanation of this analytical correction. While this correction is easy to apply in an OLS model, it becomes considerably more complicated in the dynamic panel context, and we therefore use an instrumental variable strategy for most of our estimators.

# 3 Data

To demonstrate these issues, we use data collected by the authors as part of the Learning and Educational Achievement in Punjab Schools (LEAPS) project, an ongoing survey of learning in Pakistan. The sample comprises 112 villages in 3 districts of Punjab: Attock, Faisalabad, and Rahim Yar Khan. Because the project was envisioned in part to study to dramatic rise of

> 9Item Response Theory provides the standard error for each score from the inverse Fisher information matrix after ML estimation of the IRT model. This standard error is reported in many educational datasets.

10
private schools in Pakistan, the 112 villages in these districts were chosen randomly from the list of all villages with an existing private school. As would be expected given the presence of a private school, the sample villages are generally larger, wealthier, and more educated than the average rural village. Nevertheless, at the time of the survey, more than 50 percent of the province’s population resided in such villages (Andrabi, Das and Khwaja, 2006).

The survey covers all schools within the sample village boundaries _and_ within a short walk of any village household. Including schools that opened and closed over the three rounds, 858 schools were surveyed, while three refused to cooperate. Sample schools account for over 90 percent of enrollment in the sample villages.

The first panel of children consists of 13,735 third-graders, 12,110 of which were tested in Urdu, English, and mathematics. These children were subsequently followed for two years and retested in each period. Every effort was made to track children across rounds, even when they were not promoted. In total, 12 percent and 13 percent of children dropped out or were lost between rounds one and two, and two and three, respectively.<sup>10</sup> In addition to being tested, 6,379 children—up to ten in each school—were randomly administered a survey including anthropometrics (height and weight) and detailed family characteristics such parental education and wealth, as measured by principal components analysis analysis of 20 assets. When exploring the economic interpretation of persistence, we also use a small subsample of approximately 650 children that can be matched to a detailed household survey that includes, among other things, child and parental time use and educational spending.

For our analysis, we use two subsamples of the data: all children who were tested in all three years (N=8120) and children who were tested _and_ given a detailed child survey in all three years (N=4031). Table 2 presents the characteristics of these children split by whether they attend public or private schools. The patterns across each subsample is relatively stable. Children attending privates schools are slightly younger, have fewer elder siblings, and come from wealthier and more educated households. Years of schooling, which largely captures grade retention, is lower in private schools. Children in private schools are also less likely to have a father living at home, perhaps due to a migration or remittance effect on private school attendance.

The measures of achievement are based on exams in English, Urdu (the vernacular), and mathematics. The tests were relatively long (over 40 questions per subject) and were designed to maximize the precision over a range of abilities in each grade. While a fraction of questions

> 10Attrition in private schools is 2 percentage points higher than in public schools. Children who drop out between rounds one and two have scores roughly 0.2 s.d. lower than children that don’t. Controlling for school type and drop out status, drop outs in private schools are slightly better (0.05 sd) than children in public schools, although the difference is only statistically significant for math. Given the small relative differences in attrition between public and private schools, additional corrections for attrition are unlikely to significantly our results.

11
changed over the years, the content covered remained consistent, and a significant portion of questions appeared across all years. To avoid the possibility of cheating, the tests were administered directly by our project staff and not by classroom teachers. The tests were scored and equated across years by the authors using Item Response Theory so that the scale has cardinal meaning. Preserving cardinality is important for longitudinal analysis since many other transformations, such as the percent correct score or percentile rank, are bounded artificially by the transformations that describe them. By comparison, IRT scores attempt to ensure that change in one part of the distribution is equal to a change in another, in terms of the latent trait captured by the test. Children were tested in third, fourth, and fifth grades during the winter at roughly one year intervals. Because the school year ends in the early spring, the test scores gains from third to fourth grade are largely attributable to the fourth grade school.

# 4 Results

## 4.1 Cross-sectional and Graphical Results

Before presenting our estimates of learning persistence and the implied private school effect, we provide some rough evidence for a significant private school effect using cross-sectional and graphical evidence. These results do not take advantage of the more sophisticated specifications above but nevertheless provide initial evidence that the value-added of private schools is large and significant.

### 4.1.1 Baseline estimates from cross-section data

Table 3 presents results for a cross-section regression of third grade achievement on child, household, and school characteristics. These regressions provide some initial evidence that the public-private gap is more than omitted variables and selection. Adding a comprehensive set of child and family controls reduces the estimated coefficient on private schools only slightly. Adding village fixed effects also does not change the coefficient, even though the _R_<sup>2</sup> increases substantially. Across all baseline specifications, the gap remains large: over 0.9 standard deviations in English, 0.5 standard deviations in Urdu, and 0.4 standard deviations in mathematics.

Besides the coefficient on school type, few controls are strongly associated with achievement. By far the largest other effect is for females, who outperform their male peers in English and Urdu. However, even for Urdu, where the female effect is largest, the private school effect is still nearly three times as large. Height, assets, and whether the father (and for Column 3, mother) is educated past elementary school also enter the regression as positive and significant. More elder brothers and more years of schooling (i.e. being previously retained) correlates with lower achievement. Children with a mother living at home perform worse although this result is driven 12
by an abnormal subpopulation of two percent of children with absent mothers. Overall, these results confirm mild positive selection into private schools but also suggest that, controlling a host of other observables typically not available in other datasets (such as child height and household assets) does not alter significantly the size of the private schooling coefficient.

## 4.1.2 Graphical evidence

Figure 1 plots learning levels in the tested subjects (English, mathematics, and the vernacular, Urdu) over three years. While, levels are always higher for children in private schools, there is little difference in learning gains (the gradient) between public and private schools. This illustrates why a specification that uses learning gains (i.e., assumes perfect persistence) would conclude that private schools add no greater value to learning than their public counterparts.

Many of the dynamic panel estimators that we explore identify the private school effect using children who switch schools. Figure 2 illustrates the patterns of achievement for these children. For each subject we plot two panels: the first containing children who start in public school and the second containing those who start in private school. We then graph achievement patterns for children who never switch, switch after third grade, and switch after fourth grade. For simplicity, we exclude children who switch back and forth between school types.

As the table at the bottom of the figure shows, very few children change schools. Only 48 children move from public to private schools in fourth grade, while 40 move in fifth grade. Consistent with the role of private schools serving primarily younger children, 167 children switch to public schools in fourth grade, and 160 switch in fifth grade. These numbers are roughly double the number of children available for our estimates that include controls, since only a random subset of children were surveyed regarding their family characteristics.

Even given the small number of children switching school types, Figure 2 provides preliminary evidence that the private school effect is not simply a cross-sectional phenomenon. In all three subjects, children who switch to private schools between third and fourth grade experience large achievement gains. Children switching from private schools to public schools exhibit similar achievement patterns, except reversed. Moving to a public school is associated with slower learning or even learning losses. Most gains or losses occur immediately after moving; once achievement converges to the new level, children experience parallel growth in public and private schools.

# 4.2 OLS and Dynamic Panel Value-Added Estimates

Tables 4 (English), 5 (Urdu), and 6 (mathematics) summarize our main value-added results. All estimates include the full set of controls in the child survey sample, the survey date, round (grade) dummies, and village fixed effects. For brevity, we only report the persistence and 13
private school coefficients.<sup>11</sup> We group the discussion of our results in three domains: estimates of the persistence coefficient, estimates of the private schooling coefficient, and regression diagnostics.

# 4.2.1 The persistence parameter

We immediately reject the hypothesis of perfect persistence ( _β_ = 1). Across all specifications and all subject (except M1 which imposes _β_ = 1), the estimated persistence coefficient is significantly lower than one, even in the specifications that correct for measurement error only and should be biased upward (M3 and M4). The typical lagged value-added model (M2), which assumes no omitted student heterogeneity and no measurement error, returns estimates between 0.52 and 0.58 for the persistence coefficient. Correcting only for measurement error by instrumenting using the two alternate subjects (M3), or using the analytical correction described by in the appendix (M4), increases the persistence coefficient to between 0.70 and 0.79, _consistent with significant measurement error attenuation_ . This estimate, however, remains biased upward by omitted heterogeneity.

Moving to our dynamic panel estimators, Panel B of each table gives the Arellano and Bond (1991) difference GMM estimates under the assumption that inputs are strictly exogenous (M5) or predetermined (M6). In English and Urdu, the persistence parameter falls to between 0.19 and 0.35. The estimates are (statistically) different from models that correct for measurement error only. _In other words, omitted heterogeneity in learning exists, and biases the static estimates upward_ . For mathematics, the estimated persistence coefficient is indistinguishable from zero, considerably below all the other estimates. These estimates are higher and somewhat more stable in the systems GMM approach summarized in Panel C (M7, M8).

With the addition of a conditional mean stationarity assumption (Panel D), we can more precisely estimate the persistence coefficient. In this model, we only use moments in levels to illustrate a dynamic panel estimator that improves over the lagged value-added model estimated by OLS but doesn’t require estimating a system of equations. The persistence coefficient rises substantially to between 0.39 and 0.56. This upward movement is consistent with a violation of the stationarity assumption (the fixed-effect still contributes to achievement growth) but an overall reduction in the omitted heterogeneity bias. Across the various dynamic panel models and subjects, estimates of the persistence parameter vary from 0.2 to 0.55. However the highest dynamic panel estimates come from assuming conditional mean stationary, which is likely too strong in the context of education.

> 11As discussed, time-invariant controls drop out of the differenced models. For the system and levels estimators we also assume, by necessity, that time-invariant controls are uncorrelated with the fixed effect or act as proxy variables.

14
# 4.2.2 The contribution of private schools

Assuming perfect persistence biases the private school coefficient downward. For English, the estimated private school effect in the restricted model that incorrectly assumes _β_ = 1 (Panel A, Table 4) is negative and significant. For Urdu and mathematics, the private school coefficient is small and insignificant or marginally significant (Panel A, Tables 5 and 6). By comparison, all the dynamic panel estimates are positive and statistically significant, with the exception of one of the difference GMM estimates, which is too weak to identify the private school effect with any precision.

Panel C (levels and differences GMM) illustrates the benefit of a systems approach. Adding a levels equation (Panel C, Tables 4-6), using the assumption that inputs are constantly correlated or uncorrelated with the omitted effects, reduces the standard errors for the private school coefficient while maintaining the assumption that inputs are predetermined but not strictly exogenous. Under the scenario that private school enrollment is constantly correlated with the omitted effect (M7), the private school coefficient is large: 0.19 to 0.32 standard deviations depending on the subject and statistically significant. This estimate allows for past achievement shocks to affect enrollment decisions but assumes that switching school type is uncorrelated with unobserved student heterogeneity. This is our preferred estimate.

An overarching theme in this analysis is that the persistence parameter influences the estimated private school effect but that it is rarely possible to get enough precision to distinguish estimates based on different exogeneity conditions. This is largely due to the small number of children switching between public and private schools in our sample. In Figure 3, we graph the relationship between both coefficients explicitly. Rather than estimating the persistence coefficient, we assume a specific rate and then estimate the value-added model. That is, we use _yit− βyi,t−_ 1 as the dependent variable. This provides a robustness check for any estimated effects, requires only two years of data, and eliminates the need for complicated measurement error corrections. (It assumes, however, that inputs are uncorrelated with the omitted learning heterogeneity.) As expected given the large baseline differences, the estimated private school effect strongly depends on the assumed persistence rate. Moving from the restricted valueadded model ( _β_ = 1) to the pooled cross-section model ( _β_ = 0) increases the estimated effect from negative or insignificant to large and significant. For most of the range of the persistence parameter, the private school effect is positive and significant, but pinning down the precise yearly contribution of private schooling depends on our assumptions about how children learn.

A couple of natural questions are how these estimates compare to the private-public differences reported in the cross-section and why the trajectories in Figure 1 are parallel even though the private school effect is positive. Controlling for observables suggests that after three years, children in private schools are 0.9 (English), 0.5 (Urdu), and 0.45 (mathematics) stan15
dard deviations ahead of their public school counterparts. If persistence is 0.4 and the yearly private school effect is 0.3, children’s trajectories will become parallel when that achievement gap reaches 0.5 (= 0 _._ 3 _/_ (1 _−_ 0 _._ 4)). This is roughly the gap we find in Urdu and mathematics. Any small disagreement, including the larger gap in English, may be attributable to baseline selection effects. Thus our results can consistently explain the large baseline gap in achievement, the parallel achievement trajectories in public and private schools, and the significant and ongoing positive private school effect.

## 4.2.3 Regression diagnostics

For many of the GMM estimates, Hansen’s _J_ test rejects the overidentifying restrictions implied by the model. This is troubling but not entirely unexpected. Different instruments may be identifying different local average treatment effects in the education context. For example, the portion of third grade achievement that remains correlated with fourth grade achievement may decay at a different rate than what was learned most recently. This is particularly true in an optimizing model of skill formation where parents smooth away shocks to achievement. In such a model, unexpected shocks to achievement, beyond measurement error, would fade more quickly than expected gains. Instrumenting using contemporaneous alternate subject scores will therefore more likely identify different parameters than instrumenting using previous year scores. Likewise, instrumenting using alternate lags and differenced achievement and inputs may also identify different effects. This type of heterogeneity is important and suggests that a richer model than a constant coefficient lagged value-added may be warranted.<sup>12</sup> Given the rejection of the overidentifying restrictions in some cases, the next section provides a series of robustness exercises around the estimation of the persistence parameter.

# 4.3 Robustness Checks

If our estimates are interpreted as forgetting, children lose over half of their achievement in a single year. For some subjects, such as mathematics, this fraction may be even larger. While the estimates reported may appear to be implausibly high, they match recent work on fade-out in value-added models, as well as the rapid fade-out observed in most educational interventions.

Table 7 summarizes six randomized (or quasi-randomized) interventions that followed children after the program ended. This follow-up enables estimation of both immediate and extended treatment effects. For the interventions summarized, the extended treatment effect represents test scores roughly one year after the particular program ended. For a number of

> 12Another common strategy to address potentially invalid instruments is to slowly reduce the instrument set, testing each subset, until the overidentification test is accepted or the model becomes just identified. We explored this approach but no clear story emerged. One result of note is that dropping the overidentifying inputs typically raises the the persistence coefficient slightly, to roughly 0.25 for math.

16
the interventions, the persistence coefficient is less than 0.10. In two interventions—learning incentives and grade retention—the coefficient is between 0.6 and 0.7. However, this higher level of persistence may in part be explained by the specific nature of these interventions.<sup>13</sup> Although the link between fade-out in experimental studies and the persistence parameter is not always exact, all the evidence suggests that current learning does not carry over to future learning without loss and, in fact, these losses may be substantial.

Exploring the magnitude of the _potential_ bias in a basic lagged value-added model can also give us some sense for whether our estimates are reasonable. Consider, for example, the bias in the regression _yit_ = _α_ + _βyi,t−_ 1 + _ηit_ , where we have omitted all potential inputs and corrected only for measurement error bias. Our estimates of this model suggest that the persistence coefficient is at most 0.8 to 0.9—far higher than our highest dynamic panel estimates of around 0.5. Is this discrepancy reasonable?

Aggregating all the omitted contemporaneous inputs into one variable _ηit_ implies the upward bias of the persistence coefficient is Cov( _ηit, yi,t−_ 1) _/_ Var( _yit−_ 1). If the correlation between inputs in any two periods is a constant _ρX_ , and all children in grade zero start from the same place, the persistence coefficient in a lagged value-added model for fourth grade will be biased upward to Figure 4 gives a graphical representation of this bias calculation. To read the graph, choose a true persistence coefficient, _β_ (the dotted lines), and a degree of correlation of inputs over time, _ρX_ (the horizontal axis). Given these choices, the y-axis reveals the persistence coefficient that a lagged value-added specification estimated by OLS would yield. Working with our estimates, if the true persistence effect, _β_ , is 0.4 and inputs are correlated only 0.6 over time, the (incorrectly) estimated _β_ will be 0.9. Given that the vast majority of inputs are fixed, this seems quite reasonable, and perhaps even too low.<sup>14</sup>

> 13In the case of grade retention, there is no real “post treatment” period since children always remain one grade behind after being retained. If one views grade retention as an ongoing multi-period treatment, then lasting effects can be consistent with low persistence. In the case of learning incentives, Kremer, Miguel and Thornton (2003) argue that student incentives increased effort (not just achievement) even after the program ended, leading to ongoing learning.

> 14Another way to get at the reasonableness of rapid fade out is motivated by Altonji, Elder and Taber’s (2005) assumption of equal selection on observed and unobserved variables. Absent controls and correcting for measurement error the persistence coefficient is 0.91, while the _R_<sup>2</sup> of the regression is 0.52. Adding controls raises the _R_<sup>2</sup> only modestly to 0.56 but at the same time reduces the estimated persistence coefficient to 0.74. Thus, just by explaining an additional four percent of the total variation, we reduced the persistence coefficient substantially. Assuming equal selection on observed and unobserved variables would lead to a persistence estimate _below_ our dynamic panel estimates.

17
# 4.4 Why Is Persistence So Low?

The low estimates of persistence are worrying, not only for program evaluation but for the more substantive issue of how to improve cognitive achievement. One major concern is that imperfect persistence is a psychometric testing issue, and therefore not a “true” feature of the learning dynamic. For instance, later test forms may capture fundamentally different latent traits than earlier test forms. To address this concern, we replicated our results using IRT scores based solely on a common set of items that appeared on every test form—our tests had a significant number of overlapping items in each year. Our results were similar using these scores, with the difference GMM persistence estimates in fact dropping slightly. The score equating methods used to create a single cardinal measure of learning therefore do not appear to be the driving force behind the low observed persistence.

Several other possible mechanical explanations for low persistence are also unlikely. First, artificial ceiling effects can appear like low persistence in models that use bounded scores. To address this concern, we exclusively use unbounded IRT scale scores and our exam is designed to maximize the variation over the entire range of observed abilities. Second, cheating, often driven by high stakes testing, can create artificially low estimates of persistence. Jacob and Levitt (2003), for instance, detect teacher cheating in part by looking for poor subsequent performance of students who made rapid gains. In our data, cheating is unlikely both because our test is relatively low stakes and because our project staff administered the exam directly to avoid this possibility. Third, critics of high stakes testing often argue that shallow “teaching to the test” leads to low persistence. This is also an unlikely explanation in our context; our exam is relatively low stakes, is not part of the standard educational infrastructure, and covers only subject matter that all students should know and that Pakistani parents generally demand.<sup>15</sup> We looked at a couple of other candidate explanations, but there are no “smoking guns” that could explain low persistence. Tables 8 and 9, for instance, present the results of a preliminary exercise that assesses whether low persistence arises from household and school responses to unexpected achievement shocks—an explanation that has different implications for cost-benefit analysis than simple forgetting.<sup>16</sup> Parents’ perceptions of their child’s performance reacts strongly to unexpected gains in achievement, but there is only weak evidence of substi15On average, the children tested at the end of Grade 3 could complete two-digit addition, but not subtraction or multiplication (in mathematics); recognize simple words (but not sentences) in the vernacular (Urdu); and recognize alphabets and match simple three-letter words to pictures in English.

16We examine whether inputs adjust to unexpected achievement shocks for roughly 650 children for whom we have detailed information from a survey collected at households. As a measure of the unexpected shock, we first compute the residual from a regression of fourth grade scores on third grade scores and a host of known controls. We then test whether this residual predicts changes between fourth and fifth grade in parents’ perceptions of the child’s performance, expenditure on school, and time spent helping children on homework, being tutored, doing homework, and playing.We instrument for the subject specific residuals using the alternate subject residuals to lessen measurement error attenuation.

18
tution effects. School expenditures do drop slightly as do the hours spent helping the child on his or her homework. Minutes spent playing increases, but tuition also increases. While some of these responses are in the direction of substitution, they are generally not statistically significant. Given the detailed household data we obtained, it suggests that household substitution is unlikely to be a main driver behind low persistence. This may not be particularly surprising given that very low achievement suggests that children may be below parents’ desired learning levels.

Table 9 explores the possibility that fade out captures teachers targeting poorly performing students. If teachers target poorly performing children in each classroom, persistence should be lower within schools than between schools. To test this hypothesis, we estimate a basic lagged value-added model with no controls and instrument for lagged achievement using lagged differences in alternate subjects. We estimate this model using average school scores (between school specification) and child scores measured in deviations for the school average (within school specification). If anything, the persistence coefficient is lower for the between school regressions, suggesting that within school targeting is not the primary reason for low persistence. Finally, Table 10 looks at whether heterogeneity in persistence can provide some hints about its origin(Semb and Ellis, 1994).<sup>17</sup> To obtain the most power possible, we estimate the valueadded model for specific sub-populations using the “predetermined inputs, uncorrelated effects, and conditionally stationary” based estimator (M10 of Tables 4, 5, and 6). Unfortunately, large standard errors make it difficult to find statistically different decay rates between groups. Learning in private schools seems to decay faster than learning in government schools, but the difference is not statistically significant. A similar pattern holds for richer families and children with educated parents. These results hint that learning decays faster for faster learners.

These candidate explanations do not yield a compelling story thus far; it could just be that children forget what they learned. Psychology and neuroscience provide some compelling evidence for this using laboratory experiments. Psychological research on the “curve of forgetting” dates back to Ebbinghaus’s (1885) seminal study on memorization. Rubin and Wenzel (1996) review the laboratory research spawned by this contribution. Semb and Ellis (1994) review classroom studies that test how much students remember after taking a course. Both literatures document the fragility of human memory. Cooper et al. (1996) studies the learning losses that children experience between spring and fall achievement tests. These losses are generally not as rapid as the effects we find, but the experiment is different: we estimate the depreciation with no inputs, whereas summer activities provide some stimulus, particularly for privileged

> 17To give one example, MacKenzie and White (1982) report that fade out for geographical knowledge was much higher for in-class exercises compared to field excursions (or, passive versus active learning.) Similarly, Rothstein (2008) finds heterogeneity in the long-run effects of teachers who produce equal short-run gains; Jacob, Lefgren and Sims (2008) estimate that teacher effects are only a third as persistent as achievement in general.

19
children.

# 5 Conclusion

In the absence of randomized studies, the value-added approach to estimating education production functions has gained momentum as a valid methodology for removing unobserved individual heterogeneity in assessing the contribution of specific programs or in understanding the contribution of school-level factors for learning (e.g. Boardman and Murnane, 1979; Hanushek, 1979; Todd and Wolpin, 2003; Hanushek, 2003; Doran and Izumi, 2004; McCaffrey, 2004; Gordon, Kane and Staiger, 2006). In such models, assumptions about learning persistence and unobserved heterogeneity play central roles. Our results reject both the assumption of perfect persistence required for the restricted value-added model and of no learning heterogeneity required for the lagged value-added model. Our results for Pakistan should illustrate the danger of incorrectly modeling or estimating education production functions: the restricted value-added model is fundamentally misspecified and can even yield wrong-signed estimates of a program’s impact. Underscoring the potential of affordable, mainstream, private schools in developing countries, we find that Pakistan’s private schools contribute roughly 0.25 standard deviations more to achievement each year than government schools, an effect greater than the average yearly gain between third and fourth grade.

Our estimate of persistence is consistent with recent work on teacher effects, with analytical and empirical estimates of the expected bias under OLS, and with experimental evidence of program fade out in developing and developed countries. But the economic interpretation still remains an open area of enquiry. Our context and test largely rule out mechanical explanations of low persistence. We also find little evidence that low persistence results from substitution by parents and teachers; the behavioral adjustments we are able to measure are unlikely to represent the primary reason achievement gains fade out. Simple forgetting, consistent with a large body of memory research in psychology, appears to be a likely explanation and hence a core component of education production functions, although more research is needed to provide direct evidence for it.

Our results also highlight that short evaluations, even when experimental, may yield little information about the cost-effectiveness of a program. Using the one or two year increase from a program gives an upper-bound on the longer term achievement gains. As our estimates suggest, and Table 7 confirms, we should expect program impacts to fade quickly. Calculating the internal rate of return by citing research linking test scores to earnings of young adults is therefore a doubtful proposition. The techniques described here, with three periods of data, can theoretically obtain a lower bound on cost-effectiveness by assuming exponential fade out. At the same time, the causes of fade out are equally important: if parents no longer need to hire 20
tutors or buy textbooks (the substitution interpretation of imperfect persistence), a program may be _cost_ -effective even if test scores fade out.

Moving forward, empirical estimates of education production functions may benefit from further unpacking persistence. Overall, the agenda pleads for a richer model of education and for empirical techniques for modelling the broader learning process, not simply to add nuance to our understanding of learning, but to get the most basic parameters right.

21
# A Analytical Corrections for Measurement Error

Consider the lagged value-added model where _yit_<sup>_∗_and</sup><sup>_y_</sup> _i,t_<sup>_∗_</sup> _−_ 1<sup>aretrueachievement,</sup><sup>_vit_istheerrorterm,andwehaveputasidethe</sup> possibility of omitted heterogeneity. Since achievement is a latent variable, we can only estimate it with error. Thus, we actually estimate and OLS is inconsistent because _yi,t−_ 1 is correlated with _εi,t−_ 1.

The analytic correction we apply replaces _yi,t−_ 1 with the best linear predictor ̃ where _λ_ and _ri,t−_ 1 are parameters. To see why this works, add and subtract _βyi,t−_ 1 from (11) to get where the second line follows from _yi,t_<sup>_∗_</sup> _−_ 1<sup>=</sup><sup>_yi,t−_1</sup><sup>_−εit−_1.Assumingexogeneitywithrespectto</sup> _vit_ + _εit_ , OLS is consistent if ̃ These conditions are automatically satisfied since the fitted value _yi,t−_ 1 and independent variables x _it_ are orthogonal to the residual _yi,t_<sup>_∗_</sup> _−_ 1<sup>_−y_ ̃</sup><sup>_i,t−_1bythedefinitionoftheprojection(13).</sup> The only difficulty is estimating the projection parameters _λ_ and _ri,t−_ 1 since the dependent variable _yi,t_<sup>_∗_</sup> _−_ 1<sup>isunobserved.Butitturnsoutthatwedonotneedtoobservethetruescore.</sup> The orthogonality conditions that define the projection (13) are Solving for _λ_ , we have 22
Plugging (20) into (19) and solving for _ri,t−_ 1 yields where m _x ≡_ 1 _−_ x _it_ (x<sup>_′_</sup> _it_<sup>x</sup><sup>_it_)</sup><sup>_−_1x</sup><sup>_′_</sup> _it_<sup>isanannihilatorvectorand</sup><sup>_ei,t−_1istheresidualfromaregres-</sup> sion of _yi,t−_ 1 on x _it_ . We can estimate _ri,t−_ 1 by computing _σe_<sup>2</sup> _i,t−_ 1<sup>fromtheregressionof</sup><sup>_yi,t−_1</sup> on x _it_ and taking _σε_<sup>2</sup> _i,t−_ 1<sup>fromIRT—i.e.,fromtheinverseFisherinformationmatrixof</sup><sup>_θ_ˆ.In-</sup> tuitively, _ri,t−_ 1 is the heteroscedastic reliability ratio of the score minus the variation explained by the independent variables. That is, the reliability ratio of _yi,t−_ 1 _− E_ [ _yi,t_<sup>_∗_</sup> _−_ 1<sup>_|_x</sup><sup>_it_]</sup><sup>_._</sup>

We compute _λ_ by plugging _ri,t−_ 1 into (20) to get The best predictor is This takes the familiar form of an empirical Bayes estimate that shrinks the observed score to the predicted mean. The shrinkage performs the same function as blowing up the coefficient using the reliability ratio after estimation. Here, however, our shrunken estimate provides a more efficient correction by using the full heteroscedastic error structure (Sullivan, 2001) .

Table A1 reports persistence coefficients corrected only for measurement error using the instrumental variable (using alternate subjects) and the analytical correction approach. Each cell contains the estimated coefficient on lagged achievement from a regression with no controls and the associated standard error. Where applicable, we also report the p-value for Hansen’s overidentification test statistic. This is possible for the instrumental variables estimators since we have three subject tests and three years of data.

Absent any correction (OLS), the estimated persistence coefficient ranges between 0.65 and 0.70. Instrumenting using alternate subjects raises the estimated coefficient significantly to 0.85 for English, 0.89 for mathematics, and 0.97 for Urdu. However, the overidentifying restriction is rejected at the one percent level in all three subjects. This suggests that measurement errors may be correlated across subjects at the same sitting and that this correlation may differ depending on the subject. By comparison, when we instrument for lagged achievement using double lagged scores we cannot reject the overidentifying restrictions. Unfortunately, in the context of dynamic panel methods, additional lags to address measurement error require _T_ = 4. The final line of Table A1 shows estimates based on our analytical correction around 0.9. Of course, all of these estimates remain biased upward by learning heterogeneity.

23