如何从模型固定效果中表达？

Question

我在一个对象中有lmer的固定效果，fixed.effects3，我想创建一个表达式，我可以将其用作使用curve（）绘图的公式。我试过这个：

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

但是当我尝试用curve（）绘制它时，我得到了这些错误：

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)