From 27d28f757afc8b900ef47f6c344b9699532a5fbf Mon Sep 17 00:00:00 2001
From: solliolli <olli.herrala@aalto.fi>
Date: Fri, 16 Aug 2024 10:09:57 +0300
Subject: [PATCH] Add CHD example comparison

---
 experiments/CHD.jl                         |  41 ++++
 experiments/CHD.pdf                        | Bin 0 -> 11632 bytes
 experiments/CHD_preventative_care_data.csv | 102 ++++++++++
 experiments/array_script_CHD.sh            |  10 +
 experiments/slurmjob_CHD.jl                | 209 +++++++++++++++++++++
 5 files changed, 362 insertions(+)
 create mode 100644 experiments/CHD.jl
 create mode 100644 experiments/CHD.pdf
 create mode 100644 experiments/CHD_preventative_care_data.csv
 create mode 100644 experiments/array_script_CHD.sh
 create mode 100644 experiments/slurmjob_CHD.jl

diff --git a/experiments/CHD.jl b/experiments/CHD.jl
new file mode 100644
index 00000000..5af67b95
--- /dev/null
+++ b/experiments/CHD.jl
@@ -0,0 +1,41 @@
+using DelimitedFiles, Plots, Statistics
+t1 = zeros(101)
+t2 = zeros(101)
+val1 = zeros(101)
+val2 = zeros(101)
+for i in 0:100
+    try
+        result = readdlm((@__DIR__)*"/results/CHD_"*string(i)*".csv", ',')
+        t1[i+1] = result[2]
+        val1[i+1] = result[1]
+    catch e
+        println("Missing file "*(@__DIR__)*"/results/CHD_"*string(i)*".csv")
+    end
+    try
+        result = readdlm((@__DIR__)*"/../../DecProg_0_1_0/experiments/results/CHD_"*string(i+1)*".csv", ',')
+        t2[i+1] = result[2]
+        val2[i+1] = result[1]
+        # println(val1[i+1]*1000/result[1])
+    catch e
+        t2[i+1] = 7200
+        val2[i+1] = Inf
+    end
+end
+
+#println(sort(t1))
+#println(sort(t2))
+
+# log-axis, all data
+# scatter(t1, t2, xlabel="v1.2.0", ylabel="v0.1.0", markershape=:x, label="Solution time (seconds)", xlim=[0.001,1800], ylim=[0.001,7200], xaxis=:log, yaxis=:log)
+# plot!([0.001,7200], [0.001,7200], ls=:dash, lc=:black, label=false)
+
+# log-axis, >1s
+# scatter(t1, t2, xlabel="v1.2.0", ylabel="v0.1.0", markershape=:x, label="Solution time (seconds)", xlim=[1,1800], ylim=[1,7200], xaxis=:log, yaxis=:log)
+# plot!([1,7200], [1,7200], ls=:dash, lc=:black, label=false)
+
+# "regular" axis
+scatter(t1[t2.<7200], t2[t2.<7200], xlabel="v1.2.0", ylabel="v0.1.0", markershape=:x, label="Solution time (seconds)", xlim=[0,900], ylim=[0,7200], legend=:right)
+scatter!(t1[t2.>=7200], t2[t2.>=7200], markershape=:x, label=false)
+plot!([0,7200], [0,7200], ls=:dash, lc=:black, label=false)
+
+Plots.pdf("CHD")
\ No newline at end of file
diff --git a/experiments/CHD.pdf b/experiments/CHD.pdf
new file mode 100644
index 0000000000000000000000000000000000000000..e815b8abfe35c71023cef286020b2601b37c90b3
GIT binary patch
literal 11632
zcmZX4bwHF|u(yD8r_v=Qy-O|KEg%iju!Qsi(kNY$(kvj|U4jTmcP$NqbhmW)Exz%+
z_ul<y_sp4@GxM7{bI$YZFsMn(vU79rU@)w&jZR>2({Rx^n!Uvk7RKO|f`Cn39U(N#
za!ML3e-CXO9i&ZN!8FX$f&eZ650?NpA2$~l03ZNoP=h#HxS50h5}@IwforC*fY`hP
zLs&3GMKQn*7LS?%|5-uP(+Nz&DP`(vYUgMLS72%dcA?>cGpJ*5s=B$_**L(*z(3>v
zui$a~kBNUN_=gp4qobRHD-G9wmV4CpKe_*`@((LF2B(sZg$vCa-p5r8|0Br%A3<fX
zgO#f_4L1OTQ`W}L6%65&(2!D4kOzB6gUuZ+z>i`su5h^QF}VJGEdP8SUH+T=sQ&+&
z@#r!H47d1ipVh%Gj&2ZhxYNSIG@P=I@EQ0IxN=Sz4_7%2*GJ?ZBO|8);QhlQB|*c*
z!37fF0)hZEoaUzg89hop^8L%m{fB}3QQhO5nyD+?5r;>8k3(Hmv$tUL$MqhEkCazn
z3ma2OM-Q4eTyW$0xVUI|xVQ`-fpc($+w4Na_Xyb|>)(-m9RK6vKZtSh1OEk$<_)K;
z9o&V#F#Ch`-{^Qi4-AK+fkF;HE;T$df^&aoF`>m8#><dMy!$>xvyH=rKIVBJwiT;B
zchG7+SdWGB#W*X@Q(BRTrC771?RV&Vbdd4Gof!24-S68YvG;f9orW1-e~KAXceDSp
zk^g;t_eJLY_0`Pa`+NHxS9+20tcSz3EYXJtA790$9*DvY*<I{txuLJll?~Tvmg@aU
z(AduM!~H79Z*HQS-t4jW(Isaa?~9AJL9V2~Z@1omDcy?RHqHd2_Z#M>^I#ZnDfn@%
z_vaqFWu{-;+>9HEd9+mDy<S^%>!1qHd3SJv9Nl|x`*Lt?+Bb2-dEsGT)j?N46a{mq
zNBng%>&#@H=P55K+4~p-{4?Ih!Am)X83}w?eROq4TW5W!)@r3_>PN!P)<kyKV*m5w
zqu1({+4>6KrmHMMK*O_FjZ;$6*v#h5{;%pqXV2`Xr)(vk7tL%EV>i$W`St3DGU{F3
zF|RCZv&Ma$E&@o|mbz61zkAhR*12pQC#n>$3jnY3@|mwKYvbbcei)j>wktXRo_Q40
zX7KRqt$XP6N;juoTjtgNTu?sg>J)aT>;X=a)wsM#^aQbpG0Aagjj1+y4|aS!@>NjH
zYZ<MxPLngw2Ic2ZyN3bq-6X_Em|bRqvM^HB8*NL~3#yF2KE=DX!9v`#-a=@WBGf6T
z%T87|IX6Pk!t_hR7CGtER!bUd=)}G4Zs(aXv5mEToyOv|wvWg<9HoHZ7XEr=E9kX>
zu8DY->T`U11f<`j=goN-S9DzJAcX4%wp#Ahde;+59p$DED|iKvZnDnE4eiZtD8gjX
zN)9uook@&57g-h<ac@fIAQB@BOYzfF=BVT<>Gz4d8fhb?4rX91x()rW%{oEzt;Sdu
z)sRGnW;|3mB5JUP1(O=Z&i1Z%Ah_!~*3t>9`ps?lGahFYObJAd$!!bq&<2Svo=-5R
zLg>dJ{Ve4ESmB=NB-3@}CZfF1L=FF`!r0~_=#GL){fhPSC67*S{#~8m-skWAL!YN&
zw4=x>u!eTg@z*s+R<j_xS2`|McGeTOB;;gff&Jm$#A}*d>B+IpY7V;!Dzz&n%LP)J
z=%g@ObqV9Q=4LeIAyl3DBPuj0?9xRUa!qs9KkHN~$#N!X)36*^rCJOjl=o%Or?sMz
zQ^638(E<b{K8pR@%uSWh?b3V69v@G>UmpNzbLbdlUcvjq1MH#OgQ!CWE@9(D^alEK
zYFFFxJM)XKfB+kgY(b{=q_qGT#j442lGF>GSj?QR4Nmc7BF0cgDt+E!y`S=o%H$|5
zFFt66JS~F6?@<i2_p8j@4AI7y&a7I^Io*53+V;&3%<=fz`UaI?t|w*w<RQnoP|%b6
zRRmDIO5*~uV!oEaLB}b}9txr>8@8>%$J<kuR#(sk>#b=iFe^!l8Z?FnlN<!rr!fwp
zH0fJ*JY}+s<Zh2RBpSQ<X{q;Vky{?3hS`l`VTV*Zg_fa<3wl}<n*d9ve#2Qz9r+d(
z(;d>n7}U?L@I3~*h?<;27Ira^VBPQx)6RNS;9l2+y|ygdx%9mA`Qgy?2?Unm{lICP
z!i8jZ`ugnRLqF;eC(js~vQc0xLaiokKyts0iH6F#^EChca*Q*f6>0B6!kVt8N`M}9
z_O4`8F_XY>NTBaB%m6#=6lQiJB>+pSXEGi?M3Y3ZJP+D;9e5Q=Qz#gjW*g0;Z#hzi
zOo3=~O7&e2gl-8(M@az1VdD{<1CaEwpmUzoQ!AID@f4&N?$+DUafKP;T)d1HTg6Z}
zHc&lGZ`m~ruOfe>rs%cS2HlMwC>O<Q`k$10%$Q4_849|tDL=#;eX|0tKKm#e6K&PL
zk=(_7UC7(Q(xy$kT$Y?HoN-6GNJn{xUA(8Yp6FAwG>yNW7S1_UIS$f7WU8=>j%_<Q
z`Zo!y#J5hC4Q{7Qy?_wdWvKW(a3bJ?LuvEt=@uLx7V;0=JOT;GhJDxw2<Y^*&)CnN
zhtaQv8yxCRfi_%^V_y3@^rb(Y?0&{duE@my>_d~bUUf7UfD9%cuXprv5P%s>L8s(A
z0LqZr9b^rP%;XrU`P@WEneawBKz56CESh)QVvY@h9Hm(^ufygg@9U7yTvtE*&6Ev7
zs!n%<1`fO~5DJFvNYjDsWSxo8!yH!Fs!A1T-{$a|4CiqT@^r4(uu--#G08%w3OiMK
z${RjNcsDTc^VokEARcm{&Et4!<@K^~9`5z%?>Dvm*mrGOcEcsSG9muRWd8jHJ1CaD
z1auy$P+;g2rohZ;D+KPw3UV6ePOAQdJsevG$jho5q7sgTbguxIoxtp+$new$W-)dW
zu^?c1z1|_!*N$K2$hBUQPdO87+Y6wZy*x8_9(BsYwQST##qD?>HNM*!Nm6QDKmU>l
z>skHYDiXR)>bi+1(m#^KXW5G~YPT&4@PDO&HAF%;`53?1yszMioiay3X-;bq_R{Ts
z7M*Hs_5YiLE@8B!i>_fWH5NATUN!en1i%xZ2nPXPF4iYibX<mMWi^6GX9ZF2;xXkN
z1{ApV=rs-h6aO5{8k>!3*=sVc{zL4Qtq(Fc1mP-b$+4oCy<HA$sMC(6cMT8_NiTt^
zba*&Uy^q5Vh;h~2GUC1@yP$UfXB!qPXLi@6D%VZ{TO`4)Qc78kk-Nnfh)-yLDfmNz
z9tW+xu}?WivQhqhleuwwNFKGnDpgrmY&Or&LW$CkZQo=rU>VI|DWTzpM2e)sap_u)
zx)b)aVHmi0y*ik(iT$wJ>Rv=Q084%ceGY<i`r9iAk0lb|*Sl_9Kw6V3D4s-|n!w=s
zg$|RE1sHpqU1s*c{6fW`U9Q|Q)*I(ds<{@jdWf$Zs|yWjNZ-Y`*;txTWR1WHEk=%I
zYLwxa^+Fkw?4~%Q$+8H1M9KnJ^=}BgjiWU$T!70>U)iGA*<n>C>T2s-rEK)p{W!tT
z@G_F3M32KYz6aGX(*_MEJORI?^=9cTM`Ws*P&Pm1_EMN0M&+EvWEiLTzyWPxJ1@W_
zOxrtdY}}@NrlKX(Qhi958ld(eH}^9wVW{?Zj@W>{cH6N?p%myQ)5d-ZS*_q5JA>gS
zrv60!2ulJbMMMxr86da>ZKELmGY=9@k8gVUug{!!u97SJv4c(@DJ<M&Rx5@$(con#
z8XUV<^Xbe!VX@EDUhA6zUyPLkK@1e&!IkfdVOdL{jveuexW8zZ`jz*7kSZnk>W%8&
z@04y__cHJ}vtKGIZY0hw{98`eQfNE(>YE0dV`gZPk@lc>FRkTm{QBkQy;Ra%v&Ho{
z=ijm526jC^$N4r;txAn*dWn}x0)i^f163rcXH5VL7atM<<MoaEL^%>&mt|r#_LCaW
zzD5;QFXb+j;e)M0PeXP8&y)NbVsdPUL^Y&=6bFmkX;kml;bz89P|oWg&qt(5cRcP?
z8ejf=A>&v6{!O{Ggk;OW7roN`Q};PjOOZQrb~=PW29E?WurL)~Y}ETVBGQu|00rc;
zntOpF>1URn?LV()#RIqqxZu|cJ;b}MP>L6QV%ZP6+E3}PW!_*%&4t9gDqx@&^qPAa
zi{ByFytLpHhmP)-2IZJ~gK7<3H;D#9B}NLAKLgP_J~N41C}st0lGLTB3S)hc!^LRr
zwm@Ii-?UiA{F(a@4GJAI`BX@FERiXxd}iW+{vl(*vLn^vxH<ys`BjTGVi6;Y;UL{K
zVc-iOyQxan7#d1dT<t*1;k5vt#;_fRef@`Ayy_g!NIuiC0?uYs`XA~TUy;9Q9_W2{
z4HmK5)#L_zj?!NZ%8eeuSxqQ6BN%vHAx=S835Y9VWQx(@I(Wf^%WLiX^i)?iBw8|@
z!Lz@g(AyD;QYRWbDYK!hx}ixPec`L>GXCk5Ms=GABlU*(GK*W{l|X;CiFRexz(wqU
z0@JYqp+IjK;cp$6+Gv(Cjq)<HXe)@1SrHVa%oF`n#!Gw<Oan9pin%xt>0{=b{E5c`
zLwy2{(8tBjV9?NI!W$w{|M$$PZQa4sstHKkq3ByryAE}NmrXRSX9KT9`iuD}4&e8b
zG^skN+o<BALDbM~w6I;gA35};I`W~P;YIA-F#=GB(u^m#3&aka`qVGe^Y*dq{zrac
zp;3t>0+jvn7lQZ*v89v?S$4+fb{@X45e?9u?uJF7t?wot^trBHfqRMp_iwxr0Mn3X
zV)0|yj@Pgkn1`2bt<axX60}UXwb`J`F%>ngQk{Q=)0=kjpe(wi8mkFsSb6K3Cg=Di
z9Kb|%RY(c8-0g+BZ0u4zp)7J}z}Fx#kpJcv(Dw<UFlr>PR5AMT{t#TMNjyIl>l%T(
zW26rE>6oL>wb*pi-l*rjr(s(N$`f|(iiZd2u2cnf(4Ro0mlp4^T=fT`{v7u0d^$vs
zXWAP9Blrn0^Jm(7XJnMI2*dy(6_BOM2G8en(jX4TLN79=zZ^L5neN9bik-v9!L!;V
znx5oDf#2{_UZg>@rUGFr*NM}pCZuBnv6O`rz(NRF)utAjp|Z+3kIJl*`7*hlo8VL6
zXGzQ6sj)ujH2iWHOjL8GIU=YHn%mLO0!eSs8Yb;Tf9OA1tbT+(?fQeQ#>g*Cm3gPS
zBkwd@oN=6k-u*8Wh^a6NNq?=%Se$IgCf>nMB<-w6y@|9!93TsoH8^rC@B??o@+{yJ
zJnG)`r7EE&D4Koy0_16_8baScytRwPb=6W4nid$}`CAk$)L}S%`a@vmwS_x#?w`~8
zZ-=%1r>s^bUWOIKPrJ>heOCiMRR_o%%7SG2CqOaE{>l2h$9UD0^uGQUhVAS*6vO~B
zbzL-c7iQqt8f!A{n&Qf&b~pUVY<JTd+!tKJ>H<L#4Y%@rYA}AON8+shZeZtiWI^Et
zLfK!J>%Zdvg$KvOWK4*k`$@HJq=nISSjQQj<vM+e^mGrPwVbZd0L{cy7d$(Pr_w@A
z+xrYN^NAL;e#b|Ypl<5;dB+~8w?~fuTbxGUx+n4lgN46Rd+<#VK$TWuA_`$voE^=w
zl%;GEkx2ygSPPWz-^5?}?dcN{SkU4X9;mgU1KS%^5#Sr!<M4e>J5a*-d0g(w#Fw|;
z8|IX))wINTRpl_YNV`R0`(BN*GEVEVLmzW`Novx$1O?hc(eaTIh&~Z%C00>BMrhOs
zA6=fp-Fz4qwUzJYZ|D-UvY>+|$IcgoLd4Qh8rRoNZUrRJx1HmyM>ZV#(*?>HpmPb-
z@+dS@w;i+g5+;knQYa`0QbJ$cv?-smq^q|9vQAv9lLrwGJ!YEbW3sp<<@!1*rxr*#
zdF94xWKE5w(#OEJseYK6Pz<EjyuUbO<Ocq%lNnNS(~~Q!)y7Bnu~yH?!Ir+mbv%@@
zfC)4elKXhSb_siDEvDj|B;`uAnj^(rwF|$yGV@dBQ2MkF!ME|Ro`H#17Atjq0`QwP
zMlt3_LFW2nx%<WqPQ&@W6D6tQ3$GCRc^M|+yxM*Lqe8&#eykAHFZS)|%^n}(;Ll*G
z?8XcWVJSGKYLt@HT6mCMAbA>$?%*kswg|0|>PDG);-?IU%@V0?Bs})cNvnQO6m<x}
z&<&CLztwF45*C4dez<|qot#a>a(>py)PK`zCAtigY;i7T7HApa_tSUM{rv1<(j;Ew
zZZuu7;R3N{mTH|n>g5mlggjd38tyNBm)a%5>3BpQIR#`_562$Bsexli7txuAk|Zg%
zh3I=O-3`6Oel5{fph`=***R}3ps3_Yas12t-N0k(*{}zX^O+&du3k@oCkdbM<_ch8
zW$csI!gbBw*s)rYlfh`$`nD#IERtn`B$cds6OBew_m1;eFJ-$XvV5W@i&hrrTPmyt
z`17Rs0TpgsE2A=?YUqbwd%HXQZCK-C43kR$g=+k#0gS*;=-XS{s(O;7FJ~{~U917<
zsse0vW4cr;$>8kIimW@7LyL_XpAKn=Z(61tM<p9xw~yz~wE{YEL+6iZ{O&wiZ+niT
zWoAyK^-$ERpZRbzY^FqPTb@e=fQafpc20J=qJ0|wAPVH4y|Nx^Fn5|95^!Fe+J)s-
zm0)nve*W9vlqfWy2ENYCb5LM;<%&o_rf-@8S7=U>b4nUKjsP+)AhFr)%z91;GDMa_
zpBIQktv-XK&Y-0@g}O3cmA>M;?<({XX!i5B;*4=5xeFS%g>*?)p+Ue>eKx8l-^Py$
zAGgB3#;-_xl|1rF0>>W;pBqzSGu5s~&D`Wk;%?kG_z21P?Y526oMtMJbkHqO@*=G~
z*>b@UNx1xa%AQFlcD{TP8-gSzLvWZ}cPW0_`m2`XR$esN@S18Q*7)WKn-mnN1y;vM
zrnXit;Fj^70Ke}m{KCHcIO)g}*5ct$uc^{w5;}-2dAa}G??u*XMGN&=7Z5E8nf%qN
zTw{+9i+`=1F{a6i`UYtR6>z;_bZi2IH5ic*;_T7b?cR%umFbRdLOheHwf79i6Z>s?
zXD>`3iCqYr$O%=xAKLDjs>#rFst^Np0ABfq#HyI3Xxs68EE6_p;h|8tyNL}>5vN%7
zQ7P`lsy_e*IpK8HnB|(9_4?wPN~4WZULy*lRKDj^c<6cENlsx-0qw4~n<G^R{EnM3
zYc5N$u!$na@(B;fa@{r-fcOQ{|5^Ku6l>p^IA>+fU>Wyu1s;0H&Yb6#BgLuWHD529
zp_Mb+8$KAws59JEslTWF+;iq3ivE&cKwmJ*e69v2jx0CpHO4ZjCU4^DhjWpVUGE#B
zK5D9~y^qE2IcQQ>(W$(Psgs3Z{_;FeNo;d_13O}0B}YKFhPUcC2fUr}jPF|ClU|YN
zu1rd79ZjU|`}nCI1?rmyNno_L2gq9ue}sa2<W!iPfZC-`drwqO=&w}~_)V|zuCHFn
z`#Y~qeu}v7Sh*#q$cFa7-Hj0UF>GG*;xHmWUwkSibn$-0$<Zm#i>kK{T<E{<AoPBP
zP(Jm`&pi1oP6{~4{LgO`ur4cKwboTRd59=hbalPy%GOzByXr^3=cEu$`O=Gj{qcw>
z8a$|Su}VpS&5xtd?uZDEJem9Ur5B6+0C>D}bqf!0IkC>_u%d*~ZtPV2Qls8ev9A-J
zLn{kH6bC$yvO}bU4cwk@1t1-fN?^gweYZVBlqM0oWC@(l(C<}cv*92G$GUcBk|jJ^
zT`D&bq0y_T+#q*>8ykwjr1F%&^ni6@76-w19Lp1|<jHMLT@9nOa)6Ol-G+T0wVwO@
zsgrVBu_P8mLgPF?#qv<B{Fd)yBk(BLdEywwCx5F8d!?7kiq<oAGI}v<9<6@?a`M@F
z0k><~jAtP3dx2_?PdBgVeXyK&4Cy>3?$#lxhs)9})_n%5k5RPXsV>zL$}dPuMq$+e
z%;pz&Xiz4<!2F*uZ3;y@#(4?N=Ow3DQ0aI%WS*p$vprmZ)c<4$aRT?75)wWDqsDn+
z-&t&S3jy%(d4y1CcIz`2C&HK1jZPXK;-bi&eqxtQQkb_`^kJDDa4d$joBE}C=Xu91
zd`=5Rg@k;ExHgNvV#QmiJbKiVc9>WG6~;x|h%acAezlMGyML@hv=a@#kc*AALAWQf
zjd%#0di9jgsOZeVQN+UO7Hd&V)ZbLU+HEs{f1JT+X`8ZQiWHqow0zPWrE&D<y#f|2
z;s<Kyy}HQdkg)3T)ZJahf_aGgB`HKW)Kst94GlT@1iLLtcy<B-cD6+U-__4lY!`vY
z3^SLcNyKlwyH%F5*rMF&2aV=F5I5?Id{w;-lL)zG3~YG~!GLejUh|&3KC5dB=vFQa
zF-!-XwGWTbsfBlPc7e-txB~3Jz+=)kV#Y<L?XRW0=cA|r>vt$pLXT^E5>$aw$-w(g
zec&;b0@>>wk@!%orbFPB$K%dtWnB?g3vV$IufQPCMHWrRr{Z3JK#>yrV|CwCmbRaV
zPZV5J@Q7N{uJ;jQ<gnJ-7O)tyI^h^A&gdCM#|5J|@B}<#j*KpJm!SfzU;QmR$Hmfz
z^#L^{R1dxQ^qPJQm-LhswdsiT=C4Y2a#gixZQGK=Zdz)VnwBXldl4(KG`Iwg<-Gor
z%g}2}1EtO%63f@TDbM(xQv<*FqO=bCB1(B%MFLl4=f98RH>~Fiw~&ox;rQT^hD0x5
zu|OY<B2!dp^@}~(0jCBN`3C~T$|t2lpcka;!7v~kr*ua^UzGWe0$G`(Aum73c`WTW
zfJ}=y+73ksxO?PQ{mVDzJxiO?fB{UQ)|mA(%&+y0ry^Oj#!K^i=d(4>qvo#=N%(kV
zbsbd;`$DmbkAQtlA|od7j4F>i9*!Lc5x+A0a9m$#mCS^Tw*F+05kciP^Ly3DSx+Qj
zU*w5`v43nm^WmiT4MT+7=}vGmm&3*zkd$cK8KPXb)-h?aG*;m8c=|`&O_9PeDbF@m
z&gB6i3HVvGM*D?RWiy@&%@{a*e%&Zf)gSXPpOEa>_>Hse*YXM8*Fthf2_1^l0nQeg
zdwC<Ao`ZgV_mI1Z_k|KQMr~>5gH(_<#-(mQ^=#Br7Q4Yq;Ck|dx5Sj`+eyU099A!@
zIhe#jVasx2vRpq?pX$9k*;7qO{nMpSp)VozU9h7$vG~h{(|e-pho0gkvJ;W4Ayg;i
z<&IbP*GwMp1GB=Xys+<*Ka*O!rLSY*%tWnM%THD8vlsJ``-ZKEOl#YTN15QasP5d*
z$fcgutix=%Q$L;70*&>;2O;d$jqx~xMrI{MQF+BLQAY~FRQ+EKuj7Ctu~|1q2JMq4
z^+#RqXd5(G9Lqj8>|LR!4mzzp8tV(>j^YDC1vvbI3RKC@)T&!BMn{mS2H>K*s&QF2
zzkd@ORvL}F*_WUDt3S0Z8a|8j?HHc9qg#GPwSsxr)>?5w_gcx3u-+v+$>wWRL1-ak
zZDzUGw{bsuues?g>3f5mOLF#&AkG~!mw~JG-My#_^6|ljhZk?R2r_=!kh)7&9QH*g
zd$#S{wBda_q!A07Z<0RZPgCP-Op)RL&dERTD9yQW6KGuoc{z?AMdBRf!F^_aLLPO2
zIlkAhMA+|MIr2i#)PpY+bF>h52J!CZ_%3++{I#3v<iw*Q{bP^0huGwkiQhd5YtM3K
ztS~O6TQikhC3BGtN?lBPys!RM73{?^PO1_8@hGRVW#6B_fFdTG&l1iKcY#p2+$7S)
zP3}9&c1Yp9<gezyMM5oi4!A*lBt`j@_T}xB{PUdBl?z^~nrBVYRy<^4>2-GY?Qq>G
z)}7NrIdLJX41(qdq#G_{p)HI4{L`6m<w+mirN6;#k-gZzDw|T|CiS5A7G<25_i~pu
zkM?>wP8#0cKknbnMqM)(0Y}Ha>qsakeuP0Vd)(l=&iFKx=ZWKfp$OyrEXr5<(?f$C
z;B@A_C4heU>3mw@y=+~zMGJkrZQJMzxPVvNPW0?qCY=1Mx!hg4zuP5z$+2$#D!8|Q
z(RtwHT=q6<8~o-dyk_uESWbtbxq(+zAmR74UOyxb!vlnpw}U+OyRn@@DJd868o5<(
z;UMo!;pX_kL73C1=sO4{qHo;4$k@Mz$3@{NQ{Le~<kR*?V*fVJn-6=)4SclkY3Eb=
zmht<g@`2$abE-=23nn)1Vh8x+(|o$k3n)T~+aFiLt;SHPdo+$SR<2X5%{6}z$ZWy=
zrN@8a6XzzQ<P+q>gfGlEAGUlceE<&%8T(33M7MCkuv;$-o^C&;y>s?QAK8x=W?1!Q
ze@00kr&!bbEfe%(_Y@Lj;K2wu{Pc#<-Fxl2>WssTaAp1^o7FzV!bkOCmeQWp`$4j<
z8sw}x3D?4-|8;#qHRxvWgfJV&80FsCjP;ZAQ1P$4U-xn{zi&sW-&4K?fo?a}t=*h4
zeaR-?i;6jn8;d@ipBhtFpP)QM-ygl+lE<$nf8e2zTpgyrxefpQ%bohS>+gG*@kj5*
zk~i0FD<6M+J-*U}+F^G&p)b`YuAvcuHP;45g*uI_oN<k|eKHxaXPh6{lECBCexKjH
zi><|-EC@S_ewsC>LCEBn!XVQ!%OTLM6p^BZHgFYk+kSdJ<!{@Z_hIxY3t2t_|N5x<
z;(7TpAZcPi6^#kH#KW_H);nQR!-1GzHWX(pdVzc%J;o3<@rqS$F=rkRcs!=B_6+mn
z@Pwa@14ma>y3)GN);d&uW?kc$?hR0s=M{O*+GSSc^Os&h4q^I9&5Eq8jQhXLYXVMe
zGr@`J-})(L1tZ9h<T@&^yKQIsx<}m_-}0Uy&E|F#qsCrLbZvLaGbVX)qwppUN*k_@
zXD2v9!oFeKZ<Ob65ROwN6Dxg7WCtA_v2{P29V^7X|2apsgpg=S;{9vKsW?1mk-&{P
z_%e(w1mT^N-ls0ynDjAeTpJ;}w>wNEn&h!JTn;mM(A|x9!onvqM{-gneb_Q)XdIGw
zWg~C8zSk@+qg&<Qf|I;<a>Ma{Ywv7BjhIW9jS5X=S_mjvhsu)W6`hOw+3fFdygu=p
z=GRr%tsrB`lp1w^;p(wS9<ktyjK1*Wyrn;?2zBL!T8w_BdLMOkhTe!pAf<t8HWOVl
zC<)km0=%`qlDW70n#3ru6+flg$O;@U(Ww?2;K=B=BcVL_;X-;n%i)X1uTens{2tKR
zY#L5UnbOwQdvVK79Ip(v_{{VIv6~>4wWN>0*3YDnGet#sUc+dKdDEV1#bAka5-&l{
zCXu;PUNL?BUa!QI@LEz)Na^I^QaX&c)O9zJm*L0Cl;+F#F{P%93^8~|pOs9Cb18MR
zQdJpzuq-?mtpoj-hiBx|36<)D#~B+BUTtn@7`uB^Dgm`E!@wzQ&of_r8XpsLinV@d
z@_QdjnW<07l@s{6B$eP5ualJ%*+HrYXIRdW-m{Yu?{|^+eyd>=Do{w3+Ou&FBqqU0
zYgl-*C2IE|&Zx~#k<5e7aK}<PH5VK9EN6~y&&mQqJA4JT6<UV2j*d=leCzAFpP2cO
zh+{GJe+*q-UgJlUM{l%IRdEeYec_?^Ly@+d8kR|m*%rZjP`t0~Jx~^29<)%Wv&wyt
z=$d}9j#8R<K~(Mq<A%>pfT&<gQQwva{m9!HC}BqZkhJ~dJ+|a0dI8)#=QPvnc{PI_
zuFBq)ZOGZHpqKTN+%740hZCJWW!BTSv((_^k)>5R(b?a!^PiH<Q{Bj0v`>DW?t+HD
zbuOKd)P$cDdr9a$L*B0N4;1_vwhn(q0~?f-+pd{zobcz8>5Ae}pMODY+E0@o0sWYy
zC3ybwoKKe}GZ(gg<jl2bA2JJc?Z0zt2PopuamdBzyXX&dXq)xq<x5Q4bWG}&uz%!h
zvOgN;*vnZPU<nznBXK}VHI1ylxD^t1i7?$HnnPs{g-$q?S8ODNJj>Bm8E>{!;kNr>
z)|Pb`BonnkbnPH3x=u7#=>RNzUIWeXm>d6+*3>(Vb8+i!2ihnwyhL=0RKFZx5l(o&
zjj=jKad$6b)TUzx>b_zteBMA|29J{Cdct9OqKH|E56wr9AUTOB443kGdPW05BCplh
z2LB15Prst`6GY{(GzbG=tyXyiYe$@0I9>6}jjAxxCsR%hy79b3DHAe*UQ(}+X?AKv
zSi3`TIn`_!Ns5;Vm8l89!e8NJV$c-XFD5xN!dXySt!bR_4PAkP5Q;LXJ@ze{_WYM#
zGp~k$4a9RE3$K?(6*_s@?|*JYoWDCUokg^UA-`?+$>zn3^RIWi?(}LM4$Jq;)-E;n
zI3tVs5$DgVsYLYUq__M@+6I}o3Mjr)fx;~GO?19e3TE^y&{U~!FJs4GMs}{eJi{YS
zLzJu|whH;OM)Bl4o|Gfsa#X=3S=|Tr28mM%v(q^CUg0jyaFV>KuCr!p-EM-GnWFH8
za1B@^e(Cu;d&LRXXBg0U=2Ml%LM-nchPZOFcV7X?so42O*v+YYPi@xq$Vi_c@<T0E
zJ;Lc4?<Mtqgts4vwYw0$R?9>3a<9tL$RBca9tqdjoZ8}Z-fw<IW<mau&jYX)!106H
z*<g8Tx0+a{7E3Wy{9|NY?Myg6qFWM8@{Bb;v}D+esRPrvG~qUvaqpNB#-BPjE@Por
zIY=O@8zG~-nGBW@<f84*aI4AL0xXX^9qtM){<gE;U760TOs_Euv?=x;h9DPDE8Vso
z&70RMxB-Pu1?$6+s7)-IUL={Cmr)?wK>86dh)T%2JduO#<P~f8yGACOoWdyc`$={n
zqb;<f<x5CZrnG00uOup1Tt3G%MtD(vLg?*m56{$w&VNq-T5ZrJnMYy}dxsqf{K_;n
zC_Vgz=&ZuPxM{;+xcenNZf@`|DZOz?C*@KeXR!h&Syln33C7n-A$)6DR>fx0)v9;%
zzzDaE33T)7G}q$oVmbb<3uc={?R2|_622EIn!YU*c=2mpXq#2Z_Rk^!drNa@IjMeO
zl<bxIU<YHKLOzSAFdRMshwn_aU-G63Bk+`P^3m?Cg;Zjp_23wwC{FpE^lrXs;;M^L
zfVw&+E6j2zdEulg`MH#(bz@i2EP-i~(Cb_F!)^j|oTbcUF7A*flGqOvbd%}fuSF;D
zwQXvm1Q~n$y)(j{x3epTyR2ska_^dlx}G1KRwO=`61IDeaSE%?pM)9KPwEZiC$-ih
zlL>u(bNj`M8Eh70XAq7*#Cjg=1l+=Bm(t7k91%#@FB*z}?syYi!%IP6E)P`za)(#|
zxfx}v#YVfBy(sv8BUgqZjDgC#m}fyQmcQHuUi6TQT;zCQ!E84k8p651ksGqc+!m=N
zD0PU);TDKWA$Z0n6`YRGeH4EMdzxRs8%-kYoiJgE8MMtwxWo6+bS_U(Lz|+CGG)rM
zerboNBr(w|om%@#u6`G@<%J@aR3Ldpi|@;YHI*Lym6A|-ii}tD-JTn?Yupi$)2K5Q
z&<*WjXNgN-x&f87^;~6&`@?Rh3N=pwwxus1ZtB5&vPN!Ff0unV-{it}G9r25VL4JC
zKbdyKaFpU~<#pWWpPVkJ+`S;}7b4_<T4I+O!WxGIFS1%UwmriiM)s%_`BPL5GR<$z
zS>B0YmwOSLPdn5oZN2@v73x)(PCjj@D!(;#H;fYd(7e#8`sHIzF=};>P@;(b>=}f=
zetD4{;^vsY<De$(Fvjxl$s2Qv(s*&$r2KMGd&*F?v*Bv^RHX8la3ROnY~Kq5BUHwF
zQ|`c3kE73BO1)%Du8wDxVCglV6C16&a$e-SP9@IwwJc_{Bs&@pRH4Ccm?)P$a^4Q1
zhN|D)K<bJ<@{YN};Du62o9c5C`cL&-TYdZo-)6ud$V9Nh)l5foeWFXV=c1VhU(~|`
zYK1UqoC}w6GQNP?&;20QY^+|M<&mKX*coiGJ7Hons>!6J-+R^7?5(v|tEn9%WgjZV
zT+3H8E%a5qzD9jlSSRjnoLonny%?adxOb+DAgf3&CNwDd$47fAuOLWD!B%?|Iy<%t
zg1DZ&4IrJ$vduJR@4MDJMjP;SaLTk$vQ0^1Qpt$MG~0`_sBhjJx~{YitSX%it863h
zZou8CU+uXe-7Ng#d0@r02(z@B5f^!zj5VT+#K;X}(pMpN-Az?_bp|^q4XX*G`b|xO
zac`im!Tl4K8UhspXTtZitRo#mqB>Bd1;R70)IM!!NYiS9yT=dH?o=Uaifz1A=WO_Y
zBlQ{13DegpG$p||;;$?v^mIxaR(b4LRj%DKSyfgrLWv6}U7e<T@^}rKxb1Ci8Z`(E
zoGtzQ40NEF+%@Ez5x_MU5!-GS>byvoZX4BP`)-zdV_a61CT8PgG8~Z#3HsD-FKrax
zH2dcExm42%0hekkJp<>>w?K}BhR%?0^UK;cjSN^i(C!l!QXMG&nyHICvdXx<ahbGe
z^^|Wm8>ioyYwkWBc3x_hxALiGQYqn_KfvNvBf{-c5xljO1lSu%o))_7*Q{(()y}a`
zmQUb_>rJy8Etg(Uj$O>Q-odx0ZsnXK%}>;GM|7MIe0B4F$79suKTWZJ_l<f$z?K+X
zGy;G3&OSEK((v)}0(ogHY5t{g!#imI%sBi@<K~990{?Fszrg>Y0RZs!+W*A|<OTg7
z8m|Dn=l6f{0R-S3(f`m~A*MEVV91~5Xbl@LFx+HLbw@{6n#Z1Ic+0YagQX+QV+Zx0
Zu2>gWQ;6#y-{6LVcrh3lWK?A_{vUw_2z&ql

literal 0
HcmV?d00001

diff --git a/experiments/CHD_preventative_care_data.csv b/experiments/CHD_preventative_care_data.csv
new file mode 100644
index 00000000..b280758e
--- /dev/null
+++ b/experiments/CHD_preventative_care_data.csv
@@ -0,0 +1,102 @@
+risk_levels,prior_risk_distribution,TRS_if_healthy,TRS_if_sick,GRS_if_healthy,GRS_if_sick
+0,0.11293990700000,0.39740624000000,0.04893543700000,0.02252611500000,0.00283870400000
+0.01,0.12465484400000,0.14375880600000,0.03547850500000,0.07311748200000,0.03579809800000
+0.02,0.09890158000000,0.10569106600000,0.05107539100000,0.12684417900000,0.10940458400000
+0.03,0.07045429200000,0.06429586300000,0.05294659100000,0.18990717200000,0.11918963900000
+0.04,0.04390572300000,0.05056367600000,0.05803457100000,0.11897423000000,0.13241767500000
+0.05,0.06083044000000,0.03052270900000,0.04297088900000,0.11373319200000,0.13705796900000
+0.06,0.03031388200000,0.02095879400000,0.04511271500000,0.08843425800000,0.09294875300000
+0.07,0.03966468200000,0.03275560100000,0.03669923600000,0.04883144600000,0.08072231700000
+0.08,0.02749771500000,0.03116940200000,0.02995658000000,0.03414356500000,0.06469767600000
+0.09,0.03349420200000,0.02003049500000,0.03945800100000,0.02638107900000,0.03036662900000
+0.1,0.01894239500000,0.00655804300000,0.03828682200000,0.01679492900000,0.03121447300000
+0.11,0.01655726700000,0.00880561100000,0.02568766200000,0.01527908700000,0.02907518800000
+0.12,0.02527272300000,0.00632351800000,0.03223956000000,0.02132732600000,0.01200694000000
+0.13,0.02489687600000,0.00062280900000,0.02773266400000,0.01853267700000,0.01152809400000
+0.14,0.01810201800000,0.00678644900000,0.03350329900000,0.00744960200000,0.01785592900000
+0.15,0.00000000000000,0.00833391300000,0.02556680800000,0.00742368700000,0.01787462400000
+0.16,0.01774679900000,0.00060849200000,0.01627451000000,0.01190778400000,0.01298614400000
+0.17,0.02056393700000,0.00288562700000,0.02217757000000,0.00953420300000,0.01666081400000
+0.18,0.01324023100000,0.00270084800000,0.02124150300000,0.00576241100000,0.01480715600000
+0.19,0.01531385100000,0.00589751000000,0.02479449000000,0.01144990600000,0.01115178600000
+0.2,0.01733562400000,0.00335190500000,0.02195192200000,0.00876195700000,0.01413693600000
+0.21,0.02116504900000,0.00402482500000,0.01956580500000,0.00895109400000,0.00413693600000
+0.22,0.01000000000000,0.00167715800000,0.01787645700000,0.00831543600000,0.00112293600000
+0.23,0.01601149200000,0.00307826100000,0.02319705100000,0.00539759000000,0.00000000000000
+0.24,0.00000000000000,0.00264528400000,0.01020576700000,0.00014600000000,0.00000000000000
+0.25,0.01846514000000,0.00356238400000,0.01657776100000,0.00003414600000,0.00000000000000
+0.26,0.00000000000000,0.00206032000000,0.01563507400000,0.00001677000000,0.00000000000000
+0.27,0.00694367200000,0.00181872700000,0.01312182500000,0.00002267700000,0.00000000000000
+0.28,0.00297060000000,0.00063827400000,0.01557198900000,0.00000000000000,0.00000000000000
+0.29,0.00469651100000,0.00374629400000,0.01719464600000,0.00000000000000,0.00000000000000
+0.3,0.00544617600000,0.00124890500000,0.01135377200000,0.00000000000000,0.00000000000000
+0.31,0.00123335500000,0.00149568700000,0.01212877300000,0.00000000000000,0.00000000000000
+0.32,0.00378659900000,0.00041196500000,0.01390704100000,0.00000000000000,0.00000000000000
+0.33,0.00830991800000,0.00142893800000,0.01750051600000,0.00000000000000,0.00000000000000
+0.34,0.00000000000000,0.00043471200000,0.00608215400000,0.00000000000000,0.00000000000000
+0.35,0.00939161600000,0.00056847000000,0.00733665200000,0.00000000000000,0.00000000000000
+0.36,0.00000000000000,0.00016719700000,0.00514977400000,0.00000000000000,0.00000000000000
+0.37,0.00000000000000,0.00020063600000,0.00276073400000,0.00000000000000,0.00000000000000
+0.38,0.00000000000000,0.00016719700000,0.00109835500000,0.00000000000000,0.00000000000000
+0.39,0.00161001700000,0.00253184700000,0.00279437100000,0.00000000000000,0.00000000000000
+0.4,0.00000000000000,0.00087930000000,0.00135650700000,0.00000000000000,0.00000000000000
+0.41,0.00000000000000,0.00122641700000,0.00150422700000,0.00000000000000,0.00000000000000
+0.42,0.00624891500000,0.00088858100000,0.00147217300000,0.00000000000000,0.00000000000000
+0.43,0.00000000000000,0.00013375800000,0.00251989500000,0.00000000000000,0.00000000000000
+0.44,0.00000000000000,0.00086195100000,0.00391024100000,0.00000000000000,0.00000000000000
+0.45,0.00597112500000,0.00235978300000,0.00372377600000,0.00000000000000,0.00000000000000
+0.46,0.00000000000000,0.00124009800000,0.00280890000000,0.00000000000000,0.00000000000000
+0.47,0.00176904300000,0.00013375800000,0.00000000000000,0.00000000000000,0.00000000000000
+0.48,0.00000000000000,0.00026751500000,0.00147103500000,0.00000000000000,0.00000000000000
+0.49,0.00323970300000,0.00004749700000,0.00148041900000,0.00000000000000,0.00000000000000
+0.5,0.00000000000000,0.00010031800000,0.00104513800000,0.00000000000000,0.00000000000000
+0.51,0.00000000000000,0.00190945300000,0.00116669500000,0.00000000000000,0.00000000000000
+0.52,0.00000000000000,0.00198581600000,0.00399394800000,0.00000000000000,0.00000000000000
+0.53,0.00000000000000,0.00006690000000,0.00129514400000,0.00000000000000,0.00000000000000
+0.54,0.00000000000000,0.00005890000000,0.00009660000000,0.00000000000000,0.00000000000000
+0.55,0.00430848100000,0.00029373000000,0.00000000000000,0.00000000000000,0.00000000000000
+0.56,0.00402460400000,0.00003340000000,0.00025752700000,0.00000000000000,0.00000000000000
+0.57,0.00000000000000,0.00217202900000,0.00269709800000,0.00000000000000,0.00000000000000
+0.58,0.00000000000000,0.00000000000000,0.00009690000000,0.00000000000000,0.00000000000000
+0.59,0.00000000000000,0.00205186700000,0.00107007000000,0.00000000000000,0.00000000000000
+0.6,0.00119385400000,0.00000000000000,0.00015545500000,0.00000000000000,0.00000000000000
+0.61,0.00000000000000,0.00027794900000,0.00016539900000,0.00000000000000,0.00000000000000
+0.62,0.00000000000000,0.00005370000000,0.00018467900000,0.00000000000000,0.00000000000000
+0.63,0.00346794400000,0.00003350000000,0.00003030000000,0.00000000000000,0.00000000000000
+0.64,0.00000000000000,0.00003340000000,0.00000000000000,0.00000000000000,0.00000000000000
+0.65,0.00658831300000,0.00005210000000,0.00378400900000,0.00000000000000,0.00000000000000
+0.66,0.00000000000000,0.00000000000000,0.00000000000000,0.00000000000000,0.00000000000000
+0.67,0.00000000000000,0.00004440000000,0.00000000000000,0.00000000000000,0.00000000000000
+0.68,0.00000000000000,0.00047234300000,0.00000000000000,0.00000000000000,0.00000000000000
+0.69,0.00340850900000,0.00000000000000,0.00000000000000,0.00000000000000,0.00000000000000
+0.7,0.00197460100000,0.00000000000000,0.00000000000000,0.00000000000000,0.00000000000000
+0.71,0.00000000000000,0.00003750000000,0.00000000000000,0.00000000000000,0.00000000000000
+0.72,0.00000000000000,0.00001110000000,0.00000000000000,0.00000000000000,0.00000000000000
+0.73,0.00000000000000,0.00001170000000,0.00000000000000,0.00000000000000,0.00000000000000
+0.74,0.00000000000000,0.00003340000000,0.00000000000000,0.00000000000000,0.00000000000000
+0.75,0.00386885800000,0.00000139000000,0.00000000000000,0.00000000000000,0.00000000000000
+0.76,0.00000000000000,0.00003340000000,0.00000000000000,0.00000000000000,0.00000000000000
+0.77,0.00000000000000,0.00000000000000,0.00000529000000,0.00000000000000,0.00000000000000
+0.78,0.00000000000000,0.00000000000000,0.00010941900000,0.00000000000000,0.00000000000000
+0.79,0.00000000000000,0.00000000000000,0.00000000000000,0.00000000000000,0.00000000000000
+0.8,0.00334460600000,0.00000000000000,0.00000000000000,0.00000000000000,0.00000000000000
+0.81,0.00000000000000,0.00003430000000,0.00000000000000,0.00000000000000,0.00000000000000
+0.82,0.00000000000000,0.00004220000000,0.00132149500000,0.00000000000000,0.00000000000000
+0.83,0.00000000000000,0.00005170000000,0.00128714500000,0.00000000000000,0.00000000000000
+0.84,0.00000000000000,0.00003340000000,0.00028930000000,0.00000000000000,0.00000000000000
+0.85,0.00000000000000,0.00003340000000,0.00029201000000,0.00000000000000,0.00000000000000
+0.86,0.00425751100000,0.00003340000000,0.00063885500000,0.00000000000000,0.00000000000000
+0.87,0.00000000000000,0.00002340000000,0.00000000000000,0.00000000000000,0.00000000000000
+0.88,0.00000000000000,0.00000000000000,0.00058710800000,0.00000000000000,0.00000000000000
+0.89,0.00000000000000,0.00000000000000,0.00000000000000,0.00000000000000,0.00000000000000
+0.9,0.00000000000000,0.00000658000000,0.00000000000000,0.00000000000000,0.00000000000000
+0.91,0.00000000000000,0.00000020900000,0.00000000000000,0.00000000000000,0.00000000000000
+0.92,0.00000000000000,0.00000000000000,0.00000000000000,0.00000000000000,0.00000000000000
+0.93,0.00238340700000,0.00000000000000,0.00000000000000,0.00000000000000,0.00000000000000
+0.94,0.00000000000000,0.00000000000000,0.00000000000000,0.00000000000000,0.00000000000000
+0.95,0.00290569600000,0.00000000000000,0.00000000000000,0.00000000000000,0.00000000000000
+0.96,0.00000000000000,0.00000000000000,0.00000000000000,0.00000000000000,0.00000000000000
+0.97,0.00000000000000,0.00000000000000,0.00000000000000,0.00000000000000,0.00000000000000
+0.98,0.00000000000000,0.00000000000000,0.00000000000000,0.00000000000000,0.00000000000000
+0.99,0.00038569700000,0.00000000000000,0.00000000000000,0.00000000000000,0.00000000000000
+1,0.00000000000000,0.00000000000000,0.00000000000000,0.00000000000000,0.00000000000000
\ No newline at end of file
diff --git a/experiments/array_script_CHD.sh b/experiments/array_script_CHD.sh
new file mode 100644
index 00000000..bacbd89a
--- /dev/null
+++ b/experiments/array_script_CHD.sh
@@ -0,0 +1,10 @@
+#!/bin/bash
+#SBATCH --time=0:40:00
+#SBATCH --mem=16G
+#SBATCH --cpus-per-task=8
+#SBATCH --output=results/CHD_%a.out
+#SBATCH --error=/dev/null
+#SBATCH --array=0-100
+
+module load julia
+srun julia slurmjob_CHD.jl $SLURM_ARRAY_TASK_ID
\ No newline at end of file
diff --git a/experiments/slurmjob_CHD.jl b/experiments/slurmjob_CHD.jl
new file mode 100644
index 00000000..09891357
--- /dev/null
+++ b/experiments/slurmjob_CHD.jl
@@ -0,0 +1,209 @@
+using Pkg
+Pkg.activate((@__DIR__)*"/..")
+using Logging
+using JuMP, Gurobi
+using DecisionProgramming
+using CSV, DataFrames, PrettyTables
+using DelimitedFiles
+
+
+
+# Setting subproblem specific parameters
+idx = parse(Int64, ARGS[1])
+chosen_risk_level = string(idx)*"%" 
+
+
+
+# Reading tests' technical performance data (dummy data in this case)
+data = CSV.read("CHD_preventative_care_data.csv", DataFrame)
+
+
+# Bayes posterior risk probabilities calculation function
+# prior = prior risk level for which the posterior risk distribution is calculated for,
+# t = test done
+# returns a 100x1 vector with the probabilities of getting CHD given the prior risk level and test result
+# for no test done (i.e. t = 3) returns a zero vector
+function update_risk_distribution(prior::Int64, t::Int64)
+    if t == 1 # the test is TRS
+        # P(TRS = result | sick) = P(TRS_if_sick = result) * P(sick) = P(TRS_if_sick = result) * P(prior_risk)
+        numerators = data.TRS_if_sick .* data.risk_levels[prior]
+
+        # P(TRS = result) = P(TRS_if_sick = result) * P(sick) + P(TRS_if_healthy = result) * P(healthy)
+        denominators = data.TRS_if_sick .* data.risk_levels[prior]  + data.TRS_if_healthy .* (1 - data.risk_levels[prior])
+
+        posterior_risks = numerators./denominators
+
+        # if the denominator is zero, post_risk is NaN, changing those to 0
+        for i = 1:101
+            if isnan(posterior_risks[i])
+                posterior_risks[i] = 0
+            end
+        end
+
+        return posterior_risks
+
+
+    elseif t == 2 #the test is GRS
+        numerators = (data.GRS_if_sick .* data.risk_levels[prior])
+        denominators = data.GRS_if_sick .* data.risk_levels[prior]  + data.GRS_if_healthy .* (1 .- data.risk_levels[prior])
+
+        posterior_risks =  numerators./denominators
+
+        # if the denominator is zero, post_risk is NaN, changing those to 0
+        for i = 1:101
+            if isnan(posterior_risks[i])
+                posterior_risks[i] = 0
+            end
+        end
+
+        return posterior_risks
+
+
+    else # no test performed
+        risks_unchanged = zeros(100,1)
+
+
+        return risks_unchanged
+
+    end
+end
+
+# State probabilites calculation function
+# risk_p = the resulting array from update_risk_distribution
+# t = test done
+# h = CHD or no CHD
+# returns the probability distribution in 101x1 vector for the states of the R node given the prior risk level (must be same as to function update_risk_distribution), test t and health h
+function state_probabilities(risk_p::Array{Float64}, t::Int64, h::Int64, prior::Int64)
+
+    #if no test is performed, then the probabilities of moving to states (other than the prior risk level) are 0 and to the prior risk element is 1
+    if t == 3
+        state_probabilites = zeros(101)
+        state_probabilites[prior] = 1.0
+        return state_probabilites
+    end
+
+    # return vector
+    state_probabilites = zeros(101)
+
+    # copying the probabilities of the scores for ease of readability
+    if h == 1 && t == 1    # CHD and TRS
+        p_scores = data.TRS_if_sick
+    elseif t ==1    # no CHD and TRS
+        p_scores = data.TRS_if_healthy
+    elseif h == 1 && t == 2 # CHD and GRS
+        p_scores = data.GRS_if_sick
+    else # no CHD and GRS
+        p_scores = data.GRS_if_healthy
+    end
+
+    for i = 1:101 #iterating through all risk levels 0%, 1%, ..., 99%, 100% in data.risk_levels
+        for j = 1:101 #iterates through all risk estimates in risk_p
+            #finding all risk estimates risk_p[j] within risk level i
+            # risk_level[i] <= risk_p < risk_level[i]
+            if i < 101 && data.risk_levels[i] <= risk_p[j] && risk_p[j] < data.risk_levels[i+1]
+                state_probabilites[i] += p_scores[j]
+            elseif i == 101 && data.risk_levels[i] <= risk_p[j] #special case: the highest risk level[101] = 100%
+                state_probabilites[i] += p_scores[j]
+            end
+        end
+    end
+
+    return state_probabilites
+end
+
+
+@info("Creating the influence diagram.")
+diagram = InfluenceDiagram()
+
+const H_states = ["CHD", "no CHD"]
+const T_states = ["TRS", "GRS", "no test"]
+const TD_states = ["treatment", "no treatment"]
+const R_states = [string(x) * "%" for x in [0:1:100;]]
+
+
+add_node!(diagram, ChanceNode("R0", [], R_states))
+add_node!(diagram, ChanceNode("R1", ["R0", "H", "T1"], R_states))
+add_node!(diagram, ChanceNode("R2", ["R1", "H", "T2"], R_states))
+add_node!(diagram, ChanceNode("H", ["R0"], H_states))
+
+add_node!(diagram, DecisionNode("T1", ["R0"], T_states))
+add_node!(diagram, DecisionNode("T2", ["R1"], T_states))
+add_node!(diagram, DecisionNode("TD", ["R2"], TD_states))
+
+add_node!(diagram, ValueNode("TC", ["T1", "T2"]))
+add_node!(diagram, ValueNode("HB", ["H", "TD"]))
+
+
+generate_arcs!(diagram)
+
+X_R0 = ProbabilityMatrix(diagram, "R0")
+X_R0[chosen_risk_level] = 1
+add_probabilities!(diagram, "R0", X_R0)
+
+
+X_H = ProbabilityMatrix(diagram, "H")
+X_H[:, "CHD"] = data.risk_levels
+X_H[:, "no CHD"] = 1 .- data.risk_levels
+add_probabilities!(diagram, "H", X_H)
+
+X_R = ProbabilityMatrix(diagram, "R1")
+for s_R0 = 1:101, s_H = 1:2, s_T1 = 1:3
+    X_R[s_R0, s_H,  s_T1, :] =  state_probabilities(update_risk_distribution(s_R0, s_T1), s_T1, s_H, s_R0)
+end
+add_probabilities!(diagram, "R1", X_R)
+add_probabilities!(diagram, "R2", X_R)
+
+
+cost_TRS = -0.0034645
+cost_GRS = -0.004
+forbidden = 0     #the cost of forbidden test combinations is negligible
+Y_TC = UtilityMatrix(diagram, "TC")
+Y_TC["TRS", "TRS"] = forbidden
+Y_TC["TRS", "GRS"] = cost_TRS + cost_GRS
+Y_TC["TRS", "no test"] = cost_TRS
+Y_TC["GRS", "TRS"] = cost_TRS + cost_GRS
+Y_TC["GRS", "GRS"] = forbidden
+Y_TC["GRS", "no test"] = cost_GRS
+Y_TC["no test", "TRS"] = cost_TRS
+Y_TC["no test", "GRS"] = cost_GRS
+Y_TC["no test", "no test"] = 0
+add_utilities!(diagram, "TC", Y_TC)
+
+Y_HB = UtilityMatrix(diagram, "HB")
+Y_HB["CHD", "treatment"] = 6.89713671259061
+Y_HB["CHD", "no treatment"] = 6.65436854256236
+Y_HB["no CHD", "treatment"] = 7.64528451705134
+Y_HB["no CHD", "no treatment"] = 7.70088349200034
+add_utilities!(diagram, "HB", Y_HB)
+
+generate_diagram!(diagram)
+
+
+@info("Creating the decision model.")
+model = Model()
+z = DecisionVariables(model, diagram)
+
+# Defining forbidden paths to include all those where a test is repeated twice
+forbidden_tests = ForbiddenPath(diagram, ["T1","T2"], [("TRS", "TRS"),("GRS", "GRS"),("no test", "TRS"), ("no test", "GRS")])
+fixed_R0 = FixedPath(diagram, Dict("R0" => chosen_risk_level))
+scale_factor = 10000.0
+x_s = PathCompatibilityVariables(model, diagram, z; fixed = fixed_R0, forbidden_paths = [forbidden_tests], probability_cut=false)
+
+EV = expected_value(model, diagram, x_s)
+@objective(model, Max, EV)
+
+@info("Starting the optimization process.")
+optimizer = optimizer_with_attributes(
+    () -> Gurobi.Optimizer(Gurobi.Env()),
+    "TimeLimit" => 1800,
+    "IntFeasTol"=> 1e-9,
+    "MIPGap" => 1e-6
+)
+set_optimizer(model, optimizer)
+
+optimize!(model)
+
+if termination_status(model) == OPTIMAL
+    writedlm((@__DIR__)*"/results/CHD_"*string(idx)*".csv", [objective_value(model), solve_time(model)], ',')
+end
+