src/Services/Health/GenCardioLipids.php line 18

Open in your IDE?
  1. <?php
  3. namespace App\Services\Health;
  5. class GenCardioLipids
  6. {
  8.     protected $params;
  9.     protected $sex;
  10.     protected $age;
  11.     protected $sysBloodPresure;
  12.     protected $treatmentForHypertension;
  13.     protected $smoking;
  14.     protected $diabetes;
  15.     protected $hDL;
  16.     protected $totalCholesterol;
  18.     public function __construct(
  19.         ParametersInterface $params,
  20.         int $sex,
  21.         int $age,
  22.         int $sysBloodPresure,
  23.         bool $treatmentForHypertension,
  24.         bool $smoking,
  25.         bool $diabetes,
  26.         int $hDL,
  27.         int $totalCholesterol
  28.     ) {
  29.         $this->params = $params;
  30.         $this->sex = $sex;
  31.         $this->age = $age;
  32.         $this->sysBloodPresure = $sysBloodPresure;
  33.         $this->treatmentForHypertension = $treatmentForHypertension;
  34.         $this->smoking = $smoking;
  35.         $this->diabetes = $diabetes;
  36.         $this->hDL = $hDL;
  37.         $this->totalCholesterol = $totalCholesterol;
  38.     }
  40.     public function generateCvdRiskPredictionUsingLipids()
  41.     {
  42.         $sex = $this->sex;
  43.         $age = log($this->age);
  44.         $sbp = log($this->sysBloodPresure);
  45.         $tcl = log($this->totalCholesterol);
  46.         $hdl = log($this->hDL);
  47.         $trtbp = $this->treatmentForHypertension ? 1 : 0;
  48.         $smoke = $this->smoking ? 1 : 0;
  49.         $diab = $this->diabetes ? 1 : 0;
  51.         $sumBetaX_No_SBP_Trt = $age * $this->params->getAGE() + $sbp * $this->params->getSBPNOTRT() + $tcl * $this->params->getTCL() + $hdl * $this->params->getHDL() + $smoke * $this->params->getSMOKE() + $diab * $this->params->getDIAB();
  53.         $sumBetaX_SBP_With_Trt = $age * $this->params->getAGE() + $sbp * $this->params->getSBPTRT() + $tcl * $this->params->getTCL() + $hdl * $this->params->getHDL() + $smoke * $this->params->getSMOKE() + $diab * $this->params->getDIAB();
  55.         $riskScore_No_SBP_Trt = 1 - pow(0.88936, exp($sumBetaX_No_SBP_Trt - 23.9802));
  57.         $riskScore_With_SBP_Trt = 1 - pow(0.88936, exp($sumBetaX_SBP_With_Trt - 23.9802));
  59.         $sumBetaX_No_SBP_Trt_OPTIMAL = $this->params->getAGE() * $age + $this->params->getSBPNOTRT() * log(110) + $this->params->getTCL() * log(160) + $this->params->getHDL() * log(60);
  61.         $sumBetaX_No_SBP_Trt_NORMAL = $this->params->getAGE() * $age + $this->params->getSBPNOTRT() * log(125) + $this->params->getTCL() * log(180) + $this->params->getHDL() * log(45);
  63.         $optimalRisk = 1 - pow(0.88936, exp($sumBetaX_No_SBP_Trt_OPTIMAL - 23.9802));
  65.         $normalRisk = 1 - pow(0.88936, exp($sumBetaX_No_SBP_Trt_NORMAL - 23.9802));
  67.         $riskScore_No_SBP_Trt = $riskScore_No_SBP_Trt * 100;
  68.         $riskScore_With_SBP_Trt = $riskScore_With_SBP_Trt * 100;
  69.         $sumBetaX_No_SBP_Trt_OPTIMAL = $sumBetaX_No_SBP_Trt_OPTIMAL * 100;
  70.         $sumBetaX_No_SBP_Trt_NORMAL = $sumBetaX_No_SBP_Trt_NORMAL * 100;
  71.         $optimalRisk = $optimalRisk * 100;
  72.         $normalRisk = $normalRisk * 100;
  74.         return compact(
  75.             'sumBetaX_No_SBP_Trt',
  76.             'sumBetaX_SBP_With_Trt',
  77.             'riskScore_No_SBP_Trt',
  78.             'riskScore_With_SBP_Trt',
  79.             'sumBetaX_No_SBP_Trt_OPTIMAL',
  80.             'sumBetaX_No_SBP_Trt_NORMAL',
  81.             'optimalRisk',
  82.             'normalRisk'
  83.         );
  84.     }
  86. }