Re: [R] confint for glm (general linear model)

2011-03-08 Thread cossincos
it use normal dist. instead of T dist for Z statistic, which usually corresponds to asymptotic normal dist. > -0.283+c(-1,1)*qnorm(0.975)* 0.121 [1] -0.5201556 -0.0458444 -- View this message in context: http://r.789695.n4.nabble.com/confint-for-glm-general-linear-model-tp954071p3341527.html S

Re: [R] confint for glm (general linear model)

2009-12-13 Thread Peter Dalgaard
David Winsemius wrote: On Dec 12, 2009, at 8:19 PM, casperyc wrote: for an example, counts <- c(18,17,15,20,10,20,25,13,12) outcome <- gl(3,1,9); treatment <- gl(3,3) glm.D93 <- glm(counts ~ outcome + treatment, family=poisson()) confint(glm.D93) confint.default(glm.D93) # based on asympto

Re: [R] confint for glm (general linear model)

2009-12-12 Thread David Winsemius
On Dec 12, 2009, at 8:19 PM, casperyc wrote: for an example, counts <- c(18,17,15,20,10,20,25,13,12) outcome <- gl(3,1,9); treatment <- gl(3,3) glm.D93 <- glm(counts ~ outcome + treatment, family=poisson()) confint(glm.D93) confint.default(glm.D93) # based on asymptotic normality to verify

Re: [R] confint for glm (general linear model)

2009-12-12 Thread casperyc
for an example, counts <- c(18,17,15,20,10,20,25,13,12) outcome <- gl(3,1,9); treatment <- gl(3,3) glm.D93 <- glm(counts ~ outcome + treatment, family=poisson()) confint(glm.D93) confint.default(glm.D93) # based on asymptotic normality to verify the confidence interval (confint.default(glm.D93

Re: [R] confint for glm (general linear model)

2009-12-11 Thread Peter Ehlers
I suspect that you don't know about 'profile' confidence intervals. If that's true then I can recommend the discussion in MASS (the book) in section 8.4. In a nutshell, I don't think that you want to do a profile confint calculation manually (unless typing instructions that use the function profil

Re: [R] confint for glm (general linear model)

2009-12-09 Thread brestat
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. Errorz valuePr(>|z|) >

Re: [R] confint for glm (general linear model)

2009-12-09 Thread 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

Re: [R] confint for glm (general linear model)

2009-12-09 Thread casperyc
I think the help page are exactly the same... I just want to verify the confidence interval manually. That's all I want. Thanks. casper brestat wrote: > > This functions are different. I advice you study them: > > ?confint # profile likelihood > ?confint.default # t-distribution > > Walmes

Re: [R] confint for glm (general linear model)

2009-12-09 Thread David Winsemius
On Dec 9, 2009, at 9:21 PM, casperyc wrote: does no one know this? Have you read the Posting Guide? -- View this message in context: http://n4.nabble.com/confint-for-glm-general-linear-model-tp954071p956658.html Sent from the R help mailing list archive at Nabble.com. David Winsemius, M

Re: [R] confint for glm (general linear model)

2009-12-09 Thread casperyc
does no one know this? -- View this message in context: http://n4.nabble.com/confint-for-glm-general-linear-model-tp954071p956658.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mailing list https://stat.ethz.ch/m

[R] confint for glm (general linear model)

2009-12-06 Thread casperyc
Hi, I have a glm gives summary as follows, Estimate Std. Errorz valuePr(>|z|) (Intercept) -2.03693352 1.449574526 -1.405194 0.159963578 A0.01093048 0.006446256 1.695633 0.089955471 N0.41060119 0.224860819 1.826024 0.