Hi,
Are you sure? That's not what I got:
> require(epicalc)
> ?logistic.display
> model0 <- glm(case ~ induced + spontaneous, family=binomial, data=infert)
> logistic.display(model0)
Logistic regression predicting case
crude OR(95%CI) adj. OR(95%CI) P(Wald's test)
induced (cont. var.) 1.05 (0.74,1.5) 1.52 (1.02,2.27) 0.042
spontaneous (cont. var.) 2.9 (1.97,4.26) 3.31 (2.19,5.01) < 0.001
P(LR-test)
induced (cont. var.) 0.042
spontaneous (cont. var.) < 0.001
Log-likelihood = -139.806
No. of observations = 248
AIC value = 285.612
My impression is that you did something else and you are not telling us the
full story. Here is my sessionInfo():
R version 2.14.0 Patched (2011-11-12 r57642)
Platform: x86_64-apple-darwin9.8.0/x86_64 (64-bit)
locale:
[1] US.UTF-8
attached base packages:
[1] grid splines stats graphics grDevices utils datasets
methods
[9] base
other attached packages:
[1] ggplot2_0.8.9 proto_0.3-9.2 reshape_0.8.4 plyr_1.6
[5] epicalc_2.14.0.0 nnet_7.3-1 MASS_7.3-16 survival_2.36-10
[9] maptools_0.8-10 lattice_0.20-0 foreign_0.8-47 sp_0.9-91
[13] maps_2.2-2
loaded via a namespace (and not attached):
[1] tools_2.14.0
*
*
Could you please share your sessionInfo() as well as your OS, i.e., see
http://www.r-project.org/posting-guide.html?
HTH,
Jorge.-
2011/12/13 å欣 <>
> According to the example of logistic.display:
> model0 <- glm(case ~ induced + spontaneous, family=binomial, data=infert)
> summary(model0)
> logistic.display(model0)
>
> induced: 3levels 0,1,2
> spontaneous: 3levels 0,1,2
>
> So if 0 is reference, we should get 2 OR for
> " induced1"," induced2"," spontaneous1"," spontaneous2"
>
> But the acturally OR is as the following,which is not what I expected:
>
> crude OR(95%CI) adj. OR(95%CI) P(Wald's test)
> P(LR-test)
> induced (cont. var.) 1.05 (0.74,1.5) 1.52 (1.02,2.27) 0.042
> 0.042
> spontaneous (cont. var.) 2.9 (1.97,4.26) 3.31 (2.19,5.01) < 0.001
> < 0.001
>
>
> Can anyone give me some suggestions?
>
> Many thanks!
>
>
>
>
>
>
> Hi sir:
> I follow your suggestion:
>
> result<-glm(y ~ factor(age) + factor(gender) + CD4,family = binomial)
> logistic.display(result)
>
> Error in coeff[, 1] : incorrect number of dimensions
>
>
>
>
>
>
> At 2011-12-14 01:59:36,"Jorge I Velez" <> wrote:
>
> Hi there,
>
> Try
>
> require(epicalc)
> logistic.display(result)
>
> HTH,
> Jorge
>
> On Tue, Dec 13, 2011 at 7:16 AM, å欣 <> wrote:
>
>> Hi all:
>> My data has 3 variables:
>> age(3levels : <30y=1 30-50y=2, >50y=3)
>> gender(Male=0, Female=1)
>> CD4 cell count(raw lab measurement)
>> y(1:death 0:alive)
>>
>> I perform logistic regression to find out the factors that influence y.
>>
>> result<-glm(y ~ factor(age) + factor(gender) + CD4,family = binomial)
>>
>> >From the result,I can get OR(Odds Ratio) of gender via exp(Estimate of
>> Female, since Male is regarded as reference group and has no result).But
>> how can I compute the 95%CI of OR of gender?
>>
>>
>> Thanks a lot for your help.
>>
>> My best!
>>
>>
>>
>>
>> [[alternative HTML version deleted]]
>>
>> ______________________________________________
>> [email protected] 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.
>>
>
>
>
>
>
>
[[alternative HTML version deleted]]
______________________________________________
[email protected] 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.
