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

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 asymptotic

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-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-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

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,

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: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

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|)

[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