Dekker1/dlmo-model
Archived
1
0
Fork 0
Code

Account for outliers within the groups based on consistency

Browse Source
This commit is contained in:
Jip J. Dekker 2022-01-10 13:04:58 +11:00
parent 193bcc65fa
commit d279777f80
No known key found for this signature in database
GPG Key ID: 517DF4A00618C9C3
1 changed files with 21 additions and 9 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

30
dlmo.mzn
View File

@ -7,9 +7,14 @@ float: DLMO_threshold = 13.043;
array[Patient] of opt float: dlmo; array[Patient] of opt float: dlmo;
array[Patient, Horizon] of opt float: sample; array[Patient, Horizon] of opt float: sample;
% consistency in terms of STDEV % consistency
array[Patient] of opt float: consistency; array[Patient] of opt float: consistency;
% grouping parameters
int: nr_groups;
float: success_rate;
array[_] of float: known_consistency = [deopt(consistency[pt]) | pt in Patient where occurs(consistency[pt])];
test has_dlmo(Patient: pt) = occurs(dlmo[pt]); test has_dlmo(Patient: pt) = occurs(dlmo[pt]);
% --- Derived data --- % --- Derived data ---
@ -24,34 +29,41 @@ array[Patient] of var Horizon: critical_sample = [
| pt in Patient]; | pt in Patient];
% --- Decisions --- % --- Decisions ---
enum Group = G(1..3); enum Group = G(1..nr_groups);
array[Group] of var 0.0..5.0: split; array[Group] of var min(known_consistency)..max(known_consistency): split;
array[Patient] of var Group: group; array[Patient] of var Group: group;
array[Group] of int: group_count :: output :: output_only = [count(pt in Patient)(fix(group[pt]) = grp) | grp in Group]; array[Group] of var int: group_count :: output = [grp: count(pt in Patient)(group[pt] = grp) | grp in Group];
% First round: % First round:
array[Group] of var Horizon: first_sample; % When to start sampling on the first round array[Group] of var Horizon: first_sample; % When to start sampling on the first round
array[Group] of var 2..length(Horizon): num_samples; % How many samples to take on the first round array[Group] of var 2..length(Horizon): num_samples; % How many samples to take on the first round
% --- Decision Consequences ---
% For each patient: will their DLMO be detected on the first round?
array[Patient] of var bool: first_detected = [
has_dlmo(pt) /\ (critical_sample[pt] >= first_sample[group[pt]] /\ critical_sample[pt] + 1 <= first_sample[group[pt]] + (num_samples[group[pt]] - 1))
| pt in Patient
];
% --- Constraint --- % --- Constraint ---
% Start first group with any positive % Start first group with any positive
constraint split[G(1)] = 0.0; constraint split[G(1)] = min(known_consistency);
constraint forall(grp in Group where grp != min(Group)) (split[grp] > split[enum_prev(grp)]); constraint forall(grp in Group where grp != min(Group)) (split[grp] > split[enum_prev(grp)]);
% Ensure consistency of patient is between split % Ensure consistency of patient is between split
constraint forall(pt in Patient where occurs(consistency[pt])) (deopt(consistency[pt]) >= split[group[pt]]); constraint forall(pt in Patient where occurs(consistency[pt])) (deopt(consistency[pt]) >= split[group[pt]]);
constraint forall(pt in Patient where occurs(consistency[pt])) (group[pt] == max(Group) \/ deopt(consistency[pt]) < split[enum_next(group[pt])]); constraint forall(pt in Patient where occurs(consistency[pt])) (group[pt] == max(Group) \/ deopt(consistency[pt]) < split[enum_next(group[pt])]);
% The DLMO should be detected in the first round % The DLMO should be found for patients with success_rate
constraint forall(pt in Patient where has_dlmo(pt)) ( constraint forall(grp in Group)(
critical_sample[pt] >= first_sample[group[pt]] /\ critical_sample[pt] + 1 <= first_sample[group[pt]] + (num_samples[group[pt]] - 1) sum(pt in Patient where has_dlmo(pt) /\ group[pt] == grp) (first_detected[pt]) >= (success_rate * group_count[grp])
); );
% --- Objective --- % --- Objective ---
solve minimize total_tests; solve minimize total_tests;
% Pareto Front objectives % Pareto Front objectives
var int: total_tests ::output = sum(grp in Group)(count(pt in Patient) (group[pt] == grp) * num_samples[grp]) var int: total_tests ::output = sum(grp in Group)(count(pt in Patient) (group[pt] == grp) * num_samples[grp]);
% var int: total_not_detected ::output = length([pt | pt in Patient where has_dlmo(pt)]) - count(first_detected); % var int: total_not_detected ::output = length([pt | pt in Patient where has_dlmo(pt)]) - count(first_detected);