Hello,   I have a large dataset (n= 78,000) and am trying to run a hierarchical regression, controlling for multiple variables. I am new to using SAS and I am having difficulty getting SAS to run the regression reading all four levels of the province data that I have. The class level information from the results indicates it is only reading "Ontario's" data, and not the other 3 (AB, BC, QC) provinces. Other analysis I have run, have read each level of data fine.    This is the code I am using: proc glimmix data=work.combined; CLASS approach province income enrolment urban CANNABIS_FREQ BINGE_DRINKING WEEKLY_SPENDING SMOKING_STATUS; MODEL CANNABIS_FREQ = SEX3 GRADE GETHNCL1 GETHNCC1 GETHNCN1 GETHNCA1 GETHNCH1 GETHNCM1 SMOKING_STATUS BINGE_DRINKING WEEKLY_SPENDING / SOLUTION; RANDOM CANNABIS_FREQ / SUBJECT=approach SOLUTION; run;    and I attached my results output.    Any help is greatly appreciated! ... View more