windows Vista R 2.10.1 Is it possible to get p values from gee? Summary(geemodel) does not appear to produce p values.:> 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}}
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 -- View this message in context: http://r.789695.n4.nabble.com/gee-p-values-tp2533835p2534302.html Sent from the R help mailing list archive at Nabble.com.
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
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)>>>
"Peng, C" <cpeng.usm at gmail.com> 9/10/2010 8:06 AM
>>>
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
--
View this message in context:
http://r.789695.n4.nabble.com/gee-p-values-tp2533835p2534302.html
Sent from the R help mailing list archive at Nabble.com.
______________________________________________
R-help at r-project.org 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.
Confidentiality Statement:
This email message, including any attachments, is for th...{{dropped:6}}