计算P值

Question

对于以下数据，如何使用“ t样本t检验”和“ Wilcoxon Mann-Whitney检验”来计算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 。