[R] confint for glm (general linear model)

[David Winsemius]

On Dec 9, 2009, at 9:50 PM, casperyc wrote:

>
> I think the help page are exactly the same...

I cannot tell what you maen by this.

> I just want to verify the confidence interval manually. That's all I  
> want.

Then provide some reproducible code and data. ... as the Posting Guide  
explains and provides explicit examples of the manner for most  
effectively presenting data objects.

-- 
David

>
> Thanks.
>
> casper
>
>
>
> brestat wrote:
>>
>> This functions are different. I advice you study them:
>>
>> ?confint # profile likelihood
>> ?confint.default # t-distribution
>>
>> Walmes Zeviani - Brazil
>>
>>
>>
>> casperyc wrote:
>>>
>>> Hi,
>>>
>>> I have a glm gives summary as follows,
>>>
>>>               Estimate             Std. Error        z value     
>>> Pr(>|z|)
>>> (Intercept) -2.03693352     1.449574526 -1.405194 0.159963578
>>> A            0.01093048       0.006446256  1.695633 0.089955471
>>> N            0.41060119      0.224860819  1.826024 0.067846690
>>> S           -0.20651005      0.067698863 -3.050421 0.002285206
>>>
>>> then I use confint(k.glm) to obtain a confidnece interval for the
>>> estimates.
>>>
>>>> confint(k.glm,level=0.97)
>>> Waiting for profiling to be done...
>>>                   1.5 %      98.5 %
>>> (Intercept) -5.471345995  0.94716503
>>> A           -0.002340863  0.02631582
>>> N           -0.037028592  0.95590178
>>> S           -0.365570347 -0.06573675
>>>
>>> while reading the help for 'confint', i found something like  
>>> confint.glm
>>> for general linear model.
>>> I load the MASS package by clicking on the Menu( or otherwise how  
>>> should
>>> I load the package?)
>>>
>>> then I still cant use the confint.glm command, what have I dont  
>>> wrong?
>>>
>>>
>>> How do I calculate this confidence interval for glm estimate  
>>> manually??
>>>
>>> for A, I use
>>> 0.01093048 + c(-1,1) * 0.006446256 * qt(0.985,df=77)
>>> which is a different interval i got from the  
>>> confint(k.glm,level=0.97)
>>> above.
>>>
>>> To be short, what's the right command to find the confidence  
>>> interval for
>>> glm estimats?
>>> How do I verify it manully?
>>>
>>> Thanks.
>>>
>>> casper
>>>
>>>
>>>
>>
>>
>
> -- 
> View this message in context: http://n4.nabble.com/confint-for-glm-general-linear-model-tp954071p956671.html
> Sent from the R help mailing list archive at Nabble.com.
>
> ______________________________________________
> R-help at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.

David Winsemius, MD
Heritage Laboratories
West Hartford, CT