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
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
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))
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
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
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,
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
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
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|)
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
10 matches
Mail list logo