Displaying 2 results from an estimated 2 matches for "beta_0i".
Did you mean:
beta_0
2010 Jul 21
0
Piecewise regression using lme()
....frame(y,t=1:65,group=as.factor(group))
the model I'd like to fit is:
y_t=
(beta_01+beta_11*t+error_1)*(group==1)+(beta_02+beta_12*t+error_2)*(group==2)+(beta_03+beta_13*t+error_3)*(group==3)
It looks like a mixed effects model were the fixed effect are the piecewise
linear regression parts (beta_0i+beta_1i*t) and the random effects are the
variance components error_i.
The problem is that I can't find the way the write this model correctly
using the lme() function of the package nlme. I can't remove the global
intercept and the global variance component. Here's what I tried:
#1...
2013 Feb 04
10
duda con lmer. Añadir predictor a nivel de grupos
Hola a todos.
Estoy utilizando la función lmer del paquete lme4 para ajustar un modelo
mixto.
Tengo varias variables en mi data.frame, unas son a nivel individual y
otras a nivel de comarcas. Listo algunas.
ingre_6 : Ingresos (nivel individual)
iscs_a : un indicador sintético resumen de otras variables, calculado
mediante componentes principales.
sau_com : superficie agraria útil de cada una