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

PROC GLIMMIX DATA=bmd;
  EFFECT AgeSpl = SPLINE(Age / KNOTMETHOD=PERCENTILES(3));
  CLASS Gender;
  MODEL Spnbmd = Gender|AgeSpl;
RUN;

PROC GLIMMIX DATA=bmd;
  EFFECT AgeSpl = spline(Age);
  CLASS Gender;
  MODEL Spnbmd = Gender|AgeSpl;
  ESTIMATE '10m' Intercept 1 Gender 0 1 AgeSpl [1,10] Gender*AgeSpl [1,2 10];
  ESTIMATE '15m' Intercept 1 Gender 0 1 AgeSpl [1,15] Gender*AgeSpl [1,2 15];
  ESTIMATE '20m' Intercept 1 Gender 0 1 AgeSpl [1,20] Gender*AgeSpl [1,2 20];
  ESTIMATE '10f' Intercept 1 Gender 1 0 AgeSpl [1,10] Gender*AgeSpl [1,1 10];
  ESTIMATE '15f' Intercept 1 Gender 1 0 AgeSpl [1,15] Gender*AgeSpl [1,1 15];
  ESTIMATE '20f' Intercept 1 Gender 1 0 AgeSpl [1,20] Gender*AgeSpl [1,1 20];
  ESTIMATE '10d' Gender 1 -1 Gender*AgeSpl [1,1 10] [-1,2 10];
  ESTIMATE '15d' Gender 1 -1 Gender*AgeSpl [1,1 15] [-1,2 15];
  ESTIMATE '20d' Gender 1 -1 Gender*AgeSpl [1,1 20] [-1,2 20];
RUN;