Wie kann man einen Ausdruck aus modellfixierten Effekten bilden?

Question

Ich habe die fixen Effekte von lmer in einem Objekt, fixed.effects3, und ich möchte einen Ausdruck erstellen, den ich als Formel für das Plotten mit curve() verwenden kann. Ich habe dies versucht:

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

aber wenn ich dann versuche, es mit curve() darzustellen, erhalte ich diese Fehler:

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

Kann mir jemand sagen, was ich falsch mache?

Answer 1

Hier ist ein einfaches Beispiel, bei dem ein enthaltener Datensatz verwendet wird. Die Erweiterung auf Ihren Fall sollte sehr einfach sein; wenn Sie eine globale Zuweisung verwenden, wie ich es unten getan habe, achten Sie darauf, dass Sie einen obskuren Variablennamen haben:

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)