* Percentage of body fat;

DATA DENSITY;
  INFILE 'C:\FAT\DENSITY.TXT';
  INPUT DENSITY PBFAT AGE WEIGHT HEIGHT NECK CHEST ABDOMEN HIP THIGH KNEE ANKLE BICEPS FOREARMS WRIST;
  LABEL DENSITY='Density determined from underwate weighing'
        PBFAT='Percent body fat from Siris equation'
        AGE='Age in years'
        WEIGHT='Weight in lbs'
        HEIGHT='Height in inches'
        NECK='Neck circumference in cm'
        CHEST='Chest circumference in cm'
        ABDOMEN='Adbomen 2 circumference in cm'
        HIP='Hip circumference in cm'
        THIGH='Thigh circumference in cm'
        KNEE='Knee circumference in cm'
        ANKLE='Ankle circumference in cm'
        BICEPS='Biceps (extended) circumference in cm'
        FOREARMS='Forearm cirumference in cm'
        WRIST='Wrist circumference in cm';
RUN;

PROC PRINT DATA=DENSITY;
  TITLE 'PRINT DATA';
RUN;

PROC MEANS DATA=DENSITY;
RUN;

PROC UNIVARIATE DATA=DENSITY;
RUN;

PROC CORR DATA=DENSITY;
  TITLE 'Correlation Matrix';
  VAR ALL;
RUN;

PROC REG DATA=DENSITY;
  TITLE 'Regression Analysis';
  MODEL DENSITY = PBFAT AGE WEIGHT HEIGHT NECK CHEST ABDOMEN HIP THIGH KNEE ANKLE BICEPS FOREARMS WRIST /
                SELECTION=FORWARD;
RUN;

PROC REG DATA=DENSITY;
  TITLE 'Regression Analysis';
  MODEL DENSITY = PBFAT AGE WEIGHT HEIGHT NECK CHEST ABDOMEN HIP THIGH KNEE ANKLE BICEPS FOREARMS WRIST /
                SELECTION=BACKWARD;
RUN;

PROC GPLOT DATA=DENSITY;
  PLOT DENSITY*PBFAT;
RUN;

PROC REG DATA=DENSITY;
  TITLE 'Regression Analysis';
  MODEL DENSITY = PBFAT;
  PLOT PREDICTED. *PBFAT = 'P' DENSITY*PBFAT='*' / OVERLAY;
RUN;