你好
我正在尝试通过REML生成最小二乘法。想法是为研究中的每个重复个体获得一个单一的值(即,每个个体有3个代表,我希望一个数字代表一个个体)。过去,当我在JMP 14中使用相同的模型对相同的数据进行处理时,lsmeans非常接近算术平均值。现在,在JMP 15中,所有个体的lsmeans彼此非常接近,并且与每个个体的算术平均值非常不同。版本之间有什么变化可以解释这一点吗?
在此先感谢您的帮助或建议!
因此,您在JMP 14和JMP 15的“拟合模型”启动对话框中使用相同的设置分析了相同的数据集,但是LSMean不同吗?我意识到这个问题似乎很明显,但是过去发生过很多次,两种情况下的分析都不相同。也许您可以保存该平台的脚本并将其发布在此处。该法案将大大促进可复制性。
此外,如果唯一变量是单个变量,则无需使用随机效应模型和REML估计。ANOVA模型还将产生有效的LSMeans估计。
感谢您的快速答复。令人尴尬的是,分析似乎并不相同...
但是,关于使用ANOVA生成lsmeans:当我使用ANOVA时,我得到的结果与使用REML时截然不同。两种分析中的特征(“ ADI”)和效果(“ UniqueID”)相同。关于为什么会如此的任何见解?(很抱歉,如果这应该是一个不同的讨论主题,请在这里首次发布,如果最终成为我的另一个愚蠢错误,我们也深表歉意!)
REML:
Fit Model(
Y( :ADI ),
Effects,
Random Effects( :UniqueID ),
NoBounds( 1 ),
Personality( "Standard Least Squares" ),
Method( "REML" ),
Emphasis( "Minimal Report" ),
Run(
:ADI << {Summary of Fit( 1 ), Analysis of Variance( 0 ),
Parameter Estimates( 1 ), Scaled Estimates( 0 ),
Plot Actual by Predicted( 0 ), Plot Regression( 0 ),
Plot Residual by Predicted( 0 ), Plot Studentized Residuals( 0 ),
Plot Effect Leverage( 0 ), Plot Residual by Normal Quantiles( 0 )}
),
SendToReport(
Dispatch( {"Response ADI"}, "Effect Details", OutlineBox, {Close( 0 )} ),
Dispatch(
{"Response ADI", "Effect Details", "UniqueID"},
"Least Squares Means Table",
OutlineBox,
{Close( 0 )}
)
)
);方差分析:
Oneway( Y( :ADI ), X( :UniqueID ), Means( 1 ), Mean Diamonds( 1 ) )
我很抱歉。我的意思是,如果您在Fit Model启动对话框中从:UniqueID效果中删除&Random属性,则Fit Least Squares也应给出LSMeans。我并不是说您应该使用Oneway平台。
而且,我想得更多,如果个体代表的样本是受试者群体,那么由于个体引起的变异应视为随机效应。根据将:UniqueID视为固定效果还是随机效果,LSMeans及其标准错误将有所不同。