Option "Spalten speichern" mit gemischtem ModellHallo zusammen,
Wenn wir ein gemischtes Modell mit der Fit Model-Plattform ausführen, können wir Residuen und vorhergesagte Werte speichern, aber wir haben auch die Möglichkeit, bedingte Residuen und bedingte vorhergesagte Werte zu speichern (die sowohl feste als auch zufällige Effekte berücksichtigen). Wir haben jedoch keine ähnlichen Optionen für die Spalten der studentisierten Residuen und der StdErr-Pred-Formel.
Was ist der beste Weg, um die "bedingte Überprüfung" dieser beiden Spalten zu erhalten oder möglicherweise neu zu berechnen?
Beachten Sie, dass ich nicht mit JMP Pro arbeite.
Vielen Dank im Voraus für Ihre Hilfe!
Re: Option zum Speichern von Spalten mit gemischtem ModellJa das ist richtig, ich werde das tun!
Re: Option zum Speichern von Spalten mit gemischtem ModellSie können diese Funktion in der Wunschliste der JMP-Community vorschlagen.
Re: Option zum Speichern von Spalten mit gemischtem ModellHallo GM,
Ja, Sie haben Recht, wir können die Standardfehler von zufälligen Effekten wie diesen extrahieren, aber was ich suche, sind
die Standardfehler der Residuen und der vorhergesagten Werte. Ich glaube nicht, dass es hier die gleichen Informationen gibt.
Re: Option zum Speichern von Spalten mit gemischtem ModellVielen Dank, Cameron, für diese detaillierte Antwort, ich werde sie testen!
Aber in der Tat ist es kein einfacher Weg zu bekommen, was wir wollen. Ich frage mich, warum es nicht direkt über das Menü zum Speichern von Spalten verfügbar ist.
Re: Option zum Speichern von Spalten mit gemischtem ModellWenn Sie die Zufallseffekte mithilfe der Fit Model-Plattform mit Personality = Standard Least Squares und Method = REML schätzen, können Sie die Standardfehler der Zufallseffekte extrahieren, indem Sie eine Datentabelle aus der Random Effect Prediction-Tabelle erstellen. Die Mixed Model-Persönlichkeit gibt die Standardfehler nicht mit den Random-Effekt-Koeffizienten im Bericht aus.
Re: Option zum Speichern von Spalten mit gemischtem ModellDu hast recht. Duh. Ich habe keinen Zugriff auf JMP Pro, also kann ich es nicht überprüfen, aber ich glaube nicht, dass Sie alle Komponenten erhalten können, die Sie benötigen, um die standardisierten bedingten Residuen direkt aus dem Bericht zu berechnen. Die Spalte PredSE gibt die X-Matrix als zweites Argument von Vec Quadratic () an, und der Bericht enthält die MSE. Sie müssten G und Z konstruieren.
Sie könnten Z von so etwas bekommen:
Design Nom( :RandomEffectCol << Get Values, {list of levels}) Wenn Sie mehrere zufällige Effekte haben, können Sie diese einfach miteinander verketten.
Um G zu erhalten, würde dies von Ihrer gewählten Korrelationsstruktur abhängen. Wenn Sie eine einfache Korrelationsstruktur mit allen Nullen auf den Off-Diagonalen haben, können Sie G konstruieren, indem Sie den Vektor der Varianzkomponenten mit einer Identitätsmatrix multiplizieren.
Sobald Sie X, G, Z und die MSE haben, können Sie die Varianz der Fehlermatrix folgendermaßen berechnen:
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`;
Von dort würden Sie die Fläche der diagonalen Elemente von Var_e nehmen, um die bedingten Standardfehler für Ihre bedingten Residuen zu erhalten. Ich denke nicht, dass JMP eine Funktion hat, nur die diagonalen Elemente der Matrix zurückzugeben, aber hier ist eine einfache, die den Trick ausführen wird:
getDiag = function({X},
n = ncol(X);
diags = {};
for(i=1, i<=N Col(X),i++,
Insert Into(diags,X[i,i])
)
);Das ist viel Aufwand! Entschuldigung, ich wünschte, ich wüsste einen leichteren Weg.
Re: Option zum Speichern von Spalten mit gemischtem ModellVielen Dank für Ihre Antwort, Cameron!
Wenn Sie jedoch den Std Err von Residuen gespeichert haben, denke ich, dass dies mit denen zusammenhängt, die keine zufälligen Effekte enthalten, und nicht mit den bedingten Residuen. Es scheint nicht möglich zu sein, den Conditional Std Err von Residuen zu speichern. Ich denke, wir können die Berechnung basierend auf den Hats-Werten wiederholen, aber das Problem ist das gleiche: Sie hängen wahrscheinlich mit den "Basis-Residuen" zusammen.
Re: Option zum Speichern von Spalten mit gemischtem ModellSie können die bedingten Residuen durch Division durch ihre Standardfehler studieren. Speichern Sie beide Spalten in der Tabelle und berechnen Sie dann die studentisierten Werte mit einer Spaltenformel. Ich bin mir nicht sofort sicher, wie ich Ihnen mit der Std Error Pred Formula helfen kann.