Bookmark and Share

Notice: On April 23, 2014, Statalist moved from an email list to a forum, based at statalist.org.


[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

RE: st: How is confidence interval calculated in weighted GLM in Stata?


From   "Tan, Ji (GRS, Wuhan)" <[email protected]>
To   "[email protected]" <[email protected]>
Subject   RE: st: How is confidence interval calculated in weighted GLM in Stata?
Date   Sat, 1 May 2010 09:46:16 +0800

Sorry I did not know the convention here, and Thanks for the reference, while I don't have immediate access to. Let me rewrite in plain text:

Say in a simple Poisson regression,  E(Yi)=ui=Oi*exp(b), Log likelihood function is {sum (Yi*b-ui)}, when we add weighting wi, it became {sum wi(Yi*b-ui)}, The estimation of b is {(sum wi*Yi)/(sum wi*Oi)}, and the Fisher Score algorithm gives its distribution as N(b, (sum wi^2*ui)/(sum wi*ui)^2). 

However, this calculation only agree with Stata when we have equal weights, i.e. all the wi are the same. So would anyone kindly explain this to me?

Best regards,
Ji

-----Original Message-----
From: [email protected] [mailto:[email protected]] On Behalf Of Maarten buis
Sent: Friday, April 30, 2010 8:02 PM
To: [email protected]
Subject: Re: st: How is confidence interval calculated in weighted GLM in Stata?

--- On Fri, 30/4/10, Tan, Ji (GRS, Wuhan) wrote:
> Say in a simple Poisson regression,  EYi= i=Oi*exp()
> Log likelihood function is  So  And Fisher Score
> algorithm gives: This is in line with Stata.
> 
> However,  in the weighted case, log likelihood
> function changes to , ,
> I  find Stata does not follow this numerically. So
> would anyone kindly explain this to me?

You can only sent plain text and no attachments to
the statalist, so your message got hopelessly 
garbled. So I can't answer you question directly,
I can only give you some literature recomendations:

James W. Hardin and Joseph M. Hilbe "Generalized 
Linear Models and Extensions, second edition. 
College Station, TX: Stata Press.
<http://www.stata-press.com/books/glmext.html>

Hope this helps,
Maarten

--------------------------
Maarten L. Buis
Institut fuer Soziologie
Universitaet Tuebingen
Wilhelmstrasse 36
72074 Tuebingen
Germany

http://www.maartenbuis.nl
--------------------------


      

*
*   For searches and help try:
*   http://www.stata.com/help.cgi?search
*   http://www.stata.com/support/statalist/faq
*   http://www.ats.ucla.edu/stat/stata/
Notice of Confidentiality 
This transmission contains information that may be confidential.  It has been prepared for the sole and exclusive use of the intended recipient and on the basis agreed with that person.  If you are not the intended recipient of the message (or authorized to receive it for the intended recipient), you should notify us immediately; you should delete it from your system and may not disclose its contents to anyone else.

This e-mail has come to you from the Wuhan branch of Watson Wyatt Consultancy (Shanghai) Limited, a Towers Watson company.

*
*   For searches and help try:
*   http://www.stata.com/help.cgi?search
*   http://www.stata.com/support/statalist/faq
*   http://www.ats.ucla.edu/stat/stata/


© Copyright 1996–2018 StataCorp LLC   |   Terms of use   |   Privacy   |   Contact us   |   Site index