Auto-suggest helps you quickly narrow down your search results by suggesting possible matches as you type.

Showing results for

- JMP User Community
- :
- Discussions
- :
- SAS PROC MIXED and JMP Output are Different

Topic Options

- Subscribe to RSS Feed
- Mark Topic as New
- Mark Topic as Read
- Float this Topic for Current User
- Bookmark
- Subscribe
- Printer Friendly Page

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Get Direct Link
- Email to a Friend
- Report Inappropriate Content

Sep 18, 2016 1:01 AM
(1367 views)

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:

Solution for Fixed Effects | ||||||
---|---|---|---|---|---|---|

Effect | Condition | Estimate | Standard Error | DF | t Value | Pr > |t| |

Intercept | 4.4444 | 1.0943 | 4 | 4.06 | 0.0153 | |

Condition | A | 1.1111 | 1.5476 | 4 | 0.72 | 0.5125 |

Condition | B | 0 | . | . | . | . |

Type 3 Tests of Fixed Effects | ||||
---|---|---|---|---|

Effect | Num DF | Den DF | F Value | Pr > F |

Condition | 1 | 4 | 0.52 | 0.5125 |

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:

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

6 REPLIES

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Get Direct Link
- Email to a Friend
- Report Inappropriate Content

Sep 18, 2016 5:46 AM
(1276 views)

HI laughtpp,

since pvalues are the same i would guess the analysis is the same. therefore, the difference is due to JMP estimation of parameters based on a different approach to contrasts of categories. check out the difference in the definition of contrasts between nominal and ordinal variables.

follow this link http://www.jmp.com/support/help/The_Factor_Models.shtml

and look for nominal factors. the table gives an explanation for how the parameters are treated.

Ron

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Get Direct Link
- Email to a Friend
- Report Inappropriate Content

Sep 18, 2016 1:10 PM
(1276 views)

Hey Ron,

Thank you so much. I think you are right about the way SAS and JMP deal with contrast of nominal data differently. I am wondering if you know anyway I can change JMP's default of contrast? I think SAS's results look more reasonable for me (std error are all the same in fixed effects factors in JMP), so I would like to change the JMP's default of contrast. The reason I really want both of them match is because I want to use the beta estimate value of the fixed effects factors to calculate the R^2 of model. If the beta values are not the same, the R^2 would be different. I really appreciate your help.

Sincerely,

Alicia

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Get Direct Link
- Email to a Friend
- Report Inappropriate Content

Sep 19, 2016 4:31 AM
(1276 views)

Three things to try:

1) In other cases I would change the categorical variable from nominal to ordinal modeling type. yet, in this case i think the model platform will not accept it.

2) Try to see if you can perform custom tests for thees parameters in order to "reveal" the contrasts you are looking for.

3) If it is just two categories try to see what you get by coding the categories as zeros and ones and run it as a continuous.

Ron

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Get Direct Link
- Email to a Friend
- Report Inappropriate Content

Sep 20, 2016 11:03 AM
(1276 views)

Hey Ron,

Thank you so much. I used the third method you suggested, and it works!!! Thank you so much.

Best,

Alicia

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Get Direct Link
- Email to a Friend
- Report Inappropriate Content

Sep 20, 2016 1:34 PM
(1276 views)

Hey Alicia,

I am glad it worked out for you. the last method is more of a hack or plumbers trick but it is true to the fundamentals. in the past, it was the actual way to enter categorical variables into regression models. one would build the contrasts himself by recoding the categories into n-1 new columns of zeros and ones. the program would just think they were continuous for estimation purposes.

Best,

Ron

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Get Direct Link
- Email to a Friend
- Report Inappropriate Content

Nov 3, 2016 8:59 AM
(1272 views)

Hey Ron, Alicia: An easier way still may be choosing "Indicator Parameterization Estimates" in the "Estimates" Submenu (red triangle at top of analyses after clicking on "Run") of the Fit Model platform. You can also choose to display these estimates by default (in addition to the JMP default) in the "Fit Least Squares" platform in "Preferences".

Mark