La funzione di simulazione con un clic può essere utilizzata per generare una distribuzione di statistiche in funzione di una formula di risposta che deve contenere una funzione casuale (come Random Normal). Ecco un articolo su come utilizzare la simulazione con un clic per l'analisi della potenza in Experimental Design.
Per i modelli a effetti fissi la formula può essere semplice come: Y = 3+ 3 * X + Random Normal (). Se aggiungo un fattore di gruppo casuale al modello Y, creo prima una colonna di raggruppamento (gruppo, immagine sotto) per specificare la proprietà della riga per i gruppi e un'altra colonna (G) con una formula per generare casualmente gli effetti di ciascun gruppo che sarà successivamente stimato come BLUP e la formula cambia in, diciamo: Y = 3 + 3 * X + G + casuale normale (); dove G si riferisce agli effetti di gruppo - vedere il file della tabella JMP allegato per maggiori dettagli.
Tuttavia, se la formula Y è specificata come sopra, la simulazione di BLUP - tramite la funzione di simulazione, come nella figura sotto - non sembra valutare G ad ogni simulazione perché, presumo, la generazione casuale di effetti di gruppo non appare direttamente nella formula Y; gli effetti di gruppo casuali sono generati dalla formula della colonna G e non dalla formula della colonna Y.
In questo modo, se provo a simulare con Y, l'unica rivalutazione che viene fatta ad ogni simulazione è quella del termine di errore casuale che è nella formula Y, cioè ottengo invece una simulazione del modello a effetti fissi. Come potrei quindi specificare una formula per Y che genera direttamente gli effetti di gruppo? Immagino che sia necessario un po 'di looping, ma non so come farlo.
Penso che tu abbia ragione nel dire che il problema è che la colonna G non viene rivalutata per ogni simulazione ... Questo potrebbe essere risolto includendo l'intera formula in una colonna, forse in questo modo:
3 +: X * 3 + Se ( Riga () == 1 | : group! = Ritardo (: gruppo), r = Casuale normale (), r) + Casuale normale ()
Qui la variabile "r" viene ridefinita per ciascun gruppo. Ciò presuppone che la tabella dei dati sia ordinata per "gruppo". La prima normale casuale nella formula corrisponde all'effetto di gruppo e la seconda è l'errore residuo.
Spero che sia di aiuto.
Adamo
