Come calcolare le stime REML dei componenti di varianza in base alla convalida del test biologico USP <1033>?Ciao,
Ho riscontrato un problema con il calcolo dei risultati dell'analisi REML in base ai dati e ai risultati presentati in USP <1033>. Ho difficoltà a calcolare i dati nella Tabella 4 e a ottenere gli stessi risultati contenuti nella Tabella 9. C'è una persona che può gestire questi dati e spiegare come procedere per effettuare l'analisi REML utilizzando questo esempio?

Grazie!
Modellazione statistica avanzataProgettazione degli esperimentiIngegneria della qualità e dei processiOggetto: Come calcolare le stime REML dei componenti di varianza in base alla convalida del test biologico USP <1033>?Per ulteriori informazioni, fornisco le fonti: 1033 VALIDAZIONE DEL TEST BIOLOGICO
Oggetto: Come calcolare le stime REML dei componenti di varianza in base alla convalida del test biologico USP <1033>?Ciao @Luk_s , e benvenuto. Sono in grado di riprodurre le tabelle 5 e 6. Ma non la 9. Sembra che i risultati nella tabella 9 siano un esempio che utilizza dati diversi. Ciò può essere visto dalla differenza in Var(error). Nella tabella 5 è 0,000766. Nella tabella 9 è 0,0022. Dovrebbero essere gli stessi se entrambi utilizzano gli stessi dati. Quindi, la tabella 9 non si basa sui dati mostrati nella tabella 4. È un esempio di qualcos'altro, come implicito nell'affermazione "Le stime dei componenti di varianza ottenute da tale analisi sono presentate nella tabella 9" (corsivo mio); nota che non fa riferimento ai dati nella tabella 4. E, come correttamente indicato nella tabella 5, "L'analista e il lotto dei media non sono stati inclusi a causa del piccolo numero di livelli (2 livelli) per ciascun fattore." , quindi sarebbe sconsigliato eseguire il tipo di analisi mostrato nella tabella 9.
Per quel che vale, e per quanto sconsiderato possa essere, ottengo quanto segue (Tabella 4, Livello 0.5):

Oggetto: Come calcolare le stime REML dei componenti di varianza in base alla convalida del test biologico USP <1033>?Ciao @MRB3855, grazie mille per l'attenzione al problema che ho presentato. Ero convinto di aver sbagliato qualcosa nello sviluppo dei dati per il modello. Ho anche fatto un controllo e siamo usciti con gli stessi risultati, cosa di cui sono davvero soddisfatto. Ho avuto a che fare anche con le tabelle dati precedenti, tuttavia, sono rimasto sorpreso dai risultati ottenuti in questo esempio.

Grazie!