Logit: perché i miei coefficienti esponenziati non sono uguali al mio rapporto di probabilità (per i predittori nominali)?<meta http-equiv="Content-Type" content="text/html; charset=UTF-8" />Non sono sicuro che questa domanda sia davvero appropriata per questo forum, ma ci proverò comunque 
Sto eseguendo logit in JMP. La mia comprensione è che se si eleva a potenza i coefficienti, si ottiene la variazione dei rapporti di probabilità da una variazione di 1 unità in X.
Funziona perfettamente per variabili indipendenti continue, ma non per variabili indipendenti nominali (0,1). Esempio di seguito:
Per le stime dei parametri sto ottenendo:
Intercetta stima termine -1,516621
Genere[Femmina] 0.22189387
Rapporto probabilità livello 1/livello 2
Maschio Femmina 0.6416016
Femmina Maschio 1.5585996
Il mio odds ratio per i maschi non dovrebbe essere exp(-1.516621)=0.220486 e per le femmine essere exp(0.20461503)=1.2270525?
Confuso da questo. Gradirei davvero un po' di aiuto. Grazie!
Ri: Logit: perché i miei coefficienti esponenziati non sono uguali al mio rapporto di probabilità (per i predittori nominali)?Sto avendo lo stesso problema. La discussione mi ha aiutato molto.
Sembra ancora strano che la categoria debba essere cambiata da nominale a ordinale, ma funziona!
GrazieRi: Logit: perché i miei coefficienti esponenziati non sono uguali al mio rapporto di probabilità (per i predittori nominali)?Ci sono molti modi per aggiungere i dati.Puoi incollare uno screenshot direttamente nel post del forum (cosa che faccio di seguito) oppure puoi utilizzare l'editor avanzato e allegare il file.
Spesso molte persone pubblicano semplicemente il codice della tabella dati / codice di analisi poiché ciò consente di replicare rapidamente il proprio lavoro senza dover scaricare alcun file.Puoi ottenere lo script direttamente dal menu a comparsa in qualsiasi tabella di dati:
Ri: Logit: perché i miei coefficienti esponenziati non sono uguali al mio rapporto di probabilità (per i predittori nominali)?CIAO criscico42 ,
Forse dovresti guardarlo da una prospettiva più ampia. ho l'impressione che jmp sia stato creato pensando in grande. in questo caso circa più di due categorie e quindi se una variabile è nominale (categorie più primitive che non possono essere ordinate) è opportuno confrontarle tutte con una categoria "media". in tal caso il coefficiente mancante può essere calcolato di uno meno la somma degli altri coefficienti (oppure semplicemente richiedendo dal menù le "stime espanse"). nel frattempo, per l'ordinale ha senso confrontare ogni categoria con quella sopra di essa. in questo caso il coefficiente per la categoria mancante è proprio l'intercetta (la costante).
inoltre, in base alla mia esperienza, è più rapido impostare in anticipo l '"ordinamento dei valori" nella tabella piuttosto che dichiarare manualmente i contrasti all'interno della piattaforma del modello, come è pratica comune in altri software statistici.
una volta che il fine settimana sarà finito, mi aspetto che riceverai altri pensieri da altri sul forum.
Ron
Ri: Logit: perché i miei coefficienti esponenziati non sono uguali al mio rapporto di probabilità (per i predittori nominali)?Grazie Ron!Che funzioni!
È un po' controintuitivo per me trattare il genere come una variabile ordinale. Il genere sembra il bambino poster per nominale. Questo significa che dovrei usare l'ordinale per qualsiasi attributo a due livelli (M/F, sposato/single, laureato/nessun college ecc.?). durante l'esecuzione di logit in JMP? Quando mai vorrei usare nominale?
Inoltre mi sembra molto strano che quando uso il valore nominale, i coefficienti sembrino "sbagliati" (cioè non puoi elevarli a potenza per arrivare ai rapporti di probabilità) ma il rapporto di probabilità mostra ancora lo stesso in entrambi i casi (con il valore nominale o variabile ordinale). Non mi sembra giusto dato il modo in cui pensavo che funzionasse la matematica logit.
Sembra che tu stia dicendo che questa è una caratteristica, non un bug, e non presumo di avere abbastanza conoscenze per non essere d'accordo, ma sicuramente crea confusione!
Sarei felice di condividere i dati o l'output del modello (dalla versione "sbagliata") ma non riesco a capire come incollarli qui senza una formattazione folle e incasinata.
Grazie ancora e fammi sapere se hai altri pensieri!
Ri: Logit: perché i miei coefficienti esponenziati non sono uguali al mio rapporto di probabilità (per i predittori nominali)?CIAO criscico42 ,
sarebbe più utile se potessi fornire un dato o almeno l'intero output.
nel frattempo ti suggerisco di definire il genere come ordinale e vedere che il rapporto di probabilità ha più senso per te.
se lo dichiari nominale i coefficienti non si contrappongono tra loro ma piuttosto stimano la differenza dalla media non ponderata tra le categorie.
se lo dichiari ordinale i coefficienti stimano la differenza tra categorie adiacenti (nel caso di genere tra i due) e penso che sia quello che ti aspetti di vedere.
eseguendo il modello una volta con il genere come nominale e ancora una volta come ordinale, dovresti ottenere essenzialmente gli stessi risultati in tutto tranne il coefficiente per il genere. usando l'opzione nominale dovresti ottenere un coefficiente che è la metà del coefficiente ordinale.
si tratta del modo in cui jmp costruisce i contrasti tra le categorie.
spero che questo ti aiuti.
Ron
