Thanks for the advice ballardw.    Jake ... View more

Thank you for the responses Reeza and ballardw. I figured out the problem. First my data was in wide format, so I converted it to long format in excel. Does anybody know how to do that conversion in SAS, for future reference?   Next, the ods statement was incorrect. See my correction below   data gpr126_methdiff; infile '/folders/myshortcuts/methylation/GPR126.csv' dsd firstobs=2; input pig $ sex $ gest $ post $ sow $ sire $ gpr126 location $ position $ methperc; /* proc print data=gpr126_methdiff; run; */ proc sort; by position gest; run; proc mixed data=gpr126_methdiff; by position; class pig sow gest post sex sire; model methperc = gest|post|sex; random sow(gest); random sire; ods output Tests3=Tests3; lsmeans gest|post|sex/pdiff; run;   Thanks again for the help!   Jake ... View more