[R] extracting standard errors in glm.nb

[David Winsemius]

Look more thoroughly?

You should see this about 2/3 the way down the output:

$ coefficients  : num [1:14, 1:4] 3.019 -0.475 -0.709 -0.724 -0.615 ...
   ..- attr(*, "dimnames")=List of 2
   .. ..$ : chr [1:14] "(Intercept)" "SexM" "SexF:AgeF1"  
"SexM:AgeF1" ...
   .. ..$ : chr [1:4] "Estimate" "Std. Error" "z value" "Pr(>|z|)"

You can also get the same information by entering:

coef(sumnb1)

-- 
David Winsemius

On Feb 15, 2009, at 9:40 AM, jpl wrote:

>
> Hi Dieter,
>
> I don't see anything in the str(model) output that matches the  
> standard
> errors of the coefficients given in summary(model).  Any other  
> suggestion?
> Thanks.
>
> jpl
>
>
>
> Dieter Menne wrote:
>>
>> jpl <s-jlubben1 <at> math.unl.edu> writes:
>>
>>>
>>> How do I extract the standard errors of the coefficients in a glm.nb
>>> model?
>>
>> library(MASS)
>> quine.nb1 <- glm.nb(Days ~ Sex/(Age + Eth*Lrn), data = quine)
>> sumnb1 = summary(quine.nb1)
>> sumnb1
>> str(sumnb1) # tells you about the anatomy. Check for $coefficients
>>
>> Dieter
>>
>> ______________________________________________
>> 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.
>>
>>
>
> -- 
> View this message in context: http://www.nabble.com/extracting-standard-errors-in-glm.nb-tp22019202p22023318.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.