5 messages · Dieter Menne, David Winsemius, jpl
Hi, How do I extract the standard errors of the coefficients in a glm.nb model? Thanks.
View this message in context: http://www.nabble.com/extracting-standard-errors-in-glm.nb-tp22019202p22019202.html Sent from the R help mailing list archive at Nabble.com.
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
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
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.
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.
Mea culpa. Thanks! jpl
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.
______________________________________________ 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-tp22019202p22028365.html Sent from the R help mailing list archive at Nabble.com.