Intervalli di confidenza individuali per modelli mistiCiao a tutti,
Il mio cliente dovrebbe disporre delle equazioni alla base degli intervalli di confidenza individuali (condizionati o meno) per i modelli misti calcolati da JMP.
Non ho trovato alcun riferimento nella documentazione.
Sarebbe possibile avere i dettagli che posso condividere con lei?
Grazie mille,
Firenze
Re: Intervalli di confidenza individuali per modelli mistiFirenze,
Ho utilizzato questo documento per gli intervalli di previsione marginali (vale a dire CI individuali). https://doi.org/10.1002/sim.8386
Quel documento sottolinea che JMP in precedenza utilizzava gradi di libertà non corretti per i suoi PI in SLS-REML. Ho corretto SLS-REML e implementato PI marginali in Mixed.
Le PI condizionali non sono state implementate in Mixed, perché non sono sicuro di cosa significhi (una stima condizionale è già una previsione, BLUP, quindi cos'è una previsione di una previsione??). È stata lasciata così com'è in SLS-REML. Non so quale formula sia stata usata. Sono probabilmente "ragionevoli", ma forse non del tutto corrette.
Spero che questo ti aiuti!
-Elisabetta
Re: Intervalli di confidenza individuali per modelli misti@eclaassen Ora sono perplesso. Nel documento PROC MIXED per MODEL / OUTP= c'è un'equazione generale per la varianza di previsione. Potrebbe essere usata per costruire PI condizionali JMP?
Re: Intervalli di confidenza individuali per modelli mistiRusso,
PROC MIXED non esegue intervalli di previsione (ad esempio per una nuova osservazione). OUTP fornisce BLUP (solo le solite previsioni condizionali con intervalli). OUTPM fornisce previsioni marginali con intervalli. Nessuno di questi è l'intervallo di previsione discusso nel documento e ora implementato in SLS-REML o nel nostro Mixed. Il documento discute un intervallo di previsione marginale, una nuova osservazione senza informazioni di blocco associate. Non so nemmeno cosa sarebbe un intervallo di previsione condizionale: una nuova osservazione che include un nuovo blocco? Che prevede doppiamente cosa esattamente? SLS-REML ha un'opzione per intervalli di previsione condizionali, ovvero prevede oltre un semplice BLUP. È matematicamente corretto nell'errore standard e nei dfs? Non lo so e ne dubito dato che i CI Indiv marginali avevano i DF sbagliati, sebbene l'errore standard fosse per lo più corretto prima della mia correzione di un paio di release fa.
PROC MIXED OUTP == Previsioni condizionali JMP (e intervalli), i soliti BLUP
PROC MIXED OUTPM == Previsioni marginali JMP (e intervalli)
C'è una confusione enorme dovuta al fatto che "previsione" viene utilizzata in modi diversi. JMP etichetta gli intervalli di previsione come CI Indiv, che includono un CI Indiv condizionale in SLS-REML. Ma pensiamo alle stime delle medie condizionali come previsioni grazie a BLUP. Quindi cos'è un intervallo di previsione di una previsione condizionale?
Re: Intervalli di confidenza individuali per modelli mistiAFAIK i limiti di previsione OUTP sono calcolati come previsione puntuale +/- un valore critico aggiustato con df x errore standard di previsione. Una nuova osservazione comporterebbe l'aggiunta di una nuova realizzazione di un errore residuo casuale indipendente a media zero e la varianza di previsione verrebbe aumentata in modo additivo, quindi sembra che i limiti potrebbero essere semplicemente ampliati di conseguenza.
La configurazione è adatta alla simulazione per verificare che i limiti funzionino come pubblicizzato sia per i BLUP che per le previsioni delle nuove osservazioni.
Penso che il caso di previsione individuale sia spesso importante e calcolare questi limiti sarebbe prezioso per i clienti. La terminologia è confusa, ma dovremmo essere in grado di escogitare qualcosa di ragionevole.
