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?
Don't miss out on SAS Innovate - Register now for the FREE Livestream!
Can't make it to Vegas? No problem! Watch our general sessions LIVE or on-demand starting April 17th. Hear from SAS execs, best-selling author Adam Grant, Hot Ones host Sean Evans, top tech journalist Kara Swisher, AI expert Cassie Kozyrkov, and the mind-blowing dance crew iLuminate! Plus, get access to over 20 breakout sessions.
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.