[R] extracting p value from GEE

Tu Yu-Kang yukangtu at hotmail.com
Mon Dec 8 13:48:57 CET 2003


Hi,

Many thanks for your kind help.

best regards,

Yu-Kang

>From: Emmanuel Paradis <paradis at isem.univ-montp2.fr>
>To: "Tu Yu-Kang" <yukangtu at hotmail.com>
>CC: r-help at stat.math.ethz.ch
>Subject: Re: [R] extracting p value from GEE
>Date: Fri, 05 Dec 2003 15:48:41 +0100
>
>At 11:53 04/12/2003 +0000, vous avez écrit:
>>Dear R users,
>>
>>If anyone can tell me how to extract the p values from the output 
>>of gee?
>
>They are easily computed from the output of summary(gee(...)) which 
>prints either a "z" or a "t" depending in the "family" option. z 
>follows, under the null hypothesis, a normal distribution N(0, 1), 
>you have the corresponding P-value with (for a two-tailed test):
>
>2 * (1 - pnorm(abs(z)))
>
>t follows a 'Student' distribution with df degrees of freedom given 
>by N- k - 1, where N is the number of observations, and k is the 
>number of estimated paramaters. I think, but am not definitely sure, 
>that N is counted among all clusters, and k is the number of 
>parameters in the GLM eventually included the estimated scale 
>(correlation parameters are not counted). As above, you have the 
>P-value with:
>
>2 * (1 - pnorm(abs(t), df))
>
>HTH
>
>Emmanuel Paradis
>
>
>>Many thanks in advance.
>>
>>Yu-Kang
>>
>>_________________________________________________________________
>>§K¶O¸ÕÅ¥ MSN ­^»y¾Ç²ß¡G©M¯u¤H¦Ñ®v½u¤W¾Ç­^¤å 

>>
>>______________________________________________
>>R-help at stat.math.ethz.ch mailing list
>>https://www.stat.math.ethz.ch/mailman/listinfo/r-help
>

_________________________________________________________________
¨Ó³}³} MSN eShop¡G«K§QªºÁʪ«Àô¹Ò¡BÀu½èªººô¸ô°Ó®a¡AÀH®É¦³Åå³ßªºÀu´f«P¾Pµ¥µÛ±z
¡I http://msn.com.tw/eshop




More information about the R-help mailing list