计算P值

Question

如何使用"t-sample t-test“和"Wilcoxon Mann-Whitney test”计算以下数据的P值：

died<-mydata[mydata$disoutcome=="died"]
survived<-mydata[mydata$disoutcome=="recovered"]
mean_died<-sapply(died,mean,na.rm=T)
SD_died<-apply(died,2,sd,na.rm=T)
mean_survived<-sapply(survived,mean,na.rm=T)
SD_survived<-apply(survived,2,sd,na.rm=T)

谢谢。

Answer 1

使用t.test和wilcox.test函数

Answer 2

这些函数已经可以处理数据帧中的数据。在您的示例中，看起来disoutcome是一个具有级别died和recovered的因子。如果您有另一个包含要比较的度量值的列measurement，则只需使用

t.test(data=mydata, measurement~disoutcome)

或

wilcox.test(data=mydata, measurement~disoutcome)

其中的每一个都将返回结果摘要，包括P值。或者，如果您不关心摘要，您可以只使用$pvalue提取P值。有关更多信息，请咨询?t.test或?wilcox.test。