proc format;
	value $amount
		"Non-smoker" 			= "    Non-smoker"      
		"Light (1-5)" 			= "   Light (1-5)"  
		"Moderate (6-15)" 	= "  Moderate (6-15)"
		"Heavy (16-25)" 		= " Heavy (16-25)"
		"Very Heavy (> 25)" 	= "Very Heavy (> 25)"
	;
run;

proc format;
	value $smoked
		"Non-smoker" = "None    "
		missing = "Missing"
		other = "Not none"
	;
run;

proc freq data = sashelp.heart;
	tables smoking_status;
run;

data fram;
	set sashelp.heart;
	smokin = put(smoking_Status, $smoked.);
run;


ods html style = analysis2
	body = "blah.html"
	path = "C:\blah\" (url = none)
	gpath = "C:\blah\" (url = none);
ods graphics on / imagefmt = staticmap;
proc freq data = fram order = formatted;
	format smoking_Status $amount.;
	tables smoking_Status;
	tables weight_status * smokin * status;
run;
ods graphics off;
ods html close;

ods graphics on / reset;
proc freq data = sashelp.heart order = formatted;
	format smoking_status $amount.;
	tables smoking_status / chisq ;
run;
ods graphics off;

ods graphics on / reset;
proc freq data = sashelp.heart order = formatted ;
	format smoking_status $amount.;
	tables smoking_status / chisq plots(only) = deviationPlot;
run;
ods graphics off;


ods graphics on;
proc freq data = fram;
	tables weight_status * smokin * status /measures;
run;
ods graphics off;

ods graphics on;
ods trace on;
proc ttest data = fram plots(only)=Summary;
	class sex;
	var cholesterol;
run;
ods trace off;
ods graphics off;


ods graphics on;
ods exclude QQPlot;
proc ttest data = fram;
	class sex;
	var cholesterol;
run;
ods graphics off;

* this is from SAS OnlineDoc 9.2 Copyright © 2002-2010, SAS Institute Inc., Cary, NC, USA ;


data Myeloma;                                                       
      input Time VStatus LogBUN HGB Platelet Age LogWBC Frac
            LogPBM Protein SCalc;
      label Time='Survival Time'
            VStatus='0=Alive 1=Dead';
      datalines;
    1.25  1  2.2175   9.4  1  67  3.6628  1  1.9542  12  10
    1.25  1  1.9395  12.0  1  38  3.9868  1  1.9542  20  18
    2.00  1  1.5185   9.8  1  81  3.8751  1  2.0000   2  15
    2.00  1  1.7482  11.3  0  75  3.8062  1  1.2553   0  12
    2.00  1  1.3010   5.1  0  57  3.7243  1  2.0000   3   9
    3.00  1  1.5441   6.7  1  46  4.4757  0  1.9345  12  10
    5.00  1  2.2355  10.1  1  50  4.9542  1  1.6628   4   9
    5.00  1  1.6812   6.5  1  74  3.7324  0  1.7324   5   9
    6.00  1  1.3617   9.0  1  77  3.5441  0  1.4624   1   8
    6.00  1  2.1139  10.2  0  70  3.5441  1  1.3617   1   8
    6.00  1  1.1139   9.7  1  60  3.5185  1  1.3979   0  10
    6.00  1  1.4150  10.4  1  67  3.9294  1  1.6902   0   8
    7.00  1  1.9777   9.5  1  48  3.3617  1  1.5682   5  10
    7.00  1  1.0414   5.1  0  61  3.7324  1  2.0000   1  10
    7.00  1  1.1761  11.4  1  53  3.7243  1  1.5185   1  13
    9.00  1  1.7243   8.2  1  55  3.7993  1  1.7404   0  12
   11.00  1  1.1139  14.0  1  61  3.8808  1  1.2788   0  10
   11.00  1  1.2304  12.0  1  43  3.7709  1  1.1761   1   9
   11.00  1  1.3010  13.2  1  65  3.7993  1  1.8195   1  10
   11.00  1  1.5682   7.5  1  70  3.8865  0  1.6721   0  12
   11.00  1  1.0792   9.6  1  51  3.5051  1  1.9031   0   9
   13.00  1  0.7782   5.5  0  60  3.5798  1  1.3979   2  10
   14.00  1  1.3979  14.6  1  66  3.7243  1  1.2553   2  10
   15.00  1  1.6021  10.6  1  70  3.6902  1  1.4314   0  11
   16.00  1  1.3424   9.0  1  48  3.9345  1  2.0000   0  10
   16.00  1  1.3222   8.8  1  62  3.6990  1  0.6990  17  10
   17.00  1  1.2304  10.0  1  53  3.8808  1  1.4472   4   9
   17.00  1  1.5911  11.2  1  68  3.4314  0  1.6128   1  10
   18.00  1  1.4472   7.5  1  65  3.5682  0  0.9031   7   8
   19.00  1  1.0792  14.4  1  51  3.9191  1  2.0000   6  15
   19.00  1  1.2553   7.5  0  60  3.7924  1  1.9294   5   9
   24.00  1  1.3010  14.6  1  56  4.0899  1  0.4771   0   9
   25.00  1  1.0000  12.4  1  67  3.8195  1  1.6435   0  10
   26.00  1  1.2304  11.2  1  49  3.6021  1  2.0000  27  11
   32.00  1  1.3222  10.6  1  46  3.6990  1  1.6335   1   9
   35.00  1  1.1139   7.0  0  48  3.6532  1  1.1761   4  10
   37.00  1  1.6021  11.0  1  63  3.9542  0  1.2041   7   9
   41.00  1  1.0000  10.2  1  69  3.4771  1  1.4771   6  10
   41.00  1  1.1461   5.0  1  70  3.5185  1  1.3424   0   9
   51.00  1  1.5682   7.7  0  74  3.4150  1  1.0414   4  13
   52.00  1  1.0000  10.1  1  60  3.8573  1  1.6532   4  10
   54.00  1  1.2553   9.0  1  49  3.7243  1  1.6990   2  10
   58.00  1  1.2041  12.1  1  42  3.6990  1  1.5798  22  10
   66.00  1  1.4472   6.6  1  59  3.7853  1  1.8195   0   9
   67.00  1  1.3222  12.8  1  52  3.6435  1  1.0414   1  10
   88.00  1  1.1761  10.6  1  47  3.5563  0  1.7559  21   9
   89.00  1  1.3222  14.0  1  63  3.6532  1  1.6232   1   9
   92.00  1  1.4314  11.0  1  58  4.0755  1  1.4150   4  11
    4.00  0  1.9542  10.2  1  59  4.0453  0  0.7782  12  10
    4.00  0  1.9243  10.0  1  49  3.9590  0  1.6232   0  13
    7.00  0  1.1139  12.4  1  48  3.7993  1  1.8573   0  10
    7.00  0  1.5315  10.2  1  81  3.5911  0  1.8808   0  11
    8.00  0  1.0792   9.9  1  57  3.8325  1  1.6532   0   8
   12.00  0  1.1461  11.6  1  46  3.6435  0  1.1461   0   7
   11.00  0  1.6128  14.0  1  60  3.7324  1  1.8451   3   9
   12.00  0  1.3979   8.8  1  66  3.8388  1  1.3617   0   9
   13.00  0  1.6628   4.9  0  71  3.6435  0  1.7924   0   9
   16.00  0  1.1461  13.0  1  55  3.8573  0  0.9031   0   9
   19.00  0  1.3222  13.0  1  59  3.7709  1  2.0000   1  10
   19.00  0  1.3222  10.8  1  69  3.8808  1  1.5185   0  10
   28.00  0  1.2304   7.3  1  82  3.7482  1  1.6721   0   9
   41.00  0  1.7559  12.8  1  72  3.7243  1  1.4472   1   9
   53.00  0  1.1139  12.0  1  66  3.6128  1  2.0000   1  11
   57.00  0  1.2553  12.5  1  66  3.9685  0  1.9542   0  11
   77.00  0  1.0792  14.0  1  60  3.6812  0  0.9542   0  12
   ;

