Raoul C Reulen <[email protected]>--
Works fine if you save the scalar you want to add before you run the main model:
sysuse auto
glm mpg rep78, fam(pois)
est sto linear
glm mpg if rep78<., fam(pois)
est sto consonly
lrtest linear consonly
loc ptrend=r(p)
xi: glm pr i.rep78, fam(pois) eform
est store m
estadd scalar ptrend=`ptrend'
estout m, cells("b (fmt(%9.2f)) ci(par( ( , ) )) p(fmt(%9.4f))") eform
stats(ptrend)
On 10/25/07, Raoul C Reulen <[email protected]> wrote:
> I'm running a glm model and like to do a a likelihood ratio test. However, the likelihood ratio test is based
> on a slightly different model. I have difficulty getting the likelihood ratio estimates in a table that I created
> with estout. See output below
*
* For searches and help try:
* http://www.stata.com/support/faqs/res/findit.html
* http://www.stata.com/support/statalist/faq
* http://www.ats.ucla.edu/stat/stata/
