当我意识到我的结果变量CORRECTNESS实际上是正确响应的比例,并且因此将其限制在0和1之间时,我正在拟合线性混合效应模型,我假设它可以取有限数量的可能值在我的数据中。
因此建议我使用具有泊松分布的广义线性模型。当我尝试使用泊松分布重新运行空模型时,我收到两个警告和一条错误消息。
我的数据如下所示(CORRECTNESS是一个数字变量,而ID是一个因数):
CORRECTNESS <- c(0.625, 0.375, 0.5, 1,1,0.875)
ID <- c('p01', 'p02','p03', 'p04', 'p05', 'p06')
guesses <- as.data.frame(cbind(CORRECTNESS,ID))
mdl_null <- glm(CORRECTNESS ~ (1|ID), data = guesses , family = 'poisson')
警告:
Warning message in Ops.factor(1, ID):
“‘|’ not meaningful for factors”Warning message in glm.fit(x = structure(numeric(0), .Dim = c(0L, 2L), .Dimnames = list(:
“no observations informative at iteration 1”
Warning message:
“glm.fit: algorithm did not converge”
错误消息是:
Error in glm.fit(x = structure(numeric(0), .Dim = c(0L, 2L), .Dimnames = list(: object 'fit' not found
或:
Error in model.matrix.default(mt, mf, contrasts): variable 1 has no levels
有人可以帮助我理解吗?如何运行该模型? 谢谢。