Hola a tod@s, tengo una duda que la comunidad R me puede ayudar. Estoy trabajando con gamlss, porque tengo una variable respuesta con valores entre 0 y 1 e incluidos estos. La distribución que utilizo com gamlss para este caso es "beta inflated" (Stasinopulos and Rigby 2007. Journal of Statistical Software 23(7)). El modelo que intento correr es:
m1<-gamlss(Teleosteos ~ cs(LT_max)+random(Familia/Genero,lambda=T),family=BEINF, data = datos) donde LT_max es la talla máxima de especies de peces y la variable respuesta es la proporción del consumo de peces teleósteos de cada especie. El efecto aleatorio es anidado de géneros y familias taxonómicas de especies (género dentro de familia). El problema que tengo es que los gamlss no toman el efecto aleatorio como random(Familia/Genero) o (1|Familia/Genero), sino como: m2<-gamlss(Teleosteos ~ cs(LT_max)+ random(Familia,lambda=T)+random(Familia:Genero,lambda=T),family=BEINF, data = datos) El modelo con el código en amarillo es ajustado con gamlss sin problemas. ¿Es correcto la forma randon(A)+random(A:B) para considerar un efecto aleatorio anidado como (1|A/B)? Muchas gracias de antemano y saludos cordiales a toda la comunidad, Santiago. -- Santiago A. Barbini *Laboratorio de Ictiología* Instituto de Investigaciones Marinas y Costeras (IIMyC) Consejo Nacional de Investigaciones Científicas y Técnicas (CONICET) Universidad Nacional de Mar del Plata (UNMdP) Funes 3350 - B7602AYL - Mar del Plata ARGENTINA http://www.iimyc.gob.ar/ -- [[alternative HTML version deleted]] _______________________________________________ R-help-es mailing list [email protected] https://stat.ethz.ch/mailman/listinfo/r-help-es
