Opzione Salva colonne con modello mistoCiao a tutti,
Quando eseguiamo un modello misto utilizzando la piattaforma Fit Model, possiamo salvare i residui e i valori previsti, ma abbiamo anche la possibilità di salvare i residui condizionali e i valori previsti condizionali (che tengono conto sia degli effetti fissi che casuali). Tuttavia non abbiamo opzioni simili per i residui studentizzati e le colonne StdErr Pred Formula.
Qual è il modo migliore per ottenere o magari ricalcolare la "verbo condizionale" di queste due colonne?
Si noti che non sto lavorando con JMP Pro.
Grazie in anticipo per il vostro aiuto!
Ri: opzione Salva colonne con modello mistoSì, è vero, lo farò!
Ri: opzione Salva colonne con modello mistoPuoi suggerire questa funzione nella sezione Elenco dei desideri della community JMP.
Ri: opzione Salva colonne con modello mistoCiao GM,
Sì, hai ragione, possiamo estrarre gli errori standard di effetti casuali come quello, ma quello che sto cercando sono
gli errori standard dei residui e dei valori previsti. Non penso che siano le stesse informazioni qui.
Ri: opzione Salva colonne con modello mistoGrazie Cameron per questa risposta dettagliata, la testerò!
Ma in effetti non è un modo semplice per ottenere quello che vogliamo. Mi chiedo il motivo per cui non è direttamente disponibile dal menu Salva colonne.
Ri: opzione Salva colonne con modello mistoSe si stimano gli effetti casuali usando la piattaforma Fit Model con Personality = Minimi quadrati standard e Method = REML, è possibile estrarre gli errori standard degli effetti casuali facendo una tabella di dati fuori dalla tabella Prediction Prediction - gli errori standard sono inclusi. La personalità del Modello misto non emette gli errori standard con i coefficienti Effetto casuale nel report.
Ri: opzione Salva colonne con modello mistoHai ragione. Duh. Non ho accesso a JMP Pro, quindi non posso verificarlo, ma non penso che sarete in grado di ottenere tutti i componenti necessari per calcolare i residui condizionali standardizzati direttamente dal report. La colonna PredSE fornisce la matrice X come secondo argomento di Vec Quadratic () e il report fornisce l'MSE. Dovresti costruire G e Z.
Puoi ottenere Z da qualcosa come:
Design Nom( :RandomEffectCol << Get Values, {list of levels}) Se hai più effetti casuali, potresti semplicemente concatenarli insieme.
Per ottenere G, ciò dipenderà dalla struttura di correlazione scelta. Se si dispone di una struttura di correlazione semplice con tutti gli 0 sulle diagonali disattivate, è possibile costruire G moltiplicando in multiplo il vettore di componenti di varianza con una matrice di identità.
Una volta che hai X, G, Z e MSE, puoi calcolare la varianza della matrice di errore in questo modo:
n = N Row(dt);
V = Z*G*Z` + MSE*Identity(n);
K = Identity(n) - Z*G*Z`*Inv(V);
Q = X*Inv(X`*Inv(V)*X)*X`;
Var_e = K*(V-Q)*K`;
Da lì, si prenderà la sqrt degli elementi diagonali di Var_e per ottenere gli errori standard condizionali per i residui condizionali. Non penso che JMP abbia una funzione per restituire solo gli elementi diagonali della matrice, ma ecco una semplice che farà il trucco:
getDiag = function({X},
n = ncol(X);
diags = {};
for(i=1, i<=N Col(X),i++,
Insert Into(diags,X[i,i])
)
);Questo è un grande sforzo! Scusa, vorrei sapere di più facile.
Ri: opzione Salva colonne con modello mistoGrazie per la tua risposta, Cameron!
Tuttavia, se hai salvato lo Std Err dei residui, penso che sia correlato a quelli che non includono effetti casuali e non ai residui condizionali. Sembra impossibile salvare gli errori di std condizionali dei residui. Penso che possiamo ripetere il calcolo in base ai valori di Hats ma il problema è lo stesso: probabilmente sono collegati ai "residui di base".
Ri: opzione Salva colonne con modello mistoPuoi studiare i residui condizionali dividendo per gli errori standard. Salvare entrambe le colonne nella tabella e quindi calcolare i valori studentizzati con una formula di colonna. Non sono immediatamente sicuro di come aiutarti con la Formula Std Error Pred.
