Analisi Dose-Risposta Probit<meta http-equiv="Content-Type" content="text/html; charset=UTF-8" />Ciao gruppo,
Sto cercando di riprodurre una semplice analisi dose-risposta che avevo eseguito in SPSS. Quindi in JMP sto usando il modulo "Fit model", selezionando per modelli lineari generalizzati, distribuzione binomiale e collegamento probit.
Ho tre variabili, numero di analisi eseguite (10 analisi), numero di analisi con esito positivo e concentrazione (da 0 a 36)
0 10 2
1 10 4
2 10 9
6 10 16
8 10 25
9 10 36
0 10 0
Seleziono il numero di test eseguiti come variabile "Y", il numero di test con risultato positivo come "Freq" e la concentrazione nei log come variabile su cui modellare gli effetti. Ma quando premo run ottengo questo messaggio nei risultati:
"Attenzione: Hessian non definito positivo. La singolarità è probabilmente dovuta a uno scarso adattamento del modello o a dipendenze lineari tra le covariate del modello. Fit e risultato sono di discutibile valore"
E l'analisi non viene eseguita.
Qualcuno può condividere alcuni commenti?
Grazie,
Dave
Ri: analisi Probit dose-rispostaQuesto è davvero l'ordine giusto per inserire le variabili. Grazie per averlo sottolineato, Michael.
Dave
Ri: analisi Probit dose-rispostaL'ultima riga 0 10 0 dai dati originali per qualche motivo isola l'analisi dal bisogno dell'ordine corretto. Se lo lasci fuori, l'ordine conta.
Michael
Ri: analisi Dose-Risposta ProbitSeguo con interesse questo thread.
Dai documenti JMP:
Se i tuoi dati sono riepilogati come frequenze di eventi e prove, specifica due continui
colonne in questo ordine : un conteggio del numero di successi e un conteggio del numero di
prove. In alternativa, puoi specificare il numero di errori anziché di successi.
(corsivo mio).
Quindi, se cambi l'ordine, vengono inserite le 2 variabili y, penso che vedrai più di quello che ti aspettavi, specialmente se cambi l'asse x in modo che sia logaritmico.
Detto questo, non riesco a spiegare perché l'ordine non sembra avere importanza con i tuoi dati originali.
Michael
Ri: analisi Dose-Risposta ProbitCaro DanO,
scusa se ti disturbo ancora, ma cambiando i numeri alla demo, mi sembra di ottenere risultati sbagliati. Dai un'occhiata a questo set di dati, con le stesse variabili ma numeri diversi. Per le tre concentrazioni più alte ottengo il 100% di positivi
Numero di risultati positivi Numero di test eseguiti Concentrazione
0 1140 0,01
241 1140 0,1
590 1140 0,5
963 1140 1
1134 1140 5
1140 1140 10
1140 1140 50
1140 1140 100
Eseguo il Probit come mi hai suggerito

Dopo aver corso ottengo uno strano risultato che non ha senso per me. Invece di fare percentuali di correzioni sulla Y, sembra che stia facendo qualcosa di diverso. Che cosa sto facendo di sbagliato?

Grazie mille per il tuo aiuto
Dave
Ri: analisi Probit dose-rispostaGrazie mille DanO, questa è una risposta molto chiara. È interessante che JMP sappia che sono collegati.
Per quanto riguarda la trasformazione del registro, sì, ne terrò conto (è stato effettivamente menzionato nel mio post originale, anche se piuttosto vago devo ammettere!)
Migliore,
Dave
Ri: analisi Dose-Risposta ProbitDave,
Anche il tuo "numero di test con risultato positivo" deve essere inserito come variabile "Y". Il ruolo Frequenza serve a indicare quante volte si è effettivamente verificata una riga nella tabella dati. In questo caso, il modello GLM prevede il numero di test eseguiti e il numero di risultati positivi come valori Y "separati". L'analisi li collegherà insieme.

Una cosa che non hai menzionato che viene spesso fatta per questo tipo di lavoro è prendere il registro della concentrazione. L'ho fatto nella finestra di dialogo del modello mostrata sopra. Evidenziare l'effetto del modello di concentrazione, fare clic sul pulsante di trasformazione e scegliere registro. I risultati sono pubblicati qui.

Dan