Question

我想绘制Y（二元结果）的预测概率超过观察到的x值范围（x =年龄）。我使用以下代码生成图：

（1）我计算了在我的自变量（年龄）的指定x值（xlevels = x.list）范围内Y的预测概率，并将其保存在一个对象中。

prob   <-  effect(c("age"),M, xlevels = x.list)

（2）然后我绘制该对象，自定义某些绘图外观（如轴标签，置信区间的颜色等）。

plot(prob, 
     xlab="x", 
     ylab="Pred. prob.", 
     confint=list(col="red", alpha=0.3),
     lines=list(col="red")
     rug=FALSE, main="")

这会产生一个几乎看起来像我想要的情节。但是，当尝试自定义主标题，y轴和x轴限制以及轴上的刻度时，绘图会产生，但不幸的是也会混乱（y轴的范围不是0到1，并且具有置信区间的实际线被推出图表边距。

  plot(prob, 
     xlab="x", 
     ylab="Pred. prob.", 
     confint=list(col="red", alpha=0.3),
     lines=list(col="red")
     rug=TRUE,
     axes=list(y=lim={c(0, 1, 0.1)})))

特别是，我想改变y轴以便......

（a）的范围为0-1

（b）其中以0.1为增量的刻度。

（c）我还想摆脱情节周围的方框，即只画出x和y轴。

我一直在尝试阅读？plot.eff，但遗憾的是无法让传统的论据起作用。关于如何修改代码以使其工作的任何输入都将非常感激。

Answer 1

col-md-3看起来很奇怪，你试过axes=list(y=lim={c(0, 1, 0.1)}))

吗？

也可以ylim=c(0,1)

R：用绘图绘制预测概率（效果（）） - 如何自定义绘图外观？

1 个答案: