在使用JMP对我们的数据集进行统计分析期间,我们遇到以下问题,我想与您分享一下,询问其他JMP用户是否曾经遇到过这种情况,并且知道为什么会发生以及如何解决。
Fit Model(
Y( :Measured RP% ),
Effects,
Random Effects( :Assay ),
NoBounds( 0 ),
Personality( "Standard Least Squares" ),
Method( "REML" ),
Emphasis( "Minimal Report" ),
Run(
:Measured RP% << {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 )}
),
Where( :Sample == "Accuracy 40" ),
SendToReport(
Dispatch(
{"Response Measured RP% Sample=Accuracy 40"},
"Parameter Estimates",
OutlineBox,
{Close( 0 )}
),
Dispatch(
{"Response Measured RP% Sample=Accuracy 40", "Parameter Estimates"},
"Std Error",
NumberColBox,
{Visibility( "Collapse" )}
),
Dispatch(
{"Response Measured RP% Sample=Accuracy 40", "Random Effect Predictions"
},
"Std Error",
NumberColBox,
{Visibility( "Collapse" )}
),
Dispatch(
{"Response Measured RP% Sample=Accuracy 40",
"REML Variance Component Estimates"},
"Std Error",
NumberColBox,
{Visibility( "Collapse" )}
),
Dispatch(
{"Response Measured RP% Sample=Accuracy 40",
"REML Variance Component Estimates"},
"95% Lower",
NumberColBox,
{Visibility( "Collapse" )}
),
Dispatch(
{"Response Measured RP% Sample=Accuracy 40",
"REML Variance Component Estimates"},
"Sqrt Variance Component",
NumberColBox,
{Visibility( "Visible" )}
),
Dispatch(
{"Response Measured RP% Sample=Accuracy 40",
"REML Variance Component Estimates"},
"Pct of Total",
NumberColBox,
{Visibility( "Collapse" )}
),
Dispatch(
{"Response Measured RP% Sample=Accuracy 40",
"REML Variance Component Estimates"},
"CV",
NumberColBox,
{Visibility( "Visible" )}
),
Dispatch(
{"Response Measured RP% Sample=Accuracy 40", "Effect Details", "Assay",
"Least Squares Means Table"},
"Std Error",
NumberColBox,
{Visibility( "Collapse" )}
)
)
)你好@lefkiosvanrooij ,
我能够基于DOCX文件中包含的内容重新创建使用较小的数据表集时遇到的错误。
因此,我认为这与您在“拟合模型”窗口中选择(或取消选择)“无限制方差分量”有关。 如果选中它,则CI的计算方式与未选中时不同,请参阅JMP帮助。这里。 可能是您不能使用REML方法,在这种情况下,您可以使用EMS方法来估计方差。 您还可以选中“仅估计方差组件”框。
我不是统计学家,所以我无法评论哪种方法适合您的分析,但是我相信如果未选中“无限制方差成分”,那么对于某些样本,您会得到拟合报告为0的错误,例如Accuracy 40。
希望这可以帮助!,
DS
如@DS声称,无边界估计选项确定使用REML估计方差分量的结果。 无限制的估计可能会导致负估计。 方差不能为负,但是这些结果实际上是协方差,可以为负。 此选项提供更好的假设检验和BLUP。 如果在模型中包括固定效果,则很重要。 (您的脚本没有显示固定的效果。)
您可以限制估计值并将其约束为非负数。 任何负估计都必须为零。