How to add group code in pairwise comparison in SAS Proc mixed

When there are a lot of groups to compare, lsmeans/pdiff option gives too many pairwise comparisons. What I need is a table like

enter image description here

which shows that Group 1,2 are different from Group 5,6, but not different from Group 3,4.

Is there an option in SAS proc mixed or other procedures to do this?

asked Nov 16 '18 at 1:52

Ben Z

605

add a comment |

When there are a lot of groups to compare, lsmeans/pdiff option gives too many pairwise comparisons. What I need is a table like

enter image description here

which shows that Group 1,2 are different from Group 5,6, but not different from Group 3,4.

Is there an option in SAS proc mixed or other procedures to do this?

asked Nov 16 '18 at 1:52

Ben Z

605

add a comment |

When there are a lot of groups to compare, lsmeans/pdiff option gives too many pairwise comparisons. What I need is a table like

enter image description here

which shows that Group 1,2 are different from Group 5,6, but not different from Group 3,4.

Is there an option in SAS proc mixed or other procedures to do this?

asked Nov 16 '18 at 1:52

Ben Z

605

When there are a lot of groups to compare, lsmeans/pdiff option gives too many pairwise comparisons. What I need is a table like

enter image description here

which shows that Group 1,2 are different from Group 5,6, but not different from Group 3,4.

Is there an option in SAS proc mixed or other procedures to do this?

sas

asked Nov 16 '18 at 1:52

Ben Z

605

asked Nov 16 '18 at 1:52

Ben Z

605

asked Nov 16 '18 at 1:52

Ben Z

605

asked Nov 16 '18 at 1:52

Ben Z

605

asked Nov 16 '18 at 1:52

Ben Z

605

add a comment |

2 Answers
2

active

oldest

votes

You are looking for the LINES option on the LSMEANS statement however it does not work in PROC MIXED so you will need to use PROC PLM.

proc mixed data=sashelp.class;

   class age;

   model weight=age;

   lsmeans age; *for check;

   store out=classmodel;

   run;

   quit;

proc plm restore=classmodel;

   lsmeans age / lines;

   run;

   quit;

enter image description here

answered Nov 16 '18 at 14:36

data _null_

5,487710

Awesome, that's what I have been looking for. Thank you!

– Ben Z
Nov 16 '18 at 15:53

add a comment |

Try using a custom format that is applied during the Proc.

proc format;

  value groupCat

    1 = 'A'

    2 = 'A'

    3 = 'AB'

    4 = 'AB'

    5 = 'B'

    6 = 'B'

  ;

run;



proc mixed …;

  … ;

  format group groupCat.;

  …

answered Nov 16 '18 at 12:29

Richard

9,63721329

Thanks, but that's not what I meant.

– Ben Z
Nov 16 '18 at 15:55

add a comment |

Your Answer

StackExchange.ifUsing("editor", function () {
StackExchange.using("externalEditor", function () {
StackExchange.using("snippets", function () {
StackExchange.snippets.init();
});
});
}, "code-snippets");

StackExchange.ready(function() {
var channelOptions = {
tags: "".split(" "),
id: "1"
};
initTagRenderer("".split(" "), "".split(" "), channelOptions);

StackExchange.using("externalEditor", function() {
// Have to fire editor after snippets, if snippets enabled
if (StackExchange.settings.snippets.snippetsEnabled) {
StackExchange.using("snippets", function() {
createEditor();
});
}
else {
createEditor();
}
});

function createEditor() {
StackExchange.prepareEditor({
heartbeatType: 'answer',
autoActivateHeartbeat: false,
convertImagesToLinks: true,
noModals: true,
showLowRepImageUploadWarning: true,
reputationToPostImages: 10,
bindNavPrevention: true,
postfix: "",
imageUploader: {
brandingHtml: "Powered by u003ca class="icon-imgur-white" href="https://imgur.com/"u003eu003c/au003e",
contentPolicyHtml: "User contributions licensed under u003ca href="https://creativecommons.org/licenses/by-sa/3.0/"u003ecc by-sa 3.0 with attribution requiredu003c/au003e u003ca href="https://stackoverflow.com/legal/content-policy"u003e(content policy)u003c/au003e",
allowUrls: true
},
onDemand: true,
discardSelector: ".discard-answer"
,immediatelyShowMarkdownHelp:true
});

}
});

draft saved

draft discarded

Sign up or log in

StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});

Post as a guest

Name

Required, but never shown

StackExchange.ready(
function () {
StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fstackoverflow.com%2fquestions%2f53330362%2fhow-to-add-group-code-in-pairwise-comparison-in-sas-proc-mixed%23new-answer', 'question_page');
}
);

Post as a guest

Name

Required, but never shown

2 Answers
2

active

oldest

votes

2 Answers
2

active

oldest

votes

You are looking for the LINES option on the LSMEANS statement however it does not work in PROC MIXED so you will need to use PROC PLM.

proc mixed data=sashelp.class;

   class age;

   model weight=age;

   lsmeans age; *for check;

   store out=classmodel;

   run;

   quit;

proc plm restore=classmodel;

   lsmeans age / lines;

   run;

   quit;

enter image description here

answered Nov 16 '18 at 14:36

