La correlazione tra le stime dei parametri è la correlazione che conta davvero. Questa quantità determina l'inflazione degli errori standard di queste stime. Riduce la potenza dei tuoi test. Allarga i tuoi intervalli di confidenza. La correlazione tra le colonne dei fattori nel progetto è semplicemente un "mezzo per raggiungere il fine". Eliminando questa correlazione si eliminerà, ad esempio, la correlazione tra le stime degli effetti principali.
Si noti che non esiste un limite assoluto per una correlazione inaccettabile perché il VIF tollerabile dipende dalla dimensione dell'effetto e dall'RMSE. Se hai un RMSE molto piccolo rispetto all'effetto, puoi tollerare un VIF elevato. D'altra parte, se hai un effetto piccolo rispetto all'RMSE, anche una piccola correlazione potrebbe influire negativamente sul tuo potere.
La matrice delle correlazioni delle stime NON è il valore assoluto della matrice corr(beta-hat).Ecco un esempio per dimostrare.Nota: ho aggiunto una colonna di 1 per il termine di intercetta.
X = [1 0,53 -1 -1 -0,53 -0,53 1,
1 1 1 1 1 1 1,
1 -1 1 -1 -1 1 -1,
1 -1 -1 1 1 -1 -1,
1 1 1 -1 1 -1 -1,
1 -1 1 1 -1 -1 1,
1 1 -1 1 -1 1 -1];
La "mappa dei colori sulle correlazioni" non include la colonna delle intercettazioni, ma ciò è irrilevante poiché le correlazioni sono a coppie.In ogni caso, se inserisci il progetto di cui sopra nella piattaforma DOE (meno la colonna di 1), o fai lo stesso e usi l'opzione multivariata --> correlazioni, otterrai la seguente matrice di correlazioni tra le colonne della matrice X:
[1 -0,0926 -0,0926 0,2819 0,2819 0,0926,
-0,0926 1 -0,1667 0,0926 0,0926 0,1667,
-0,0926 -0,1667 1 0,0926 0,0926 0,1667,
0,2819 0,0926 0,0926 1 -0,2819 -0,0926,
0,2819 0,0926 0,0926 -0,2819 1 -0,0926,
0,0926 0,1667 0,1667 -0,0926 -0,0926 1]
Se esamini la documentazione di JMP relativa alle correlazioni delle stime, scoprirai che è definita come
corr(beta-hat) := V_inv*(X'X)_inv*V_inv dove V:=sqrt(diag(X'X)_inv)
Utilizzando la matrice X come sopra, si ottiene:
X_t = Trasposizione ( X ) ;
V = mq ( diag ( Inverso ( X_t*X ))) ;
corr_beta_hat = Rotondo ( Inv ( V ) *Inverso ( X_t*X ) *Inv ( V ) ,4 ) ;
[1 -0,3101 -0,3154 -0,3154 0,3101 0,3101 0,3154,
-0,3101 1 0,3101 0,3101 -0,5 -0,5 -0,3101,
-0,3154 0,3101 1 0,3154 -0,3101 -0,3101 -0,3154,
-0,3154 0,3101 0,3154 1 -0,3101 -0,3101 -0,3154,
0,3101 -0,5 -0,3101 -0,3101 1 0,5 0,3101,
0,3101 -0,5 -0,3101 -0,3101 0,5 1 0,3101,
0,3154 -0,3101 -0,3154 -0,3154 0,3101 0,3101 1]
Questo è esattamente ciò che ottieni quando usi il comando "adatta modello" in JMP.Si noti che i valori in questa matrice non dipendono dai valori della risposta o dell'MSE.
Quindi la mia domanda è questa: se si cercasse di decidere se la quantità di confusione tra i termini del modello (o le stime) fosse accettabile, si dovrebbero rivedere le correlazioni tra le colonne della matrice di progettazione o si dovrebbero rivedere le correlazioni tra i beta-hat ?Perché uno potrebbe essere migliore dell'altro?Esiste una spiegazione intuitiva che descrive la relazione tra i due?
Finalmente sono riuscito a trovare una buona spiegazione per la creazione della mappa dei colori sulle correlazioni.In realtà l'ho trovato nella documentazione di JMP sull'adattamento dei modelli lineari.Fitting Linear Models.pdf > Capitolo 3 Report e opzioni dei minimi quadrati standard > Correlazione delle stime, pagina 189 nella documentazione di JMP 11.
Utilizzando queste equazioni, con varianza = 1, la mappa dei colori sul grafico di correlazione è solo il valore assoluto di Corr(Beta-hat)