Formatting sglot to display group data from model predictions

[cm2005]

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!

[PGStats]

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.

[DanH_sas]

What do you get if you drop the "CLASS = ID" variable?