안녕하세요,
REML을 통해 최소 제곱 평균을 생성하려고합니다. 아이디어는 내 연구에서 각 복제 된 개인에 대해 단일 값을 얻는 것입니다 (즉, 개인 당 3 명의 담당자가 있고 개인을 나타내는 하나의 숫자를 원합니다). 과거에 JMP 14에서 (동일한 데이터에 동일한 모델로)이 작업을 수행했을 때 lsmeans는 산술 평균에 매우 가깝습니다. 이제 JMP 15에서 모든 개인에 대한 lsmeans는 서로 매우 가깝고 개인당 산술 평균과 매우 다릅니다. 이것을 설명하는 버전들 사이에 어떤 변화가 있습니까?
도움이나 제안에 미리 감사드립니다!
따라서 JMP 14와 JMP 15 모두에서 모델 맞춤 실행 대화 상자에서 동일한 설정으로 동일한 데이터 세트를 분석했지만 LSMeans가 동일하지 않습니까? 나는이 질문이 명백해 보이지만 과거에 두 경우 모두 분석이 동일하지 않은 것으로 여러 번 발생했음을 알고있다. 플랫폼의 스크립트를 저장하고 여기에 게시 할 수 있습니다. 그 행동은 재현성에 먼 길을 갈 것입니다.
또한 변수가 개별 인 경우 랜덤 효과 모델 및 REML 추정을 사용할 필요가 없습니다. 분산 분석 모형은 유효한 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 속성을 제거하면 최소 제곱 맞춤에 LSMeans도 제공됩니다. 나는 당신이 Oneway 플랫폼을 사용해야한다는 것을 의미하지 않았습니다.
그리고 더 생각해 보면, 개인으로 인한 변이는 개체가 대상체의 표본을 나타내는 경우 무작위 효과로 취급되어야합니다. LSMeans와 표준 오류는 : UniqueID를 고정 효과 또는 임의 효과로 취급하는지에 따라 달라집니다.
