Modelli lineari generalizzati senza intercettaCIAO,
Ho bisogno del vostro aiuto. Voglio fare un'analisi con GLM, ma senza intercettare. Non riesco a scegliere l'opzione "Nessuna intercettazione".
Voglio ottenere gli stessi risultati come lì: http://support.sas.com/kb/24/992.html
Non ho idea di come ottenere l'opzione "No Intercept" in GLM, distribuzione=Binomiale, Link=Logit
Sarò molto grato per il tuo aiuto
Re: Modelli lineari generalizzati senza intercettaCIAO @mhans ,
È possibile in JMP standard, ma richiede una leggera ristrutturazione della tabella. Ciò di cui hai bisogno è una variabile di risultato nominale di Vittoria o Perdita, e quello che hai è una tabella riepilogativa delle frequenze di tali eventi. Ho allegato qui una tabella così ampliata, con lo script logistico nominale salvato (discusso di seguito più avanti). Nota la forma di queste tabelle, nell'immagine sotto c'è la tabella originale con la prima riga selezionata, in basso c'è la tabella espansa, con le 13 righe selezionate che corrispondono a quella riga riepilogativa. Nota che ce ne sono 13, ovvero i 13 scontri tra Mil e Det, e ci sono 7 vittorie e 6 sconfitte, che riflettono la vittoria = 7 nella tabella in alto. 
Il modello è semplice da impostare: Analizza > Adatta modello, con Risultato come Y, tutti i team come effetti del modello e Nessuna intercettazione selezionata. 
Ciò restituisce il seguente output, che vedrai corrispondere:
Fit Model(
Y( :Outcome ),
Effects( :mil, :det, :tor, :new, :bos, :cle, :bal ),
No Intercept( 1 ),
Choose High Target( 1 ),
Personality( "Nominal Logistic" ),
Run( Likelihood Ratio Tests( 1 ), Wald Tests( 0 ) ))
Il valore aggiunto della regressione generalizzata qui è che puoi specificare direttamente la struttura vincita/totale quando definisci la variabile y, che gestisce l'espansione in background.
@giuliano
Re: Modelli lineari generalizzati senza intercettaCIAO @giuliano ,
wow, grazie! Sto vedendo che è possibile in JMP. Ho solo la versione "normale" di JMP. È possibile eseguire questo modello con la versione normale (come regressione logistica senza intercetta o GLM)?
In caso negativo, potete inviarmi l'output nel modulo pdf.-form?
Tanto di cappello per questa soluzione.
Re: Modelli lineari generalizzati senza intercettaCIAO @mhans ,
Capirò meglio cosa stai cercando di fare adesso. Per ottenere una corrispondenza con questi risultati (modello Bradley-Terry) sarà necessario utilizzare la piattaforma di regressione generalizzata, disponibile solo in JMP Pro. Ho allegato una versione JMP della tabella dati utilizzata nell'esempio SAS e ho salvato lo script nella tabella. Ecco i passaggi per riprodurlo:
1. Avviare Analizza > Stima modello.
2. Specificare Regressione generalizzata come Personalità e Distribuzione come Binomiale.
3. Immettere sia Vincita che Totale come Y, Colonne e immettere tutti i predittori come effetti del modello.
4. Seleziona "Nessuna intercettazione" in basso.
La finestra di dialogo Adatta modello ora dovrebbe assomigliare a questa:
5. Fare clic su Esegui, quindi nella finestra di dialogo risultante specificare Regressione logistica come metodo di stima e fare clic su Vai:
Ciò restituirà il seguente output, che puoi confermare corrisponde ai risultati SAS.
Ecco il JSL per questa analisi, anch'esso salvato nella tabella allegata.
Fit Model(
Y( :win, :total ),
Effects( :mil, :det, :tor, :new, :bos, :cle, :bal ),
No Intercept( 1 ),
Personality( "Generalized Regression" ),
Generalized Distribution( "Binomial" ),
Run(
Fit( Estimation Method( Logistic Regression ), Validation Method( None ) )
)
);
Spero che aiuti!
@giuliano
Re: Modelli lineari generalizzati senza intercettaCIAO @giuliano ,
grazie per la tua risposta. Mi sono stancato di usare la personalità logistica nominale. La variabile di risposta y dovrebbe essere nominale, la mia variabile y era metrica (tra 0 e 1, una percentuale di vincite). Ho creato una nuova struttura per i miei dati con la variabile y: 1 o -1 e ho copiato i dati nel modulo
1 -1 0 0 0 0 0 7 1
1 0 -1 0 0 0 0 9 1
(...)
0 0 0 0 0 1 -1 6 1
(new part of data)
-1 1 0 0 0 0 0 6 -1
-1 0 1 0 0 0 0 4 -1
(...)
0 0 0 0 0 -1 1 7 -1
I risultati non sono gli stessi di SAS.Non ho idea di come ripararlo...
Re: Modelli lineari generalizzati senza intercettaCIAO @mhans ,
non mi piace il Modello lineare generalizzato personalità, il Logistica nominale La personalità in Stima modello può essere eseguita senza intercettazione selezionando la casella "Nessuna intercettazione" nella finestra di dialogo Stima modello. Questa personalità restituirà risultati equivalenti alla personalità del modello lineare generalizzato con una distribuzione binomiale dell'errore e una funzione di collegamento logit, presupponendo che tu non abbia selezionato le stime Firth Bias-Adjusted (con questa casella selezionata le stime del modello differiranno). Forse questo funzionerà per la tua situazione?
@giuliano
