Estoy construyendo un modelo lineal mixto generalizado en SAS. Desafortunadamente, esta es la primera vez que trabajo con GLMM y no estoy muy familiarizado con SAS. He estado usando proc GLIMMIX. Tengo datos de 4 sitios, cada uno especificado con un código de 4 letras en la variable sitio. Cada sitio fue muestreado 11 veces. Se utilizó un diseño cruzado para el tratamiento, por lo que 2 sitios fueron tratados durante las primeras 5 semanas del estudio, y los otros 2 sitios fueron tratados durante las últimas 6 semanas del estudio. He incorporado esto en el modelo con una variable llamada order con un valor 1 para el tratamiento 1º y 2 para el tratamiento 2º, y una variable llamada tratamiento que lee sí cuando el tratamiento estaba presente y no cuando el tratamiento estaba ausente. Mi variable de respuesta debería tener una distribución normal. Mi código actual es:
proc glimmix data=data;
class Tratamiento Sitio Orden;
modelo Respuesta = Tratamiento Orden Sitio;
aleatorio Sitio;
aleatorio Orden;
run; ¿Es esto demasiado simplista? Mi objetivo es averiguar qué efecto tuvieron el tratamiento, el orden y el sitio en la variable de respuesta y determinar qué variables son necesarias en el modelo. Sé que tendré que ejecutar el modelo varias veces con diferentes combinaciones de variables y comparar los valores de AIC para determinar qué variables son necesarias (¿a menos que haya una forma de hacer esto en SAS?). ¿Necesito cambiar alguno de los otros parámetros predeterminados? ¿Necesito hacer de Sitio y/o Orden una variable residual aleatoria en su lugar? ¡Gracias! Se agradece cualquier consejo.
