Dear All,
I am trying to compare the results of SAS PROC Mixed and JMP mixed model. I thought the results would be the same, but they are not. I hope JMP can give the same results of SAS PROC MIXED. Could anyone help me with it?
My SAS code as following:
PROC MIXED DATA = DATA;
CLASS Sub Condition Stimuli;
MODEL Y = Condition / SOLUTION;
RANDOM Sub;
RANDOM Stimuli(Condition);
run;
The results for the fixed effect are:
|
|
| 4.4444 | 1.0943 | 4 | 4.06 | 0.0153 |
---|
A | 1.1111 | 1.5476 | 4 | 0.72 | 0.5125 |
---|
B | 0 | . | . | . | . |
---|
Then, I used JMP's interface to generate the following code:
Fit Model(
Y( :Y ),
Effects( :Condition ),
Random Effects( :Sub, :Stimuli[:Condition] ),
NoBounds( 0 ),
Personality( "Mixed Model" ),
Run( Repeated Effects Covariance Parameter Estimates( 0 ) )
)
I got:
![12760_pastedImage_1.png 12760_pastedImage_1.png](https://community.jmp.com/t5/image/serverpage/image-id/3722i6EF590B08B5130C2/image-size/medium?v=v2&px=400)
![12761_pastedImage_5.png 12761_pastedImage_5.png](https://community.jmp.com/t5/image/serverpage/image-id/3723i2345ED37C43073AD/image-size/medium?v=v2&px=400)
It seems like Type 3 test of fixed effects are the same in JMP and SAS, but my fixed effects parameter estimates are totally different. I am wondering if there is anyone know how could i code the JMP to make it is the same results as SAS.
Thank you all so much,
Alicia