怎么调用假设检验算出的pvalue

Gujun

做完假设检验后，想进行下一步直接调用p.value，不知道该怎么调用~~~求大神解答，谢谢！

求调用p值

sillyboy。。。

比如说你的  结果存在result里，即 result <- survdiff(.......你那一大堆)
然后  str(result)  一下，看看P值存的变量叫什么(一般就是pvalue)
然后    result$pvalue   (具体那个变量叫什么，美元符后面就跟什么)

englianhu

1. > ##试一试
   
2. > library('survival')
   
3. > fit <- survdiff(Surv(time,status)~x,data=aml)
   
4. > fit
   
5. Call:
   
6. survdiff(formula = Surv(time, status) ~ x, data = aml)
   
7. 
   
8.                  N Observed Expected (O-E)^2/E (O-E)^2/V
   
9. x=Maintained    11        7    10.69      1.27       3.4
   
10. x=Nonmaintained 12       11     7.31      1.86       3.4
    
11. 
    
12. Chisq= 3.4  on 1 degrees of freedom, p= 0.0653
    
13. > 1 - pchisq(fit$chisq, length(fit$n) - 1)
    
14. [1] 0.06533932
    
15. >
    
16. > fit1 <- survdiff(Surv(tstart,status)~treat+sex,data=cgd)
    
17. > fit1
    
18. Call:
    
19. survdiff(formula = Surv(tstart, status) ~ treat + sex, data = cgd)
    
20. 
    
21.                             N Observed Expected (O-E)^2/E (O-E)^2/V
    
22. treat=placebo, sex=male   100       48     42.1     0.829     2.202
    
23. treat=placebo, sex=female  20        8      6.2     0.520     0.679
    
24. treat=rIFN-g, sex=male     68       17     22.7     1.434     2.389
    
25. treat=rIFN-g, sex=female   15        3      5.0     0.799     0.997
    
26. 
    
27. Chisq= 4.2  on 3 degrees of freedom, p= 0.24
    
28. > 1 - pchisq(fit1$chisq, length(fit1$n) - 1)
    
29. [1] 0.2400191

复制代码

Gujun

sillyboy。。。 发表于 2015-8-10 16:14
比如说你的  结果存在result里，即 result

谢谢！
可能是survdiff函数的原因，str()之后，显示没有p值。

Gujun

englianhu 发表于 2015-8-11 08:09

谢谢，这种方法可行！