Я прочитал несколько вопросов по этой теме, но не смог найти ответ.
У меня 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 *
Сначала , я попытался построить кривую, используя 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 Но в этом случае я не рассматривал реальную модель, а моделировал переменные отдельно, что, вероятно, является неправильным способом показать данные.
Секунда , я пытался использовать visreg, но он не работает с "poly(x,3)" в модели: «Ошибка в poly (x, 3): ' Степень «должно быть меньше, чем количество уникальных точек». Если я использую "I(x^3)", это работает, но немного меняет модель.
Третий , я пытаюсь использовать plogis, но я не знаю, как мне вставить оценки параметров в функцию, а затем в график ggplot.