如何选择效果图的唯一面板我想使用R中的“效果”包？

问题描述：

我在R中使用Fox et al的效果包来显示我的logit结果的预测概率。如何选择效果图的唯一面板我想使用R中的“效果”包？

而且他是我做了什么：

logit1 <- glm(vote00 ~ age+ I(age^2) + hsdeg + coldeg, family=binomial(link="logit"), data=data)

我的情节看起来像有以下几点：

plot(Effect(focal.predictors = c("hsdeg","coldeg", "age"), 
     mod=logit1, 
     xlevels=list(age=seq(20,90,1), hsdeg=c(0,1), coldeg=c(0,1))), 
     type = "response", 
     ylab="probability of voting", 
     main="Predicted Probabilities of Age and Education")

Plot to be modified

但是，我想删除f第一个和效果显示的第三个面板，因为第二个和第四个面板是相关的。并建议如何做到这一点？谢谢。

答

一种方法是侵入生成的效果对象并修改变量的级别。

eff <- Effect(focal.predictors = c("hsdeg","coldeg", "age"), 
       mod=logit1, 
       xlevels=list(age=seq(20,90,1), hsdeg=c(0,1), coldeg=c(0,1))) 
eff$variables$hsdeg$is.factor <- TRUE 
eff$variables$hsdeg$levels <- c(1) 
plot(eff, 
    type = "response", 
    ylab="probability of voting", 
    main="Predicted Probabilities of Age and Education")

我想你的建议，但似乎情节不能被修改。我猜这是因为我正在用两种教育（高中和大学）假人来绘制Pr（投票）与年龄的关系。我仍然试图找出一种方法来删除（hsdeg = 0; coldeg = 1）和（hsdeg = 0; coldeg = 0）的面板。任何进一步的建议表示赞赏。 –

您可以共享您的模型或底层数据吗？这样，我们可以使用您的示例而不是样本数据进行重现。如果你愿意分享你的模型，最简单的方法是在你的R控制台中运行'dput（logit1）'并将其粘贴为问题的一部分。 –

我刚刚编辑了我的原始答案，以包含数据链接。再次感谢！ –

如何选择效果图的唯一面板我想使用R中的“效果”包？

相关推荐