如何選擇效果圖的唯一面板我想使用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第一個和效果顯示的第三個面板，因爲第二個和第四個面板是相關的。並建議如何做到這一點？謝謝。

來源

2017-06-22 Kenny Chen

一種方法是侵入生成的效果對象並修改變量的級別。

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")

來源

2017-06-22 15:48:38

我想你的建議，但似乎情節不能被修改。我猜這是因爲我正在用兩種教育（高中和大學）假人來繪製Pr（投票）與年齡的關係。我仍然試圖找出一種方法來刪除（hsdeg = 0; coldeg = 1）和（hsdeg = 0; coldeg = 0）的面板。任何進一步的建議表示讚賞。 –

您可以共享您的模型或底層數據嗎？這樣，我們可以使用您的示例而不是樣本數據進行重現。如果你願意分享你的模型，最簡單的方法是在你的R控制檯中運行'dput（logit1）'並將其粘貼爲問題的一部分。 –

我剛剛編輯了我的原始答案，以包含數據鏈接。再次感謝！ –

如何選擇效果圖的唯一面板我想使用R中的「效果」包？

回答

相關問題