I want to display estimates from an HLM model for two different groups (variables: out= my outcome, ID = case identifyer, age = measure of time, diagnosis = grouping variable).
proc mixed;
class = ID;
model out = age diagnosis / solution outp=pred;
random= intercept /SUBJECT = ID TYP = UN GCORR;
run;
proc sgplot data=pred;
series y=pred x=age /group = diagnosis;
run;
Resulting plot has lines for individual cases and I'm trying to create a graph with one line per group. Added bonus if you can help with overlaying raw data within same graph and/or adding confidence intervals. Thanks in advance!
Add
proc sort data=pred; by diagnosis age; run;
before the call to proc sgplot.
Add
scatter x=age y=out / group=diagnosis;
before run; to get the raw data points.
What do you get if you drop the "CLASS = ID" variable?
Registration is open! SAS is returning to Vegas for an AI and analytics experience like no other! Whether you're an executive, manager, end user or SAS partner, SAS Innovate is designed for everyone on your team. Register for just $495 by 12/31/2023.
If you are interested in speaking, there is still time to submit a session idea. More details are posted on the website.
Learn how use the CAT functions in SAS to join values from multiple variables into a single value.
Find more tutorials on the SAS Users YouTube channel.