使用混合模型保存列選項大家好,
當我們使用擬合模型平台運行混合模型時,我們可以保存殘差和預測值,但我們也可以保存條件殘差和條件預測值(考慮固定和隨機效應)。但是,我們沒有針對學生化殘差和StdErr Pred Formula列的類似選項。
獲得或重新計算這兩列的“條件版權”的最佳方法是什麼?
請注意,我不使用JMP Pro。
在此先感謝您的幫助!
Re:使用混合模型保存列選項是的,我會這樣做的!
Re:使用混合模型保存列選項您可以在JMP社區的“願望清單”部分中建議此功能。
Re:使用混合模型保存列選項你好GM,
是的,你是對的,我們可以提取像這樣的隨機效應的標準錯誤,但我正在尋找的是
殘差和預測值的標準誤差。我不認為這是相同的信息。
Re:使用混合模型保存列選項謝謝Cameron這個詳細的答案,我會測試一下!
但事實上,獲得我們想要的東西並不容易。我想知道為什麼它不能直接從保存列菜單中獲得。
Re:使用混合模型保存列選項如果使用具有Personality =標準最小二乘法和方法= REML的擬合模型平台估計隨機效應,則可以通過從隨機效應預測表中生成數據表來提取隨機效應的標準誤差 - 包括標準誤差。 混合模型個性不會在報告上輸出帶有隨機效應係數的標準誤差。
Re:使用混合模型保存列選項你是對的。 咄。 我無法訪問JMP Pro,因此我無法查看,但我認為您無法直接從報告中獲取計算標準化條件殘差所需的所有組件。 PredSE列將X矩陣作為Vec Quadratic()的第二個參數,報告為您提供MSE。 你必須構建G和Z.
你可以從以下的東西得到Z:
Design Nom( :RandomEffectCol << Get Values, {list of levels}) 如果你有多個隨機效果,你可以將它們連接在一起。
要獲得G,這取決於您選擇的相關結構。 如果你有一個簡單的相關結構,在對角線上有所有0,那麼你可以通過將方差分量的矢量與一個單位矩陣相乘來構造G.
一旦你有了X,G,Z和MSE,你就可以像這樣計算誤差矩陣的方差:
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`;
從那裡,您將採用Var_e的對角元素的sqrt來獲得條件殘差的條件標準誤差。 我不認為JMP有一個函數只返回矩陣的對角線元素,但這裡有一個簡單的方法可以解決這個問題:
getDiag = function({X},
n = ncol(X);
diags = {};
for(i=1, i<=N Col(X),i++,
Insert Into(diags,X[i,i])
)
);這是一個很大的努力!對不起,我希望我知道一個更簡單的方法。
Re:使用混合模型保存列選項謝謝你的回答卡梅隆!
但是如果你保存了殘差的Std Err,我認為它與那些不包含隨機效應而不是條件殘差的那些相關。似乎無法保存殘差的條件標準誤差。我認為我們可以根據帽子值重做計算,但問題是相同的:它們可能與“基本殘差”有關。
Re:使用混合模型保存列選項您可以通過除以標準誤差來學習條件殘差。 將兩列保存到表中,然後使用列公式計算學生化值。 我不能立即確定如何幫助您使用Std Error Pred Formula。