Hi PG, Thanks for your response. However, the article I am referencing provides different calculations of skewness at each level - which is what I would assume. However, when I do as you say above I get identical output when I put L1 versus L2 - I assume I am just changing the name of the dataset, right? Am I supposed to put the OUTP statement somewhere different in the proc mixed statement to get the residual output for level2? Thanks again, Naomi
... View more