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を報告する近似のエラーが発生すると思います。 、Accuracy40のように。
お役に立てれば!、
DS
なので@DS主張されているように、無制限の推定オプションは、REMLを使用して分散成分を推定した結果を決定します。 無制限の見積もりは、負の見積もりになる可能性があります。 分散を負にすることはできませんが、これらの結果は実際には共分散であり、負になる可能性があります。 このオプションは、より良い仮説検定とBLUPを提供します。 モデルに固定効果を含めるかどうかは重要です。 (スクリプトは固定効果を示していません。)
推定値を制限し、非負になるように制約することができます。 負の見積もりは強制的にゼロになります。