OPTIONS LINESIZE=75 PAGESIZE=54 NODATE PAGENO=1;

DATA CREDIT;
 B0=-1; B1=.5; B2=.05; B3=.25;
DO ID=1 TO 200;
 AGE=CEIL(1.25*NORMAL(18490382)+23);
 GPT=3.25+ .55*NORMAL(842098); GPT = CEIL(10*GPT)/10;
 HRS=20+7*NORMAL(84208); HRS=CEIL(HRS); IF HRS<10 THEN HRS=0; HRS=5*MOD(HRS,5);
 RAN=UNIFORM(240801);
 IF RAN<.15 THEN MAJOR='SCI'; IF RAN>.15 AND RAN<.35 THEN MAJOR='BUS';
 IF RAN>.35 AND RAN<.75 THEN MAJOR='SOC'; IF RAN>.75 THEN MAJOR='HUM';
 RAN=UNIFORM(8208091); IF RAN>.6 THEN SEX='MALE  ' ; ELSE SEX='FEMALE';

LOGIT = B0 + B1*GPT + B3*HRS + .5*NORMAL(428021);
IF MAJOR='SCI' OR MAJOR='BUS' THEN LOGIT=LOGIT+.75;
Y =(LOGIT/(1+LOGIT)); IF Y > .75 THEN RISK = 'GOOD'; ELSE RISK = 'BAD';
OUTPUT; END;
RUN;

DATA CREDIT; SET CREDIT;
   IF AGE<23 AND UNIFORM(843290)<.25 THEN DELETE;
   IF GPT>4.0 THEN DELETE;
   DROP B0-B3 RAN LOGIT Y;

PROC PRINT;
 VAR SEX MAJOR AGE GPT HRS RISK;
              /*
DATA CREDIT; SET CREDIT;
 IF SEX='MALE' THEN NSEX=1; ELSE NSEX=0;
 IF MAJOR = 'SCI' THEN DUM_MAJ1=1; ELSE DUM_MAJ1=0;
 IF MAJOR = 'SOC' THEN DUM_MAJ2=1; ELSE DUM_MAJ2=0;
 IF MAJOR= 'HUM' THEN DUM_MAJ3=1; ELSE DUM_MAJ3=0;


PROC LOGISTIC DATA=CREDIT;
 MODEL RISK = NSEX DUM_MAJ1 DUM_MAJ2 DUM_MAJ3 GPT AGE HRS/CTABLE;
 OUTPUT OUT=PDICTS PREDICTED = PHAT;
 RUN;

DATA; SET PDICTS;
IF PHAT>.5 THEN PREDICT='BAD '; ELSE PREDICT='GOOD';

IF PHAT=. THEN DELETE;

 PTAME=PHAT; PWILD=1-PHAT;

PROC PRINT;
RUN;
*/
RUN;