*****************************************; ** 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;