Question

我以前用lme4成功地使用了Effects包，但现在我使用glm的一个非常简单的例子。在拟合二项式glm然后使用Effect I绘制交互时，会收到错误消息"Error in Analyze.model(focal.predictors, mod, xlevels, default.levels,:the following predictor is not in the model: Freq*Gp"。下面是一个简单的例子，说明了我的问题。

感谢您的帮助。

Gp<-c(rep("A",20),rep("B",20))
Freq<-c(rep(1,10),rep(2,10),rep(1,10),rep(2,10))
Resp<-sample(c(0,1),40, replace=T)
data<-data.frame(Gp)
data$Freq<-Freq
data$Resp<-Resp

m2<-glm(Resp~Freq*Gp, data,family=binomial)
eff<-Effect("Freq*Gp",m2)

Answer 1

从帮助部分：效果的focal.predictors要求预测变量的字符向量。对于交互，只需提供两种主效应的向量。

eff <- Effect(c('Freq', 'Gp'), m2)

给

Freq*Gp effect

              Gp
Freq          A   B
  1   0.4000000 0.3
  1.2 0.4196106 0.3
  1.4 0.4394784 0.3
  1.6 0.4595421 0.3
  1.8 0.4797378 0.3
  2   0.5000000 0.3

效果包：使用GLM二项模型时出错

1 个答案: