Test di sfericità non visualizzato nei risultati del modelloHo due domande sull'esecuzione di un test di sfericità per misure ripetute ANOVA se le persone sono disposte a sentirmi fuori?
Ecco una descrizione del mio esperimento e dei dati: il mio esperimento ha coinvolto test su diversi media filtranti per rimuovere il fosforo dall'acqua. Ho un totale di 12 scatole (i miei "soggetti") che rappresentano 4 media filtranti in triplice copia. Nella mia tabella di dati JMP, ogni riga rappresenta una scatola, e per le mie colonne, ho una colonna "ID" (fondamentalmente 1-12 per le mie scatole filtro), una colonna che è il media filtrante chiamato "Trattamento" (etichettato "BF" , "NF", "NS" e "WS" e ripetuti 3 volte per ogni replica) e le colonne successive sono concentrazioni di fosforo per ciascun punto temporale (totale di 11 punti temporali, non equamente distanziati). Ho anche incluso il mio file di dati con questo post.
Quindi, le mie domande:
1.) Ho eseguito un test a modelli misti come una misura ANOVA ripetuta sui miei dati, ma un articolo di JMP dice che devo testare la sfericità per vedere se ho bisogno di aggiustare il valore F nel caso in cui violino il test di sfericità. Per testare la sfericità, devo eseguire il test multivariato. Tuttavia, mi mancano dati per determinati punti temporali e JMP butta fuori l'intera materia se mancano dati con l'approccio multivariato, motivo per cui sono andato con l'approccio del modello misto (ho solo 12 soggetti). Come posso eseguire un test di sfericità quindi, che sembra richiedere l'approccio multivariato?
2.) Sto anche cercando di eseguire misure ripetute ANOVA usando l'approccio multivariato su un sottoinsieme dei miei dati almeno per quel sottoinsieme di sfericità. Ho provato a eseguire un modello MANOVA con i miei dati, ma il test di sfericità non si è presentato nei miei risultati. Se riduco il numero di punti tempo che ho (ho 11 punti tempo), che sono i miei valori "y", il test di sfericità appare, ma, a seconda di quanti punti tempo includo, il criterio di Mauchly è significativo o non significativo. Non ho molta familiarità con le statistiche in generale, quindi mi chiedevo se qualcuno sa cosa sta succedendo?
Qualcuno ha idea del perché il test di sfericità si comporti come ho descritto? Ci sono altre opzioni per testare la sfericità che non richiedono un approccio multivariato?
Grazie mille per qualsiasi input! Non sono affatto forte nelle statistiche, quindi qualsiasi aiuto sarebbe apprezzato!
Stefania
Re: test di sfericità non visualizzato nei risultati del modelloSto riscontrando un problema simile con il test di sfericità non mostrato in JMP. Ho la sensazione che sia perché i miei dati non dicono molto o a causa di problemi con le dimensioni del campione, ma purtroppo sono troppo verde per essere sicuro. Posso controllare la casella di correlazione parziale per vedere che l'assunzione di sfericità è comunque violata (senza il test)? Ho allegato il mio set di dati. Sarà grande per qualsiasi aiuto.
Re: test di sfericità non visualizzato nei risultati del modelloGrazie mille per la tua risposta dettagliata! Apprezzo molto il tuo tempo per spiegare i concetti e fornire esempi reali basati sui miei dati!
Stefania
Re: test di sfericità non visualizzato nei risultati del modelloStefania,
Ho visto degli zeri nei tuoi dati, ma nessun valore mancante. L'errore che ho ricevuto era "insufficienti gradi di errore di errore". Con 2 df per ogni gruppo di trattamento che ti dà 8 df per ogni periodo di tempo (ripetere). Quindi mi aspettavo solo 4 periodi di tempo per il test. Per essere onesti o per mostrare ciò che si vede nei dati, usare 4, 12, 72 e 120 mostrerà il problema della sfericità.
In realtà, il tuo set di dati è piuttosto interessante. Non penso che sia necessario un test statistico per notare un'importante differenza nei risultati dei materiali. I dati mostrano un modello di crescita.
Utilizzando l'ampio set di dati che era stato allegato, creare la matrice di covarianza usando l'analisi multivariata e selezionare Covarianza. Vedrete che la struttura della covarianza (e la correlazione) sono abbastanza differenti per i periodi iniziali rispetto a quelli successivi, in cui i valori misurati non cambiano molto. Anche in questo caso dovresti esaminare i dati e la covarianza per vedere che la regola della sfericità viene violata senza un test. A mio parere, un valore p non ha valore qui.
Quindi il prossimo passo è scegliere un test diverso. Dovresti almeno modellare un componente di serie temporali per le tue misure ripetute. (Vedi suggerimento sotto).
Measured Values Valori misurati 
Weibull Growth Models Modelli di crescita di Weibull 
Predicted vs. Actual Previsione vs. Effettivo
Se vuoi quantificare la differenza, ti suggerisco di modellare la crescita per ciascun gruppo di trattamento. Le immagini sopra mostrano i valori misurati, i risultati di Fit Curves usando una curva di Weibull Growth e la perforfmance. Nota che il modello Weibull Growth non è in JMP13, ma è disponibile in JMP14. I valori previsti sono stati salvati nella tabella dati e tracciati rispetto a quelli effettivi con la linea y = x. Gli script per entrambi vengono salvati nella versione impilata della tabella (allegata). Lo script Curve di adattamento potrebbe non funzionare, a seconda della versione di JMP.
Dovresti studiare l'instabilità del trattamento con BF. C'è qualcosa nel materiale, l'ambiente come la posizione della scatola o la temperatura che potrebbe influenzare la precisione delle letture per le concentrazioni di fosforo.
