PROC IMPORT OUT= WORK.bmd
            DATAFILE= "Z:\764-F11\data\bmd.txt" 
            DBMS=TAB REPLACE;
     GETNAMES=YES;
     DATAROW=2; 
RUN;

PROC SORT DATA=bmd;
  BY Gender;
RUN;

PROC GAM DATA=bmd;
  CLASS Gender;
  BY Gender;
  MODEL Spnbmd = spline(Age) / METHOD=GCV;
RUN;

PROC GAM DATA=bmd;
  CLASS Gender;
  MODEL Spnbmd = PARAM(Gender) SPLINE(Age);
RUN;

ODS GRAPHICS ON;
PROC GAM DATA=bmd PLOTS=COMPONENTS(CLM ADDITIVE);
  CLASS Gender;
  BY Gender;
  MODEL Spnbmd = spline(Age) / METHOD=GCV;
RUN;
ODS GRAPHICS OFF;