From 82186a87f8e13c7d4189ecf5de18ee81708205dc Mon Sep 17 00:00:00 2001 From: Fernando Benitez Date: Fri, 4 Nov 2022 00:26:14 +0000 Subject: [PATCH] lastest updates --- aspics/kernels/ramp_ua.cl | 10 +++++----- 1 file changed, 5 insertions(+), 5 deletions(-) diff --git a/aspics/kernels/ramp_ua.cl b/aspics/kernels/ramp_ua.cl index 3566f29..397b713 100644 --- a/aspics/kernels/ramp_ua.cl +++ b/aspics/kernels/ramp_ua.cl @@ -195,17 +195,17 @@ float get_mortality_prob_for_age(ushort age, ushort sex, int origin, ushort cvd, //printf("Age_ID: %d, Origin: %d, Initial BMI: %f, Odds: %f, Probability: %f\n", age, originNew, new_bmi, oddBMI, probaOrigin); - //new_bmi = 26.7; - //new_bmi = 40.0; + //new_bmi = 26.7; //Scenario 2 + //new_bmi = 40.0; //Scenario 3 if (new_bmi <= 0.0){ // Negatives are missing values from SPC, treat as neutral case (odds of 1) oddBMI = 1.0; } else if (new_bmi < lower_new_bmi){ ///// if New_BMI is lower than threshold, then set to threshold ////// oddBMI = params->bmi_multipliers[originNew * 4] + params->bmi_multipliers[originNew * 4 + 1] * lower_new_bmi + params->bmi_multipliers[originNew * 4 + 2] * pown(lower_new_bmi,2) + params->bmi_multipliers[originNew * 4 + 3] * pown(lower_new_bmi,3); - //} else if (new_bmi > 26.7){ + //} else if (new_bmi > 26.7){ //Scenario 4 // oddBMI = params->bmi_multipliers[originNew * 4] + params->bmi_multipliers[originNew * 4 + 1] * 0.99 * new_bmi + params->bmi_multipliers[originNew * 4 + 2] * pown(0.99 * new_bmi,2) + params->bmi_multipliers[originNew * 4 + 3] * pown(0.99 * new_bmi,3); - //} else if (new_bmi > 27.7){ + //} else if (new_bmi > 27.7){ //Scenario 5 // oddBMI = params->bmi_multipliers[originNew * 4] + params->bmi_multipliers[originNew * 4 + 1] * (new_bmi - 1) + params->bmi_multipliers[originNew * 4 + 2] * pown(new_bmi -1,2) + params->bmi_multipliers[originNew * 4 + 3] * pown(new_bmi - 1,3); - //} else if (new_bmi > 50){ //oddBMI is calculated using 50 for [!!! BMI VARIABLE !!!]} + //} else if (new_bmi > 50){ //oddBMI is calculated using 50 for [!!! BMI VARIABLE !!!]} //Scenario 6 // oddBMI = params->bmi_multipliers[originNew * 4] + params->bmi_multipliers[originNew * 4 + 1] * 50 + params->bmi_multipliers[originNew * 4 + 2] * pown(50,2) + params->bmi_multipliers[originNew * 4 + 3] * pown(50,3); }else { oddBMI = params->bmi_multipliers[originNew * 4] + params->bmi_multipliers[originNew * 4 + 1] * new_bmi + params->bmi_multipliers[originNew * 4 + 2] * pown(new_bmi,2) + params->bmi_multipliers[originNew * 4 + 3] * pown(new_bmi,3);