Estoy trabajando en lmer
en R. Tengo dudas si hay formas de poner un término quadratico (como factor fijo) forzando al origen esta relación. Actualmente estoy usando (poly(x,degree=2,raw=TRUE))
, pero esta no es forzado a 0:0. Es posible realizar esto?
El modelo es el siguiente:
model<-lmer(Y~a+b+(poly(x,degree=2,raw=TRUE))+(1|factor_aleatorio), data=data)
Y ~ x + I(x^2) - 1