La fonction de simulation en un clic peut être utilisée pour générer une distribution de statistiques en fonction d'une formule de réponse qui doit contenir une fonction aléatoire (comme Random Normal). Voici un article sur la façon d'utiliser la simulation en un clic pour l'analyse de puissance dans la conception expérimentale.
Pour les modèles à effets fixes, la formule peut être aussi simple que: Y = 3+ 3 * X + Normal aléatoire (). Si vous ajoutez un facteur de groupe aléatoire au modèle Y, je crée d'abord une colonne de regroupement (groupe, image ci-dessous) pour spécifier la propriété des lignes aux groupes et une autre colonne (G) avec une formule pour générer de manière aléatoire les effets de chaque groupe qui seront plus tard estimé comme BLUPs, et la formule devient, disons: Y = 3 + 3 * X + G + Normal aléatoire (); où G fait référence aux effets de groupe - voir le fichier de table JMP joint pour plus de détails.
Cependant, si la formule Y est spécifiée comme ci-dessus, la simulation des BLUP - via la fonction de simulation, comme dans l'image ci-dessous - ne semble pas évaluer G à chaque simulation car, je présume, la génération aléatoire d'effets de groupe n'apparaît pas directement dans la formule Y; les effets de groupe aléatoires sont générés par la formule de la colonne G et non par la formule de la colonne Y.
De cette façon, si j'essaie de simuler avec Y, la seule réévaluation qui est effectuée à chaque simulation est celle du terme d'erreur aléatoire qui est dans la formule Y, c'est-à-dire que j'obtiens à la place une simulation de modèle à effets fixes. Comment pourrais-je alors spécifier une formule pour Y qui génère directement les effets de groupe? Je suppose qu'une boucle est nécessaire mais je ne sais pas comment le faire.
Je pense que vous avez raison de dire que le problème est que la colonne G n'est pas réévaluée pour chaque simulation ... Cela pourrait être résolu en incluant la formule entière dans une colonne, peut-être comme ceci:
3 +: X * 3 + Si ( Rangée () == 1 | : groupe! = Décalage (: groupe), r = Normal aléatoire (), r) + Normal aléatoire ()
Ici, la variable "r" est redéfinie pour chaque groupe. Cela suppose que la table de données est triée par "groupe". La première normale aléatoire dans la formule correspond à l'effet de groupe, et la seconde est l'erreur résiduelle.
J'espère que cela t'aides.
Adam
