Content deleted Content added
m fix red link |
move section on Poisson from Quasi-maximum likelihood estimate |
||
Line 18:
By the usual argument for [[M-estimator]]s (details in Wooldridge <ref name= "Woolridge" />), the asymptotic variance of {{radic|''N''}} ''(θ<sub>MLE</sub>- θ<sub>0</sub>) is A<sup>−1</sup> BA<sup>−1</sup>'' where ''A<sup>−1</sup> = E[ ∑<sub>t</sub>∇<sup>2</sup><sub>θ</sub> logf<sub>t</sub> (y<sub>it</sub>│x<sub>it</sub> ; θ)]<sup>−1</sup> and B=E[( ∑<sub>t</sub>∇<sub>θ</sub> logf<sub>t</sub> (y<sub>it</sub>│x<sub>it</sub> ; θ) ) ( ∑<sub>t</sub>∇<sub>θ</sub> logf<sub>t</sub> (y<sub>it</sub>│x<sub>it</sub>; θ ) )<sup>T</sup>]''. If the joint conditional density of y<sub>i</sub> given x<sub>i</sub> is correctly specified, the above formula for asymptotic variance simplifies because information equality says ''B=A''. Yet, except for special circumstances, the [[joint probability distribution#Joint density function or mass function|joint density]] modeled by partial MLE is not correct. Therefore, for valid inference, the above formula for asymptotic variance should be used. For information equality to hold, one sufficient condition is that scores of the densities for each time period are uncorrelated. In dynamically complete models, the condition holds and thus simplified asymptotic variance is valid.<ref name= "Woolridge" />
== Example: pooled QMLE for Poisson models==
Pooled QMLE is a technique that allows estimating parameters when [[panel data]] is available with Poisson outcomes. For instance, one might have information on the number of patents files by a number of different firms over time. Pooled QMLE does not necessarily contain [[unobserved effects]] (which can be either [[random effects]] or [[fixed effects]]), and the estimation method is mainly proposed for these purposes. The computational requirements are less stringent, especially compared to [[fixed-effect Poisson model]]s, but the trade off is the possibly strong assumption of no [[unobserved heterogeneity]]. Pooled refers to pooling the data over the different time periods 'T'', while QMLE refers to the quasi-maximum likelihood technique.
The [[Poisson distribution]] of <math>y_i</math> given <math>x_i</math> is specified as follows:<ref name="CameronTrivedi">Cameron, C. A. and P. K. Trivedi (2015) Count Panel Data, Oxford Handbook of Panel Data, ed. by B. Baltagi, Oxford University Press, pp. 233–256</ref>
: <math>f(y_i \mid x_i ) = \frac{e^{-\mu_i} \mu_i^{y_i}}{y_i!}</math>
the starting point for Poisson pooled QMLE is the conditional mean assumption. Specifically, we assume that for some <math>b_0</math> in a compact parameter space '''B''', the conditional mean is given by<ref name="CameronTrivedi"/>
: <math> \operatorname E[y_t \mid x_t]=m(x_t, b_0) = \mu_t \text{ for } t= 1,\ldots, T. </math>
The compact parameter space condition is imposed to enable the use of [[M-estimator|M-estimation techniques]], while the conditional mean reflects the fact that the population mean of a Poisson process is the parameter of interest. In this particular case, the parameter governing the Poisson process is allowed to vary with respect to the vector <math>x_{t}\centerdot</math>.<ref name="CameronTrivedi"/> The function ''m'' can, in principle, change over time even though it is often specified as static over time.<ref name="Woolbridge2002">Wooldridge, J. (2002): Econometric Analysis of Cross Section and Panel Data, MIT Press, Cambridge, Mass.</ref> Note that only the conditional mean function is specified, and we will get consistent estimates of <math>b_{0}</math> as long as this mean condition is correctly specified. This leads to the following first order condition, which represents the quasi-log likelihood for the pooled Poisson estimation:<ref name="CameronTrivedi"/><!--not sure if the correct ref is this one or the previous one-->
: <math>\ell_i(b)=\sum[y_{it} \log(m(x_{it},b))-m(x_{it},b)]=0</math>
A popular choice is <math>m=(x_t,b_0)=\exp(x_t b_0)</math>, as Poisson processes are defined over the positive real line.<ref name="Woolbridge2002"/> This reduces the conditional moment to an exponential index function, where <math>x_t b_0</math> is the linear index and exp is the link function.<ref>McCullagh, P. and J. A. Nelder (1989): Generalized Linear Models, CRC Monographs on Statistics and Applied Probability (Book 37), 2nd Edition, Chapman and Hall, London.</ref>
==References==
| |||