Sto usando effetti annidati nella regressione generalizzata e scopri che la parametrizzazione sembra strana, forse sbagliata? Nelle stime dei parametri mi aspetterei una riga per ogni livello del fattore di nidificazione, con solo il nome di quel livello in []. Invece usa la differenza tra quel livello e il livello di riferimento. Per le ultime righe si ottengono voci con [livello di riferimento - livello di riferimento]. Non ha senso per me. Per illustrare, esegui il codice allegato e vedrai le righe nella tabella delle stime dei parametri, incluse le seguenti:
ter m | Stima | Errore Std | Wald ChiSquare | Prob> ChiSquare | Inferiore al 95% | Superiore al 95% |
... | | | | | | |
Regione [South-Southwest]: Stato [South Carolina-Wyoming] | -3,902821 | 115,44,404 mila | 0.0011429 | 0,9730 | -230,169 | 222,36,334 mila |
| ... | |
Regione [Southwest-Southwest]: Stato [New Mexico-Wyoming] | 1547.9147 | 121,83,319 mila | 161,42,203 mila | <.0001 * | 1309.1261 | 1786.7034 |
Regione [Southwest-Southwest]: Stato [Texas-Wyoming] | 1854.3527 | 115,99,596 mila | 255,56,355 mila | <.0001 * | 1627.0048 | 2081.700 6 |
la sottrazione della categoria di riferimento per Stato va bene per me, è come se la parametrizzazione di SAS (so che è diversa dalle altre piattaforme JMP), ma per la Regione, non ha senso? Come lo interpreto?
Grazie
Gunter
dt = Open( "$SAMPLE_DATA/Seasonal Flu.jmp" );
Fit Model(
Y( :Flu Cases ),
Effects( :Region, :State[:Region] ),
Personality( "Generalized Regression" ),
Generalized Distribution( "Normal" ),
Run(
Fit(
Estimation Method( Standard Least Squares ),
Validation Method( None ),
Normal Quantile Plot( 1 ),
Active Parameter Estimates( 1 )
)
),
SendToReport( Dispatch( {}, "Model Launch", OutlineBox, {Close( 0 )} ) )
);
