options ls=80 nodate nonumber;

* 6451  HW 1 problem;
* reads hw1_CRP.csv  from the Desktop;
 
data a;
  infile "C:\Documents and Settings\Administrator\Desktop\hw1_CRP.csv"
     DSD dlm="," firstobs=2;
  input id treatment$ sex$ age WT_kg BMI Base_CRP Base_Chol 
     Final_CRP Final_Chol;

data b;
  set a;
  CRP_change = Final_CRP - Base_CRP;
  Chol_change = Final_Chol - Base_Chol;

proc glm data=b;
  class treatment;
  model age WT_kg BMI Base_CRP Base_Chol 
        CRP_change Chol_change = treatment;
  lsmeans treatment / pdiff stderr;
  means treatment;

proc freq data=b;
  tables sex * treatment / norow nopercent chisq;  
  
proc sort data=b;
  by treatment;
  
proc means n mean stderr t probt data=b;
  var CRP_change Chol_change;
  by treatment;

proc ttest ci=none data=b;
  paired Final_CRP * Base_CRP  Final_Chol * Base_Chol;
  by treatment;

  run;  quit;