data _null_

5,487710

Awesome, that's what I have been looking for. Thank you!

– Ben Z
Nov 16 '18 at 15:53

add a comment |

You are looking for the LINES option on the LSMEANS statement however it does not work in PROC MIXED so you will need to use PROC PLM.

proc mixed data=sashelp.class;

   class age;

   model weight=age;

   lsmeans age; *for check;

   store out=classmodel;

   run;

   quit;

proc plm restore=classmodel;

   lsmeans age / lines;

   run;

   quit;

enter image description here

answered Nov 16 '18 at 14:36

data _null_

5,487710

Awesome, that's what I have been looking for. Thank you!

– Ben Z
Nov 16 '18 at 15:53

add a comment |

You are looking for the LINES option on the LSMEANS statement however it does not work in PROC MIXED so you will need to use PROC PLM.

proc mixed data=sashelp.class;

   class age;

   model weight=age;

   lsmeans age; *for check;

   store out=classmodel;

   run;

   quit;

proc plm restore=classmodel;

   lsmeans age / lines;

   run;

   quit;

enter image description here

answered Nov 16 '18 at 14:36

data _null_

5,487710

You are looking for the LINES option on the LSMEANS statement however it does not work in PROC MIXED so you will need to use PROC PLM.

proc mixed data=sashelp.class;

   class age;

   model weight=age;

   lsmeans age; *for check;

   store out=classmodel;

   run;

   quit;

proc plm restore=classmodel;

   lsmeans age / lines;

   run;

   quit;

enter image description here

answered Nov 16 '18 at 14:36

data _null_

5,487710

answered Nov 16 '18 at 14:36

data _null_

5,487710

answered Nov 16 '18 at 14:36

data _null_

5,487710

answered Nov 16 '18 at 14:36

data _null_

5,487710

Awesome, that's what I have been looking for. Thank you!

– Ben Z
Nov 16 '18 at 15:53

add a comment |

Awesome, that's what I have been looking for. Thank you!

– Ben Z
Nov 16 '18 at 15:53

Awesome, that's what I have been looking for. Thank you!

– Ben Z
Nov 16 '18 at 15:53

add a comment |

Try using a custom format that is applied during the Proc.

proc format;

  value groupCat

    1 = 'A'

    2 = 'A'

    3 = 'AB'

    4 = 'AB'

    5 = 'B'

    6 = 'B'

  ;

run;



proc mixed …;

  … ;

  format group groupCat.;

  …

answered Nov 16 '18 at 12:29

Richard

9,63721329

Thanks, but that's not what I meant.

– Ben Z
Nov 16 '18 at 15:55

add a comment |

Try using a custom format that is applied during the Proc.

proc format;

  value groupCat

    1 = 'A'

    2 = 'A'

    3 = 'AB'

    4 = 'AB'

    5 = 'B'

    6 = 'B'

  ;

run;



proc mixed …;

  … ;

  format group groupCat.;

  …

answered Nov 16 '18 at 12:29

Richard

9,63721329

Thanks, but that's not what I meant.

– Ben Z
Nov 16 '18 at 15:55

add a comment |

Try using a custom format that is applied during the Proc.

proc format;

  value groupCat

    1 = 'A'

    2 = 'A'

    3 = 'AB'

    4 = 'AB'

    5 = 'B'

    6 = 'B'

  ;

run;



proc mixed …;

  … ;

  format group groupCat.;

  …

answered Nov 16 '18 at 12:29

Richard

9,63721329

Try using a custom format that is applied during the Proc.

proc format;

  value groupCat

    1 = 'A'

    2 = 'A'

    3 = 'AB'

    4 = 'AB'

    5 = 'B'

    6 = 'B'

  ;

run;



proc mixed …;

  … ;

  format group groupCat.;

  …

answered Nov 16 '18 at 12:29

Richard

9,63721329

answered Nov 16 '18 at 12:29

Richard

9,63721329

answered Nov 16 '18 at 12:29

Richard

9,63721329

answered Nov 16 '18 at 12:29

Richard

9,63721329

Thanks, but that's not what I meant.

– Ben Z
Nov 16 '18 at 15:55

add a comment |

Thanks, but that's not what I meant.

– Ben Z
Nov 16 '18 at 15:55

Thanks, but that's not what I meant.

– Ben Z
Nov 16 '18 at 15:55

add a comment |

draft saved

draft discarded

Thanks for contributing an answer to Stack Overflow!

Please be sure to answer the question. Provide details and share your research!

But avoid …

Asking for help, clarification, or responding to other answers.

Making statements based on opinion; back them up with references or personal experience.

To learn more, see our tips on writing great answers.

draft saved

draft discarded

Sign up or log in

StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});

Post as a guest

Name

Required, but never shown

Post as a guest

Name

Required, but never shown

Sign up or log in

StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});

Post as a guest

Name

Required, but never shown

Sign up or log in

StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});

Post as a guest

Name

Required, but never shown

Sign up or log in

StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});

Post as a guest

Name

Required, but never shown

Name

Required, but never shown

Name

Required, but never shown

This page is only for reference, If you need detailed information, please check here

搜尋此網誌

Vfrdtyky