モデル固定効果から式を作成するには？

Question

私はlmerの固定効果をfixed.effects3というオブジェクトに持っています。私はcurve（）でプロットするための式として使用できる式を作りたいと思います。私は、曲線（とそれをプロットしてみたとき、私は、これらのエラーを取得する）、その後

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)