如何从模型固定效果中制作表情？

Question

我在对象fixed.effects3中有来自lmer的固定效果，并且我想创建一个表达式，我可以将其用作使用曲线()进行绘图的公式。我试过这个：

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)