DISCUSSION PAPER I I I - ECONOMEmIC X T H O D S FOR AGRICULTURAL SUPPLY FERTILIZER VSE AND CROP RESPONSE by I Xchael J. Hartley Ecc~~cmis t I I I I I - I I I I I Development Research Department i Economics and Research Staff World aank . * I ! I '3 -i C ac l I .I II I ! I i The views presented here are those o f the author, and they should not be interpreted as reflecting those of the World Bank ,I I I I I METHODS FOR AGRICULTURAL SUPPLY UNDER UNCERTAINTY: *~ FERTILIZER USE AND CROP RESPONSE I I ~ by I Michael J. Hartley I , E c o n o m i s t P u b l i c E c o n o m i c s D i v i s i o n I~ D e v e i o p m e n t R e s e a r c h D e p a r t m e n t I I I V o r l d Bank a , e I * ~orthcornin~,Journal of Mathematical Analysis and Applications. July, L The paper develops general econometric methods for the treatmeht I of modelling decisionkaking under conditions of uncertainty and risk. - 4 I problem proposed to the author by John Holsen, Chief Economist, South ~ s j a - I Regional Office -- the problem of modelling f e r t i l i z e r use and crop respknse ~ i n the context of irrigated and/or rainfed water supply when t h e i r is unier- 8 I ; tainty and r i s k with respect t o weather conditions over the crop-cycle ahd I I the subsequent price a t marketing -- suggested the present approach. A s a by-product, we also derive the stochastic duality theory associated w i t h the present primal problem. ~ have included government extension service advice to farmers on cultivation practices; the introduction and distribution of High Yielding v a r i e t i e s ( W s ) to replace T ~ d i t i o u a lLocal Varieties (TLVs) of seed, etc.; the provis$on of I public worlcs and other complementary inputs to f e r t i l i z e r (such a s i r r i g a t i o n I I * and v a t e r c o n t r o l systems, pesticides, etc.); subsidies t o the output p$ice of I crops u t i l i z i n g f e r t i l i z e r; provision of credit to farmers t o r f e r t i l i z $ r . I purchases via government credit institutions and/or subsidies t o i n t e r e t i rates; etc. I . I I n the case where a l l other aspects regarding the use of f e r t j l i z e r I I are l e f t t o the voluntary decisions of farmers, assessment of the costs!and I benefits associated with any c ~ l l e c t i o nof government programs requires ;an explicit behavioral-model of the farmer's response to the resulting 'stllucture I 'of inceatives." i~ We shall consider these issues within the context of a s t a t i c I I neoclassical mod-: of individual farmer behavior. We shall adopt the maintained hypothesis that tho bbjective of each farm unit is to maximize the conditional expectation of profit with respect to the levels of all v a r a b l e 4 inputs (including f e r t i l i z e r ) , subject t o the constraint of a given technology for a given crop, and given prices 'of a l l variable inputs a d output. ~~ I refer to t h i s a s the Standard Neoclassical b d e l (SNM) of agrLcultura1 production sod input.>emand-ncudn the derived concepts of output supply, .s I cost and profit funceons, obtained a s a consequence of duality theory. I I a I I Econometric'methods appropriate for the estimation of the par+eters I I ,:ontained within the production function of a SNM f o r a specific crop w i l l be I treated in section 2. We shall extend the present state-of-the- art w i t 4 I I respect to t h i s problem by exhibiting econometric and computational methbds I I behaved ,"but o t h e r w i ~ cunrestricted parametric specifications of the I production function required under any maintained hypothesis. Our apprca I w i l l be couched entirely i n terms of the so- called "primal problem," invo ving C a profit equation (containing the production function) the derived syItem I I of input denand functions (obtained, say, under maximization of the I I conditional expectation of profits), rather than postulating a functional form for either the p r o f i t function or cost function, to be estimated i n con junction vit h the derived input demand system-the so-called "dual" to the I primal problemL/ We believe this to be a more useful approach since i information to guide the choice of functional forms is more readily availeble for the production function.?! Further, once th; parameters i n the produ function of the primal problem. have been estimated, the implied profit, l and output supply density functions can be derived and their expectat1ons;and I other moments evaluated. This has not been accepted current practice- except i n s p e c i a l i ~ e d 4 . cases- because of the apparent necessity of obtaining explicit closed-fo analytic solutions to the system of marginal productivity equations i n p l i j d by . e C a par=tic~;lar choice of a production function. Tradition has required this/ i n o r d e r that: (a) t h e d e r i v e d functinna!frln Fr)t tho I +?ST:+ ? ~ n n n J~ n - ~ , ? : * ~ r I (including the resulting homogeneity conditions %- prices) is consistent Qth I .. #he postulated production function, and (b) t o s&cify the s e t of "cross- - i I I I i I -I/ I See, e.g., Sidhu and Baanante I19791 for a recent example of the use oif a 1 translog profit function. -2/ Such information is available from the agronomic l i t e r a t u r e on cultivation practices for each crop- often obtained from experimental t r i a l s . 1 I equation" restrictions, functionally relating every demand equation ~ a r + e t e r I to the set of original parameters i n the production function. By exhibi~ting I an algorithm, i n vhich the solution values to 3well-behaved optinizat&on I 2roble.m specified i n the m+intained hypothesis automatically s a t i s f y botb (a) 1 l and (b), we eliminate a major constraint on the set of "econometrically y . I estimable" parametric forms for a productiou function and a s ~ o c i a t e dinp t demanZ functions. This obviates the necessity of proceeding t o the b I "econometrics of duality" approach- since the associated profit, cost an I supply density functions can then be derived d i r e c t l y from within the I I stochastic structure of the p r i d problem and, hence, t h e i r expectation+ and I I . I variances can be numericallz evaluated. I 1 I The econometric approach of section 2 w i l l then be extended to! jt several cases which a r e of particular relevance to the analysis of f e r t i i z e r use and crop response. In section 3 we consider the comparison of irri ed and rainfed crops, and t r e a t the typical situation i n which the availabiiity I of irrigation and water-management systems is a binary, exogenously-dete event- say, as a consequence of the differential access of farmers to pu 4 f a c i l i t i e s . In the situation where water-managesent f a c i l i t i e s a r e avai able, I . * C w e assume that t h i s permits the farmer t o control the level of soil-mois/ure I accessible to the crop. This not only provides the farmer wfth Rn addft4onal control variable. but also eliminate= much of the uncertainty with r'9 p e i t t o w - I 1 unpredictable weathe; conditions. Since water is a complementary i n k t , both E - - w I of these affect the ;sage of and response ta f e r t i l i z e r . In contrast, firmer* without access to such f a c i l i t i e s must t r e a t water (or soil-moisture) as la P s t a t e variable, beyond their present control. Here, farmers maximize ex ected profits with control over one less variable input than those wfth access t o I I i r r i g a t i o n f a c i l i t i e s . This, by appeal to ,the 'theory of the second-best'," ~ c e t c r i s paribus, w i l l reduce the expected levels of p r o f i t s and, i n pract/ice, I also f e r t i l i z e r wage. As an i l l u s t r a t i o n of these issuer, w e consider the I irrigated case, i n which the actual levels of water a r-e not observed- thu - I representing a l a t e n t variable input. This is contrasted with the rainfe! I - case, i n which the amount of r a i n f a l l is reported but uncontrolled. - ~ By I I pooling data on both types of farmers cultivating a given crop t o e s t i m a t ~the parameters within t h e i r common production functions, we may, i n principle, quantitatively assess the benefits of watermanagement public works projects I and t h e i r c r o p s p e c i f i c effect on f e r t i l i z e r use and yields. I I Our discussioi of the tainfed case i n section 3 extends naturaliy to I the problem of decision-making under uncertainty. Thus, i n section 4, we^ - I ~ consider a general econometric approach, requiting the additional I spettffication of the (subjective) probability density function ,f a l l s t a I variables to which the farmer attaches uncertainty. This 'endogenizesw t4e sources of uncertainty for all 'Ptate variables within the scope of the model i and, by "integrating out' their combined effect, we can obtain ML estimat s of the original (production function) parameters of interest, i n addition to the ~ I parameters i n the marginal distribution of these "new" l a t e n t efidogenous variables. There are two principal sources of uncertainty i n the. economics 1of I crop cultivation. The f i r s t arises i n the case'of rainfed crops. Here, input I) decisions over the crop's planting t o harvest lffe-span must invariably b taken prior to harvest. This, i n conjunction with the fact that weather is, apart f;;m seasonal trends, a random s t a t e of nature implies that, under dhe maintained hypothesis, f e r t i l i z e r use and crop yields w i l l , i n general, "fY with a l l of the parameters of the r a i n f a l l distribution, and not just {he with - - I expected value of r a i n f a l l , a s typical. practice suggests. To the extedt that there a r e systematic climatic differellces a t d i f f e r e n t locations, and racional behavior on the part of farmers, our methods permit modeling (he a ~ I range of cross- sectional f e r t i l i z e r use and crop responses to such I uncertainties. h e may a l s o t r e a t uncertainty with respect t o the price of dlutput a t harvest-time along similar lines. Indeed, it is precisely t h i s pro which occasioned the seminal contribution of Nerlove [I9581 ro the I on uncertain future prices, whers farmers a r e assumed to form "adaptivel expectations" involving the notion of a "normal" (in the sense of long-run) I formulation^ price. Tile dynamics of supply response resulting from such is 1 now a common theme within the post-Nerlovian agricultural supply l i t e r a t u r e - see, again, Askari and Cummings (19771. The adaptivcexpectatiolla hypothesis is, however, but one of many currently popular expectations formulation I- Other time- series specifications include the partial- adjustment model; autoregressive integrated moving-average (ARIMA) processes f o r prices; 1 rational and so-called "pseudo" rational price expectations models; etcl-see 4 i Nerlove, ~ r e t h i rand Carvalho (1979, ch: XIII]. A l l such choices, hove er, boil down to making an assumption regarding the ~ a r a m e t r i cform of the , 'ixpectation of 'the price a t harvest- time, based upon information availa Ie to - I, I *The - farmer a t the t i m e a decision regarding variable inputs must be pad?. I &nce, by adding a normal random error component, these candidates t a l l within I If our approach t c estimation- requiring the probability density function future prices. b It is also important to note that, by integrating out the effe t of a l l l a t e n t endogenous variables, ve obtain the marginal distribution of ~~ I p r o f i t s end all obsened decision variables. Once t h i s is obtained, we &an define the conditional density of profits, given the values of a l l d e c i s o n i variables. I In general, &moments of t h i s conditional p.d.f. are requi ed i n order t o completely characterize the consequences for decision making unjer I uncertainty. There is, however, no reason t o assume, a priori, t h a t f a d e r s are "neutral with respect to risk." To the extent that a t t i t u d e s toward r i s k l a r e an important common behavioral characteristic, and presuming that sudh risks vary with respect to time and place, our approach may easily be 1 generalized to consider, instead, an objective function containing the 1 conditional expectation and -conditional variance of profits a s arguments within a specified parametric form of the fanner's u t i l i t y function. ~ h j s I permits a direct econometric treatment of the consequences of r i s k and I uncertainty on f e r t i l i z e r and other input use. I We should also stress, to the extent that: (1) the standard treatment of decision-making under tmdertainty has been largely couched terms of theories regarding the expected values ~ .- of the uncertain variables, and (2) economic units a r e not -indiffeient t o risk, the standard approac the problem w i l l therefore provide an incomplete explanation. - I In additfop, I since, except i n highly .apecializeJ cases, the structure of a neoclassicall L I model represents a n o n l i k a r transformation of the s t a t e variables into tbe 't I decision variables (incluaing output and prof it) it follows that the enptcted , value of (say) profit can not i n general be calculated solely from infonn tion - a bated on the expectations of the uncertain s t a t e variables--as i n the cask of I "certainty equivalence," see Theil (19611. This also suggests that theor e s of so-called rational expectations- see, e.g., Lucas and Sargent (19811-are I unlikely t o provide a complete treatment of behavior under uncertaint , i Rather, we need, i n addition, rational. variance models, rational t h i r 1- moment model., etc. This suggests a convergence t o our-approach for dealing ivith uncertainty- in which we seek t o parametrize only the I probability density function of the uncertain state variables, and estimate I I its parameters based only on information available to the fanner up td the I time that decisions are required, ~ I 2. The Standard Seoclassical Model of F e r t i l i z e r Use I and Response for a Given Crop: I I In t h i s section we. present the standard neoclassicai model a f ~~ I ' f e r t i l i z e r use m d output response for a given crop. Kc consider the maintained hypothesis of (expected) prof1t maxfmization for the indivqdual I farm unit i n vhich che prices of all vatiable inputs and output a r e gdven and I known t o the farmer with certainty at the t i m e of decision making. We shall I consider a situation i n which the land-use decision has already been dade, 1/ I The problem, here, is t o decide upon the optinal levels of the decisiohr I variables-a. set of variable inputs, including fertilizer-dgiven knovlCdge of - the production technology for the crop and given knowledge 02 a l l pric other ecate variables which serve to disting+sh the individual farm u its, n L In subsequent sections we shall extend t h i s p d a l to incorporate khe vbrious ? I- ! other features noted i n section 1. I I 1 ~ -1/ Extension of the model to a multi-crop setting in which a land-use decision is also made will not be treated i n the present paper and will be - discussed i n a sequel. - (a) Notation 1 I I The general situation requires each Earm unit to choose b 'tween t ..., 04 -1, c o p U e ahall consider the (conditional) question optiaal f e r t i l i z e r usage given that-the fasm unit, i, has decided to allcca portion of i z s available land, aij' to the arbitrary crop j. Let i observation index for the farm unit, 1-l,Z,. ..,N, where N denotss t?e d z e of a random sample of such units. Let s ( i ) denote the farm unit denote the day and ~(1) denoto the year of observation. Thus, permite the random sample to constitute a pure crosa-section, a pure time- series, or, more importantly, a mixed tCmeseries/crosssection data sample. Further, the analysis can easily be restructured f o r the case of str t i f i e d C random samples associated with each yesr of obserration, by the intrbduction of suitable sampling weights attached t6 each sanple observatioa. I Foz any variable, (say) z, we shall use the shorthand notation, xi, t o refer to an o b s e r ~ a t i o non farm unit, s(i), and, where relevant t b crop j, the notation, x etc. Further, within year ~ ( i ) and for crop, j/i),we i j' , . shall let tj(i) denote :he date of harvest; tj(l) - L (1)denote t h e d a t e of j I planting and R (1) denote the length (in days) of the cr3p cycle. Fo any j a observation, 1, and crop, j, we define the following variables:- 1 F = F K, --vector of f e r t f l t z c r i q ~ v (tC~: J Y ! ~ , ~ . ? ~ :7: ~ 5 r - e -1 . *; I 4 I t9 I= element vector of N,'P,K nutrient levels), xO = KO --vector of other variable inputs, I -1il j ? i , -1/ Different varieties of a given crop, to the extent that t h e i r prdductioa functions d i f f e r , should be viewed as different 'crops.' The nodation here (particularly the subscript j) is also applicable to the seduel, i n which c r o r c h o i c e vill also be modelled. - = area of land under cultivation, "ij a = L vector of fixed non-land inputs and other noatprice state 3- I variables, = output, - = price of output, ~ P i j, vector of f e r t i l i z e r prices O = K0 vector of other variable input prices , lij J Crj = t o t a l fixed costs and = total profit . We now turn to the Standard Neoclassical Model. ~ ~ (b) The Deterministic Standard ~ e o c l a s s i c a lModel A deterministic conditional neoclassical model of f e r t i ' izer usage i and output response, given that a positive area of land has been e,llocated t o -crop I j, assumes that farmers maximize profits with reipect to a l l variable -- input Levels associated with the area of land, aij. t The technolo y associated 1 with crop j, here taken as invariant over the sample period, may e summal-ized by a production function, I where x0 d 9 is an -ij -= I-ij XXIJxoO] is a K 1-vector-- of variable input levels, an -ij - -1 M vector of parameters of interest. Expsted profit associated d i t h crop j . j- I may then be defined F' 0' C vhere q, s [q-ij q 1 is the Kj-vector of a l l variable input p r i c e s . The I I A, -ij maintained hypothesis is that farmers maximize the function, @ w i t h respect j' I t o x for given q?aluesof aij, zij, pi, and zij. %b,f Under the customa y -il I assumption of continuous d i f f e r e n t i a b i l i t y f o r f with respect to t h i s 1 leads to the familiar marginal productivity conditions, l 1 o r , i n the alternative o m , exhibiting the role of relative prices, where 31, * 1 is the vector of normalized or relative prices-- 9'j U (19781, using' the crop price, pi,, a s aukeraire. Under the conditio implicit function theorem, a solution, , for the optimal G v e l s exists, and, at l e a s t locally, may be calculated via L 1 *and a,quation system, where xf denotes a K 3 - I -.. The major ?roblem with the generality of t h i s theory- at b a s t in terms of current practice i n both mathematical economics and ecdnometric I empirical work-Is that, even f o r "well-behaved" production functionb ( i .e., I those satisfying the typical neoclassicai second-order conditions fok a I I mm-mum for any 9 which solves (2.3b)), the implicit function theoken is -ij I I only an existence-theorem, and does not necessarily permit closed-fob analytic solutions for the K input demand functions 1 --. i?. Indeed, onby f o r I c e r t a i n specific choices f o r the production function, 'J' is it poss ble to analytically solve (2.3 b) uniquely for the optimal x* functions, as, ) -j ~ e.g., i n the Cobb-Dooglas and quadratic cases. However, given a particular w 11- i behaved functional form for f j, and given specific values for tte vafiables, 4tj, aid, E ~&~t h e, parameters, 9 a suitable algorithm, progr -j' computer, can e a s i l y calculate numerical solution values f o r the opt$mal i n p u t ' levels, x* , such that I I -ij ~ Our purpose i n the remaf~derof t h i s section w i l l be t o inkroduca r I stochnstie version of the Standard ~ e d c l n s e i c a lModel and to e*hibit I computational ~ e t h o d sfor the calculation of maximum likelihood est3 a t e s of the parameters of the p;oduction function contained within the prof i k equation ~ & I and system of input demznd equations regardless of whether o r not th4 choice ir I of the production f unct fon permi I marginal productivity equations of the form (2.4). Our approach w i l d be to employ, instead, an i t e r a t i v e procedure, such a s the avido on-~letcheJ - P ~ ~ ~ ~ I I a (1963, 1964, 19681 (DFP) algorithm, to solve t'.o-.(expected) p r o f i t rn ximi- zation problem, (2.5), for the optimal x* values, given values of t +j + parameters, 9 and the s t a t e variables,Sj,aij and zij, associate -1 ' each individual sample member. This set of "inner" maximization pro lems w i l l 9 then be embedded within an "outer" maximization problem, i n which t h log- likelihood function of the parameters is maximized given the data, 1 l z i j qij. pijs aijs aj}* -and the solution values, By i terating i between the "inner" and "outer" maximization problems such that the ikelihood d of the observed sample of data under the maintained hypothesis is a onotone increasing sequence, convergence t o the maximum likelihood estimates of the parameters can be establishod- see Hartley 11981b] l! This simple divice- replacement of the derived analytic solution for the optimal input dmand 9 functions by an unconstrained optimization algorithm- not only penni s ti estimation of zgpply response model? f o r 3well-behaved choice of dhe production function (and automatically imposes a l l cross- equation ree!t r i c t i o n s o 4 on the parameters of the demand system), but also permits extension a rich array of alternative maintained' hypotheses encompassing problems of r i s k and I uncertainty. (c) A Stochastic Standard Neoclassical Hodel: e I t We now turn to a stochastic specification of the Standard - tieoclaesical Model to implement t h i s approach We shall assume that 'f choose input levels t o maximize the condition& expectation, = !pr o f i t s , n , given the vector of decision var9ableas, x ij -ij' -1/ I In keceral, convergence to a local maximum of the likelihood func ion is a l l that can be guaranteed --Hartley see [1981b]. - , I~ I I uncertainty regarding the level of p r o f i t s , xij, t h a t w i l l o b t a i n 2 I ti I follows that the conditional p.d.f. of xij given -i j' is given by x whereas the marginal p.d.f. of x -4is 1 (2. lb) Hence, the joint p.d.f. of x and x is the (K + 1)- variate normal I1 i J -iJ 1 I density, It should be stressed that the foregoing analysis is upon the fact that farmers have devoted a positive land area, a to crqp I j. We therefore define the index set of relevant observt.tions, I containiug N such o b ~ e r v aions r 2 The conditional 1 given the decision to cultivate crop j, is thus defined 5y . , * 1 P I -1/ See Zellner, Kmenta and Dr'eze (19661 for a Cobb-Douglas version of t h s specification. -2 1 . Observations with a = 0 should be excluded from the conditional estimation of 9 anadz In the sequel, however, we exhibit vhich u t i l i z e d e i o i o b a t i o n that a particular farmer has cultivate crop j. bhere C is the implied symmetric covariance matrix of the joint p.d.fl - 1 ~ of eij a d lij'i.e., consistent with (2.8) and (2.10)' i.e., II ..aloj ..a.llj .. . "'%lCjj .. bKjOj :a*Kjlj..* bK jKjj 1 The problem, in general, is to provide econometric methods a d caputational algorithms to permit the maximization of the log-likelihdod function L*(8 C ) with respect to 0 and C without requiring a closed-form j - j 1' --j 1 analytic solution for the function, x*, of (2.4), -j which solves the mar productivity equations (2.3b)' but, rather, requires only that numeric 1-4- solu!Aons to the optimization problem (2.5). (d) The Algorithm I,= We are now in a position to state an algorithm by which the imum n n Likelihood (ML) estdmates, 8 and -j y be calculate I I n 8 19: :vecr.1 t i . 1 5 ) .. ~ - m - - - d -J I 3 denote &he R Z H l)(K +2) parameter vector, containing the) Mj- * - ' 1 1+z(K1 + j rector of structural parameters in the production function, 8 and the(I .E e -1 ' additional parameters in the lower triangle of the symmetric covarianc matrix, C1 'where vecC denotes the functionally independent elements -1 in C organized into a vector, 1 ' -4 w CI( r( \d n --*--+--- 9) rlscec 325 5 U r( t cecum 5 0, m uo4rl&d u rlo 9) 0 a~m U cd (r d (r 4 K 0 "a" " 2 G (d u CJ rl Ua, LC rl 0 cd rl U 9) (r a 4 0 2 50 *I *. 9) M U !? ' e tii may be used a s an estimate of the (asmptotic) covariance matrix of a I It is hportant to stress-that, in the context of the present 1 problem, the algorithm applies whether o r not the explicit analytic solu ions - t o &he"inner" maximization p r o b l a s i n ( I ) of Step (n.1.r) can be calcu ated i for any a(nsr). If the konner is the case, the solution values, x -J be calculated d i r e c t l y using the explicit expressions for the optimal da$mnd functions, x*(nsr), of (2.4). Otherwise, we must employ an i t e r a t i v e -3 procedure (such as the DFP algorithm), with the use of the actual decisic!m vectors', x -iJ' a s the i n i t i a l vector t o commence the algorithm for the "i marbization problem, (2.5) -L! I n the present case it w i l l also be no - ..-. - that the solution vectors, {& 1, are independent . l d -iJ the parameters, 8 Hence, i n Step (n.l.r), since x*(','I x*(n~o) for -1 -1) --Ij values r = M + 1, M + 2,...,R (corresponding to the perturbations of J J J the C elements), the calculations i n ( i ) may be by-passed for r > M dn J J' I other problems (see the discussion of r i s k and uncertainty i n section 4 1 belov), we shall confront situation. in ~ h i c h ~ t helements of E e e n t r r tt/e J optimal decision vectors, and nere, evaluation of x*(npr) for a l l r val es - -11 Y ,,-- . _ . - - J C - , Y C . L . '2 In other contexts we ahall confront cases i n which the log- - - l * likelihood function, L*(8 ,C f, and/or the conditional expectation of I I .I . 1-J 1- I I i -1/ I f the actual input vector, sl,is not "too far*' from the optimal ve tor, , then we should obtain -t h e global maximum for (2.5) even i.n cases r e f j permits multiple local maxima as solutions to (2.3b). I p r o f i t s , E[n ( x ] of (2.6), themselves, can only be evaluated by 5 i j -ij numerical aethods. Our algorithm applies, mutatis mutandis, t o thesd situations. Further, i n problems involving decisiolrrnaking under w e may even postulate a more general objective function--e.g., a - function of 3 specific parametric form i n which both the conditional expectation and the conditional variance (risk) of profits (or order conditional moments) enter as arguments. There can all be treaked within the confines of ouz approach for general awell-behavad" choice$ of functional forms. ~ l i u s t r a t i o n sof such cases w i l l be given i n subsebuant . sections (e) Stochastic Duality Theory: We shall cortlplete t h i s section with a discussioa of the theory of the p r o f i t , oueput supply and cost density functions for th representative micro farm unit. However, rather than postulating func;tional - forms for the (expected) profit and/or cost functions, and deriving t4e input I d.emand and output supply functions corresponding to them under perfec l y I competitive profit-maximization o r cost-minimization, we shall, methods for any s e t of s t a t e variables, i j ~ j ~ P iand j qij. We begin with the marginal p.d.f. for proflts ~learld,frcm %if 8 the joint p.d.f. of x and x given i n ( 2 . 1 1 ~ ) weohave ~ i j +j with marginal expectation, and marginal variance, ., In general) even though h is a notmal p.d. f in the absence of clo:? expressions for E[sr ] and V[x ] e~aluationof (2.23), , (2.24a) and ( .24b), iJ I 1) f w i l l require the use of numerical integration (see, e.g., Haber [1970)/ using I the estimates, 8 and C This w i l l provide point estimates of the de -j 1 expectation and variance of p r o f i t s for any choice of the s t a t e variab es, t I ! a i j ~ ~ i j . P i j ad Ldjm This provides an implicit representation of the p r o f i t densi y I I function h and expectation function, E[x 1, treated previously i n the I I x, i J i context of a deterministic SIN. Indeed, a s Lau [I9781 has shown, the se of the Legendre transformation, t o obtain a closed Pcrm expression for th{ form of the profit function corresponding to a s i v e n prodcction function un er d p r o f i t maximization, La, i n practice, applicable t o only z iimited set of choices for the ~ r i g i n a lproduction function. . I f such analytic expressions * I are infeasible, a flexible functional form (say, the trans-log) is o f t e selected- f o t tye -rqff.Y z u n c i 4 , 7 ~ ,an? ~ C ~ ~ T - ?~- CC IO ~Z: ~ ~ > - _ > L O : I::'.I ? ' - -9 resulting input demand functions- see, e.g., Sidhu and Baanante - approac& however, obviates t h e necessity for such e can numerically evaluate (2.2.3) , (2.24 a) and (2.24 b) for well-behavkd choice of the production function, fj, i n (2.1). - 22- I We next consider the output supply, yi j, defined by (2.1) Cud note & that provided the conditions of the implicit function theorem are s a i s f i e d ..,K (i.e., suppose af /ax +0 for some k-1,. ), then we are assur j i j k J existence of an inverse function of the form with k=l (say): where x' [xijl zij2], with Jacobi.ur, -13 Thus, the marginal p.d.f. of yij is given by and can be evaluated by numerical methods. It follows that the of the quantity of crop j eapplied is given by the function of ~ I v l t h variance, vhich require further numerical integration for evaluation as Q ljandid. -?- Finally, we consider the stochastic version of the ex e c t e cost function. L e t Cij=GjE~~+ ci be the t o t a l cost and ~:onsidert h ;l transf ormatlon, Since the conditional p.d.f. of n given x haa already been def - i j -ij (2.11), w e have Hence the marginal p.d.f for Cij is and the general analysis follows similar l i n e s for evaluation of E A A any choice of the stat= variable values 8.t 0 and C etc. -3 1' i In short, &exploiting the distribution theory (and ass ciated numerical integrations) i n the above primal problem xe have abtaindd the I marginal p.d.f.'s cf all duality concepts uhich are internally con I the stochastic version of the primal SNM-:hereby obviating the nadd for the - -- so-called profit and cost- function econometric approaches- .see, e.a/., ~ a u 4 [I9781 acd McFadden [1978], unless we a r e i n the peculiar s i t u a t i o of having prior informscion regarding the specification of these l a t t e r cor;cefpts, inotead of the production function, i t s e l f . - '3 1 I - 3. I r r i ~ a t 2 dversus Pninfed Crops: - I I I I! - r The f i r s t extension of the Standard IIeoclas'.ical Model t o be considezed is the problem of modelling f e r t i l i z e r use and response crop i n one cf two situations: Either J +? (1) the farmer has access t o i r r i g a t i o n (and water manageme t) systems t o control water inputs, but the amounts of water availabl to the crop are not recorded, or (2) the farmer has no control over water inputs ( i n the abs nce of t i r r i g a t i o n .ad water management systerns) and, instead, r e l i e s upon , * the actual r a i n f a l l , which is recorded. C ' The problem is how to estimate the parameters i n the production unction and determine the demand for f e r t i l i z e r , etc., i n the case of data f crop, j, f o r farmers i n e i t h e r of these circumstances. I To simplify our analysis we shall assume that the presence or I absenre of i r r i g a t i o n and water Panagement systems is exogenous1 jr I determined. Thus we map defineathe binary variable, -- - i 1 i f irrigaticn and water management system a r e available to unit i on crop j 0 othervise. 1 Suppose we now mite the production function f o r the particular top j as 4 = f ( X ' 0 1, Y i j j -ijBuijBa: jB%jl-j .# ;-.. . . - e -e.13L,J a , ~ - i a o - ~ z c r~bz;;,.G ;atc. ,npuL. Iil case { L ) i n which '9 Z 6' = 1, the variable wis is an unobserved decision variable i n he contra1 of i j C the farmer, and, hence, i. to be treated as an additional variab e i n u t ;thus , augmenting the X j-vector, 3j. On the other hand, i n case (2), +ere = 0, the variable, w , plays the role of an observed s t a t e ariable, and iJ hence augments the Lj-vector, z,j. We ccnsider each case i n tur (a) Case (1): Irrigated Crops I n case (1) we have the model consisting of: c,ij=5 1 where denotes the normalized price of water (yos %,ij i f uater is a free public good) and, for any 9 the values of x and -f ' -ij w('Ir are the joint solutions to the problem of -imi~ing the cond .il expectation of profit, . with respect to both xt and w Let denote the marginal p.d.f. of the j- il -u* j errors i n the complete set of input demand equations, defined by P I E acd l e t g denote the El:,j defined as (1) gclu, j irij -j (3-5) ( 1 2- 1 - * - (1)) -j 1 zj 5" - I - - - - "K+ Z - d v-- 5 LI 0 cU w r'l r'l r* t! a t.~ 3 u Q) m at+. .r) 3 cu d d g 5 5 Ip 3* cI u 0 rn .m ld 9) * dm C, w4 d.r) C, 3* I 3 n L) .r) w r'l a * h -4 d +In n 3* Ids" 3&> a rl � ! 2 o ld a "r 0 wd 1 5 u4 rd KT .r) r .. u g a z w 0 d g rn a ti4 .? n 3a 3" .r) 4 .r) 4 .r) -.r) V r'l P. .'s 2 ld r( d.0 0 rna 9) bl I Wdrn a a P) a a! a at a U 8 513 H U 3-$ 4 dm rd cl m 0 P w at 0 .r) tt .r) 4 1 u b. cU !i 0, U n 5 3" w r'l a0 n .r) -l sz a d 4 a 49- N 9)cdWd 0 ld rlma 9) a UUd 3 +&-- S U 0 4 3 5 0) P) r( 4 0 a I n general, the 'incomplete' data sample (due to a l a t e n t dbpendent variable, wi ) means that certain of '.he parameters i n 9 and ~j( l ) ,w 11 be C j -j inestimable vithout further information. Though it may not be possib e tc i. decide analytically exactly which parameters t h i s w i l l affect, t h i s d es not A -of e preclude use of our algorithm. Feasible i n i t i a l values f o r all th complete data parameters, 9(') and c ( ~ ) ( ~ )w,i l l still have t o be cho -3 j i t e r a t i o n may proceed. A l l parameters, , i n i t e r e t i o n n such t =j (L*(n*r) - is zero, ~ 1 . 2 , . , are clearly inestimable, a n d thus may .. j 1 be deleted from ( i ) of Step (n.1.r). For the remainder of the giaramet/ers, functional dependencies may still exist- in which case the log- likeli ood h surface will have a "ridge" at the global maximum, represenring the 4 rameter sibspace within which all parameter points are observationally equivallent. 1/ The algorithm then w i l l simply converge to one such solution on the rddge, - Second, we may employ a Moore-Panrose generalized inverse o$ H(1) I j (given i n (2.20) and (2.22)) i f it is of i n i e r e s t to determine, f o r a 1 I particular functional form of f J i n (3.1), exactly how many independent parameters a r e present. This w i l l be revealed -rank of H(') and can be calculsted by means of the Singular Value L j Decomposition algorithm--see Goluh r19681. 7'-"r d , ' - "0nvn! Zho ' z ~ ; l z ; : j zi.23 .. ;.-J 3 , ,:5 1-3 ~e-;dLfi LC L~ie -3 decision variables, thooagh kcown t o be of i=poetsnco, do no; happen to be - t * present in the data sample, it may be empting to, instead, apply the ode1 of m section 2, treating the x vector as the complete set of decision var able. -il f h As is well-known, t h i s r e s u l t s i n a specification error in employing t e -1/ For a more complete discussion of identification problems caused b unobserved dependent variable i n such contexts, see Hartley [1981a I . I production function, (2.1), when i n fact (3.1) is know to hold. $us I I 1 1 I the and estimates w i l l not be consistent. Rather, to correct for such a -1 1 1 , the specification error, one should endure the additional computati I I by estimating the 'observationally relevant" parameters i n i,of (31.1) and I Fourth, i f wij is also observed in case (1)--as, e.g., wh n farm- t level soil- moisture measurements are available t o construct an ind mount of water accessible to the area, we a r e back to the correct model i n section 2, but using the density, h(1) q j .=ij~zij~~ij (b) Case (2): Rainfed Crops We now consider a model in which wij is a state variable, and hence 4 - beyond the control of the farm unit; but where 'wi is observed, a s n index of l I data on r a i n f a l l over the cropcycle. In t h i s situation, provided ji, -is know t o the farmer at the tine the amer must decide upon h i s optdmal f e r t i l i z e r and ofher variable input levels, xP and xu then we -ij 4 3 ' i with "complete data" (as i n section 2), but with a structure given y: vij, now enters as an argument i n the input demand functions, normalized price, %,ij. Further, while the production crop j are identical--and given by equation (3.1) i n both cases--the demand functions f o r f e r t i l i z e r and other inputs, x('I* -j versus xj-(2)*, w i l l d i,f f e r depending upon whether or not i r r i g a t i o n is present. Thus, i n case (2) we proceed a s i n section 2.c, but with an - additional s t a t e variable, vi Here, defines the conditional expectation of profit, as opposed to x(I)* of (3.3) in i J case (1). As before, i f and then h(2) (2) "'"I - - (2)* -g6;1, j(mj Xi, kij - zij (2)*) . I 0 . u ) ~ , ~ , j ( * i j s ~ i9j 8u,ij(zij-zij - )- I Thus, the joint p.d.f. of E ( ~ a) d u ( ~ )is given by +j - - ( 2 ) ( ~ ( 2 ) (2)) N ( o , z ~ ~ ) )ri2) ; =j i j *xij ,--* 4 and the log- likelihood function f o r cage (2) is given by w I . =*(2) (2) I (2) J ( i j * " iErj ( 1 - 6 i j ) - b g X , x , j i j B i j - and the marginal p.d.f. of n Is h(2) (xi,) = ... (2 iJ h j ( ~ ~ ~ ~ f ~ , ) d x s j (c) Pooled Samples: I In cases i n which the data sample, x i , ,pi , , a , , bll, contains farm units which employ crop j under both irrigated and raiqfed - condition..we may po-1 the complete sample to improve the precision (of our I estimbtes, obtained by treating the two subsamples within czses (1) '~ separately. The improved precision occurs due to the common paramete s, 8 ic -Y i n each of the log- likelihoods, (3.8) and (3.14), since the function, (3.1) applies t o each case. In t h i s case the p o l e d log- likelihood is simply Apart from improving the efficiency of the estimation of 8 , pooling -j "salvage" certain parameters i n 8, of case (1) which would have othe P1 i 'I inestimable due t o the f a i l u r e to report wi Thus, provided the Sam le j I contains both irrigated and rainf ed farmers cultivating a given crop, the - i1 reported data for r a i n f a l l i n the l a t t e r may be pooled with data for rrigated . . I farmers t o identify the parameters even though w is not observed. ndeed, Ij a 1 - T C - 11 7 ~ ~ 1 ~ ~ ! ~ ~ : : p ~~f 2r3?-7L2rri ;deri: i : i - 3 > i 1 - 1 : y . '9 . - 4. F e r t i l i z e r Usage and Response under Uncertainty: * e Y 0 The fact that the length of the production period 1 harvest is of significaat duration for most crops creates various typ s of i 1 additional problems involving decision-making under uncertainty. Inp t decisions with respect to f e r t i l i z e r use must often be made w e l l i n a vance of I I - 31 - I 4 the harvest time. As a result, variations i n the pattern of r a i n f a l l nd other weather conditions, vhich occur a s random events i n the interven ng i time-span, create uncertainty with respect to both the level of output and the price a t vhich the crop map be sold, and, thus, influence decisions reg)arding f e z t i l i z e r use. W e shall now i l l u s t r a t e econometric methods to handle (such types of decision-makin~unde~uncertaiatpsituations. (a) Case (3): Rainfed Crops with Uncertain Rainfall 1 I 1 I To i l l u s t r a t e our approach to t h i s c l a s s of problems, we shal I reconsider the problem of deciding upon the optimal levels of f e r t i l i z e (and I other variable inputs) for the case of rainfed cro-ps. In case (2) of sgction 3, we assumed that rainfall, vij was the only source of orater for crops i n the I absence of i r r i g a t i o n and water-management systems. EIowever,- the.major prablem v i t h the use of the model, (3.10a) and (3.10b), and associated 1 coaditional log- likelihood, L*(2) (0 ,Z(2) ) of (3.15), is that the valu 1 . - j 1 the s t a t e variable, wl denoting the r a i n f a l l index over the crop cycl j' cannot be assumed to be imom to the decision-maker at the t i m e of fertfllizer I application, etc. , e Suppose that r a i n f a l l conditions i n the t i m e period under 1 C"v,":~"7: !3T; ,-2 +z_C> _ - -- . - -h i - l " * , - , d f i T ,- p a -" ' J C i- - . J - - 7 : ' ' I i . J 't - trend and a sizeable random white-noise component with location- specific - ' parameters. Also, sup&se that, by utilizillg a sufficiently lengthy serkes of w I prior historical data on nonthly o r even daily r a i n f a l l , i n relation to ihe I I -expected time--span of the crop, j, the probability distribution, ( ), can be estimated for an appropriate rainfall index, vij. ~v,ij(wij;-iJ Here, denotes a vector of parameters associated v i t h the t i m e , t j ( i ) Jij place, P ( i ) of observation i on crop j> Finally, suppose that the lamer I has knowledge of the historical rainfall patterns a t his location and, indeed, I I enploys % ij(wij;&j) F 0 t o guide h i s decisions on zij and sij. 8 In a particular example, agronomic/meteorological knowledge - . used to define a suitable rainfall- index, w il' associated with crop j location L*(i), a s w e l l a s a simple stochastic model f o r the r a i n f a l l * *P bution a t P ( i ) f o r the typical tlmrspan (from planting, - t i ) ( ) t o e hamest, tj(i)) associated with crop j. For present purposes we assum that such a p.d.f., i %,ij ("ijiAij ), has been defined and that its paranet rs can be estimated via MI, methods from a .extraneoixs sample of r a i n f a l l data f o r 1 I each location, L*(i), and crop j. L e t aj denote the vector of a l l s t a t e variables observed a t t h e I t i m e of decision-making, I In case (2) of section 3 we described a neoclassical model f c r muimira(tion of the conditional expectation of profits, given the set of decision the set of state variables, %j, lij' - and the r a i n f a l l index, w ~ ~ for -4.8n (21,d ~ [ n ~ ~W ilj , alj ~1 , ~ir, however, , )* knowledge of w1j and may be denoted by ni (2 s j w i of . $ffect by taking the marginal expectation of x ij= present situation, ~ - 1/ In practice, 15 farm level data are not collected, we would have to rely upon data from neighboring weather reporting stations. I I I I Thus, i n comparing (2)* of (3.11) with n of (4.2). we see that th iJ Q consequence of uncertainty is to replace the latent variablle, wij, by the parameters, I$ i n its p.d.f., and to change the form of the o jective - -ij' % , i d S function to one requiring numerical itegration for its evaluation. 1 . Let x(3)* be the solurioa to the problem, -iJ which, by the implicit function theorem, vill represent the optimal declision- . vector aa a function of the form, (3>* *(3>* zij (4.4) -ij (2i,;2j&&,)* Comparison with x(2)* of (3.10b), the optimal decisfon vector when -1 j known, exhibits again the replacement of. w i Jby the parameter vector, and the difference i n th6 functional form of the demand functions--x 8 versus x(3)* m -J Given x ( ~ ) * ,.the observed decision vector, zij, can now be -i J =. represented by the model, (3)* (3) .(3)* (3) = X "ij -iJ +4j -j ( 4 j s j*Aij) + '2 4~ * ( 4 . 5 ) -(2) v ~ e r eu ( ~ )follovs the p.d .f., of ( 3 . ~ a ) , i.e., -ij gu,j' Thus, the marginal p.d.f. of x is simply *j +- I and comparison with (3.U) reveals that uncertainty with respect to chamsea the mean of x from 3(2)* to x(3)*, while preseming the ran -ij -il covariance matrix, z!~)*, as i n (3.2a) of case (2). Consider next the conditional p.d..f. of profits, zij, given that is determined through the model (4.5 j in ignorance, ex-ante, of the ex 3 j I - post value of the rainfall index, wij, which subsequently obtains as 4 randon . I In the case where the resulting w value is i j ..- .. ,at harvest- tine), then, conditional upon both x- -and the p.d.f. of case (2), ~ vhere g(2) is defined in (3.2 b); n(2)* (w ) is defined by. (3.11) , and is a cI:, l i j ij function of the actual value, wij; and x(3)* is defined by (4.4), and depends 4 -if I conditional joint p.d.f. of kjand E ~ ~ , is whereas the unconditional p.d.f. is defined by ~ * s alo aa 5 5 d Q-4 Ur, r'l n W WWW0 Mdd ea4 3 guu M cd "4ci a cd d cd td #4 cd a td PI 5 d cd I degenerate point-mass distribution a t w13 ' whence n(3)* of (4.2) c o ~ v e r ~ e s i J 1 The algorithm of. section 2(c) applies, mutatis mutandis, do this 1 I I ,' case. The major difference is the use of nunarical integration to evaluate I I I the objective function, (3)* of (4.2j, associated with each Inciivid/ual, which , I I is then maximized with respect to x i n (I) of Step (n.1.r). I 4 3 I It should also be evident that this approach readily gene the treament of other sources of uncertain? with respect to any o variables la s An obvious extedsiorr would be to treat prices of 4 3 pij, a t harvest t i m e as uncertain a t the Lime of planting and devel to generate the p.d.f. of pij based on historical data on prices and other . 1/ I _determinants up to the t i m e of planting- (b) Risk and Uncertaintp: ~* We now generalize our discussion to the case where, in the presence I of uncertainty, farmers are not assumed to be neutral to risk. - - We dpfine r i s k ~ I as the conditional variance of profits, given the decision variables1,x and +IS the vector of s t a t e variables, s of (&.I), where wlj is now an el6 ent 4 3 I A of z Suppose, upon reordering, that ge partition s as -ij- -1 j * where s ( l ) is h o r n with eerfalnty a t the t i m e 05 decision making, but -il uncertainty surrounds the elements of a(2 . (2 ...--- Let g (2) ( E ;~l j )~d e k t e the -ij 9 -9 j i -11 Such uncertainties can be eliminated i f farmers contract t (I) - k t time, E i s for delivery of their crop' a t time, t j ( i ) , a t the currently pJevailin$ futures market prices. ~ I I . i 4 11' g u I-ln f;l s u m -a 4 4 3 I LC A KT cd 0 ku cd * Q Q)4 d 4a w al LCLC- r, 14 -m 3 ili Ill w s Y B w rrj PO 5 LI 0 rw 4'), Thus x(4)* depends on the state variables known wi:h certainty, s the -ij 1 parameters, 8 i n the production function; the parameters, P , defining th -j ' covariance structure of &e errors; the parameters, l,,a s w e l l as other J possible determinants, Ln the distribution of the s t a t e which uncertainty is attached; and, f i n a l l y , the parameters,%, i n the u t i l i w function, defining the trade-off between the expectation and r i s k of p r o f i t s , Hence, i n the case where s ( ~ )is l a t e n t we ha-~bthe mod&* -il I I I = x(4)* + -= X U.,=., 4 s ! . r (4 2Obir l i j -ij Eij *~ 1 -A, -j j j j ' l j 9 + where bij)= N(o,q) (4.21a) j ~ u , j ( ~iui ij - - ) = N ( d t l - % ' jb - -j d ~*-l ) , j --j a (4.21b) j +i -j -A". Thus, the joint p.d.f. of r and x is given by ij -ij with log- likelihood function, I f , on the other hand, s ( ~ )is obaerved -ij (4.20a) is replaced by I I p . / d . f . where fi is defined by (2.4) ~ i i hthe notation change, and the is j i to be inserted into (4.23). I Thus, apart from these changes, the algorithm cf section 2 uiay be I - usad to estimate the parameter vector, a I vecZ2' 02'], where ljis y a s s d -1 [ti-J n to have been estimated from aa extrarieous sample--providing a gendral a p p r d for the econometrPc treatment of risk i n the presence of uncettai ty. 5. Conclusions: I This paper has outlined a general approach to the agricultural productgon, input demand and output sugply i n the coltext of-- and uncerrainty. As a by-product, we have.also indicated.(section 2e) how tge , 4, econometrics of the dual concepts of profit and cost functions ma be treatd in the context of a stochastic model. Our purpose has been to utfllize the issues of f e r t i l i z e r use, irrigation systems and crop response as a vehicle m I illustrate these ideas. It should, however, be clear that the approach is . I rn applicable to a vast Array of other agricultural (and non-agriculkral) I problems. I Implenentation of these methods on the Bangladesh [I981 data has. ]i I a t present, to await clearance by the Internationgl Pertilizer ~evklopment Ic I Center and the Government of Bangladesh. It is hgpedBhowever. t h ~ in a t I subsequent paper we w i l l be able to report on such erpirlcal resul References Askari, Hossein and J.T. Cunmings 119771, Agricultural SupplyiResponse: A Survey of the Econometric Evidence, Praeger, New Yor*. I I Bangladesh F e r t i l i z e r Study [1981], Agricultural Production, $ e r t i l i z e r E f s e and 4 1979f80, Bangladesh, Joint Report of the Bangladesh i)griculturalih Research Council and the International F e r t i l i z e r Deqelopment Cinter, I preliminary Draft (December, 1981). I I I Davidson, W.C. [1968], "Variance Algorithm for Minimization," Computer Journal, 10, 406-410. II Fletcher, R. and M.J.D. Powell [1963], "A Rapidly convergent descent ~NetacPa for Minfaizaticn ,"Computer Journal, 6, 163-68. I Haber, S. [1970], "Numerical Integration of 24ultiple Integrals," SUIM M e w , 12. Hartley, Michael, J. [1981a], "Neoclassical Econometrics, part One: GeParrl Considerations," World Bank mineo, December, 1981, lj39. Rartley, Michael J [1981b], "Neoclassical Econometrics, Part ~ . I o MottwatZma : and Synthesis," World Benk mimeo, December, 1981, 1-91. Lau, Larry [1978], "Application of R o f i t Functions," i n M. Fubs and D. McFadden, eds., Production Economics: A Dual Apprbach to T k a q ~ and Applications, Vol. 1, Amsterdam: North-Holland. Lucas, Robert E., Jr. a d Thomas J. Sargent [I9811 , Ecoaometric Practice, Vols. 1 and 2, Minnesota Press. i I Hesadden, Daniel, M. [1978], "Estimation Techniques for the labt tic it^ of Substitution and Other Production Parametefs." i n M. Fuss and D. McFadden (eds.) ,- Production Economics: dual ~ p ~ k o a cth W r y o and A ~ ~ l i c a t i o n Vs o l . , 1 , A m s t e r d a m : N o r t k H o l l a n d . I 9 ~ b l o v e ,Marc [I9581 , - Responst! to Price, Baltimore: Johns Hoptins University Press. L ~ Q l o v e , Marc, David Grether and Jose Luis Carvalho [1979], Economic Time Series, New York: Academic Ress. Powell, M.J.D. [1964], "An Efficient Method for Finding the Function of Several Variables without ~ a l c u l a t i -n g Computer Journal, 7, 155-162. Sidhu, S.S. and C.A. Bzanante [1979], "Farm-level Mexican Wheat Varieties i n the Indian Agricultural Economics (August), 61(3), 455-62. Theil, Henri [I9611, Economic Forecasts and Policy, 2nd. $d., Amsterdam: - North-Holland. I I Zellner h o l d , Jan Ruenta and Jacques Drgza [1966], '~pecificationand , Estimation of Cobb-Douglas Production Function Hqdels," Zconornetrica (October), 34, 784-795. i I