使用混合模型保存列选项大家好,
当我们使用拟合模型平台运行混合模型时,我们可以保存残差和预测值,但我们也可以保存条件残差和条件预测值(考虑固定和随机效应)。但是,我们没有针对学生化残差和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。