La función de simulación de un clic se puede usar para generar una distribución de estadísticas en función de una fórmula de respuesta que debe contener una función aleatoria (como Random Normal). Aquí hay un artículo sobre cómo usar la simulación de un clic para el análisis de potencia en Diseño experimental.
Para los modelos de efectos fijos, la fórmula puede ser tan simple como: Y = 3+ 3 * X + Random Normal (). Si agrego un factor de grupo aleatorio al modelo Y, primero creo una columna de agrupación (grupo, foto a continuación) para especificar la propiedad de la fila a los grupos y otra columna (G) con una fórmula para generar aleatoriamente los efectos de cada grupo que luego serán estimado como BLUP, y la fórmula cambia a, digamos: Y = 3 + 3 * X + G + Random Normal (); donde G se refiere a los efectos de grupo; consulte el archivo de tabla JMP adjunto para obtener más detalles.
Sin embargo, si la fórmula Y se especifica como anteriormente, la simulación de BLUP, a través de la función de simulación, como en la imagen a continuación, no parece evaluar G en cada simulación porque, supongo, la generación aleatoria de efectos grupales no aparece directamente en la fórmula Y; Los efectos de grupo aleatorio son generados por la fórmula de la columna G, no por la fórmula de la columna Y.
De esta manera, si trato de simular con Y, la única reevaluación que se realiza en cada simulación es la del término de error aleatorio que se encuentra en la fórmula Y, es decir, obtengo una simulación de modelo de efectos fijos. ¿Cómo podría especificar una fórmula para Y que genere directamente los efectos del grupo? Supongo que se necesitan algunos bucles, pero no sé cómo hacerlo.
Creo que tiene razón en que el problema es que la columna G no se vuelve a evaluar para cada simulación ... Esto podría remediarse incluyendo la fórmula completa en una columna, tal vez así:
3 +: X * 3 + Si ( Fila () == 1 El | : grupo! = Retraso (: grupo), r = Normal aleatorio (), r) + Normal aleatorio ()
Aquí la variable "r" se redefine para cada grupo. Esto supone que la tabla de datos está ordenada por "grupo". La primera normal aleatoria en la fórmula corresponde al efecto de grupo, y la segunda es el error residual.
Espero que esto ayude.
Adán
