計算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 。