Your original data must have looked something like the following:
sib.data <- data.frame(sib=rep(c(0,1,0,1), c(159,4,162,37)),
sex=rep(c(0,0,1,1), c(159, 4, 162, 37)))
as that gives the 2x2 table you showed (with 'Response' -> 'sib'):
> table(sib.data)
sex
sib
> -Original Message-
> From: r-help-boun...@r-project.org [mailto:r-help-bounces@r-
> project.org] On Behalf Of Jerome Myers
> Sent: Wednesday, January 18, 2012 7:27 AM
> To: r-help@r-project.org
> Subject: [R] confint function in MASS package for logistic regression
> analysis
>
> I have
Yes, the results from confint() are much more accurate than yours and
SPSS's. (As Bill Venables once said in a similar circumstance: this is
not the place to report bugs in SPSS.)
Hint: the word 'profile' appears all over the place on the help pages.
confint() uses profile likelihood methods.
On Jan 18, 2012, at 9:27 AM, Jerome Myers wrote:
> I have the following binary data set:
> Sex
> Response 0 1
>0 159 162
>1 4 37
> My commands
> library(MASS)
> sib.glm=glm(sib~sex,family=binomial,data=sib.data)
>
4 matches
Mail list logo
