*****************************************;
** SAS Scoring Code for PROC HPGenselect;
*****************************************;

length I_hf_category $ 12;
label I_hf_category = 'Into: hf_category' ;
label U_hf_category = 'Unnormalized Into: hf_category' ;

label P_hf_category1 = 'Predicted: hf_category=1' ;
label P_hf_category2 = 'Predicted: hf_category=2' ;

drop _LMR_BAD;
_LMR_BAD=0;

*** Check interval variables for missing values;
if nmiss(male,index_dx_out,age,dx_defibrillator,hosp_chf,rx_ace,rx_antagonist,
        rx_bblocker,rx_digoxin,rx_loop_diuretic,rx_nitrates,rx_thiazide,
        dx_afib,dx_anemia,dx_cabg,dx_cardiomyopathy,dx_copd,dx_depression,
        dx_htn_nephropathy,dx_hyperlipidemia,dx_hypertension,dx_hypotension,
        dx_mi,dx_obesity,dx_oth_dysrhythmia,dx_psychosis,dx_rheumatic_heart,
        dx_sleep_apnea,dx_stable_angina,dx_valve_disorder,hf_systolic,
        hf_diastolic,hf_left,hf_unspecified) then do;
   _LMR_BAD=1;
   goto _SKIP_000;
end;

*** Compute Linear Predictors;
drop _LP0;
_LP0 = 0;

*** Effect: male;
_LP0 = _LP0 + (0.3236513639157) * male;
*** Effect: index_dx_out;
_LP0 = _LP0 + (-0.18719146708481) * index_dx_out;
*** Effect: age;
_LP0 = _LP0 + (-0.00574721311529) * age;
*** Effect: dx_defibrillator;
_LP0 = _LP0 + (0.27503224445615) * dx_defibrillator;
*** Effect: hosp_chf;
_LP0 = _LP0 + (0.34628912993253) * hosp_chf;
*** Effect: rx_ace;
_LP0 = _LP0 + (0.22174836881711) * rx_ace;
*** Effect: rx_antagonist;
_LP0 = _LP0 + (0.16600764344235) * rx_antagonist;
*** Effect: rx_bblocker;
_LP0 = _LP0 + (0.08725681577962) * rx_bblocker;
*** Effect: rx_digoxin;
_LP0 = _LP0 + (0.1632243359033) * rx_digoxin;
*** Effect: rx_loop_diuretic;
_LP0 = _LP0 + (0.08425128361259) * rx_loop_diuretic;
*** Effect: rx_nitrates;
_LP0 = _LP0 + (0.12922471804603) * rx_nitrates;
*** Effect: rx_thiazide;
_LP0 = _LP0 + (-0.16081892568955) * rx_thiazide;
*** Effect: dx_afib;
_LP0 = _LP0 + (-0.00226707256102) * dx_afib;
*** Effect: dx_anemia;
_LP0 = _LP0 + (-0.16535349893014) * dx_anemia;
*** Effect: dx_cabg;
_LP0 = _LP0 + (-0.04017489921599) * dx_cabg;
*** Effect: dx_cardiomyopathy;
_LP0 = _LP0 + (1.41511291168383) * dx_cardiomyopathy;
*** Effect: dx_copd;
_LP0 = _LP0 + (-0.03702295023012) * dx_copd;
*** Effect: dx_depression;
_LP0 = _LP0 + (-0.03382925166781) * dx_depression;
*** Effect: dx_htn_nephropathy;
_LP0 = _LP0 + (-0.03383047561884) * dx_htn_nephropathy;
*** Effect: dx_hyperlipidemia;
_LP0 = _LP0 + (-0.00180475696809) * dx_hyperlipidemia;
*** Effect: dx_hypertension;
_LP0 = _LP0 + (-0.09853867451809) * dx_hypertension;
*** Effect: dx_hypotension;
_LP0 = _LP0 + (-0.01728150097827) * dx_hypotension;
*** Effect: dx_mi;
_LP0 = _LP0 + (0.65177847636946) * dx_mi;
*** Effect: dx_obesity;
_LP0 = _LP0 + (-0.14195611374389) * dx_obesity;
*** Effect: dx_oth_dysrhythmia;
_LP0 = _LP0 + (0.11665163824424) * dx_oth_dysrhythmia;
*** Effect: dx_psychosis;
_LP0 = _LP0 + (-0.06819845439405) * dx_psychosis;
*** Effect: dx_rheumatic_heart;
_LP0 = _LP0 + (-0.07388907047358) * dx_rheumatic_heart;
*** Effect: dx_sleep_apnea;
_LP0 = _LP0 + (-0.03556032318895) * dx_sleep_apnea;
*** Effect: dx_stable_angina;
_LP0 = _LP0 + (-0.0156569460021) * dx_stable_angina;
*** Effect: dx_valve_disorder;
_LP0 = _LP0 + (-0.1636839757241) * dx_valve_disorder;
*** Effect: hf_systolic;
_LP0 = _LP0 + (0.75495402521349) * hf_systolic;
*** Effect: hf_diastolic;
_LP0 = _LP0 + (-0.95085645879015) * hf_diastolic;
*** Effect: hf_left;
_LP0 = _LP0 + (0.76641501222714) * hf_left;
*** Effect: hf_unspecified;
_LP0 = _LP0 + (-0.57722085981031) * hf_unspecified;

*** Predicted values;
drop _MAXP _IY _P0 _P1;
_TEMP = -1.37218706653502  + _LP0;
if (_TEMP < 0) then do;
   _TEMP = exp(_TEMP);
   _P0 = _TEMP / (1 + _TEMP);
end;
else _P0 = 1 / (1 + exp(-_TEMP));
_P1 = 1.0 - _P0;
P_hf_category1 = _P0;
_MAXP = _P0;
_IY = 1;
P_hf_category2 = _P1;
if (_P1 >  _MAXP + 1E-8) then do;
   _MAXP = _P1;
   _IY = 2;
end;
select( _IY );
   when (1) do;
      I_hf_category = '1' ;
      U_hf_category = 1;
   end;
   when (2) do;
      I_hf_category = '2' ;
      U_hf_category = 2;
   end;
   otherwise do;
      I_hf_category = '';
      U_hf_category = .;
   end;
end;
_SKIP_000:
if _LMR_BAD = 1 then do;
I_hf_category = '';
U_hf_category = .;
P_hf_category1 = .;
P_hf_category2 = .;
end;
drop _TEMP;