我阅读了有关此主题的一些问题,但找不到答案。
我有30个样本,我的模型是:
y ~ poly(x,3) + z
结果是:
Estimate SE Z P
(Intercept) -2.851 0.234 -12.173 < 0.0001
poly(x)1 -0.784 1.036 -0.758 0.449
poly(x)2 1.937 0.845 2.293 0.022 *
poly(x)3 2.754 0.768 3.587 0.0003 **
z 0.342 0.105 3.268 0.001 *
First,我尝试使用geom_smooth
绘制曲线。
a)geom_smooth(method = "betareg", se=FALSE, formula = y ~ poly(x,3))
b)geom_smooth(method = "betareg", se=FALSE, formula = y ~ z)
然后,我得到了这些图形:https://i.stack.imgur.com/PSKNG.png但是,在这种情况下,我没有考虑实际模型,而是分别对变量建模,这可能是显示数据的错误方法。
Second,我尝试使用visreg
,但是在模型中的"poly(x,3)"
上不起作用:“ poly(x,3)中的错误:'度'必须小于唯一点数”。如果使用"I(x^3)"
,它可以工作,但是会稍微改变模型。
第三,我正在尝试使用plogis
,但我不知道如何在函数中然后在ggplot图中插入参数估计值。
我认为gridExtra可以为您提供帮助。因为您自己创建了单独的ggplot对象(每个要获取的图形一个)。这是您可以开始使用的一些示例。