[R-sig-ME] gee p values

John Sorkin jsorkin at grecc.umaryland.edu
Fri Sep 10 23:22:54 CEST 2010


I am sending this Email to r-sig-mixed-models at r-project.org having previously sent it to r-help and it was suggested that I follow up on this listserver. 
windows Vista
R 2.10.1

Is it possible to get p values from gee? Summary(geemodel) does not appear to produce p values (see code at end of this Email message)

I received the following response from Peng:
There are two z-scores reported in the summary: Naive z and Robust z.

pvalue=2*min(pnorm(z-score), 1-pnorm(z-score))   # two-sided test

I replied to Peng as follows:
> Peng,
> If the answer were as simple as you suggest, I would expect that gee would
automatically produce the p
> values. Since gee does not produce the values, I fear that the computation may
be more complex, or perhaps
> computing p values from gee may be controversial. Do you know which, if either
of my speculations is true?
> Thank you,
> John

Ben Bolker responded: 
  May be worth following up on r-sig-mixed-models .  My guess would be
that if you're willing to treat your data set as 'large' (e.g. your
guess is that the 'residual degrees of freedom', whatever that may
mean, are > 40 ), then you could go ahead and use the naive translation
from Z-score to p-value; otherwise it probably devolves to the
usual 'effective residual degrees of freedom for complex multilevel/
smoothing models' can of worms.













> fit4<- gee(y~time, id=Subject, data=data.frame(data))
Beginning Cgee S-function, @(#) geeformula.q 4.13 98/01/27
running glm to get initial regression estimate
(Intercept)        time 
  1.1215614   0.8504413 
> summary(fit4)

 GEE:  GENERALIZED LINEAR MODELS FOR DEPENDENT DATA
 gee S-function, version 4.13 modified 98/01/27 (1998) 

Model:
 Link:                      Identity 
 Variance to Mean Relation: Gaussian 
 Correlation Structure:     Independent 

Call:
gee(formula = y ~ time, id = Subject, data = data.frame(data))

Summary of Residuals:
          Min            1Q        Median            3Q           Max 
-2.5224390768 -1.4384989365 -0.0006304408  1.4385426203  2.5229173416 


Coefficients:
             Estimate Naive S.E.  Naive z Robust S.E.  Robust z
(Intercept) 1.1215614  0.8023886 1.397778  0.31918831  3.513792
time        0.8504413  0.0993967 8.556031  0.03851821 22.078938

Estimated Scale Parameter:  2.642821
Number of Iterations:  1

Working Correlation
     [,1] [,2] [,3]
[1,]    1    0    0
[2,]    0    1    0
[3,]    0    0    1



John David Sorkin M.D., Ph.D.
Chief, Biostatistics and Informatics
University of Maryland School of Medicine Division of Gerontology
Baltimore VA Medical Center
10 North Greene Street
GRECC (BT/18/GR)
Baltimore, MD 21201-1524
(Phone) 410-605-7119
(Fax) 410-605-7913 (Please call phone number above prior to faxing)

Confidentiality Statement:
This email message, including any attachments, is for th...{{dropped:6}}




More information about the R-sig-mixed-models mailing list