小编use*_*003的帖子

我重新运行从此前公布的数据Kaplan-Meier生存曲线,使用在出版物中使用的准确数据集(夏邦杰等2008 - 中环尾狐猴(狐猴咔嗒)的近交衰退:遗传多样性预测寄生,免疫活性和幸存者).该出版物使用LIFETEST在SAS Version 9中运行曲线,以分析由遗传杂合性和动物性别构成的死亡年龄(n = 64).她报告的卡方值为6.31,p值为0.012; 但是,当我在R中运行曲线时,我得到Chi平方值0.9和ap值0.821.谁能解释一下这个？

使用的R代码:年龄是死亡时间,mort是审查代码,性别是性别的阶层,而ho2是描述要比较的两个群体的因素.

> survdiff(Surv(age, mort1)~ho2+sex,data=mariekmsurv1)
Call:
survdiff(formula = Surv(age, mort1) ~ ho2 + sex, data = mariekmsurv1)

              N Observed Expected (O-E)^2/E (O-E)^2/V
ho2=1, sex=F 18        3     3.23    0.0166    0.0215
ho2=1, sex=M 12        3     2.35    0.1776    0.2140
ho2=2, sex=F 17        5     3.92    0.3004    0.4189
ho2=2, sex=M 17        4     5.50    0.4088    0.6621

Chisq= 0.9  on 3 degrees of freedom, p= 0.821 


> str(mariekmsurv1)
'data.frame':   64 obs. of  6 variables:
 $ id   : Factor w/ 65 …