src/Services/Health/CvdRiskPredictionLipids.php line 24

Open in your IDE?
  1. <?php
  3. namespace App\Services\Health;
  5. class CvdRiskPredictionLipids
  6. {
  8.     private $M_AGE = 3.06117;
  9.     private $M_SBP_NO_TRT = 1.93303;
  10.     private $M_SBP_TRT = 1.99881;
  11.     private $M_TCL = 1.1237;
  12.     private $M_HDL = -0.93263;
  13.     private $M_SMOKE = 0.65451;
  14.     private $M_DIAB = 0.57367;
  16.     private $F_AGE = 2.32888;
  17.     private $F_SBP_NO_TRT = 2.76157;
  18.     private $F_SBP_TRT = 2.82263;
  19.     private $F_TCL = 1.20904;
  20.     private $F_HDL = -0.70833;
  21.     private $F_SMOKE = 0.52873;
  22.     private $F_DIAB = 0.69154;
  24.     public function calcGenCvdRiskPredictionLipids(
  25.         int $iSex,
  26.         int $iAge,
  27.         int $iSysBloodPress,
  28.         bool $bTreatmentForHypertension,
  29.         bool $bSmoking,
  30.         bool $bDiabetes,
  31.         int $iHDL,
  32.         int $iTotalCholesterol,
  33.         bool $debug = false
  34.     ) {
  35.         if ($debug) {
  36.             echo 'iSex: ' . $iSex . '<br/>';
  37.             echo 'iAge: ' . $iAge . '<br/>';
  38.             echo 'iSysBloodPress: ' . $iSysBloodPress . '<br/>';
  39.             echo 'bTreatmentForHypertension: ' . (int)$bTreatmentForHypertension . '<br/>';
  40.             echo 'bSmoking: ' . (int)$bSmoking . '<br/>';
  41.             echo 'bDiabetes: ' . (int)$bDiabetes . '<br/>';
  42.             echo 'iHDL: ' . $iHDL . '<br/>';
  43.             echo 'iTotalCholesterol: ' . $iTotalCholesterol . '<br/><br/>';
  44.         }
  46.         $risk = new PredictedRisk($iSex, $bTreatmentForHypertension);
  48.         $sex = $iSex;
  49.         $age = log($iAge);
  50.         $sbp = log($iSysBloodPress);
  51.         $tcl = log($iTotalCholesterol);
  52.         $hdl = log($iHDL);
  53.         $trtbp = $bTreatmentForHypertension ? 1 : 0;
  54.         $smoke = $bSmoking ? 1 : 0;
  55.         $diab = $bDiabetes ? 1 : 0;
  57.         if ($debug) {
  58.             echo 'sex: ' . $sex . '<br />';
  59.             echo 'age: ' . $age . '<br />';
  60.             echo 'sbp: ' . $sbp . '<br />';
  61.             echo 'tcl: ' . $tcl . '<br />';
  62.             echo 'hdl: ' . $hdl . '<br />';
  63.             echo 'trtbp: ' . $trtbp . '<br />';
  64.             echo 'smoke: ' . $smoke . '<br />';
  65.             echo 'diab: ' . $diab . '<br /><br />';
  66.         }
  68.         if ($iSex == GenderType::FEMALE) {
  69.             $sumBetaX_No_SBP_Trt = $age * $this->F_AGE + $sbp * $this->F_SBP_NO_TRT + $tcl * $this->F_TCL + $hdl * $this->F_HDL + $smoke * $this->F_SMOKE + $diab * $this->F_DIAB;
  71.             $sumBetaX_SBP_With_Trt = $age * $this->F_AGE + $sbp * $this->F_SBP_TRT + $tcl * $this->F_TCL + $hdl * $this->F_HDL + $smoke * $this->F_SMOKE + $diab * $this->F_DIAB;
  73.             $riskScore_No_SBP_Trt = 1 - pow(0.95012, exp($sumBetaX_No_SBP_Trt - 26.1931));
  75.             $riskScore_With_SBP_Trt = 1 - pow(0.95012, exp($sumBetaX_SBP_With_Trt - 26.1931));
  76.         } elseif ($iSex == GenderType::MALE) {
  77.             $sumBetaX_No_SBP_Trt = $age * $this->M_AGE + $sbp * $this->M_SBP_NO_TRT + $tcl * $this->M_TCL + $hdl * $this->M_HDL + $smoke * $this->M_SMOKE + $diab * $this->M_DIAB;
  79.             $sumBetaX_SBP_With_Trt = $age * $this->M_AGE + $sbp * $this->M_SBP_TRT + $tcl * $this->M_TCL + $hdl * $this->M_HDL + $smoke * $this->M_SMOKE + $diab * $this->M_DIAB;
  81.             $riskScore_No_SBP_Trt = 1 - pow(0.88936, exp($sumBetaX_No_SBP_Trt - 23.9802));
  83.             $riskScore_With_SBP_Trt = 1 - pow(0.88936, exp($sumBetaX_SBP_With_Trt - 23.9802));
  84.         } else {
  85.             return null;  // ERROR invalid gender
  86.         }
  88.         if ($debug) {
  89.             echo 'sumBetaX_No_SBP_Trt: ' . $sumBetaX_No_SBP_Trt . '<br/>';
  90.             echo 'sumBetaX_SBP_With_Trt: ' . $sumBetaX_SBP_With_Trt . '<br/>';
  91.             echo 'riskScore_No_SBP_Trt: ' . $riskScore_No_SBP_Trt . '<br/>';
  92.             echo 'riskScore_With_SBP_Trt: ' . $riskScore_With_SBP_Trt . '<br/><br />';
  93.         }
  95.         $risk->riskScore_No_SBP_Trt = $riskScore_No_SBP_Trt;
  96.         $risk->riskScore_With_SBP_Trt = $riskScore_With_SBP_Trt;
  98.         if ($iSex == GenderType::FEMALE) {
  99.             $sumBetaX_No_SBP_Trt_OPTIMAL = $this->F_AGE * $age + $this->F_SBP_NO_TRT * log(110) + $this->F_TCL * log(160) + $this->F_HDL * log(60);
  101.             $sumBetaX_No_SBP_Trt_NORMAL = $this->F_AGE * $age + $this->F_SBP_NO_TRT * log(125) + $this->F_TCL * log(180) + $this->F_HDL * log(45);
  102.         } elseif ($iSex == GenderType::MALE) {
  103.             $sumBetaX_No_SBP_Trt_OPTIMAL = $this->M_AGE * $age + $this->M_SBP_NO_TRT * log(110) + $this->M_TCL * log(160) + $this->M_HDL * log(60);
  105.             $sumBetaX_No_SBP_Trt_NORMAL = $this->M_AGE * $age + $this->M_SBP_NO_TRT * log(125) + $this->M_TCL * log(180) + $this->M_HDL * log(45);
  106.         } else {
  107.             return null;  // ERROR invalid gender
  108.         }
  110.         if ($debug) {
  111.             echo 'sumBetaX_No_SBP_Trt_OPTIMAL: ' . $sumBetaX_No_SBP_Trt_OPTIMAL * 100 . '<br/>';
  112.             echo 'sumBetaX_No_SBP_Trt_NORMAL: ' . $sumBetaX_No_SBP_Trt_NORMAL * 100 . '<br/><br/>';
  113.         }
  115.         if ($iSex == GenderType::FEMALE) {
  116.             $optimalRisk = 1 - pow(0.95012, exp($sumBetaX_No_SBP_Trt_OPTIMAL - 26.1931));
  118.             $normalRisk = 1 - pow(0.95012, exp($sumBetaX_No_SBP_Trt_NORMAL - 26.1931));
  119.         } elseif ($iSex == GenderType::MALE) {
  120.             $optimalRisk = 1 - pow(0.88936, exp($sumBetaX_No_SBP_Trt_OPTIMAL - 23.9802));
  122.             $normalRisk = 1 - pow(0.88936, exp($sumBetaX_No_SBP_Trt_NORMAL - 23.9802));
  123.         } else {
  124.             return null;  // ERROR invalid gender
  125.         }
  127.         $risk->optimalRisk = $optimalRisk;
  128.         $risk->normalRisk = $normalRisk;
  130.         if ($iSex == GenderType::FEMALE) {
  131.             if ($bTreatmentForHypertension) {
  132.                 $consti_num = exp(-($this->F_SBP_TRT * $sbp + $this->F_TCL * $tcl + $this->F_HDL * $hdl + $this->F_SMOKE * $smoke + $this->F_DIAB * $diab - 26.1931) / 2.32888);
  133.                 $consti_denom = pow((-log(0.95012)), (1 / 2.32888));
  134.                 $consti = $consti_num * 1 / $consti_denom;
  135.                 $expo = 1 * (1 / 2.32888);
  136.                 $term = pow(-log(1 - $riskScore_With_SBP_Trt), $expo);
  137.             } else {
  138.                 $consti_num = exp(-($this->F_SBP_NO_TRT * $sbp + $this->F_TCL * $tcl + $this->F_HDL * $hdl + $this->F_SMOKE * $smoke + $this->F_DIAB * $diab - 26.1931) / 2.32888);
  139.                 $consti_denom = pow((-log(0.95012)), (1 / 2.32888));
  140.                 $consti = $consti_num * 1 / $consti_denom;
  141.                 $expo = 1 * (1 / 2.32888);
  142.                 $term = pow(-log(1 - $riskScore_No_SBP_Trt), $expo);
  143.             }
  145.             $heart_age = $consti * $term;
  146.         } elseif ($iSex == GenderType::MALE) {
  147.             if ($bTreatmentForHypertension) {
  148.                 $consti_num = exp(-($this->M_SBP_NO_TRT * $sbp + $this->M_TCL * $tcl + $this->M_HDL * $hdl + $this->M_SMOKE * $smoke + $this->M_DIAB * $diab - 23.9802) / 3.06117);
  149.                 $consti_denom = pow((-log(0.88936)), (1 / 3.06117));
  150.                 $consti = $consti_num * 1 / $consti_denom;
  151.                 $expo = 1 * (1 / 3.06117);
  152.                 $term = pow(-log(1 - $riskScore_With_SBP_Trt), $expo);
  153.             } else {
  154.                 $consti_num = exp(-($this->M_SBP_NO_TRT * $sbp + $this->M_TCL * $tcl + $this->M_HDL * $hdl + $this->M_SMOKE * $smoke + $this->M_DIAB * $diab - 23.9802) / 3.06117);
  155.                 $consti_denom = pow((-log(0.88936)), (1 / 3.06117));
  156.                 $consti = $consti_num * 1 / $consti_denom;
  157.                 $expo = 1 * (1 / 3.06117);
  158.                 $term = pow(-log(1 - $riskScore_No_SBP_Trt), $expo);
  159.             }
  161.             $heart_age = $consti * $term;
  162.         } else {
  163.             return null;  // ERROR invalid gender
  164.         }
  166.         if ($debug) {
  167.             echo 'consti_num: ' . $consti_num . '<br/>';
  168.             echo 'consti_denom: ' . $consti_denom . '<br/>';
  169.             echo 'consti: ' . $consti . '<br/>';
  170.             echo 'expo: ' . $expo . '<br/>';
  171.             echo 'term: ' . $term . '<br/>';
  172.         }
  174.         $risk->heart_age = round($heart_age);
  176.         return $risk->toArray();
  177.     }
  179. }