Как сделать выражение из модели с фиксированными эффектами?

Question

У меня есть фиксированные эффекты от lmer в объекте fixed.effects3, и я хочу создать выражение, которое можно использовать в качестве формулы для построения графика с помощью кривой (). Я попробовал это:

math.growth <- expression(eval(fixed.effects3[1]) + x*eval(fixed.effects3[2]) + x^2*eval(fixed.effects3[3]))

но затем, когда я пытаюсь построить его с кривой (), я получаю следующие ошибки:

curve(math.growth, from=-2, to=3)
Error in eval(expr, envir, enclos) : 
could not find function "math.growth"

Может кто-нибудь сказать мне, что я делаю не так?

Answer 1

Вот простой пример использования включенного набора данных, расширение его на ваш случай должно быть очень простым;если вы используете глобальное присваивание, как я делал ниже, убедитесь, что у вас есть скрытое имя переменной:

fe<<-fixef(lmer(Reaction ~ Days + (1|Subject) + (0+Days|Subject), sleepstudy))
fef<-function(x) fe[1] + fe[2]*x
curve(fef, from = 1, to = 10)