/* SAS code to analyze data from Benign Breast 
   Disease 1-3 Matched Case-Control Study.
   This code is stored as benignbreast.sas */

/* These data come from  Hosmer and Lemeshow 
   (2000) Applied Logistic Regression: Second 
   Edition.  These data are copyrighted by 
   John Wiley & Sons Inc. and must be 
   acknowledged and used accordingly.


Variable  Description	Codes/Values	Name
1		Stratum			1 - 50			STRATUM
2		Observation 
        within Stratum	1 = Case,		OBS
						2 - 4 = Control 
3		Age at Interview  Years			AGMT
4		Final Diagnosis   1 = Case, 
						  0 = Control	FNDX
5		Highest Grade 
		in School		5 - 20			GRADE
6		Degree			0 = None		DEG
						1 = High School
						2 = Jr. College
						3 = College
						4 = Masters
						5 = Doctoral
7		Regular Medical 
			Check-ups	1 = Yes, 
						0 = No			CHK
8		Age at First 
          Pregnancy		 Years			AGP1
9		Age at Menarche	 Years			AGMN
10	No. of Stillbirths,  0 - 7			NLV
		Miscarriages etc. 
11		Number of Live 
          Births		 0 - 11			LIV
12		Weight of the 
           Subject 		Pounds			WT
		At Interview
13		Age at Last 
        Menstrual Period  Years			AGLP
14		Marital Status	1 = Married		MST
						2 = Divorced
						3 = Separated
						4 = Widowed
						5 = Never Married

*/

data set1;
  infile "c:\documents and settings\kkoehler.IASTATE\my documents\courses\st565\data\bbdm13.dat";
  input stratum obs agmt fndx grade deg chk 
        agp1 agmn nlv liv wt aglp mst;
		nvmr = 1; if(mst <5) then nvmr=0;
		if(fndx=0) then fndx=2;
  run;

proc phreg data=set1;
  model fndx = chk agmn wt nvmr / ties=discrete rl;
  strata stratum;
  run;