这是我的第一篇帖子,大家好!我希望它不会让人感到高兴。
对于我的一些数据,我有计算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)