binom.test

时间:2016-11-21 15:09:12

标签: r tidyr

这是我的第一篇帖子,大家好!我希望它不会让人感到高兴。

对于我的一些数据,我有计算Cochran Q-test的问题。 我收到的错误是:

  

binom.test中的错误(长度(符号[符号> 0]),长度(符号[符号!= 0]),:    ' N'必须是正整数> =' x'

此代码有助于查看正在发生的事情(我的数据的一部分):

Data <- c(0, 0, 0, 1, 0, 0, 1, 1, 0, 1, 0, 0, 1, 1, 0, 0, 1, 1, 0, 1, 
      0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
      0, 0, 1, 0, 0, 0, 1, 1, 0, 1)

Data <- as.data.frame(matrix(data, 5, 10))

library(RVAideMemoire)
library(tidyr)
tData <- as.data.frame(t(Data))
colnames(tData) <- c(rep(1:5))
tData$id <- factor(row.names(tData))
tData %>% tidyr::gather(tree, count, c(1:5)) -> tData2
xy <- cochran.qtest(count ~ id | tree, data = tData2, alpha = 0.05)

但如果我将数据分组,那么它就可以了。问题在哪里?

data <- c(0, 0, 0, 1, 0, 0, 1, 1, 0, 1, 0, 0, 1, 1, 0, 0, 1, 1, 0, 1, 
      0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
      0, 0, 1, 0, 0, 0, 1, 1, 0, 1)

data <- as.data.frame(matrix(data, 5, 10))

library(RVAideMemoire)
library(tidyr)
tData <- as.data.frame(t(Data[1:4,]))
colnames(tData) <- c(rep(1:4))
tData$id <- factor(row.names(tData))
tData %>% tidyr::gather(tree, count, c(1:4)) -> tData2
xy <- cochran.qtest(count ~ id | tree, data = tData2, alpha = 0.05)

0 个答案:

没有答案