options linesize=80 pagesize=60; data cfkids; infile 'NewCFkids-SAS.data'; input id fev1 age female pseudoA f508 panc age0 ageL; run; data cfkids; set cfkids; f508_1 = 0; if f508=1 then f508_1 = 1; f508_2 = 0; if f508=2 then f508_2 = 1; run; proc mixed data=cfkids method=reml; class id; model fev1 = age0 ageL female f508_1 f508_2 female*ageL f508_1*ageL f508_2*ageL / s; repeated / type=cs subject=id; run; proc mixed data=cfkids method=reml; class id; model fev1 = age0 ageL female f508_1 f508_2 female*ageL f508_1*ageL f508_2*ageL / s; random intercept / type=un subject=id g; run; proc mixed data=cfkids method=reml; class id; model fev1 = age0 ageL female f508_1 f508_2 female*ageL f508_1*ageL f508_2*ageL / s; random intercept ageL / type=un subject=id g; run; proc mixed data=cfkids method=reml; class id; model fev1 = age0 ageL female f508_1 f508_2 female*ageL f508_1*ageL f508_2*ageL / s; random intercept / type=un subject=id g; repeated / type=sp(pow)(ageL) subject=id; run; proc mixed data=cfkids method=reml; class id; model fev1 = age0 ageL female f508_1 f508_2 female*ageL f508_1*ageL f508_2*ageL / s; random intercept / type=un subject=id g; repeated / type=sp(pow)(ageL) subject=id local; run; proc mixed data=cfkids method=reml; class id; model fev1 = age0 ageL female f508_1 f508_2 female*ageL f508_1*ageL f508_2*ageL / s; random intercept ageL / type=un subject=id g; repeated / type=sp(pow)(ageL) subject=id; run; proc mixed data=cfkids method=reml; class id; model fev1 = age0 ageL female f508_1 f508_2 female*ageL f508_1*ageL f508_2*ageL / s; random intercept ageL / type=un subject=id g; repeated / type=sp(pow)(ageL) subject=id local; run;