Hey all I have a question which I cannot seem to solve. I am working with brain imaging data. I have two "treatment" conditions (before and after) as well as several brain "regions" (n=18). In each region I would like to examine the effect from my treatment. The test is a repeated meassure design where the same "subjects" undergo both treatment and non-treatment condition. First I am interested in the two way ANOVA for "condition" "Region" and "condition*Region". This seems to be working just great! However secondly I want to do post-hoc testing to see which brain regions are significantly affected by my treatment. am therefore only interested in how my treatment affects each indivual brain region (the interaction between "treatment*region") not if different brain regions differ from one another! In the code below SAS only provides me with output of the "treatment*region" interaction. However as far as I can tell SAS still has done all the other multiple comparisons between brain regions. Since I have 18 brain regions this means that I end up with no significant results.. Here is my code: PROC sort data=DataPet; by Region subjects Condition; proc print data=DataPet; run; proc mixed data=DataPet; class subjects Region condition; model Activity=condition Region condition*Region; lsmeans Region*condition/diff adjust=Bon;  ods output diffs=temp; Repeated / Subject=subjects; run; data DataPet2; set temp; if Region=_Region; proc print data=DataPet2; var Region _Region condition _condition estimate stderr df tvalue probt; run; I would be very happy if anyone could tell me if I am doing the multiple comparison correctly. And if not, how I can modify my code to get the post-hoc output I am interested in, namely only the treatment*region interaction. Thank you very much in advance! ... View more