ods trace on;
ods graphics on;   
proc phreg data=Myeloma plots(cl)=(survival CUMHAZ);
	model Time*VStatus(0)=LogBUN HGB ;
run;
ods graphics off;
ods trace off;

proc format;
	value thing
	0 = " Paper"
	1 = "China";
run;

data plates;
	set myeloma;
	plate = put(platelet, thing.);
run;

ods graphics on;   
proc lifetest data=plates /*plots = (survival(atrisk=0 to 100 by 20 test) )*/;
	time Time*VStatus(0);
	strata Plate;
	survival plots = (survival(atrisk=0 to 100 by 20 test) );
run;
ods graphics off;


ods trace on;
ods graphics on  / imagefmt = staticmap tipmax=1100;
proc reg data = fram (where=(status="Dead" and sex="Male")) plots(only) = cooksd;
	model diastolic = systolic weight cholesterol;
run;
ods graphics off;
ods trace off;

title;

proc sgplot data = fram;
	histogram cholesterol;
run;

proc sgplot data = fram;
	histogram cholesterol;
	density cholesterol/ type = kernal;
	density cholesterol/ type = normal;
run;




proc sgplot data = fram;
	vbox systolic / category = weight_status;
run;



proc format;
	value $weight
	"Underweight" = "1. Underweighted"
	"Normal" = "2. Normal" 
	"Overweight" = "3. Overweight"
;
run;

proc sgplot data = fram;
	*detail on a reference line;
	refline 120 / axis = y name = "sys" legendlabel = "Normal Systolic BP";
	* draw the line;
	keylegend "sys"/location = inside position = topleft;

	* tell it to use ordered values on the X axis;
	format weight_status $weight.;
	xaxis discreteorder = formatted;

	* draw the box;
	vbox systolic / category = weight_status ;
run;

proc sgplot data = fram;
	scatter x = height y = weight /group = sex;
run;


proc sort data = fram;
	by sex; 
run;

proc sgplot data = fram;
	reg x = height y = weight /group = sex;
run;

proc sort data = fram;
	by sex smokin; 
run;


proc format;
	value $weight
	"Underweight" = "1. Underweighted"
	"Normal" = "2. Normal" 
	"Overweight" = "3. Overweight"
;
run;

data test;
	set sashelp.heart;
	w = put(weight_status, $weight.);
run;


proc sgpanel data = test;
	panelby sex w / noVarName; 
	scatter x = systolic y = diastolic ;
run;


proc sgpanel data = test;
	panelby sex w /layout = lattice columns= 2 rows = 3 noVarName; 
	scatter x = systolic y = diastolic ;
run;



proc sgpanel data = sashelp.heart;
	panelby sex weight_status /layout = lattice columns= 2 rows = 3 noVarName; 

	format weight_status $weight.;

	rowaxis grid DISCRETEORDER= FORMATTED ; 
	scatter x = systolic y = diastolic ;
run;



/*
proc sgplot data = sashelp.heart;	
run;
*/