如何在ddply中運行fisher.test？

d <- data.frame(topic=c("a","a","a","b","c","c"), year=c(2001,2002,2004,2003,2011,2012), 
       I=c(3,2,4,3,0,1), II=c(2,1,2,3,4,0), III=c(0,0,1,2,3,0)) 
library(plyr) 
chip <- ddply(df.agg[,-2], "topic", function(x){ 
    round(fisher.test(x[,-1])$p.value, 3) 
    }) 
#Error in fisher.test(x[, -1]) : 'x' must have at least 2 rows and columns

如何在ddply內製作fisher.test？我想要一行的主題（如b）的值爲NA，但其他行報告的值爲p。如何在ddply中運行fisher.test？

來源

2017-04-16 Z. Zhang

是'df.agg [ - 2]'應該是'D'或' d [ - 2]'？ – r2evans

是的，我很抱歉的錯誤。 –

錯誤消息非常多地解釋它。您至少需要兩行和兩列才能工作，但只有一行「topic」是「b」。 – AkselA

您可以在基礎R與split和sapply做到這一點：

sapply(split(d, d$topic), function(i) 

    if (nrow(i) == 1) { 

    NA 

    } else { 

    round(fisher.test(i[,3:5])$p.value, 3) 

})

結果：

a  b  c 
1.000 NA 0.125

來源

2017-04-16 22:31:27 ulfelder

如何在ddply中運行fisher.test？

回答

相關問題