diff --git a/data/covid19_DTM/interim/QALY_model/long_COVID/average_QALY_losses.csv b/data/covid19_DTM/interim/QALY_model/long_COVID/average_QALY_losses.csv index 1f7eab0de..daae269d8 100644 --- a/data/covid19_DTM/interim/QALY_model/long_COVID/average_QALY_losses.csv +++ b/data/covid19_DTM/interim/QALY_model/long_COVID/average_QALY_losses.csv @@ -1,429 +1,429 @@ hospitalisation,age,mean,sd,lower,upper -Non-hospitalised,0,0.7935148886177248,0.31417909127761423,0.2423022087164734,1.3867414434929877 -Non-hospitalised,1,0.7827494130434783,0.3097146420948966,0.2395149785780798,1.3677860389502539 -Non-hospitalised,2,0.7712002244030746,0.30490267873782745,0.23650705811652592,1.347298894566754 -Non-hospitalised,3,0.7596813652244611,0.30010216731982525,0.23350584122365797,1.326856441352946 -Non-hospitalised,4,0.7482408870740573,0.29533393440503697,0.23052449642670778,1.3065493551390754 -Non-hospitalised,5,0.7368675175767143,0.2905926542348832,0.22755963232128923,1.2863535774735702 -Non-hospitalised,6,0.7255945296806691,0.28589278645854105,0.2246203862763359,1.2663321215080803 -Non-hospitalised,7,0.7144123144004324,0.28122985282734186,0.22170391132718031,1.2464649032685233 -Non-hospitalised,8,0.7033327079834,0.2766089071509234,0.2188133968142906,1.226773792050471 -Non-hospitalised,9,0.6923647074251721,0.2720337764746854,0.21595125733362705,1.2072753252261421 -Non-hospitalised,10,0.6815038211960737,0.2675023396904403,0.2131161481716525,1.1879599911054044 -Non-hospitalised,11,0.6707528947783638,0.26301574723985033,0.21030879671849867,1.1688306106519284 -Non-hospitalised,12,0.6601287137832386,0.2585813406230745,0.20753384436973676,1.1499161257235309 -Non-hospitalised,13,0.6496123391488879,0.25419066177624083,0.20478594896876634,1.1311850714457812 -Non-hospitalised,14,0.6392338755946715,0.24985697790217817,0.2020734968526022,1.1126953691559838 -Non-hospitalised,15,0.628960516394513,0.24556572673826393,0.1993872954019011,1.0943834207481642 -Non-hospitalised,16,0.6188301943300375,0.2413336593680989,0.19673792964129624,1.0763223656060572 -Non-hospitalised,17,0.6088119705630577,0.23714707108834698,0.19411675412788063,1.0584524247803198 -Non-hospitalised,18,0.5989231785197543,0.23301361891831826,0.19152861124355527,1.0408070596451031 -Non-hospitalised,19,0.5891840019609387,0.2289422159078727,0.18897913172373068,1.0234251716732503 -Non-hospitalised,20,0.5795425406404991,0.22490994590460262,0.18645386787650484,1.0062069474968398 -Non-hospitalised,21,0.5700199194502871,0.2209261665871037,0.18395873526921233,0.9891932882753016 -Non-hospitalised,22,0.560647002691227,0.21700451760302575,0.1815023433977218,0.9724436822314451 -Non-hospitalised,23,0.5513684681113125,0.21312063509026077,0.17906934283457887,0.955852156707679 -Non-hospitalised,24,0.542205012493777,0.2092836980624765,0.17666553314345293,0.9394588175861959 -Non-hospitalised,25,0.5331560027208048,0.20549348185606553,0.17429078457226863,0.923262823311068 -Non-hospitalised,26,0.5242070694194199,0.2017437291664715,0.17194116763412412,0.9072371166675576 -Non-hospitalised,27,0.5153909323713911,0.19804889040910592,0.1696257954504701,0.8914446501736304 -Non-hospitalised,28,0.5066689120771534,0.19439203743811367,0.16733401803759956,0.8758120052238206 -Non-hospitalised,29,0.4980560652033406,0.19077986410099462,0.16507006241649913,0.8603684234400621 -Non-hospitalised,30,0.48954792287719284,0.18721048966738799,0.16283275543145329,0.8451058898036843 -Non-hospitalised,31,0.48115176279359045,0.18368718150431668,0.16062416356752304,0.8300387283163306 -Non-hospitalised,32,0.47283797983823694,0.18019707257967416,0.15843620075840215,0.815111283166894 -Non-hospitalised,33,0.4646133028640916,0.17674321396277662,0.15627079859037468,0.8003369707819383 -Non-hospitalised,34,0.4564867842864103,0.17332965721442697,0.1541305158309651,0.7857335141536175 -Non-hospitalised,35,0.44843640084216596,0.16994689597496765,0.15200938347156073,0.7712598911259065 -Non-hospitalised,36,0.44047745227001206,0.16660170663701399,0.14991167285382384,0.7569456200394192 -Non-hospitalised,37,0.4325923770959779,0.16328654830229639,0.1478326517618652,0.7427582157257817 -Non-hospitalised,38,0.42479723650052387,0.1600085098872848,0.14577678554822232,0.7287285131533786 -Non-hospitalised,39,0.41707195277658343,0.1567589712364413,0.14373866393724696,0.7148193733115009 -Non-hospitalised,40,0.40943830924538743,0.1535474815097419,0.1417242961863202,0.7010722313541701 -Non-hospitalised,41,0.4018688617998613,0.15036224123277953,0.13972627456185963,0.6874362210876448 -Non-hospitalised,42,0.39434623879298997,0.14719585352625505,0.13773996034173427,0.6738795580051353 -Non-hospitalised,43,0.38691289457926703,0.14406677745828086,0.13577695786345886,0.660482115872698 -Non-hospitalised,44,0.3795152715093309,0.14095198082564428,0.1338228058153735,0.6471445907432892 -Non-hospitalised,45,0.37218106320238015,0.1378635273148763,0.13188508951675093,0.6339192134579709 -Non-hospitalised,46,0.36493160152523535,0.13481068315311648,0.12996963379233958,0.6208460675629366 -Non-hospitalised,47,0.35770584467693645,0.13176726596873148,0.1280600035240858,0.6078123604633976 -Non-hospitalised,48,0.350538460416561,0.128748290868091,0.1261656341769568,0.5948829823531331 -Non-hospitalised,49,0.34339896025101013,0.1257407593627082,0.12427837171123783,0.5820020471338104 -Non-hospitalised,50,0.33633086678332025,0.12276342788056624,0.12240999681011563,0.5692504862448811 -Non-hospitalised,51,0.3292830118786437,0.11979444512894105,0.12054679934796154,0.5565342999657885 -Non-hospitalised,52,0.32230208498696855,0.11685387944369699,0.11870138088073971,0.5439400095668674 -Non-hospitalised,53,0.31534623100648784,0.11392389543991802,0.11686255381322944,0.5313909476620361 -Non-hospitalised,54,0.30843654909027024,0.11101358376213329,0.11503602862754432,0.5189263425461095 -Non-hospitalised,55,0.3015685769068029,0.10812109781206652,0.11322065175310768,0.5065383414269583 -Non-hospitalised,56,0.2947103185646623,0.10523286646130685,0.11140791516372041,0.4941687007563716 -Non-hospitalised,57,0.2878796610194109,0.10235658781527415,0.10960265460985635,0.48185063235886616 -Non-hospitalised,58,0.281078714483823,0.09949321916481546,0.10780547603854382,0.4695883448193632 -Non-hospitalised,59,0.2743185246686917,0.09664750468010344,0.10601935961956188,0.4574022969453085 -Non-hospitalised,60,0.2675449826959003,0.09379651598408459,0.10422993228962933,0.4451941173212019 -Non-hospitalised,61,0.26085610365366424,0.0909819048529854,0.10246331653053661,0.4331427078682745 -Non-hospitalised,62,0.2542172640878498,0.088188983003333,0.10071030375636537,0.4211850270620155 -Non-hospitalised,63,0.2476623948624547,0.08543213287186598,0.09897991170612314,0.4093827973148247 -Non-hospitalised,64,0.24113485305062402,0.08268737069427638,0.09725710653920333,0.3976331422338611 -Non-hospitalised,65,0.23464605523374754,0.07995955610507398,0.09554493683268384,0.38595676715040733 -Non-hospitalised,66,0.22819820054051512,0.07724963793823228,0.09384399864742402,0.37435693562442185 -Non-hospitalised,67,0.22171127343412578,0.07452384640021663,0.09213313172079426,0.3626901056038824 -Non-hospitalised,68,0.21532235306006137,0.07184010347162018,0.09044863360836018,0.3512046363294416 -Non-hospitalised,69,0.20880271248230606,0.06910192268213117,0.08873004931400347,0.33948706771036413 -Non-hospitalised,70,0.2024116298253288,0.06641871657644091,0.08704593094089927,0.3280062732692569 -Non-hospitalised,71,0.19607635176058785,0.06375978375180816,0.0853770377989497,0.316630642353807 -Non-hospitalised,72,0.1897299795088367,0.06109696398931985,0.08370572581517342,0.3052395714186561 -Non-hospitalised,73,0.1834680812510294,0.0584705079420584,0.0820572017006151,0.2940053562653685 -Non-hospitalised,74,0.17723473424628533,0.05585688370131766,0.08041671744131272,0.2828271909309509 -Non-hospitalised,75,0.17109255761631192,0.05328242519255187,0.0788007552611922,0.2718176755572011 -Non-hospitalised,76,0.16495031121335932,0.05070879662537961,0.07718529338084695,0.2608126885650739 -Non-hospitalised,77,0.15872635231898707,0.04810178002485507,0.07554888860688104,0.24966581824540907 -Non-hospitalised,78,0.15276392009606496,0.04560547042183685,0.07398175664433843,0.2389931573301006 -Non-hospitalised,79,0.14667311781237452,0.043056281149969323,0.07238142957806416,0.22809368640775943 -Non-hospitalised,80,0.14076698744134142,0.04058554276818634,0.07083011853690484,0.2175295951878685 -Non-hospitalised,81,0.1349938153989904,0.03817155006923483,0.0693141954703488,0.207207967565043 -Non-hospitalised,82,0.12923108227153324,0.03576302600750779,0.0678015002822642,0.19690923509692723 -Non-hospitalised,83,0.12363054361234647,0.03342355461279164,0.06633180288474917,0.18690482781413256 -Non-hospitalised,84,0.11828395096128308,0.03119150057296133,0.06492908485990648,0.17735833025144296 -Non-hospitalised,85,0.11301285810714262,0.0289922809966988,0.06341548736416927,0.1679503164982513 -Non-hospitalised,86,0.10804355725144464,0.026920479952600874,0.061870929404384975,0.15908472964646467 -Non-hospitalised,87,0.10331065387783417,0.02494876678558532,0.060398683873747826,0.15064416942451891 -Non-hospitalised,88,0.09894204047042399,0.023130449094346577,0.059043579882771,0.14285634028494815 -Non-hospitalised,89,0.09486232549578025,0.021434018507333882,0.05777747743634109,0.13558608339449563 -Non-hospitalised,90,0.09090642347274541,0.019790783880827577,0.05654905725772251,0.12853856861368118 -Non-hospitalised,91,0.08745567935323478,0.018359307167724277,0.05547649889161276,0.12239333147447279 -Non-hospitalised,92,0.08428190916290065,0.017044530582410105,0.05448932467408543,0.116743019734608 -Non-hospitalised,93,0.0810221461407603,0.015696050946339683,0.053475106242703116,0.11094022417338502 -Non-hospitalised,94,0.07803624561500065,0.014463270736131318,0.05254541162429662,0.1056240421490883 -Non-hospitalised,95,0.07547168361778242,0.013406894799561975,0.0520623275658384,0.10105956153602945 -Non-hospitalised,96,0.0732158168989003,0.012479980565553977,0.05125475863913113,0.09704573928308177 -Non-hospitalised,97,0.07121992174987073,0.011662126420764942,0.050539496838064885,0.09349558738536061 -Non-hospitalised,98,0.0694951062481589,0.010957518873527404,0.049920565908253146,0.09042875852937317 -Non-hospitalised,99,0.0676621558481355,0.010211005131325116,0.04926764636460227,0.08694604233024915 -Non-hospitalised,100,0.06660967501050723,0.00978464969843044,0.04888019548180275,0.08516312927937361 -Non-hospitalised,101,0.06444900691655336,0.008911326539880656,0.04824326244189084,0.08134413602632591 -Non-hospitalised,102,0.0632482299581572,0.008430091001136241,0.047857918548256675,0.07906973277082777 -Non-hospitalised,103,0.06089595018634276,0.00749521963573015,0.047100342826921646,0.07490610407480318 -Non-hospitalised,104,0.05781585450864537,0.006304275173193031,0.045725965163107,0.06979073120140537 -Non-hospitalised,105,0.05433204300154784,0.005054278460054339,0.04482553774948081,0.06380253286448412 -Non-hospitalised,106,0.04259856995048379,0.0025554534387112883,0.037669669526621076,0.04764855470370469 -Cohort,0,2.0388817915749544,0.8992141967774195,0.5285961407065378,3.980240087870274 -Cohort,1,2.0134230834527584,0.8883216919041147,0.5225216223479228,3.9322309279392007 -Cohort,2,1.9847195763351742,0.8756827870389561,0.5157198246577558,3.8768241463635977 -Cohort,3,1.9560070244021366,0.8630201205795083,0.5089184342877962,3.8213261295778276 -Cohort,4,1.9274548056137717,0.8504205393786162,0.5021560053814846,3.766107506046169 -Cohort,5,1.8989978158261847,0.8378461386377747,0.4954183257906335,3.711010166476132 -Cohort,6,1.8707553516469388,0.8253584528832262,0.4887324952467594,3.656296395324024 -Cohort,7,1.842675764306794,0.8129277856156978,0.4820871802419533,3.6018422099383716 -Cohort,8,1.8147965033329083,0.8005725484476655,0.4754910027153992,3.547726921975398 -Cohort,9,1.787145249718417,0.7883061863989813,0.4689503559995824,3.4940086154921612 -Cohort,10,1.7596967302087998,0.7761139755008981,0.4624597114479176,3.4406261647003626 -Cohort,11,1.732456702830841,0.7639980377048134,0.45602047992448974,3.3875895342564264 -Cohort,12,1.7054880237544654,0.7519912666631426,0.44964687716760604,3.3350385403375893 -Cohort,13,1.6787110425008647,0.7400503435440826,0.443321083650812,3.2827909406191416 -Cohort,14,1.65224360399696,0.7282378843244827,0.43706963817329364,3.2321802367690022 -Cohort,15,1.6259527390496797,0.7164823314726737,0.43086269459502313,3.181975425298696 -Cohort,16,1.5999884309039099,0.7048636748505103,0.42473399738692313,3.132324424983575 -Cohort,17,1.5742271773860805,0.6933157150742842,0.4186557876081998,3.082913170431326 -Cohort,18,1.5487379041299882,0.6818752474043206,0.41264356561644,3.0339166404283153 -Cohort,19,1.52360027094271,0.670584751516135,0.4067152573709569,2.985536382187579 -Cohort,20,1.4986115597519816,0.6593362413693363,0.4008252077736613,2.9372619278643395 -Cohort,21,1.4738569227228404,0.6481753981066387,0.3949925422884715,2.8893107554112394 -Cohort,22,1.449458179172876,0.63716725179694,0.389244672715591,2.8419910859717934 -Cohort,23,1.4252032840060007,0.6261994514956624,0.3835337601693089,2.79477363572516 -Cohort,24,1.4011759842255296,0.6153170123398394,0.37787860788799466,2.747872652779641 -Cohort,25,1.3773774162012642,0.6045210705947947,0.372279398876615,2.7012951734454296 -Cohort,26,1.3537569921977228,0.5937855130305209,0.3667246271042562,2.6549202165548405 -Cohort,27,1.330442779360257,0.5831785972170478,0.3612431316453844,2.6090705798522125 -Cohort,28,1.307292980449148,0.5726262127639976,0.35580278716922636,2.5633999347766183 -Cohort,29,1.28436914387444,0.5621615775241053,0.350417400454356,2.518066435217357 -Cohort,30,1.261659370027181,0.551779222725452,0.3450841783337321,2.473046541656386 -Cohort,31,1.23919523067138,0.5414964710274143,0.3398101677362464,2.428423394114581 -Cohort,32,1.2168733540213308,0.5312601795771764,0.33457183629213977,2.3839508720098106 -Cohort,33,1.1947245850290604,0.5210875217239584,0.3293760381100992,2.3397119236717665 -Cohort,34,1.1727872918589524,0.5109994381535137,0.32423135349487575,2.2958066465001026 -Cohort,35,1.1509882197543881,0.5009589386952332,0.31912102269615916,2.2520653398529835 -Cohort,36,1.1293880375007925,0.49099852392734006,0.31405869163902583,2.208641827241521 -Cohort,37,1.10793049833849,0.48109021395063456,0.30903143970080216,2.165408751228977 -Cohort,38,1.0866778293389552,0.47126711690283324,0.30405329274979614,2.1225222219145503 -Cohort,39,1.0655657303503345,0.46149725479863085,0.2991094799211967,2.0798366548310137 -Cohort,40,1.04467490388135,0.4518229592895164,0.29421826757361597,2.03755020256099 -Cohort,41,1.0239164884910839,0.4421998476223247,0.2893592666991986,1.995460302284042 -Cohort,42,1.0032381935507533,0.4326025936633972,0.2845203816633372,1.9534532760666588 -Cohort,43,0.9827888610733586,0.4231078065833193,0.2797354836678798,1.911884289260335 -Cohort,44,0.9623947121578841,0.4136286255022229,0.2749647081219969,1.8703566980807 -Cohort,45,0.9421540655754211,0.4042158828082801,0.2702304010127748,1.8291066905776687 -Cohort,46,0.9221412285408568,0.3949077345137071,0.26554948321257377,1.7883109585695522 -Cohort,47,0.90216185597366,0.385607785576177,0.2608772804226187,1.7475310920479352 -Cohort,48,0.882334292132964,0.37637635483446563,0.25624079671797495,1.7070453621029646 -Cohort,49,0.8625657831801716,0.36716831449218124,0.25161860230331257,1.6666506542618256 -Cohort,50,0.8429997933271356,0.35805581191325264,0.24704355291097338,1.6266775475558584 -Cohort,51,0.8234786674703262,0.3489617193077022,0.24247927605981223,1.586777852655728 -Cohort,52,0.8041534436979831,0.3399614247273939,0.23796043159954386,1.5472956747805133 -Cohort,53,0.7848973039454056,0.33099333363912453,0.23345770579373523,1.507954070692058 -Cohort,54,0.7657797396972051,0.3220923711846878,0.22898701492150123,1.468913087611602 -Cohort,55,0.7467903076451918,0.3132541112351842,0.2245458639751767,1.4301545073178001 -Cohort,56,0.7278355539958338,0.30443390314304786,0.22011255906853136,1.3914794520703906 -Cohort,57,0.70897396693391,0.2956610504331953,0.21570050662101664,1.3530220460985172 -Cohort,58,0.6902152306393208,0.28694094074925197,0.2113118584463951,1.314808343019425 -Cohort,59,0.6715949804628504,0.2782913232731402,0.20695479465643865,1.276919394454409 -Cohort,60,0.6529563747594994,0.26963752637995103,0.2025928736960714,1.2390232832441421 -Cohort,61,0.6345903092168211,0.26111950039867,0.1982934973241661,1.2017461947989823 -Cohort,62,0.6163955122010181,0.2526888219799046,0.19403316071932997,1.1648723927669504 -Cohort,63,0.5984706763273917,0.24439252871907244,0.1898347857706342,1.1286112655564149 -Cohort,64,0.5806525916311895,0.23615313427260107,0.18566042423793358,1.0926189413604745 -Cohort,65,0.5629755938665605,0.22798724420337071,0.1815180198471186,1.056970021129552 -Cohort,66,0.545446982989119,0.2198985269915308,0.17740924338887865,1.021681395843492 -Cohort,67,0.5278422889398746,0.2117817666743538,0.1732817336736266,0.9862896624841411 -Cohort,68,0.5105497163950123,0.20381973836853745,0.16922594192412682,0.9516021241298777 -Cohort,69,0.4929301838506308,0.19571348432604257,0.1650927183820824,0.9163035837918675 -Cohort,70,0.4757098534049711,0.1878030953209463,0.16105148701148245,0.8818912237531586 -Cohort,71,0.45868455423154764,0.17999287372088915,0.1570546238436607,0.8479434287580013 -Cohort,72,0.4416703571371421,0.1721974798898497,0.1530591176146059,0.8140865179732711 -Cohort,73,0.4249300876127187,0.16453894592301782,0.14912642758510272,0.7808545837452481 -Cohort,74,0.4083100875092608,0.15694602834723526,0.14522061856932106,0.7479355998814706 -Cohort,75,0.3919798832524217,0.1494968627402586,0.14138141096853707,0.7156696546812567 -Cohort,76,0.37569195525618304,0.1420774295873372,0.1375508259286675,0.6835593491030069 -Cohort,77,0.3592289058066633,0.13458859134289253,0.13367781972717088,0.6511746176913152 -Cohort,78,0.343509791260289,0.12745131129948703,0.12997802959012805,0.6203421972852867 -Cohort,79,0.3274914788294062,0.1201883642608726,0.1262066812123317,0.5889901745347922 -Cohort,80,0.3120068273716592,0.11318002411140045,0.12255934717463676,0.5587649473328322 -Cohort,81,0.29691452784868966,0.10636116536735063,0.11879478804504452,0.5293808544966433 -Cohort,82,0.2818898604472989,0.09958430382663175,0.11488829438202917,0.5001981066756073 -Cohort,83,0.26732992829914204,0.09302937170292258,0.11109665634701975,0.4719904124459059 -Cohort,84,0.253469816555173,0.08680175232929532,0.10748151687804182,0.44520717018197187 -Cohort,85,0.23984012005589433,0.0806892726060211,0.10392152300256963,0.4189295535671334 -Cohort,86,0.2270252880434281,0.07495447252501987,0.10056926623995394,0.3942830081141291 -Cohort,87,0.2148500284592758,0.06951762169122021,0.09737988183901913,0.3709187960180863 -Cohort,88,0.20363899780774777,0.06452295260821458,0.09443892978343728,0.34944263470617937 -Cohort,89,0.1931919138320061,0.059879627501292026,0.09169490341914917,0.32945984076764745 -Cohort,90,0.1830810247281454,0.05539674604750873,0.08903631594809414,0.309605232002492 -Cohort,91,0.17427988013690718,0.05150582925847875,0.08671894371183428,0.2920951791637816 -Cohort,92,0.16619881829461683,0.047943431040790366,0.08458886906549214,0.2760332714344278 -Cohort,93,0.15790941208349174,0.04430008793797959,0.08240238269539632,0.25956942464864324 -Cohort,94,0.15032862700715366,0.04098085717699488,0.08040056523643792,0.2445271888495746 -Cohort,95,0.1438279093043865,0.03814662164357434,0.07870824235268185,0.2316405012560124 -Cohort,96,0.13811708389006946,0.035667793013521513,0.07769236272148212,0.22032918043453312 -Cohort,97,0.13306987501337675,0.03348738444845809,0.07754401594075466,0.21034005826593719 -Cohort,98,0.12871242065672317,0.03161457213092548,0.07741316000007889,0.20169342933163853 -Cohort,99,0.12408255000249598,0.02963532449046922,0.07680689294277858,0.19251126605173982 -Cohort,100,0.12142994449731342,0.02850972547340673,0.07614048736323602,0.1872601947877508 -Cohort,101,0.11597101071787298,0.026205693168365813,0.07478969317975462,0.1767571662749944 -Cohort,102,0.11294156852830967,0.024942114920838813,0.07404012326847498,0.17100162514167888 -Cohort,103,0.1069928537902417,0.022493403391879805,0.07226963787464409,0.1591262770087007 -Cohort,104,0.09919085502183704,0.019391897247520804,0.06928474483606596,0.14353838413445247 -Cohort,105,0.09038069302967518,0.01614804541063871,0.06496914874808764,0.12531369735717118 -Cohort,106,0.063693491517932,0.009205741086214874,0.048807087035780905,0.0826411756936002 -ICU,0,3.275433421012084,1.0422986814322854,1.3067409069197538,5.421511283325188 -ICU,1,3.233787088118184,1.0298527951467984,1.2904461068261395,5.355999064535664 -ICU,2,3.1868891307012395,1.0153084953113007,1.2721531662102865,5.2802124931694685 -ICU,3,3.1399798397170517,1.0007312951480236,1.2538589922493373,5.204291094989306 -ICU,4,3.09333401448425,0.9862245495383478,1.2356688444629402,5.128747808347906 -ICU,5,3.0468466765019087,0.9717420612057341,1.217543229299834,5.053361898972334 -ICU,6,3.0007113318448306,0.9573571438734285,1.1995562061251843,4.978496516485738 -ICU,7,2.9548446429921795,0.9430336758743435,1.1816763541059825,4.9039787316921695 -ICU,8,2.909307483531888,0.928793373361041,1.1639271082571585,4.829917989676072 -ICU,9,2.8641448716172935,0.9146520523076512,1.1463258239099419,4.756394374868528 -ICU,10,2.8193160493345744,0.9005917527901586,1.1288571540894088,4.683322579190145 -ICU,11,2.774830472169128,0.8866147899811885,1.1115248394416417,4.6107160720165545 -ICU,12,2.730790027572505,0.8727604724902306,1.0943678010996958,4.538768635265736 -ICU,13,2.687065823606213,0.8589766221687234,1.0773370311598496,4.467227076217114 -ICU,14,2.6438487537480073,0.8453383044454391,1.060505317361138,4.396459175145917 -ICU,15,2.6009235366644,0.8317594971230907,1.0437906861470498,4.326046144706561 -ICU,16,2.558533081084677,0.8183361890694237,1.0272857245175742,4.256457316152915 -ICU,17,2.516477411192994,0.8049888094817452,1.0109142325909786,4.187305249370656 -ICU,18,2.4748681196102633,0.7917615813271468,0.9947187460798419,4.118805950473033 -ICU,19,2.433834222800422,0.778705484332107,0.9787484451742399,4.051208572225277 -ICU,20,2.393047341443447,0.7656909479458537,0.9628781048793219,3.983880820741737 -ICU,21,2.3526453446335323,0.7527728587520718,0.947160235958398,3.9170905013687003 -ICU,22,2.3128255314441843,0.7400292777101872,0.9316700466406824,3.8512186361470593 -ICU,23,2.2732443151771835,0.7273255400625898,0.9162764119997671,3.7856072226534345 -ICU,24,2.234037282525568,0.7147157552519863,0.9010309656243166,3.720519907339323 -ICU,25,2.1952061843338706,0.7022014139709354,0.8859342859550051,3.6559632563295135 -ICU,26,2.1566689377171953,0.6897513835527506,0.8709549338316566,3.591784119833823 -ICU,27,2.1186329874269623,0.6774475758259939,0.8561720111599336,3.5283821705829026 -ICU,28,2.0808684110233466,0.6652014087378991,0.8414976076797891,3.465322905449102 -ICU,29,2.0434748408175296,0.6530528454045613,0.8269696499933459,3.4028005908481314 -ICU,30,2.006432863578063,0.640995534563172,0.8125805971529351,3.340782762990689 -ICU,31,1.9697935433724547,0.6290503845165011,0.7983498443196122,3.279370530252542 -ICU,32,1.9333891291004068,0.617154094383938,0.7842131156787004,3.2182517386446947 -ICU,33,1.8972694727853505,0.6053274515106322,0.7701893088877889,3.1575264163412045 -ICU,34,1.8614966512510955,0.5935957054167642,0.7563020220840622,3.0973167567964524 -ICU,35,1.825951686914155,0.5819149645839314,0.7425055664246447,3.0374049645307375 -ICU,36,1.790732832414406,0.570324271924973,0.7288373877617351,2.9779810307439347 -ICU,37,1.7557486713170465,0.55879053562486,0.7152623179668843,2.91887985793659 -ICU,38,1.721100015309572,0.5473534730453642,0.7018188139770635,2.8602951119678344 -ICU,39,1.6866823589680655,0.535975265776457,0.6884666720549523,2.802038062716486 -ICU,40,1.6526265370219468,0.524706547493263,0.6752558946700913,2.7443567570990797 -ICU,41,1.6187881545407787,0.5134947797935351,0.6621309616078426,2.6869893836323304 -ICU,42,1.5850821247156535,0.5023101873854224,0.6490590384628463,2.629785876751285 -ICU,43,1.5517499690148973,0.491244040948725,0.6361326510407825,2.573196262802215 -ICU,44,1.5185093177519797,0.4801935018955296,0.6232432381054523,2.516708171460736 -ICU,45,1.485519676136697,0.46921918621817377,0.6104518715099517,2.4606200158237135 -ICU,46,1.452901633187132,0.45836652634588043,0.5978047617778415,2.4051560161731635 -ICU,47,1.4203392918988544,0.4475215875254002,0.5851803465263112,2.3497469130513102 -ICU,48,1.3880247734155111,0.4367560786718274,0.5726523150432281,2.2947475378778655 -ICU,49,1.355807181374196,0.426016885473685,0.560162468725077,2.2398906377333687 -ICU,50,1.3239195604932408,0.4153895305648442,0.5478003384040416,2.1856014667389436 -ICU,51,1.2921054929670044,0.40478312575463476,0.5354670904438384,2.131423633186151 -ICU,52,1.2606104099983937,0.3942868875092313,0.5232570985198185,2.0778019982926392 -ICU,53,1.2292279818738887,0.3838283457073078,0.5110907676567884,2.024371714365454 -ICU,54,1.1980710927749136,0.3734488863831615,0.49901146734246965,1.9713387091689396 -ICU,55,1.1671226539656867,0.36314346116280144,0.48701250719728983,1.9186762600303353 -ICU,56,1.1362305079557125,0.352859722767216,0.47503507856478056,1.866119275806007 -ICU,57,1.1054896771229643,0.3426323652523603,0.4631157010067975,1.8138406567331462 -ICU,58,1.0749158128551488,0.3324678895259795,0.4512603073760109,1.761871815753623 -ICU,59,1.0445668238884978,0.3223872737316145,0.43949116786151643,1.7103176134432398 -ICU,60,1.0141873082192914,0.3123030421263013,0.427709530095791,1.6587344967292792 -ICU,61,0.9842507369152343,0.3023794694648013,0.41609823792170325,1.6079527855110332 -ICU,62,0.9545922285985832,0.29255976330005573,0.4045935696229097,1.5576850521962349 -ICU,63,0.9253724832486121,0.28289900459115946,0.3932576494915183,1.5082107853712596 -ICU,64,0.8963256963689032,0.2733064671825561,0.3819876703819985,1.4590695247477483 -ICU,65,0.8675077390375846,0.2638016367688965,0.3708051937962681,1.4103597757410713 -ICU,66,0.8389304676325362,0.25438882619166814,0.35971477332081314,1.3621032270714215 -ICU,67,0.810228106781539,0.2449451842273646,0.3485747311196668,1.3136734957722387 -ICU,68,0.782033112837798,0.23568422204875752,0.3376299146205452,1.266158124021099 -ICU,69,0.7533040110319792,0.2262570749160371,0.32647683467453725,1.217776758079099 -ICU,70,0.7252241268923965,0.21706059584919785,0.31557387722947944,1.170554351932331 -ICU,71,0.6974607652908559,0.20798301888439594,0.3047921801443294,1.1239210816547427 -ICU,72,0.6697140844398151,0.19892487008836807,0.29401547507164716,1.0773680593103339 -ICU,73,0.6424125398835316,0.1900282318763242,0.2834099000357649,1.0316225470465537 -ICU,74,0.6153056344533102,0.18121003892078194,0.27287831755276026,0.9862595213446891 -ICU,75,0.5886698220115667,0.1725610794856738,0.26252801975096496,0.9417448838245412 -ICU,76,0.5621014794181617,0.16394860696637917,0.2522023783159424,0.8973977729179602 -ICU,77,0.5352459579591098,0.15525736952188923,0.24176362434929727,0.8526251110875377 -ICU,78,0.5096023631730455,0.1469763572163862,0.2317939254452188,0.8099403955746522 -ICU,79,0.48346917299498315,0.13855093386079564,0.22163247633750083,0.7664918392146394 -ICU,80,0.4582050844567361,0.1304225242431886,0.211807110474894,0.7245506745665453 -ICU,81,0.4335797048635797,0.12251512456953399,0.2022284589119513,0.6837269691756505 -ICU,82,0.4090632523392228,0.11465723876875225,0.19269066420741715,0.6431367351432371 -ICU,83,0.38530378667159576,0.10705728014987795,0.18344572906542236,0.6038547114792229 -ICU,84,0.36268514271917546,0.09983702564122204,0.17464309493848365,0.5665108884600095 -ICU,85,0.34044135578923657,0.0927498512260305,0.16598500183081974,0.529831769897134 -ICU,86,0.3195265107850452,0.08609980303351633,0.1578427361997582,0.49538951621707 -ICU,87,0.29965468929370603,0.07979388905920894,0.1501052705601531,0.4627044924361333 -ICU,88,0.2813560873346568,0.07399904068739432,0.1429791442453401,0.4326428593485637 -ICU,89,0.2643039612793491,0.06860951344543324,0.13633741145242034,0.40465881116099905 -ICU,90,0.24780008166810003,0.06340319817580072,0.12990840463504744,0.3775999587299512 -ICU,91,0.23343415097219627,0.05888124580831036,0.12431118732719752,0.354070771953404 -ICU,92,0.22024360887826472,0.054737631186920106,0.11917118593451224,0.33248470741001046 -ICU,93,0.20671258454569222,0.050495303478891,0.1138981512437593,0.3098399229315839 -ICU,94,0.19433859188431016,0.046625249528626594,0.10907530058381415,0.2890669263220663 -ICU,95,0.1837285698761847,0.04331569407059615,0.1047849892795412,0.2712575121830722 -ICU,96,0.17440915433351775,0.04041635617417341,0.10075156731628923,0.25561622500137343 -ICU,97,0.1661747090699182,0.037861461919284904,0.09719003830120573,0.24179743272827448 -ICU,98,0.15906842025853898,0.035662807415821736,0.09411435312387309,0.22987335209895202 -ICU,99,0.15152149225198994,0.03333407318730086,0.09084685516598558,0.21721066842193773 -ICU,100,0.1472028402001833,0.03200702516852799,0.08897402202927242,0.20996716005425145 -ICU,101,0.1383213720634288,0.02928308295536382,0.0851255090729158,0.19506833404668852 -ICU,102,0.13340564623586215,0.027784293525836744,0.0829913316347619,0.18701473239133207 -ICU,103,0.12378784721192897,0.0248675613774783,0.07986996556033245,0.17138533698432665 -ICU,104,0.11132918022875264,0.021145161724653045,0.07345449361399496,0.15117620018523592 -ICU,105,0.09773199451387127,0.017214157978992305,0.06713233148647271,0.12926681939238052 -ICU,106,0.06338673047498927,0.00897350431789788,0.04768255774619171,0.08022856831025232 -Non-hospitalised (no AD),0,0.06018970662432619,0.0029729533081749764,0.05458845828834137,0.06621683529190213 -Non-hospitalised (no AD),1,0.06018887377052802,0.0029729533081749734,0.05458762543454323,0.06621600243810398 -Non-hospitalised (no AD),2,0.06018761788132908,0.002972953308174975,0.05458636954534426,0.06621474654890504 -Non-hospitalised (no AD),3,0.06018606558792035,0.0029729533081749734,0.05458481725193558,0.06621319425549631 -Non-hospitalised (no AD),4,0.060184270466377525,0.0029729533081749742,0.05458302213039274,0.06621139913395348 -Non-hospitalised (no AD),5,0.06018226478015189,0.002972953308174975,0.05458101644416707,0.06620939344772783 -Non-hospitalised (no AD),6,0.06018007074524419,0.002972953308174975,0.05457882240925936,0.06620719941282015 -Non-hospitalised (no AD),7,0.06017770486662893,0.0029729533081749742,0.05457645653064412,0.0662048335342049 -Non-hospitalised (no AD),8,0.06017518003040379,0.0029729533081749742,0.05457393169441901,0.06620230869797973 -Non-hospitalised (no AD),9,0.06017250665743683,0.002972953308174973,0.05457125832145206,0.06619963532501279 -Non-hospitalised (no AD),10,0.06016969339957393,0.0029729533081749742,0.05456844506358912,0.0661968220671499 -Non-hospitalised (no AD),11,0.06016674758828379,0.0029729533081749755,0.05456549925229897,0.06619387625585975 -Non-hospitalised (no AD),12,0.060163675538624714,0.0029729533081749747,0.054562427202639936,0.06619080420620067 -Non-hospitalised (no AD),13,0.06016048276351087,0.0029729533081749755,0.05455923442752605,0.06618761143108683 -Non-hospitalised (no AD),14,0.06015717412968474,0.002972953308174973,0.05455592579369996,0.06618430279726069 -Non-hospitalised (no AD),15,0.060153753974322585,0.0029729533081749755,0.05455250563833781,0.06618088264189854 -Non-hospitalised (no AD),16,0.06015022619419019,0.002972953308174976,0.05454897785820538,0.06617735486176615 -Non-hospitalised (no AD),17,0.06014659431513313,0.0029729533081749742,0.054545345979148306,0.06617372298270911 -Non-hospitalised (no AD),18,0.06014286154714716,0.002972953308174976,0.05454161321116234,0.06616999021472311 -Non-hospitalised (no AD),19,0.060139030828658596,0.002972953308174975,0.05453778249267379,0.06616615949623456 -Non-hospitalised (no AD),20,0.06013510486258882,0.002972953308174975,0.054533856526604044,0.06616223353016477 -Non-hospitalised (no AD),21,0.060131086146061996,0.002972953308174973,0.05452983781007721,0.06615821481363794 -Non-hospitalised (no AD),22,0.06012697699512739,0.0029729533081749734,0.05452572865914258,0.06615410566270334 -Non-hospitalised (no AD),23,0.0601227795655218,0.002972953308174973,0.054521531229537006,0.06614990823309774 -Non-hospitalised (no AD),24,0.06011849587025124,0.0029729533081749742,0.05451724753426645,0.0661456245378272 -Non-hospitalised (no AD),25,0.06011412779459189,0.002972953308174974,0.054512879458607104,0.06614125646216783 -Non-hospitalised (no AD),26,0.06010967710897696,0.0029729533081749764,0.05450842877299215,0.06613680577655293 -Non-hospitalised (no AD),27,0.0601051454801377,0.0029729533081749734,0.054503897144152914,0.06613227414771365 -Non-hospitalised (no AD),28,0.06010053448079093,0.0029729533081749742,0.05449928614480614,0.06612766314836688 -Non-hospitalised (no AD),29,0.060095845598107794,0.002972953308174973,0.054494597262123,0.06612297426568374 -Non-hospitalised (no AD),30,0.06009108024115403,0.002972953308174973,0.054489831905169255,0.06611820890872998 -Non-hospitalised (no AD),31,0.06008623974745676,0.0029729533081749725,0.05448499141147197,0.06611336841503272 -Non-hospitalised (no AD),32,0.060081325388824956,0.002972953308174977,0.054480077052840144,0.06610845405640091 -Non-hospitalised (no AD),33,0.060076338376529544,0.0029729533081749725,0.05447509004054478,0.0661034670441055 -Non-hospitalised (no AD),34,0.06007127986593051,0.002972953308174974,0.05447003152994574,0.06609840853350644 -Non-hospitalised (no AD),35,0.06006615096062477,0.002972953308174975,0.05446490262463997,0.06609327962820073 -Non-hospitalised (no AD),36,0.060060952716176584,0.0029729533081749755,0.05445970438019177,0.06608808138375252 -Non-hospitalised (no AD),37,0.06005568614348299,0.0029729533081749764,0.054454437807498215,0.06608281481105895 -Non-hospitalised (no AD),38,0.0600503522118191,0.0029729533081749794,0.05444910387583428,0.06607748087939506 -Non-hospitalised (no AD),39,0.06004495185160074,0.0029729533081749725,0.054443703515615983,0.06607208051917669 -Non-hospitalised (no AD),40,0.060039485956898134,0.002972953308174975,0.054438237620913314,0.0660666146244741 -Non-hospitalised (no AD),41,0.06003395538772788,0.0029729533081749755,0.05443270705174306,0.06606108405530385 -Non-hospitalised (no AD),42,0.060028360972148025,0.002972953308174977,0.0544271126361632,0.06605548963972398 -Non-hospitalised (no AD),43,0.06002270350817776,0.0029729533081749742,0.054421455172192976,0.0660498321757537 -Non-hospitalised (no AD),44,0.060016983765559556,0.0029729533081749742,0.05441573542957477,0.0660441124331355 -Non-hospitalised (no AD),45,0.06001120248738046,0.002972953308174976,0.054409954151395636,0.06603833115495643 -Non-hospitalised (no AD),46,0.060005360391566474,0.0029729533081749734,0.05440411205558169,0.06603248905914244 -Non-hospitalised (no AD),47,0.05999945817226218,0.002972953308174975,0.0543982098362774,0.06602658683983813 -Non-hospitalised (no AD),48,0.05999349650110758,0.0029729533081749742,0.054392248165122795,0.06602062516868354 -Non-hospitalised (no AD),49,0.0599874760284208,0.002972953308174974,0.05438622769243602,0.06601460469599676 -Non-hospitalised (no AD),50,0.05998139738429579,0.002972953308174975,0.054380149048310976,0.06600852605187174 -Non-hospitalised (no AD),51,0.05997526117962301,0.0029729533081749716,0.05437401284363822,0.06600238984719899 -Non-hospitalised (no AD),52,0.059969068007039106,0.0029729533081749725,0.054367819671054314,0.06599619667461504 -Non-hospitalised (no AD),53,0.059962818441812976,0.002972953308174975,0.05436157010582817,0.06598994710938892 -Non-hospitalised (no AD),54,0.0599565130426722,0.0029729533081749734,0.05435526470668741,0.06598364171024815 -Non-hospitalised (no AD),55,0.05995015235257605,0.0029729533081749734,0.05434890401659127,0.065977281020152 -Non-hospitalised (no AD),56,0.05994373689943891,0.0029729533081749734,0.054342488563454125,0.06597086556701486 -Non-hospitalised (no AD),57,0.0599372671968079,0.0029729533081749742,0.054336018860823114,0.06596439586438387 -Non-hospitalised (no AD),58,0.0599307437444992,0.0029729533081749725,0.05432949540851441,0.06595787241207517 -Non-hospitalised (no AD),59,0.059924167029195,0.002972953308174975,0.054322918693210175,0.06595129569677095 -Non-hospitalised (no AD),60,0.05991753752500556,0.0029729533081749742,0.05431628918902077,0.0659446661925815 -Non-hospitalised (no AD),61,0.05991085569399765,0.002972953308174973,0.054309607358012865,0.06593798436157361 -Non-hospitalised (no AD),62,0.0599041219866929,0.0029729533081749695,0.054302873650708135,0.06593125065426884 -Non-hospitalised (no AD),63,0.05989733684253742,0.0029729533081749725,0.05429608850655263,0.06592446551011337 -Non-hospitalised (no AD),64,0.059890500690345856,0.002972953308174974,0.054289252354361064,0.06591762935792182 -Non-hospitalised (no AD),65,0.059883613948719924,0.002972953308174976,0.054282365612735105,0.06591074261629588 -Non-hospitalised (no AD),66,0.05987667702644561,0.002972953308174974,0.05427542869046083,0.06590380569402159 -Non-hospitalised (no AD),67,0.059869690322867805,0.002972953308174976,0.05426844198688298,0.06589681899044375 -Non-hospitalised (no AD),68,0.05986265422824638,0.0029729533081749747,0.0542614058922616,0.06588978289582233 -Non-hospitalised (no AD),69,0.05985556912409298,0.0029729533081749734,0.054254320788108205,0.06588269779166894 -Non-hospitalised (no AD),70,0.05984843538349119,0.0029729533081749742,0.05424718704750637,0.06587556405106713 -Non-hospitalised (no AD),71,0.059841253371400366,0.002972953308174977,0.05424000503541553,0.06586838203897633 -Non-hospitalised (no AD),72,0.05983402344494412,0.0029729533081749725,0.05423277510895936,0.06586115211252007 -Non-hospitalised (no AD),73,0.059826745953685354,0.0029729533081749742,0.05422549761770056,0.06585387462126131 -Non-hospitalised (no AD),74,0.05981942123988732,0.002972953308174975,0.05421817290390251,0.06584654990746329 -Non-hospitalised (no AD),75,0.059812049638762674,0.002972953308174977,0.05421080130277787,0.06583917830633865 -Non-hospitalised (no AD),76,0.059804631478710825,0.0029729533081749742,0.05420338314272604,0.06583176014628678 -Non-hospitalised (no AD),77,0.059797167081543964,0.0029729533081749716,0.05419591874555917,0.06582429574911992 -Non-hospitalised (no AD),78,0.059789656762703064,0.0029729533081749708,0.05418840842671828,0.065816785430279 -Non-hospitalised (no AD),79,0.059782100831463666,0.002972953308174941,0.054180852495478944,0.06580922949903956 -Non-hospitalised (no AD),80,0.05977449959113101,0.002972953308174721,0.05417325125514671,0.06580162825870646 -Non-hospitalised (no AD),81,0.05976685333921586,0.0029729533081732334,0.05416560500323431,0.06579398200678828 -Non-hospitalised (no AD),82,0.0597591623675301,0.002972953308163265,0.05415791403156734,0.06578629103508232 -Non-hospitalised (no AD),83,0.059751426961871974,0.0029729533081019432,0.05415017862602478,0.06577855562929986 -Non-hospitalised (no AD),84,0.05974364739973688,0.002972953307764127,0.05414239906452613,0.06577077606647992 -Non-hospitalised (no AD),85,0.05973582393718804,0.002972953305958905,0.05413457560537844,0.0657629526002713 -Non-hospitalised (no AD),86,0.05972795675695905,0.002972953297509701,0.05412670844106837,0.06575508540291308 -Non-hospitalised (no AD),87,0.05972004575945764,0.002972953261262115,0.05411879751185992,0.0657471743319262 -Non-hospitalised (no AD),88,0.059712089969221845,0.0029729531265658047,0.05411084197540122,0.06573921826861784 -Non-hospitalised (no AD),89,0.05970408573562444,0.00297295267262083,0.05410283859706729,0.06573121311472857 -Non-hospitalised (no AD),90,0.059696020591985074,0.002972951187685736,0.054094776251147704,0.06572314496065011 -Non-hospitalised (no AD),91,0.059687869941636065,0.0029729472985144012,0.05408663292826477,0.06571498642570497 -Non-hospitalised (no AD),92,0.05967957475106672,0.0029729377495929625,0.05407835572851964,0.06570667187641314 -Non-hospitalised (no AD),93,0.05967094591317844,0.002972912202116854,0.05406977502383197,0.06569799124560728 -Non-hospitalised (no AD),94,0.05966167172626026,0.0029728538342118083,0.054060610806057986,0.06568859872824886 -Non-hospitalised (no AD),95,0.059651363725680875,0.002972741917394383,0.05405051366445364,0.06567806383643075 -Non-hospitalised (no AD),96,0.05963942326959273,0.0029725446274290663,0.054038944916218044,0.06566572341038598 -Non-hospitalised (no AD),97,0.059625116371827866,0.002972223212821139,0.05402524358566234,0.06565076490225365 -Non-hospitalised (no AD),98,0.05960792694467467,0.0029717508566120036,0.0540089441100678,0.06563261785775693 -Non-hospitalised (no AD),99,0.059583799282003476,0.00297091279412954,0.05398639541470578,0.06560679117398124 -Non-hospitalised (no AD),100,0.05956190791472424,0.002970197144817233,0.05396585237991085,0.06558344895630543 -Non-hospitalised (no AD),101,0.05950867054635202,0.0029678294277307865,0.05391707595324063,0.0655254114669595 -Non-hospitalised (no AD),102,0.0594606866939578,0.002965747219059815,0.05387301512507332,0.06547320631060528 -Non-hospitalised (no AD),103,0.059321798493173546,0.00295891368937816,0.05374700176381414,0.06532046435620446 -Non-hospitalised (no AD),104,0.05894714438734247,0.0029398748113338815,0.05340821821258321,0.06490721234594495 -Non-hospitalised (no AD),105,0.05797359909115071,0.0028901491209288948,0.052528359534297867,0.06383285714690984 -Non-hospitalised (no AD),106,0.04749747381282853,0.002360538574786213,0.04305005692010718,0.05228304164394481 +Non-hospitalised,0,0.7579775022370616,0.28888480884894013,0.2193688877419769,1.316145814375487 +Non-hospitalised,1,0.7477034333094617,0.28478853976382107,0.21682077747815068,1.2983694284227385 +Non-hospitalised,2,0.7366895040665494,0.28037014528335724,0.2140989685950534,1.2794770922285077 +Non-hospitalised,3,0.7257050484029705,0.2759621067446168,0.21138492324329466,1.2605930945618269 +Non-hospitalised,4,0.7147956045555821,0.27158366599007816,0.20868955820823992,1.2412222307793834 +Non-hospitalised,5,0.7039506351997731,0.26722984281022816,0.20601054899958082,1.2219620138777083 +Non-hospitalised,6,0.6932016579847579,0.26291400014348404,0.2033554302601771,1.2028552634750616 +Non-hospitalised,7,0.682539666144833,0.2586319572372425,0.20072218223291327,1.1838978117023649 +Non-hospitalised,8,0.6719758975932972,0.2543883715175261,0.1981135294173865,1.165112502176015 +Non-hospitalised,9,0.6615188979870361,0.2501867674868043,0.1955315494908017,1.1465147544584238 +Non-hospitalised,10,0.6511644612958093,0.24602516081193182,0.1929753077998541,1.1280964581436752 +Non-hospitalised,11,0.6409153046936713,0.2419046019334466,0.1904454881229699,1.1098623795952467 +Non-hospitalised,12,0.6307873102945568,0.23783188079292433,0.1879458658841829,1.0918416861950389 +Non-hospitalised,13,0.620762600172943,0.23379914894611495,0.18547225929266378,1.0740011460980938 +Non-hospitalised,14,0.6108696326129516,0.22981869117309373,0.18303139381871145,1.0563932562809712 +Non-hospitalised,15,0.6010774130289339,0.2258770075209276,0.18061598127997092,1.0389606320928868 +Non-hospitalised,16,0.5914217924660746,0.22198961341455853,0.17823448255874394,1.0217694763936627 +Non-hospitalised,17,0.5818735299128102,0.21814380778980655,0.17588000787368338,1.0048089495537218 +Non-hospitalised,18,0.5724490035359577,0.2143466826896506,0.17355642022900877,0.9880260060764708 +Non-hospitalised,19,0.5631672892339695,0.21060649316925117,0.17126822294375948,0.9714465121043273 +Non-hospitalised,20,0.5539793113519541,0.20690201208119638,0.1690038192925125,0.9550757317548536 +Non-hospitalised,21,0.544905025224905,0.20324191031761182,0.1667679071317435,0.938904258503473 +Non-hospitalised,22,0.5359736034799935,0.19963882246008272,0.1645673724472347,0.9229859532762247 +Non-hospitalised,23,0.5271327113579214,0.19607018905392973,0.16238981891634474,0.9072244751799253 +Non-hospitalised,24,0.5184018985443212,0.19254451689103314,0.16023984641303712,0.8916560115549436 +Non-hospitalised,25,0.5097805429238822,0.18906160449074122,0.15811728259353403,0.8762794648822291 +Non-hospitalised,26,0.5012550232036221,0.1856156669637362,0.1560188691127499,0.8610424826856522 +Non-hospitalised,27,0.49285627342787464,0.18222008883341223,0.15395191518997148,0.8460320047258585 +Non-hospitalised,28,0.48454766471643373,0.17885920974150213,0.15190769810166838,0.831182446341692 +Non-hospitalised,29,0.47634341596210134,0.17553923535437474,0.14988955894076803,0.8165192838348921 +Non-hospitalised,30,0.46823926769391644,0.1722584353489049,0.1478964515288434,0.802034902740136 +Non-hospitalised,31,0.4602420815356352,0.16901984277300808,0.14592997525468096,0.7877416020963444 +Non-hospitalised,32,0.4523237917835855,0.165811566115682,0.14398340944056517,0.7735891443674671 +Non-hospitalised,33,0.44449073312156157,0.1626364415469709,0.14205821866718896,0.7595888888054135 +Non-hospitalised,34,0.43675144652844344,0.15949822863555002,0.14015640010792452,0.7457561422335065 +Non-hospitalised,35,0.4290850327487835,0.1563881504417132,0.13827292300239438,0.732053505295075 +Non-hospitalised,36,0.4215059524699665,0.1533124860558209,0.13641119863388385,0.7185068823828118 +Non-hospitalised,37,0.4139975280051523,0.15026427848980836,0.13456719689951618,0.7050864347421745 +Non-hospitalised,38,0.40657494858779597,0.14725009329509983,0.13274451144872976,0.6918193677652814 +Non-hospitalised,39,0.39921914724812296,0.1442619780567181,0.13093853468671962,0.6786715680427039 +Non-hospitalised,40,0.39195074663529095,0.14130876853425164,0.12915417226733777,0.6656799624296605 +Non-hospitalised,41,0.3847436902138929,0.13837957714787613,0.12738513435227894,0.6527979296636032 +Non-hospitalised,42,0.3775814617717799,0.13546758836318984,0.12562740572281753,0.6399959294331123 +Non-hospitalised,43,0.3705043031809652,0.13258986244696605,0.12389061949374042,0.6273459951716536 +Non-hospitalised,44,0.3634613672401434,0.12972514658840556,0.12216250381043307,0.614757148651876 +Non-hospitalised,45,0.3564788975072765,0.12688459344262704,0.12044933212398218,0.6022763695362174 +Non-hospitalised,46,0.34957711730766944,0.12407676609536347,0.11875594909899276,0.5899398563187692 +Non-hospitalised,47,0.34269805287085436,0.12127751425710254,0.11706833609437149,0.5776438901711631 +Non-hospitalised,48,0.33587458433469075,0.11850070776398026,0.11539438702117699,0.5654473178608566 +Non-hospitalised,49,0.3290777311387782,0.1157343678724922,0.11372706319649638,0.5532972723314487 +Non-hospitalised,50,0.32234880115409503,0.11299580786180928,0.11207631855092276,0.541253219879782 +Non-hospitalised,51,0.31563916822630705,0.11026488564797803,0.11043035704040921,0.5292432364875093 +Non-hospitalised,52,0.30899316250455794,0.10756011843659058,0.10879988144634836,0.5173476634002866 +Non-hospitalised,53,0.30237099620310565,0.10486507143215004,0.10717522443606076,0.5054947779902725 +Non-hospitalised,54,0.29579269577287304,0.10218813598437827,0.10556121229735921,0.4937209305742565 +Non-hospitalised,55,0.2892540025263445,0.09952761865120444,0.10395678839327044,0.48201858263561714 +Non-hospitalised,56,0.2827244813540917,0.09687102303300218,0.10235453279976771,0.4703388016152199 +Non-hospitalised,57,0.27622111259272575,0.09422545420162237,0.1007585407715554,0.4587744180558228 +Non-hospitalised,58,0.26974588457155785,0.09159180315364762,0.09916927201556423,0.44725198075293293 +Non-hospitalised,59,0.26330928638515194,0.0889744476748782,0.09758925979882208,0.4357347777367092 +Non-hospitalised,60,0.256859841437498,0.08635227812413165,0.09600594560592046,0.42413574535596815 +Non-hospitalised,61,0.2504907612094512,0.08376366153778453,0.09444202169834702,0.41268991659440357 +Non-hospitalised,62,0.24416911224274704,0.08119507279376072,0.0928894600474171,0.401354451092161 +Non-hospitalised,63,0.23792717466624075,0.07865975559288817,0.09135613377038791,0.39019244096508326 +Non-hospitalised,64,0.23171105564687158,0.07613563142662043,0.0898288878365925,0.37907807177709024 +Non-hospitalised,65,0.2255316107756768,0.07362717883397091,0.08831036193808389,0.3680309003198074 +Non-hospitalised,66,0.21939092823077505,0.07113527497916644,0.08680105906226089,0.3570547192326378 +Non-hospitalised,67,0.2132128425574979,0.06862884848196119,0.08528234377441793,0.34601303984297044 +Non-hospitalised,68,0.20712781985860768,0.06616120633733806,0.08378611154357446,0.3351398213613452 +Non-hospitalised,69,0.20091812214148408,0.06364357679906806,0.08225906954436123,0.32404499287934885 +Non-hospitalised,70,0.19483056381786065,0.06117663278844665,0.08076161527704322,0.31317079162960954 +Non-hospitalised,71,0.1887958927753787,0.05873212710850711,0.0792768053273009,0.30239311434627664 +Non-hospitalised,72,0.18275040845184173,0.05628416143468042,0.0777890272156154,0.29159797778316654 +Non-hospitalised,73,0.17678511283368767,0.05386976074161861,0.07632058461614102,0.28094822012086595 +Non-hospitalised,74,0.1708467563336155,0.05146728498879681,0.07485842584254938,0.27034856286688036 +Non-hospitalised,75,0.1649949815385959,0.049100952419130046,0.07341718798424395,0.2599056092423906 +Non-hospitalised,76,0.15914289111742766,0.046735517349299045,0.07197554926774397,0.24949993558256997 +Non-hospitalised,77,0.15321270578534552,0.0443395340165216,0.07051713403033635,0.2389798399686725 +Non-hospitalised,78,0.1475314001937873,0.04204547299609932,0.06911959450116371,0.2288387146724407 +Non-hospitalised,79,0.1417275507986476,0.03970296539156606,0.06769169046226524,0.21846762264196323 +Non-hospitalised,80,0.13609939557078268,0.03743273338673544,0.0663065725004043,0.20840915012032063 +Non-hospitalised,81,0.1305976883398977,0.035214826689420616,0.06495218726611511,0.19856256957654497 +Non-hospitalised,82,0.1251056983507864,0.033002135099350696,0.06378370812965053,0.18876440877485512 +Non-hospitalised,83,0.1197680416262691,0.030853094659190936,0.06277057511589144,0.17922759660187532 +Non-hospitalised,84,0.11467218097683135,0.028802951369442817,0.061568942177260534,0.17012308367382062 +Non-hospitalised,85,0.10964808105141391,0.026783199442698233,0.05995941759949787,0.16115277101378003 +Non-hospitalised,86,0.1049114293434975,0.024880722528729738,0.05885885530789553,0.15273441185168599 +Non-hospitalised,87,0.10039993195561946,0.023070417916999594,0.057548005893160066,0.14471806295489745 +Non-hospitalised,88,0.09623551907108181,0.02140122588422868,0.056432998022881015,0.13731485793059994 +Non-hospitalised,89,0.09234636354057262,0.0198442046080194,0.055674445161078195,0.13041284104958725 +Non-hospitalised,90,0.0885751249772346,0.01833631298482043,0.0547847385874496,0.12373609902518219 +Non-hospitalised,91,0.08528533832254193,0.01702304601854587,0.054119299949936646,0.11788644246920728 +Non-hospitalised,92,0.08225951588827927,0.015817142628995817,0.05323987408756229,0.11251049654162558 +Non-hospitalised,93,0.0791516527589101,0.014580682316630015,0.0521884618956791,0.10699690672235655 +Non-hospitalised,94,0.07630480624984716,0.013450709149785706,0.051473583895182314,0.10194723376004973 +Non-hospitalised,95,0.07385959084438588,0.012482807289815536,0.050867733719886406,0.09786804798980114 +Non-hospitalised,96,0.07170863595689878,0.011633874643254993,0.05033424572411254,0.09423590551043558 +Non-hospitalised,97,0.0698054999833298,0.010885154983464948,0.049905260180351345,0.09094213699946166 +Non-hospitalised,98,0.06816078031451914,0.010240404309269886,0.0494470983454829,0.08820177439792518 +Non-hospitalised,99,0.06641292039045298,0.009557652381147797,0.04900698303658158,0.08532920260685599 +Non-hospitalised,100,0.06540916394629813,0.009167915482336326,0.048752368033411514,0.08352801068140461 +Non-hospitalised,101,0.06334867469205561,0.008370064307601555,0.04804166926414752,0.0798501291142984 +Non-hospitalised,102,0.06220336091913156,0.007930767621041498,0.047667861828669125,0.07777921705800625 +Non-hospitalised,103,0.05995954321844564,0.007078175407791552,0.04709923693486728,0.0739085001910159 +Non-hospitalised,104,0.05701963201194508,0.005994018448802171,0.0461707333533145,0.0686904730545087 +Non-hospitalised,105,0.053687240116727625,0.004858904360080288,0.0449245417655166,0.06327339722844395 +Non-hospitalised,106,0.04229270248742703,0.002574467740409443,0.03746403911184274,0.04707923350772385 +Hospitalised (no IC),0,2.087018731180328,0.8355287537222525,0.6196568248365983,3.872310138079448 +Hospitalised (no IC),1,2.060780612192024,0.8253060058141501,0.6097272807842798,3.825114785890012 +Hospitalised (no IC),2,2.031277113099352,0.8135002858747475,0.5998341402362193,3.77091923932976 +Hospitalised (no IC),3,2.001769019927607,0.8016756427114288,0.5900131421876905,3.7166500520360035 +Hospitalised (no IC),4,1.9724277450699799,0.7899113424180421,0.5802781903650331,3.6626606783905404 +Hospitalised (no IC),5,1.943188357005789,0.778173343946216,0.5706398593368343,3.6088031883085017 +Hospitalised (no IC),6,1.914171503465436,0.7665177651681706,0.5611066766176642,3.555327476568615 +Hospitalised (no IC),7,1.885325579796752,0.754917884409732,0.5516855552746821,3.502117297437643 +Hospitalised (no IC),8,1.8566886389887596,0.7433905940679534,0.5423822466670906,3.449248756337475 +Hospitalised (no IC),9,1.828288838731811,0.7319482480336674,0.5332015496388556,3.3967776870865043 +Hospitalised (no IC),10,1.8001009501904606,0.7205776027999837,0.5241473096087821,3.3446468460992036 +Hospitalised (no IC),11,1.7721309532061358,0.7092806783563669,0.5152226755264581,3.292866203474581 +Hospitalised (no IC),12,1.744442311445971,0.6980874141397269,0.506430488358056,3.2415686677321074 +Hospitalised (no IC),13,1.7169548951909883,0.6869585326058286,0.4977722945589065,3.1905818783488127 +Hospitalised (no IC),14,1.6897875146512362,0.675950935218145,0.4892504041931382,3.1401573190069514 +Hospitalised (no IC),15,1.6628062476780132,0.6649996558924876,0.4808647701138403,3.090009252488404 +Hospitalised (no IC),16,1.6361622584149742,0.6541773633577953,0.47262736374539016,3.0404583920135364 +Hospitalised (no IC),17,1.6097310965519556,0.6434239176535428,0.464539767971817,2.9912401266978765 +Hospitalised (no IC),18,1.5835822133350645,0.6327727346183597,0.4565891041487128,2.942502000923794 +Hospitalised (no IC),19,1.5577958984043692,0.6222624099723034,0.4491609767754555,2.8944142880294614 +Hospitalised (no IC),20,1.5321677328025558,0.6117947568918718,0.44255892950561526,2.8465445305459243 +Hospitalised (no IC),21,1.5067834824122477,0.6014112451659255,0.436021020564045,2.7990754972996306 +Hospitalised (no IC),22,1.481765934996827,0.5911709684422354,0.4295780933274556,2.7522675899023077 +Hospitalised (no IC),23,1.4569011279006647,0.5809716879618743,0.4231763267579589,2.705670184901106 +Hospitalised (no IC),24,1.4322734024523718,0.5708542677772133,0.41683688138570163,2.65946317055697 +Hospitalised (no IC),25,1.4078837750578483,0.5608196638886596,0.41055996828261965,2.6136505082242016 +Hospitalised (no IC),26,1.383681039367701,0.55084401808583,0.40435207037928,2.568126645568307 +Hospitalised (no IC),27,1.3597943201643348,0.5409893855573276,0.3982276216460874,2.5231648884680404 +Hospitalised (no IC),28,1.3360803104998562,0.531188193631261,0.39214710298038963,2.4784667167889998 +Hospitalised (no IC),29,1.3126009757016115,0.5214705901287269,0.3861264699570591,2.4341645954672653 +Hospitalised (no IC),30,1.2893441220552473,0.5118314914876694,0.38016259804165037,2.3902355136950266 +Hospitalised (no IC),31,1.2663414743045702,0.5022865862743535,0.37426366051555604,2.346748164599211 +Hospitalised (no IC),32,1.2434883708457622,0.4927873112640851,0.368402749569622,2.3034872782746034 +Hospitalised (no IC),33,1.2208157554110353,0.4833491932524356,0.36258784071246414,2.2605206054264495 +Hospitalised (no IC),34,1.1983622117603445,0.4739912156767355,0.3568288721369636,2.2179312813950367 +Hospitalised (no IC),35,1.1760534310291924,0.4646794959120517,0.35110674913436857,2.175568500727711 +Hospitalised (no IC),36,1.1539505572521715,0.45544357107997385,0.34543721136509875,2.133562072949147 +Hospitalised (no IC),37,1.1319964470208448,0.44625776144053286,0.33980558103625313,2.09179729653422 +Hospitalised (no IC),38,1.110253865060001,0.4371521847011149,0.33422801156948256,2.0504066997772754 +Hospitalised (no IC),39,1.0886574867242422,0.42809749256814794,0.3286877236329269,2.0092591625299883 +Hospitalised (no IC),40,1.0672888406830994,0.41913226476513904,0.3232056909023285,1.9685249669862261 +Hospitalised (no IC),41,1.0460576911287611,0.4102157930984378,0.31775873235438423,1.9280224024555477 +Hospitalised (no IC),42,1.0249107749473754,0.401324752594652,0.3123331762952994,1.887646586751364 +Hospitalised (no IC),43,1.0039987647717759,0.3925291338959649,0.306967762769883,1.8477077646745297 +Hospitalised (no IC),44,0.9831452101373256,0.3837492775309622,0.3016171611905276,1.8078503934403403 +Hospitalised (no IC),45,0.9624495927216736,0.3750315980537774,0.2963069529986637,1.7682799943681007 +Hospitalised (no IC),46,0.9419871527898442,0.36641096560619607,0.2910564824556342,1.7291512242960516 +Hospitalised (no IC),47,0.9215604047593074,0.35779888809396065,0.28581502529379677,1.690068382273699 +Hospitalised (no IC),48,0.9012892718407453,0.34925054262683564,0.2806134071672718,1.6512766187004375 +Hospitalised (no IC),49,0.8810793208420926,0.34072439489166423,0.27542738533499606,1.6125893505932063 +Hospitalised (no IC),50,0.8610761304642439,0.3322865562436247,0.2702943663561956,1.5743012597567825 +Hospitalised (no IC),51,0.8411192782242843,0.3238660927752646,0.2651731606238077,1.5360941133806039 +Hospitalised (no IC),52,0.8213621448262437,0.31553215237777116,0.2601031714959018,1.4982766955092723 +Hospitalised (no IC),53,0.8016755938216429,0.3072280314697086,0.2550512486179578,1.4605941844458499 +Hospitalised (no IC),54,0.782130144599559,0.29898573406388407,0.2500355119857224,1.423189291449689 +Hospitalised (no IC),55,0.7627150291116368,0.2908011034980203,0.24505320682578072,1.3860427359392833 +Hospitalised (no IC),56,0.7433349321523569,0.28263294180747733,0.24007987544815212,1.348968667938597 +Hospitalised (no IC),57,0.7240492249802069,0.27450811101764166,0.2351307721861077,1.3120869875469878 +Hospitalised (no IC),58,0.7048676262865854,0.2664314828888155,0.23020840203846094,1.2754189762450372 +Hospitalised (no IC),59,0.6858263302089472,0.2584193387774968,0.22532206316981304,1.239037440496251 +Hospitalised (no IC),60,0.6667653156570366,0.2504027473831721,0.22043069796306833,1.2026311590036227 +Hospitalised (no IC),61,0.647981069986097,0.24251070607084305,0.21561039994319364,1.1667812893040812 +Hospitalised (no IC),62,0.6293703114621048,0.23469853973245602,0.2108346667111553,1.131286352893701 +Hospitalised (no IC),63,0.6110337216510954,0.227009659727801,0.2061293322476435,1.0963423882878647 +Hospitalised (no IC),64,0.5928047373322914,0.21937250903763894,0.20145166338500314,1.0616261024301614 +Hospitalised (no IC),65,0.5747183362327472,0.21180238131494541,0.19681063874813537,1.0272063513168779 +Hospitalised (no IC),66,0.5567819264623843,0.2043026370869334,0.1922081622097958,0.9930981594893127 +Hospitalised (no IC),67,0.5387661442717696,0.19677593728500922,0.1875854075598622,0.9588604245404967 +Hospitalised (no IC),68,0.5210674915234719,0.18939126712530557,0.1830440802101382,0.9252581938522488 +Hospitalised (no IC),69,0.5030328276008676,0.1818719722144452,0.17841667335259667,0.8910371486728439 +Hospitalised (no IC),70,0.4854042133335717,0.1745327242729214,0.17389350459527442,0.8576235594831952 +Hospitalised (no IC),71,0.4679730305662187,0.16728499490352464,0.16942106823103725,0.8246161700955988 +Hospitalised (no IC),72,0.45055116358579805,0.16004972506544232,0.1649511292114434,0.7916557818398635 +Hospitalised (no IC),73,0.43340742302845475,0.15293997612624077,0.16055261442811802,0.7592557265826726 +Hospitalised (no IC),74,0.4163846491135123,0.14588974861806864,0.15618522662475515,0.7271159804244212 +Hospitalised (no IC),75,0.39965635913211994,0.13897152230857557,0.15189344938799337,0.6955660153869548 +Hospitalised (no IC),76,0.3829692377669683,0.13207957240051024,0.14760730347925874,0.6641244762363535 +Hospitalised (no IC),77,0.36610059762099956,0.1251218518362462,0.14325887708378326,0.6323711113201335 +Hospitalised (no IC),78,0.34999161764852416,0.11848913566846765,0.13912075357812248,0.6020857691469539 +Hospitalised (no IC),79,0.3335740189064352,0.11173844030441245,0.1349092372079168,0.5712488573880863 +Hospitalised (no IC),80,0.3177009722175541,0.10522294617570918,0.1306843050449378,0.541469949759239 +Hospitalised (no IC),81,0.30222791015361383,0.09888232160898826,0.12633296439369374,0.512473611995631 +Hospitalised (no IC),82,0.28682213468897577,0.0925796075379275,0.12248237977106481,0.48363308535080424 +Hospitalised (no IC),83,0.2718907651047155,0.0864821542027147,0.11867870582390505,0.45571161899103846 +Hospitalised (no IC),84,0.2576750696098668,0.08068814993448013,0.11402637041690752,0.4291578512514244 +Hospitalised (no IC),85,0.24369393194024508,0.07500046942459014,0.11026743500703,0.40306799627571005 +Hospitalised (no IC),86,0.2305469300278332,0.06966349990099863,0.10670223738802424,0.3785604559006088 +Hospitalised (no IC),87,0.21805456355343156,0.06460330685556441,0.10330244496773991,0.3552956414354733 +Hospitalised (no IC),88,0.20655018232615233,0.0599543161096465,0.09978856331095264,0.33389114490929467 +Hospitalised (no IC),89,0.19582860615397119,0.0556322169512496,0.09650914118607751,0.3139600950781917 +Hospitalised (no IC),90,0.1854510838258567,0.0514595406554769,0.09282270289560952,0.29468305588487725 +Hospitalised (no IC),91,0.17641694387477327,0.04783800673330082,0.09000117710024141,0.2779156404327129 +Hospitalised (no IC),92,0.16812128025224116,0.04452257214341102,0.08749080420463408,0.2625244357724133 +Hospitalised (no IC),93,0.1596111764094863,0.041132405681867275,0.08492156487478632,0.24673167773436236 +Hospitalised (no IC),94,0.15182796774790905,0.03804453167744162,0.08282710893157876,0.2322972538133905 +Hospitalised (no IC),95,0.14515319867678425,0.03540859079128553,0.08054605165359227,0.21992684957476483 +Hospitalised (no IC),96,0.13928914900409395,0.03310397676725444,0.07914398028870873,0.2090654750757989 +Hospitalised (no IC),97,0.13410630308828872,0.03107761244425992,0.0775643999057094,0.19947142064052467 +Hospitalised (no IC),98,0.1296316134672758,0.029337874784803215,0.07677394412485124,0.1909547001427929 +Hospitalised (no IC),99,0.12487715800019886,0.02750027947005221,0.07468548905638764,0.18184927039270554 +Hospitalised (no IC),100,0.12215311171907861,0.026455719103147486,0.07383054873569476,0.17663853959824566 +Hospitalised (no IC),101,0.11654754554237899,0.024319303665214512,0.07208859665779742,0.1668190722468932 +Hospitalised (no IC),102,0.1134369158612906,0.023148670020970374,0.07153357513204882,0.16132329685230765 +Hospitalised (no IC),103,0.10732995531755667,0.02088299844504821,0.0694795084041144,0.15045645232850197 +Hospitalised (no IC),104,0.09932530075169964,0.018020598391351887,0.06594702528210614,0.13699314645360258 +Hospitalised (no IC),105,0.09030289989049653,0.01503964071785684,0.061341974704148955,0.12183091177321118 +Hospitalised (no IC),106,0.06326751007289154,0.00870309958408291,0.045680982200924786,0.08090663658748803 +Hospitalised (IC),0,3.308592057603336,1.1548226219230324,1.2210710621369467,5.655427317087269 +Hospitalised (IC),1,3.2663264869485866,1.1410288678038767,1.2052395558745834,5.586046937581679 +Hospitalised (IC),2,3.218827366757618,1.1249113762165985,1.18777922355743,5.506334397720329 +Hospitalised (IC),3,3.1713225849004627,1.1087575847262385,1.1703359553537445,5.426511337334654 +Hospitalised (IC),4,3.1240871828102486,1.0926819706548376,1.1519681901220005,5.347098999005797 +Hospitalised (IC),5,3.07701726893083,1.0766333622727324,1.1337088318346704,5.267878749818243 +Hospitalised (IC),6,3.0303064137857088,1.060692966003955,1.115617951121311,5.189219155798216 +Hospitalised (IC),7,2.9838720240783716,1.0448207702991175,1.097685459776494,5.110948454770126 +Hospitalised (IC),8,2.937775212899061,1.0290408254048342,1.0799288296556824,5.033178803535145 +Hospitalised (IC),9,2.8920612121061477,1.0133706463967926,1.062361491589472,4.9559924681278025 +Hospitalised (IC),10,2.8466896526535685,0.9977903382904114,1.0449789444514117,4.87930485425777 +Hospitalised (IC),11,2.801670181929064,0.9823024629799918,1.0277860100608882,4.803130586089893 +Hospitalised (IC),12,2.757104591326834,0.9669505558171574,1.0108055841855594,4.727665737821604 +Hospitalised (IC),13,2.7128644578319725,0.9516768238868253,0.9940135051069107,4.652655882489308 +Hospitalised (IC),14,2.6691402840757115,0.9365644077476294,0.9774498537745808,4.578472111080424 +Hospitalised (IC),15,2.6257173887662253,0.9215180260539074,0.9610718140539272,4.504692711916309 +Hospitalised (IC),16,2.582838128733253,0.9066439846358635,0.9456124283149843,4.431790956095628 +Hospitalised (IC),17,2.5403029818516156,0.8918542066995151,0.9304159254210256,4.359376340471618 +Hospitalised (IC),18,2.498223294040428,0.8771975950962715,0.9154173607986327,4.287666599196455 +Hospitalised (IC),19,2.4567277936691663,0.8627306037054858,0.9006317175518734,4.216912999115838 +Hospitalised (IC),20,2.4154887045374602,0.8483097273137071,0.8859517970892106,4.14647740489633 +Hospitalised (IC),21,2.374643515983859,0.83399574320142,0.8714223162013193,4.076629550479781 +Hospitalised (IC),22,2.334389093738792,0.8198751019835117,0.8571075407528248,4.007753703197473 +Hospitalised (IC),23,2.2943823154648473,0.8057986458044214,0.8430165009704915,3.9391849857404697 +Hospitalised (IC),24,2.2547583825754556,0.7918262952611279,0.8291441906265511,3.8711888859061863 +Hospitalised (IC),25,2.2155188824733894,0.7779596929748775,0.815402074460824,3.803771304912884 +Hospitalised (IC),26,2.1765816461069245,0.7641643448336783,0.8017606051693651,3.736776575987223 +Hospitalised (IC),27,2.138153709450433,0.750530976131771,0.788294823528928,3.670607992573832 +Hospitalised (IC),28,2.1000051861411135,0.7369614552764641,0.7749218050830414,3.6048252036217328 +Hospitalised (IC),29,2.062235399824444,0.7235000416386524,0.7616776585186618,3.53962371666763 +Hospitalised (IC),30,2.024824734982071,0.7101396913506144,0.748555506104926,3.4749696636143876 +Hospitalised (IC),31,1.987824021819878,0.6969035599096643,0.7355738998116591,3.410964456810353 +Hospitalised (IC),32,1.9510653076570506,0.6837215101111744,0.7226724819640684,3.3472906361961208 +Hospitalised (IC),33,1.9145981553664058,0.6706165607394868,0.709869403539209,3.2840482458031675 +Hospitalised (IC),34,1.8784843943603147,0.6576166777185406,0.6971871969418557,3.221360006849333 +Hospitalised (IC),35,1.8426047352957053,0.6446732213810761,0.6845831602317195,3.1590035909645127 +Hospitalised (IC),36,1.807057203479059,0.6318294418331395,0.672092856446845,3.0971705587957845 +Hospitalised (IC),37,1.7717500377395017,0.6190486628726521,0.6596835422684801,3.0356918503574213 +Hospitalised (IC),38,1.736783876944957,0.6063748860163686,0.6473916695974623,2.974763015896534 +Hospitalised (IC),39,1.702053823111533,0.5937661959894015,0.6351798074450742,2.914190796036033 +Hospitalised (IC),40,1.6676906426369895,0.5812786957918306,0.6230951662348215,2.854226389799176 +Hospitalised (IC),41,1.633549449890007,0.5688541556571561,0.6110859969900847,2.7946019526330295 +Hospitalised (IC),42,1.5995446717939619,0.5564595657969292,0.5991219309766196,2.7351629645470337 +Hospitalised (IC),43,1.56591807511279,0.5441960719922238,0.5872898858132645,2.6763669553608582 +Hospitalised (IC),44,1.5323863552882049,0.5319496879796168,0.5754886531576882,2.6176898423933124 +Hospitalised (IC),45,1.4991091208133338,0.5197875918032312,0.5637756712288604,2.559434737623661 +Hospitalised (IC),46,1.4662070392681807,0.5077602641419805,0.5521943235468497,2.5018297105931366 +Hospitalised (IC),47,1.433362947084149,0.4957414840877329,0.5406315022411604,2.4442915574530373 +Hospitalised (IC),48,1.4007693024613048,0.48381073954064285,0.529156254223033,2.3871817534895348 +Hospitalised (IC),49,1.3682743940563,0.47190916560028323,0.5177146898850735,2.3302253903553933 +Hospitalised (IC),50,1.336111925416903,0.46013155631919817,0.5063904105136838,2.2738568831248007 +Hospitalised (IC),51,1.3040242020014539,0.4483771808497391,0.49509177438721175,2.2176073104424585 +Hospitalised (IC),52,1.272257494131067,0.43674492371306506,0.48390673734880213,2.161931821304605 +Hospitalised (IC),53,1.2406043301313507,0.4251544686964859,0.4727616303126062,2.1064549767377363 +Hospitalised (IC),54,1.2091779184038993,0.41365168679360853,0.46169694781098075,2.051387158069572 +Hospitalised (IC),55,1.1779609098799582,0.4022309848735063,0.45070669353925663,1.9967000243473723 +Hospitalised (IC),56,1.1468001150297336,0.39083435113286347,0.4397366588612484,1.9421198303991316 +Hospitalised (IC),57,1.1157908517589263,0.37950023749664696,0.42882091804589145,1.8878233216847944 +Hospitalised (IC),58,1.0849486817618104,0.3682358493403069,0.41796517209815814,1.8338419235068808 +Hospitalised (IC),59,1.0543317201107085,0.3570644347248451,0.40719017213162234,1.7802829507820985 +Hospitalised (IC),60,1.0236827632967043,0.3458890528751278,0.3964049684510573,1.7266880216278853 +Hospitalised (IC),61,0.9934782455069854,0.3348917502023983,0.3857784100438588,1.6739132468658398 +Hospitalised (IC),62,0.9635521759295751,0.3240095862998383,0.3752517503717646,1.6216618583263487 +Hospitalised (IC),63,0.9340663749902914,0.3133035960131659,0.36488222978262935,1.5702225983972644 +Hospitalised (IC),64,0.904753107856513,0.30267323579535044,0.35457522881515346,1.5191193275857715 +Hospitalised (IC),65,0.8756685737227717,0.29214009490663095,0.34435068922231876,1.4684534883568632 +Hospitalised (IC),66,0.8468246641675513,0.2817089477903794,0.3342128674740355,1.4182472304040936 +Hospitalised (IC),67,0.8178525974098281,0.27124365331347705,0.3240317731026528,1.3678510466680809 +Hospitalised (IC),68,0.7893898248170794,0.26098080469994356,0.3140323082922808,1.3183915273707196 +Hospitalised (IC),69,0.7603861529318191,0.25053380682160353,0.30384444155918167,1.2680217958425661 +Hospitalised (IC),70,0.7320347213173718,0.24034241614055107,0.2938886698254629,1.2188419534597876 +Hospitalised (IC),71,0.7040001051870692,0.23028277178406092,0.28404675627641174,1.1702611635010793 +Hospitalised (IC),72,0.6759797794743608,0.22024463289476431,0.2742122794075626,1.1217506080240447 +Hospitalised (IC),73,0.6484060398893,0.21038543343320532,0.2645373156092906,1.074066023592803 +Hospitalised (IC),74,0.621026142606597,0.20061311244784216,0.25493296078513,1.0267657125579135 +Hospitalised (IC),75,0.5941191810997175,0.19102826164125591,0.24549724002568912,0.9803346641684145 +Hospitalised (IC),76,0.5672777194115306,0.1814837605350893,0.23608702100689877,0.934064301000441 +Hospitalised (IC),77,0.5401435120742584,0.17185186639576683,0.22651412673249752,0.8873360943064469 +Hospitalised (IC),78,0.5142305765316966,0.1626744457712193,0.21719177863501127,0.8427697633875028 +Hospitalised (IC),79,0.48782042168790457,0.1533368281436647,0.20768214309347233,0.7973927365310327 +Hospitalised (IC),80,0.4622855867806246,0.14432817293687628,0.1984772958772208,0.7535739745266931 +Hospitalised (IC),81,0.43739356535242174,0.13556422465151854,0.18949467319387983,0.7109079999189644 +Hospitalised (IC),82,0.4126091065866786,0.12685487714922083,0.18054211432671777,0.6684723845196332 +Hospitalised (IC),83,0.38858728298959533,0.11843106943660636,0.1718558884207599,0.6273903044940745 +Hospitalised (IC),84,0.3657163965542116,0.11042774230576392,0.1635771293066991,0.588321833208783 +Hospitalised (IC),85,0.3432223548344473,0.1025714878183624,0.1554271727130011,0.5499369202675022 +Hospitalised (IC),86,0.32207004482785584,0.09519928610748024,0.14775569522107984,0.5138811744039565 +Hospitalised (IC),87,0.30197071853786117,0.08820802373921127,0.1404594550709031,0.4796547576358213 +Hospitalised (IC),88,0.28346093253444316,0.08178275626893898,0.13369883201020002,0.44816624104049135 +Hospitalised (IC),89,0.26621059834256267,0.07580624663259213,0.12743484586786014,0.41884623819102174 +Hospitalised (IC),90,0.24951367232855515,0.07003215575119748,0.12133974230611622,0.3904890029818509 +Hospitalised (IC),91,0.2349785685047653,0.0650162940333103,0.11602908846995015,0.3658245788420378 +Hospitalised (IC),92,0.22163185265611635,0.06041932231224316,0.11117497855120703,0.34319241668476214 +Hospitalised (IC),93,0.2079399260684369,0.05571190020375144,0.10582834885938389,0.3199874209096867 +Hospitalised (IC),94,0.19541803435420252,0.05141650271212634,0.10076751971223387,0.29917232183940246 +Hospitalised (IC),95,0.18468057716837516,0.04774218550306203,0.09601356832618561,0.2813477946116477 +Hospitalised (IC),96,0.1752487774926033,0.04452233868891465,0.09251567906628803,0.26570429809609886 +Hospitalised (IC),97,0.16691469037862716,0.04168410236872879,0.0888494977646423,0.25189202963033147 +Hospitalised (IC),98,0.15972214714505678,0.03924078565298665,0.08629729886743093,0.23998027776046027 +Hospitalised (IC),99,0.1520834622593604,0.03665193186055516,0.0836213266078191,0.2270368294711182 +Hospitalised (IC),100,0.14771205524503142,0.03517607603691496,0.08189336267369471,0.21947167914550647 +Hospitalised (IC),101,0.13872233887769017,0.03214531665769532,0.07882042432026225,0.20390906393517144 +Hospitalised (IC),102,0.13374649001255917,0.030476669268054007,0.07713252836382381,0.19530206815662518 +Hospitalised (IC),103,0.12401128323747163,0.02722690097574827,0.07456492508542426,0.17876333752847096 +Hospitalised (IC),104,0.11140118449465622,0.023073036018184306,0.0695204317750224,0.1579900176867259 +Hospitalised (IC),105,0.09764286135808492,0.018675487484194887,0.06442269761506825,0.13551641017470126 +Hospitalised (IC),106,0.06301307845401817,0.009442466564163607,0.044935199120540587,0.08177075484206882 +Non-hospitalised (no AD),0,0.060090815673273054,0.0030609520073688616,0.054152914707442594,0.06603870837643472 +Non-hospitalised (no AD),1,0.060089981350366875,0.0030609520073688616,0.05415208038453641,0.06603787405352858 +Non-hospitalised (no AD),2,0.060088723684577586,0.003060952007368861,0.05415082271874712,0.06603661638773928 +Non-hospitalised (no AD),3,0.06008716949364365,0.0030609520073688625,0.054149268527813176,0.06603506219680531 +Non-hospitalised (no AD),4,0.0600853724165159,0.0030609520073688625,0.05414747145068542,0.06603326511967755 +Non-hospitalised (no AD),5,0.06008336474954803,0.0030609520073688616,0.054145463783717554,0.06603125745270967 +Non-hospitalised (no AD),6,0.06008116872894839,0.0030609520073688616,0.05414326776311791,0.06602906143211003 +Non-hospitalised (no AD),7,0.060078800873271915,0.0030609520073688616,0.054140899907441434,0.0660266935764336 +Non-hospitalised (no AD),8,0.0600762740782928,0.0030609520073688625,0.05413837311246233,0.06602416678145447 +Non-hospitalised (no AD),9,0.06007359877206418,0.0030609520073688608,0.0541356978062337,0.06602149147522586 +Non-hospitalised (no AD),10,0.06007078361193263,0.003060952007368862,0.05413288264610216,0.06601867631509427 +Non-hospitalised (no AD),11,0.06006783593367814,0.0030609520073688616,0.054129934967847654,0.0660157286368398 +Non-hospitalised (no AD),12,0.06006476205580132,0.0030609520073688616,0.05412686108997085,0.06601265475896297 +Non-hospitalised (no AD),13,0.06006156749400594,0.003060952007368861,0.05412366652817548,0.06600946019716762 +Non-hospitalised (no AD),14,0.06005825711732272,0.003060952007368862,0.05412035615149224,0.06600614982048437 +Non-hospitalised (no AD),15,0.06005483526482364,0.0030609520073688608,0.05411693429899318,0.06600272796798533 +Non-hospitalised (no AD),16,0.06005130583485792,0.003060952007368862,0.05411340486902745,0.06599919853801961 +Non-hospitalised (no AD),17,0.060047672354602594,0.003060952007368861,0.05410977138877212,0.06599556505776423 +Non-hospitalised (no AD),18,0.06004393803517904,0.003060952007368861,0.05410603706934857,0.06599183073834074 +Non-hospitalised (no AD),19,0.060040105815969594,0.0030609520073688608,0.054102204850139134,0.06598799851913124 +Non-hospitalised (no AD),20,0.06003617840071023,0.0030609520073688616,0.05409827743487977,0.06598407110387192 +Non-hospitalised (no AD),21,0.060032158287220874,0.0030609520073688616,0.05409425732139041,0.06598005099038252 +Non-hospitalised (no AD),22,0.06002804779214616,0.0030609520073688616,0.0540901468263157,0.06597594049530782 +Non-hospitalised (no AD),23,0.06002384907173272,0.003060952007368861,0.05408594810590224,0.06597174177489436 +Non-hospitalised (no AD),24,0.06001956413942321,0.0030609520073688616,0.05408166317359274,0.06596745684258491 +Non-hospitalised (no AD),25,0.060015194880867116,0.0030609520073688608,0.05407729391503664,0.06596308758402879 +Non-hospitalised (no AD),26,0.060010743066816453,0.0030609520073688616,0.05407284210098598,0.06595863576997811 +Non-hospitalised (no AD),27,0.06000621036427375,0.003060952007368862,0.05406830939844327,0.06595410306743542 +Non-hospitalised (no AD),28,0.06000159834618566,0.0030609520073688616,0.05406369738035518,0.06594949104934733 +Non-hospitalised (no AD),29,0.05999690849991671,0.0030609520073688603,0.05405900753408623,0.0659448012030784 +Non-hospitalised (no AD),30,0.0599921422346943,0.0030609520073688616,0.05405424126886384,0.06594003493785601 +Non-hospitalised (no AD),31,0.05998730088817876,0.0030609520073688608,0.05404939992234829,0.06593519359134044 +Non-hospitalised (no AD),32,0.05998238573228755,0.0030609520073688616,0.05404448476645706,0.0659302784354492 +Non-hospitalised (no AD),33,0.05997739797837801,0.0030609520073688608,0.05403949701254754,0.0659252906815397 +Non-hospitalised (no AD),34,0.05997233878187686,0.0030609520073688616,0.05403443781604638,0.06592023148503855 +Non-hospitalised (no AD),35,0.05996720924643023,0.0030609520073688603,0.054029308280599764,0.06591510194959191 +Non-hospitalised (no AD),36,0.05996201042763603,0.0030609520073688608,0.05402410946180554,0.06590990313079771 +Non-hospitalised (no AD),37,0.05995674333641101,0.0030609520073688556,0.054018842370580536,0.0659046360395727 +Non-hospitalised (no AD),38,0.05995140894203747,0.003060952007368864,0.05401350797620698,0.06589930164519911 +Non-hospitalised (no AD),39,0.059946008174927345,0.0030609520073688616,0.05400810720909687,0.06589390087808902 +Non-hospitalised (no AD),40,0.05994054192913685,0.003060952007368861,0.054002640963306384,0.06588843463229854 +Non-hospitalised (no AD),41,0.05993501106465942,0.0030609520073688608,0.05399711009882895,0.06588290376782108 +Non-hospitalised (no AD),42,0.05992941640952201,0.003060952007368861,0.05399151544369153,0.06587730911268369 +Non-hospitalised (no AD),43,0.059923758761705345,0.0030609520073688625,0.053985857795874864,0.065871651464867 +Non-hospitalised (no AD),44,0.05991803889090676,0.0030609520073688608,0.05398013792507627,0.06586593159406845 +Non-hospitalised (no AD),45,0.05991225754016234,0.003060952007368861,0.05397435657433186,0.06586015024332405 +Non-hospitalised (no AD),46,0.05990641542734161,0.003060952007368862,0.053968514461511144,0.06585430813050326 +Non-hospitalised (no AD),47,0.059900513246528106,0.0030609520073688608,0.05396261228069764,0.06584840594968981 +Non-hospitalised (no AD),48,0.05989455166929618,0.0030609520073688608,0.05395665070346569,0.06584244437245787 +Non-hospitalised (no AD),49,0.05988853134589451,0.0030609520073688625,0.05395063038006403,0.06583642404905617 +Non-hospitalised (no AD),50,0.05988245290634403,0.003060952007368862,0.05394455194051355,0.06583034560950568 +Non-hospitalised (no AD),51,0.05987631696145891,0.0030609520073688616,0.05393841599562842,0.06582420966462056 +Non-hospitalised (no AD),52,0.05987012410379684,0.0030609520073688603,0.05393222313796636,0.06581801680695855 +Non-hospitalised (no AD),53,0.059863874908545014,0.003060952007368861,0.05392597394271455,0.06581176761170666 +Non-hospitalised (no AD),54,0.05985756993434705,0.0030609520073688616,0.05391966896851657,0.06580546263750871 +Non-hospitalised (no AD),55,0.05985120972407613,0.0030609520073688616,0.05391330875824565,0.06579910242723781 +Non-hospitalised (no AD),56,0.05984479480555864,0.0030609520073688616,0.05390689383972817,0.06579268750872029 +Non-hospitalised (no AD),57,0.05983832569225224,0.0030609520073688616,0.05390042472642176,0.06578621839541389 +Non-hospitalised (no AD),58,0.05983180288388197,0.0030609520073688616,0.05389390191805149,0.06577969558704365 +Non-hospitalised (no AD),59,0.05982522686703769,0.0030609520073688616,0.0538873259012072,0.06577311957019935 +Non-hospitalised (no AD),60,0.0598185981157361,0.003060952007368861,0.05388069714990562,0.0657664908188978 +Non-hospitalised (no AD),61,0.059811917091949295,0.0030609520073688608,0.053874016126118834,0.06575980979511098 +Non-hospitalised (no AD),62,0.05980518424610354,0.0030609520073688586,0.053867283280273064,0.06575307694926517 +Non-hospitalised (no AD),63,0.05979840001754861,0.003060952007368863,0.05386049905171813,0.06574629272071024 +Non-hospitalised (no AD),64,0.05979156483500218,0.003060952007368862,0.053853663869171696,0.06573945753816383 +Non-hospitalised (no AD),65,0.05978467911696856,0.0030609520073688616,0.05384677815113808,0.06573257182013025 +Non-hospitalised (no AD),66,0.05977774327213559,0.0030609520073688625,0.05383984230630512,0.06572563597529725 +Non-hospitalised (no AD),67,0.05977075769974987,0.0030609520073688625,0.0538328567339194,0.06571865040291157 +Non-hospitalised (no AD),68,0.059763722789972364,0.0030609520073688616,0.05382582182414188,0.06571161549313403 +Non-hospitalised (no AD),69,0.059756638924215864,0.0030609520073688616,0.05381873795838539,0.06570453162737755 +Non-hospitalised (no AD),70,0.059749506475464655,0.003060952007368861,0.053811605509634174,0.06569739917862635 +Non-hospitalised (no AD),71,0.059742325808578696,0.0030609520073688616,0.053804424842748215,0.06569021851174038 +Non-hospitalised (no AD),72,0.059735097280582246,0.003060952007368862,0.05379719631475177,0.0656829899837439 +Non-hospitalised (no AD),73,0.05972782124093873,0.003060952007368862,0.05378992027510826,0.06567571394410038 +Non-hospitalised (no AD),74,0.05972049803181186,0.0030609520073688608,0.05378259706598139,0.06566839073497356 +Non-hospitalised (no AD),75,0.05971312798831495,0.0030609520073688616,0.05377522702248448,0.06566102069147664 +Non-hospitalised (no AD),76,0.05970571143874809,0.0030609520073688616,0.053767810472917614,0.06565360414190974 +Non-hospitalised (no AD),77,0.05969824870482418,0.0030609520073688608,0.05376034773899371,0.06564614140798586 +Non-hospitalised (no AD),78,0.05969074010188524,0.003060952007368858,0.053752839136054775,0.06563863280504689 +Non-hospitalised (no AD),79,0.059683185939107916,0.0030609520073688265,0.053745284973277505,0.06563107864226951 +Non-hospitalised (no AD),80,0.059675586519698835,0.003060952007368601,0.05373768555386886,0.06562347922286002 +Non-hospitalised (no AD),81,0.05966794214107046,0.0030609520073670666,0.053730041175243465,0.06561583484422867 +Non-hospitalised (no AD),82,0.05966025309493671,0.0030609520073568057,0.053722352129129625,0.06560814579807495 +Non-hospitalised (no AD),83,0.059652519666999394,0.0030609520072936683,0.0537146187013148,0.06560041237001497 +Non-hospitalised (no AD),84,0.05964474213466399,0.003060952006945853,0.05370684116965411,0.0655926348370037 +Non-hospitalised (no AD),85,0.05963692075393597,0.0030609520050871933,0.053699019792531666,0.06558481345266402 +Non-hospitalised (no AD),86,0.0596290557076245,0.0030609519963878994,0.05369115476309586,0.0655769483894485 +Non-hospitalised (no AD),87,0.05962114689674435,0.003060951959067391,0.05368324602461326,0.06556903950604895 +Non-hospitalised (no AD),88,0.05961319334808661,0.003060951820384107,0.05367529274498545,0.06556108568790862 +Non-hospitalised (no AD),89,0.05960519141826775,0.003060951353002472,0.053667291721834126,0.06555308284989654 +Non-hospitalised (no AD),90,0.05959712866416511,0.003060949824113668,0.05365923193360298,0.06554501712493174 +Non-hospitalised (no AD),91,0.05958898053463163,0.0030609458198238206,0.05365109157193926,0.06553686121445745 +Non-hospitalised (no AD),92,0.0595806881041859,0.0030609359882566582,0.053642818213622726,0.06552854967978977 +Non-hospitalised (no AD),93,0.0595720626090189,0.003060909684581541,0.053634243744612786,0.06551987307260371 +Non-hospitalised (no AD),94,0.05956279290583925,0.003060849589000841,0.05362509062006809,0.06551048659462162 +Non-hospitalised (no AD),95,0.05955249121784922,0.0030607343594730705,0.05361501246433945,0.06549996099823127 +Non-hospitalised (no AD),96,0.0595405599591667,0.0030605312297731168,0.05360347525431849,0.06548763502781946 +Non-hospitalised (no AD),97,0.059526266429100264,0.003060200301371463,0.053589823687951926,0.06547269845381742 +Non-hospitalised (no AD),98,0.05950909542852074,0.0030597139635340524,0.0535735961278143,0.06545458242526378 +Non-hospitalised (no AD),99,0.05948499837999192,0.003058851094606883,0.053551172947488464,0.0654288086919051 +Non-hospitalised (no AD),100,0.05946313360449935,0.003058114262249683,0.05353073754349697,0.06540551213970049 +Non-hospitalised (no AD),101,0.05940997771502873,0.0030556764613151872,0.053482310712335224,0.0653476192343096 +Non-hospitalised (no AD),102,0.05936206589542377,0.0030535326197752206,0.05343855770303505,0.06529554160634776 +Non-hospitalised (no AD),103,0.05922340749027942,0.0030464968192796425,0.053313547955712685,0.0651432115767128 +Non-hospitalised (no AD),104,0.05884938845325926,0.0030268943950772093,0.052977555406309454,0.06473110206470405 +Non-hospitalised (no AD),105,0.05787749740401882,0.0029756968363162917,0.05210498183505715,0.06365972641543047 +Non-hospitalised (no AD),106,0.04741897410018828,0.002430409946310288,0.04270425347313574,0.05214162821255264 diff --git a/notebooks/analysis/woldmuyn_long_COVID.py b/notebooks/analysis/woldmuyn_long_COVID.py index a621ba470..0f8d0e004 100644 --- a/notebooks/analysis/woldmuyn_long_COVID.py +++ b/notebooks/analysis/woldmuyn_long_COVID.py @@ -16,6 +16,12 @@ import matplotlib.cm as cm import os +# From Color Universal Design (CUD): https://jfly.uni-koeln.de/color/ +colorscale_okabe_ito = {"orange" : "#E69F00", "light_blue" : "#56B4E9", + "green" : "#009E73", "yellow" : "#F0E442", + "blue" : "#0072B2", "red" : "#D55E00", + "pink" : "#CC79A7", "black" : "#000000"} + if __name__ == '__main__': ################################ ## Define simulation settings ## @@ -73,9 +79,40 @@ result_folder = '../../results/covid19_DTM/analysis/QALY/long_COVID' states = ['QALY_NH', 'QALY_C', 'QALY_ICU','QALY_D'] - titles = ['Non-hospitalised', 'Cohort', 'ICU','Deaths'] + titles = ['Non-hospitalised', 'Non-hospitalised (no IC)', 'Non-hospitalised (IC)','Deaths'] colors = ['green','yellow','red','black'] + age_groups=['0-12','12-18','18-25','25-35','35-45','45-55','55-65','65-75','75-85','85+'] + + fig,axes=plt.subplots(nrows=1, ncols=2, figsize=(8.3,0.25*11.7)) + for ax,scenario,out in zip(axes, ['no_AD','AD'], [out_no_AD,out_AD]): + # group data + data = [out['QALY_D'].mean(dim="draws").sum(dim='doses').isel(date=-1).values/initN*100000, + out['QALY_NH'].mean(dim="draws").sum(dim='doses').isel(date=-1).values/initN*100000, + out['QALY_C'].mean(dim="draws").sum(dim='doses').isel(date=-1).values/initN*100000, + out['QALY_ICU'].mean(dim="draws").sum(dim='doses').isel(date=-1).values/initN*100000 + ] + # settings + colors = ['grey', colorscale_okabe_ito['green'], colorscale_okabe_ito['orange'], colorscale_okabe_ito['red']] + hatches = ['....','////','\\\\\\','||||'] + labels = ['Deaths', 'Non-hospitalised', 'Hospitalised (no IC)', 'Hospitalised (IC)'] + # plot data + bottom=np.zeros(len(age_groups)) + for d, color, hatch, label in zip(data, colors, hatches, labels): + p = ax.bar(age_groups, d, color=color, hatch=hatch, label=label, bottom=bottom, linewidth=0.25, alpha=0.7) + bottom += d + ax.grid(False) + ax.tick_params(axis='both', which='major', labelsize=10) + ax.tick_params(axis='x', which='major', rotation=30) + ax.set_ylim([0,7000]) + axes[0].legend(loc=2, framealpha=1, fontsize=10) + axes[0].set_ylabel('QALYs lost per 100K inhab.', size=10) + + plt.tight_layout() + plt.show() + fig.savefig('QALY_losses_per_age_group.pdf') + plt.close() + for scenario,out in zip(['no_AD','AD'],[out_no_AD,out_AD]): # With confidence interval @@ -95,7 +132,7 @@ ax.grid(False) plt.subplots_adjust(hspace=0.5) - fig.savefig(os.path.join(abs_dir,result_folder,f'QALY_losses_{scenario}.png')) + fig.savefig(os.path.join(abs_dir,result_folder,f'QALY_losses_{scenario}.pdf')) # QALYS per age group Palette=cm.get_cmap('tab10_r', initN.size).colors @@ -113,4 +150,4 @@ axs[0].legend(fancybox=True, frameon=True, framealpha=1, fontsize=15,title='Age Group', loc="upper left", bbox_to_anchor=(1,1)) plt.subplots_adjust(hspace=0.5) - fig.savefig(os.path.join(abs_dir,result_folder,f'QALY_losses_per_age_group_{scenario}.png'), dpi=600) \ No newline at end of file + fig.savefig(os.path.join(abs_dir,result_folder,f'QALY_losses_per_age_group_{scenario}.pdf')) \ No newline at end of file diff --git a/notebooks/calibration/plot_fit_BASE-COVID19_SEIQRD_hybrid_vacc.py b/notebooks/calibration/plot_fit_BASE-COVID19_SEIQRD_hybrid_vacc.py index 216ea1ab3..a1d286160 100644 --- a/notebooks/calibration/plot_fit_BASE-COVID19_SEIQRD_hybrid_vacc.py +++ b/notebooks/calibration/plot_fit_BASE-COVID19_SEIQRD_hybrid_vacc.py @@ -68,7 +68,7 @@ ################################ # Start and end of simulation -end_sim = datetime(2022,1,1) +end_sim = datetime(2021,7,1) # Confidence level used to visualise model fit conf_int = 0.05 @@ -146,7 +146,7 @@ print('2) Visualizing fit') -fig,(ax1,ax2,ax3,ax4,ax5) = plt.subplots(nrows=5,ncols=1,figsize=(8.3,11.7),sharex=True) +fig,(ax1,ax2,ax3,ax4,ax5) = plt.subplots(nrows=5,ncols=1,figsize=(8.3,0.75*11.7),sharex=True) # Plot mildly sick ax1.plot(df_2plot['M_in','mean'], color='blue', linewidth=1.5) @@ -154,8 +154,9 @@ ax1.scatter(df_cases[start_calibration:end_sim].index,df_cases[start_calibration:end_sim], color='black', alpha=0.20, linestyle='None', facecolors='black', s=10) ax1 = _apply_tick_locator(ax1) ax1.set_xlim(start_sim,end_sim) -ax1.set_ylabel('Incidence\nMild cases (-)', fontsize=13) +ax1.set_ylabel('Incidence\nMild cases (-)', fontsize=10) ax1.get_yaxis().set_label_coords(-0.1,0.5) +ax1.tick_params(axis='both', which='major', labelsize=10) ax1.grid(False) # Plot hospitalizations ax2.plot(df_2plot['H_in','mean'], color='blue', linewidth=1.5) @@ -164,24 +165,27 @@ ax2.scatter(df_hosp[pd.to_datetime(end_calibration)+timedelta(days=1):end_sim].index,df_hosp['H_in'][pd.to_datetime(end_calibration)+timedelta(days=1):end_sim], color='black', alpha=0.2, linestyle='None', facecolors='black', s=10) ax2 = _apply_tick_locator(ax2) ax2.set_xlim(start_sim,end_sim) -ax2.set_ylabel('Incidence\nHospital (-)', fontsize=13) +ax2.set_ylabel('Incidence\nHospital (-)', fontsize=10) ax2.get_yaxis().set_label_coords(-0.1,0.5) +ax2.tick_params(axis='both', which='major', labelsize=10) ax2.grid(False) # Plot hospital total ax3.plot(simtime, df_2plot['H_tot', 'mean'], color='blue', linewidth=1.5) ax3.fill_between(simtime, df_2plot['H_tot', 'lower'], df_2plot['H_tot', 'upper'], alpha=0.20, color = 'blue') ax3.scatter(df_hosp[start_calibration:end_sim].index,df_hosp['H_tot'][start_calibration:end_sim], color='black', alpha=0.2, linestyle='None', facecolors='black', s=10) ax3 = _apply_tick_locator(ax3) -ax3.set_ylabel('Load\nHospital (-)', fontsize=13) +ax3.set_ylabel('Load\nHospital (-)', fontsize=10) ax3.get_yaxis().set_label_coords(-0.1,0.5) +ax3.tick_params(axis='both', which='major', labelsize=10) ax3.grid(False) # Plot ICU ax4.plot(simtime, df_2plot['ICU_R', 'mean']+df_2plot['ICU_D', 'mean']+df_2plot['C_icurec', 'mean'], color='blue', linewidth=1.5) ax4.fill_between(simtime, df_2plot['ICU_R', 'lower']+df_2plot['ICU_D', 'lower']+df_2plot['C_icurec', 'lower'], df_2plot['ICU_R', 'upper']+df_2plot['ICU_D', 'upper']+df_2plot['C_icurec', 'upper'], alpha=0.20, color = 'blue') ax4.scatter(df_hosp[start_calibration:end_sim].index,df_hosp['ICU_tot'][start_calibration:end_sim], color='black', alpha=0.2, linestyle='None', facecolors='black', s=10) ax4 = _apply_tick_locator(ax4) -ax4.set_ylabel('Load\nIntensive Care (-)', fontsize=13) +ax4.set_ylabel('Load\nIntensive Care (-)', fontsize=10) ax4.get_yaxis().set_label_coords(-0.1,0.5) +ax4.tick_params(axis='both', which='major', labelsize=10) ax4.grid(False) # Plot fraction of immunes ax5.plot(df_2plot['R','mean'][start_calibration:'2021-03-01']/sum(initN)*100, color='blue', linewidth=1.5) @@ -190,12 +194,13 @@ ax5.errorbar(x=df_sero_herzog.index,y=df_sero_herzog['rel','mean'].values*100,yerr=yerr, fmt='x', color='black', elinewidth=1, capsize=5) yerr = np.array([df_sero_sciensano['rel','mean']*100 - df_sero_sciensano['rel','LL']*100, df_sero_sciensano['rel','UL']*100 - df_sero_sciensano['rel','mean']*100 ]) ax5.errorbar(x=df_sero_sciensano.index,y=df_sero_sciensano['rel','mean']*100,yerr=yerr, fmt='^', color='black', elinewidth=1, capsize=5) -ax5.legend(['model (mean)', 'model (95% CI)', 'Herzog et al. 2020', 'Sciensano'], loc='upper right', fontsize=13) +ax5.legend(['model (mean)', 'model (95% CI)', 'Herzog et al. 2020', 'Sciensano'], loc=2, fontsize=8) ax5.axvline(x=pd.Timestamp('2020-12-27'), linewidth=1.5, linestyle='--', color='black') ax5 = _apply_tick_locator(ax5) +ax5.tick_params(axis='both', which='major', labelsize=10) ax5.set_xlim(start_sim,end_sim) ax5.set_ylim(0,35) -ax5.set_ylabel('Seroprelevance (%)', fontsize=13) +ax5.set_ylabel('Seroprelevance (%)', fontsize=10) ax5.get_yaxis().set_label_coords(-0.1,0.5) ax5.grid(False) diff --git a/notebooks/preprocessing/woldmuyn_long_covid_average_QALY_loss.py b/notebooks/preprocessing/woldmuyn_long_covid_average_QALY_loss.py index 0f6862f4d..6502ca2e5 100644 --- a/notebooks/preprocessing/woldmuyn_long_covid_average_QALY_loss.py +++ b/notebooks/preprocessing/woldmuyn_long_covid_average_QALY_loss.py @@ -36,10 +36,11 @@ from covid19_DTM.visualization.output import _apply_tick_locator from covid19_DTM.models.QALY import life_table_QALY_model, lost_QALYs_hospital_care -import copy import emcee from tqdm import tqdm +import sys + ############### ## Load data ## ############### @@ -54,8 +55,8 @@ # --------------- # severity_groups = ['Mild','Moderate','Severe-Critical'] -hospitalisation_groups = ['Non-hospitalised','Cohort','ICU'] -color_dict = {'Mild':'g','Non-hospitalised':'g','Non-hospitalised (no AD)':'g','Moderate':'y','Cohort':'y','Severe-Critical':'r','ICU':'r'} +hospitalisation_groups = ['Non-hospitalised','Hospitalised (no IC)','Hospitalised (IC)'] +color_dict = {'Mild':'g','Non-hospitalised':'g','Non-hospitalised (no AD)':'g','Moderate':'y','Hospitalised (no IC)':'y','Severe-Critical':'r','Hospitalised (IC)':'r'} # raw prevalence data per severity group prevalence_data_per_severity_group = pd.read_csv(os.path.join(abs_dir,rel_dir,'Long_COVID_prevalence.csv'),index_col=[0,1]) @@ -91,7 +92,7 @@ sd_QoL_decrease_non_hospitalised = 0.33/np.sqrt(1146) #(0.66-0.64)/1.96 QoL_difference_data = pd.DataFrame(data=np.array([[mean_QoL_decrease_non_hospitalised,mean_QoL_decrease_hospitalised,mean_QoL_decrease_hospitalised], [sd_QoL_decrease_non_hospitalised,sd_QoL_decrease_hospitalised,sd_QoL_decrease_hospitalised]]).transpose(), - columns=['mean','sd'],index=['Non-hospitalised','Cohort','ICU']) + columns=['mean','sd'],index=['Non-hospitalised','Hospitalised (no IC)','Hospitalised (IC)']) # ------- # # results # @@ -157,11 +158,15 @@ def WSSE(theta,x,y): prevalence_func = lambda t,tau, p_AD: p_AD + (1-p_AD)*np.exp(-t/tau) -for scenario,(fig,ax) in zip(('AD','no_AD'),(plt.subplots(),plt.subplots())): - for hospitalisation in hospitalisation_groups: +fig,axes=plt.subplots(nrows=1, ncols=2, sharey=True, figsize=(8.3,0.3*11.7)) + +for ax, scenario, in zip(axes, ('AD','no_AD')): + markers = ['o', 's', '^'] + linestyles = ['-', '--', '-.'] + for (hospitalisation, marker, linestyle) in zip(hospitalisation_groups, markers, linestyles): x = prevalence_data_per_hospitalisation_group.loc[hospitalisation].index.values y = prevalence_data_per_hospitalisation_group.loc[hospitalisation].values.squeeze() - ax.plot(x,y,color_dict[hospitalisation]+'o',label=f'{hospitalisation} data') + ax.scatter(x,100*y,label=f'{hospitalisation} data', color='black', marker=marker) if hospitalisation =='Non-hospitalised' and scenario == 'no_AD': tau = taus['Non-hospitalised (no AD)'] @@ -169,15 +174,17 @@ def WSSE(theta,x,y): else: tau = taus[hospitalisation] p_AD = p_ADs[hospitalisation] - ax.plot(time,prevalence_func(time,tau,p_AD),color_dict[hospitalisation]+'--',alpha=0.7, - label=f'{hospitalisation} fit\n'rf'($\tau$:{tau:.2f},'' $f_{AD}$'f':{p_AD:.2f})') + ax.plot(time,100*prevalence_func(time,tau,p_AD),color='black',alpha=0.8, linewidth=1.5, + label=f'{hospitalisation} fit\n'rf'($\tau$:{tau:.2f},'' $f_{AD}$'f':{p_AD:.2f})', linestyle=linestyle) - ax.set_xlabel('Months after infection') - ax.set_ylabel('Prevalence') - ax.set_title('Prevalence of long-COVID symptoms') - ax.legend(prop={'size': 12}) + ax.set_xlabel('Months after infection', size=10) + ax.legend(prop={'size': 8}) ax.grid(False) - fig.savefig(os.path.join(abs_dir,fig_result_folder,f'prevalence_first_fit_{scenario}')) + ax.tick_params(axis='both', which='major', labelsize=10) + +axes[0].set_ylabel('Symptom prevalence (%)', size=10) +plt.tight_layout() +fig.savefig(os.path.join(abs_dir,fig_result_folder,f'prevalence_first_fit.pdf')) print('\n(2.2) MCMC to estimate f_AD\n') # objective functions for MCMC @@ -216,63 +223,60 @@ def log_probability(theta, tau, x, y, bounds): p_AD_summary = pd.DataFrame(index=hospitalisation_groups,columns=['mean','sd','lower','upper'],dtype='float64') for hospitalisation in hospitalisation_groups: - + # slice data x = prevalence_data_per_hospitalisation_group.loc[hospitalisation].index.values y = prevalence_data_per_hospitalisation_group.loc[hospitalisation].values.squeeze() - + # set parameters tau = taus[hospitalisation] p_AD = p_ADs[hospitalisation] - + # setup sampler nwalkers = 32 ndim = 1 pos = p_AD + p_AD*1e-1 * np.random.randn(nwalkers, ndim) - bounds = (0,1) sampler = emcee.EnsembleSampler( nwalkers, ndim, log_probability, args=(tau,x,y,bounds) ) samplers.update({hospitalisation:sampler}) - sampler.run_mcmc(pos, 20000, progress=True) - - flat_samples = sampler.get_chain(discard=1000, thin=30, flat=True) + # run sampler + sampler.run_mcmc(pos, 25000, progress=True) + # extract chains + flat_samples = sampler.get_chain(discard=5000, thin=100, flat=True) + # print results p_AD_summary['mean'][hospitalisation] = np.mean(flat_samples,axis=0) p_AD_summary['sd'][hospitalisation] = np.std(flat_samples,axis=0) p_AD_summary['lower'][hospitalisation] = np.quantile(flat_samples,0.025,axis=0) p_AD_summary['upper'][hospitalisation] = np.quantile(flat_samples,0.975,axis=0) # visualise MCMC results -fig,axs = plt.subplots(2,2,figsize=(10,10),sharex=True,sharey=True) -axs = axs.reshape(-1) -for ax,hospitalisation in zip(axs[:-1],hospitalisation_groups): +fig,axes = plt.subplots(nrows=1,ncols=3,figsize=(8.3,0.25*11.7),sharey=True) + +for ax,hospitalisation in zip(axes,hospitalisation_groups): x = prevalence_data_per_hospitalisation_group.loc[hospitalisation].index.values y = prevalence_data_per_hospitalisation_group.loc[hospitalisation].values.squeeze() - ax.plot(x,y,color_dict[hospitalisation]+'o',label=f'{hospitalisation} data') - axs[-1].plot(x,y,color_dict[hospitalisation]+'o',label=f'{hospitalisation} data') + ax.scatter(x,100*y,color='black',marker='o',label=f'{hospitalisation} data', alpha=0.8) tau = taus[hospitalisation] prevalences = [] - for n in range(200): + for n in range(1000): prevalences.append(prevalence_func(time,tau,np.random.normal(p_AD_summary.loc[hospitalisation]['mean'], p_AD_summary.loc[hospitalisation]['sd']))) mean = np.mean(prevalences,axis=0) lower = np.quantile(prevalences,0.025,axis=0) upper = np.quantile(prevalences,0.975,axis=0) - ax.plot(time,mean,color_dict[hospitalisation]+'--',alpha=0.7) - axs[-1].plot(time,mean,color_dict[hospitalisation]+'--',alpha=0.7,label=f'{hospitalisation} mean fit') - ax.fill_between(time,lower, upper,alpha=0.2, color=color_dict[hospitalisation]) - ax.set_title(hospitalisation) - ax.grid(False) + ax.plot(time,100*mean,color='black',linestyle='--',linewidth=1.5,alpha=1) + ax.plot(time,100*lower,color='black',linestyle='-',linewidth=1,alpha=1) + ax.plot(time,100*upper,color='black',linestyle='-',linewidth=1,alpha=1) -axs[-1].set_xlabel('Months after infection') -axs[-2].set_xlabel('Months after infection') -axs[0].set_ylabel('Prevalence') -axs[-2].set_ylabel('Prevalence') -axs[-1].legend(prop={'size': 12}) -axs[-1].grid(False) + ax.grid(False) + ax.set_xlabel('Months after infection',size=10) + ax.set_title(hospitalisation, size=10) + ax.tick_params(axis='both', which='major', labelsize=10) -fig.suptitle('Prevalence of long-COVID symptoms') -fig.savefig(os.path.join(abs_dir,fig_result_folder,'prevalence_MCMC_fit.png')) +axes[0].set_ylabel('Symptom prevalence (%)', size=10) +plt.tight_layout() +fig.savefig(os.path.join(abs_dir,fig_result_folder,'prevalence_MCMC_fit.pdf')) ######### ## QoL ## @@ -283,18 +287,20 @@ def log_probability(theta, tau, x, y, bounds): QoL_Belgium_func = Life_table.QoL_Belgium_func # visualise fit -fig,ax = plt.subplots() +fig,ax = plt.subplots(figsize=(0.75*8.3,0.25*11.7)) for index in QoL_Belgium.index: left = index.left right = index.right w = right-left ax.bar(left+w/2,QoL_Belgium[index],w-1,color='grey',alpha=0.5,label='data') -ax.plot(QoL_Belgium_func(LE_table.index.values),label='fit',color='b') -ax.set_xlabel('Age') -ax.set_ylabel('QoL Belgium') +ax.plot(QoL_Belgium_func(LE_table.index.values),label='fit',color='black') +ax.set_xlabel('Age (years)',size=10) +ax.set_ylabel('QoL (-)',size=10) ax.set_ylim([0.5, 0.9]) ax.grid(False) -fig.savefig(os.path.join(abs_dir,fig_result_folder,'QoL_Belgium_fit.png')) +ax.tick_params(axis='both', which='major', labelsize=10) +plt.tight_layout() +fig.savefig(os.path.join(abs_dir,fig_result_folder,'QoL_Belgium_fit.pdf')) ####################### ## Average QALY loss ## @@ -309,7 +315,7 @@ def QALY_loss_func(t,tau,p_AD,age,QoL_after): beta = QoL_Belgium_func(age+t/12)-QoL_after return prevalence_func(t,tau,p_AD) * max(0,beta) -draws = 200 +draws = 1000 # Pre-allocate new multi index series with index=hospitalisation,age,draw multi_index = pd.MultiIndex.from_product([hospitalisation_groups+['Non-hospitalised (no AD)'],np.arange(draws),LE_table.index.values],names=['hospitalisation','draw','age']) @@ -337,7 +343,7 @@ def QALY_loss_func(t,tau,p_AD,age,QoL_after): QoL_after = QoL_Belgium_func(age)-beta # integrate QALY_loss_func from 0 to LE QALY_loss = quad(QALY_loss_func,0,LE,args=(tau,p_AD,age,QoL_after))[0]/12 - average_QALY_losses[idx] = QALY_loss + average_QALY_losses.iloc[idx] = QALY_loss # save result to dataframe def get_lower(x): @@ -360,29 +366,20 @@ def get_sd(x): average_QALY_losses_summary.to_csv(os.path.join(abs_dir,data_result_folder,'average_QALY_losses.csv')) # Visualise results -fig,axs = plt.subplots(2,2,sharex=True,sharey=True,figsize=(10,10)) -axs = axs.reshape(-1) -for ax,hospitalisation in zip(axs,['Non-hospitalised (no AD)']+hospitalisation_groups): +fig,axes = plt.subplots(nrows=1,ncols=3,figsize=(8.3,0.25*11.7),sharey=True) + +for ax,hospitalisation in zip(axes,hospitalisation_groups): mean = average_QALY_losses_summary.loc[hospitalisation]['mean'] lower = average_QALY_losses_summary.loc[hospitalisation]['lower'] upper = average_QALY_losses_summary.loc[hospitalisation]['upper'] - ax.plot(LE_table.index.values,mean,color_dict[hospitalisation],label=f'{hospitalisation}') - ax.fill_between(LE_table.index.values,lower,upper,alpha=0.20, color = color_dict[hospitalisation]) - ax.set_title(hospitalisation) + ax.plot(LE_table.index.values,mean,color='black',linewidth=1.5, linestyle='--') + ax.plot(LE_table.index.values,lower,color='black',linewidth=1, linestyle='-') + ax.plot(LE_table.index.values,upper,color='black',linewidth=1, linestyle='-') + ax.set_title(hospitalisation, size=10) ax.grid(False) + ax.set_xlabel('Age of infection (years)', size=10) + ax.tick_params(axis='both', which='major', labelsize=10) -axs[-1].set_xlabel('Age when infected') -axs[-2].set_xlabel('Age when infected') -axs[0].set_ylabel('Average QALY loss') -axs[-2].set_ylabel('Average QALY loss') - -fig.suptitle('Average QALY loss related to long-COVID') -fig.savefig(os.path.join(abs_dir,fig_result_folder,'average_QALY_losses.png')) - -# QALY losses due COVID death -QALY_D_per_age = Life_table.compute_QALY_D_x() -fig,ax = plt.subplots(figsize=(12,4)) -ax.plot(QALY_D_per_age,'b') -ax.set(xlabel='age', ylabel=r'$QALY_D$') -ax.grid(False) -fig.savefig(os.path.join(abs_dir,fig_result_folder,'QALY_D.png')) \ No newline at end of file +axes[0].set_ylabel('Average QALY loss', size=10) +plt.tight_layout() +fig.savefig(os.path.join(abs_dir,fig_result_folder,'average_QALY_losses_per_age.pdf')) diff --git a/results/covid19_DTM/preprocessing/QALY_model/long_COVID/QALY_D.png b/results/covid19_DTM/preprocessing/QALY_model/long_COVID/QALY_D.png deleted file mode 100644 index 745d8d993..000000000 Binary files a/results/covid19_DTM/preprocessing/QALY_model/long_COVID/QALY_D.png and /dev/null differ diff --git a/results/covid19_DTM/preprocessing/QALY_model/long_COVID/QoL_Belgium_fit.pdf b/results/covid19_DTM/preprocessing/QALY_model/long_COVID/QoL_Belgium_fit.pdf new file mode 100644 index 000000000..d8fcd9ee7 Binary files /dev/null and b/results/covid19_DTM/preprocessing/QALY_model/long_COVID/QoL_Belgium_fit.pdf differ diff --git a/results/covid19_DTM/preprocessing/QALY_model/long_COVID/QoL_Belgium_fit.png b/results/covid19_DTM/preprocessing/QALY_model/long_COVID/QoL_Belgium_fit.png deleted file mode 100644 index dc762c1fa..000000000 Binary files a/results/covid19_DTM/preprocessing/QALY_model/long_COVID/QoL_Belgium_fit.png and /dev/null differ diff --git a/results/covid19_DTM/preprocessing/QALY_model/long_COVID/QoL_decrease.pdf b/results/covid19_DTM/preprocessing/QALY_model/long_COVID/QoL_decrease.pdf new file mode 100644 index 000000000..677c9c466 Binary files /dev/null and b/results/covid19_DTM/preprocessing/QALY_model/long_COVID/QoL_decrease.pdf differ diff --git a/results/covid19_DTM/preprocessing/QALY_model/long_COVID/average_QALY_losses.png b/results/covid19_DTM/preprocessing/QALY_model/long_COVID/average_QALY_losses.png deleted file mode 100644 index 5961a5ab7..000000000 Binary files a/results/covid19_DTM/preprocessing/QALY_model/long_COVID/average_QALY_losses.png and /dev/null differ diff --git a/results/covid19_DTM/preprocessing/QALY_model/long_COVID/average_QALY_losses_per_age.pdf b/results/covid19_DTM/preprocessing/QALY_model/long_COVID/average_QALY_losses_per_age.pdf new file mode 100644 index 000000000..c62cbd61a Binary files /dev/null and b/results/covid19_DTM/preprocessing/QALY_model/long_COVID/average_QALY_losses_per_age.pdf differ diff --git a/results/covid19_DTM/preprocessing/QALY_model/long_COVID/prevalence_MCMC_fit.pdf b/results/covid19_DTM/preprocessing/QALY_model/long_COVID/prevalence_MCMC_fit.pdf new file mode 100644 index 000000000..2cc2186a7 Binary files /dev/null and b/results/covid19_DTM/preprocessing/QALY_model/long_COVID/prevalence_MCMC_fit.pdf differ diff --git a/results/covid19_DTM/preprocessing/QALY_model/long_COVID/prevalence_MCMC_fit.png b/results/covid19_DTM/preprocessing/QALY_model/long_COVID/prevalence_MCMC_fit.png deleted file mode 100644 index 55390b085..000000000 Binary files a/results/covid19_DTM/preprocessing/QALY_model/long_COVID/prevalence_MCMC_fit.png and /dev/null differ diff --git a/results/covid19_DTM/preprocessing/QALY_model/long_COVID/prevalence_first_fit.pdf b/results/covid19_DTM/preprocessing/QALY_model/long_COVID/prevalence_first_fit.pdf new file mode 100644 index 000000000..61e9a5e13 Binary files /dev/null and b/results/covid19_DTM/preprocessing/QALY_model/long_COVID/prevalence_first_fit.pdf differ diff --git a/results/covid19_DTM/preprocessing/QALY_model/long_COVID/prevalence_first_fit_AD.png b/results/covid19_DTM/preprocessing/QALY_model/long_COVID/prevalence_first_fit_AD.png deleted file mode 100644 index f563a49a2..000000000 Binary files a/results/covid19_DTM/preprocessing/QALY_model/long_COVID/prevalence_first_fit_AD.png and /dev/null differ diff --git a/results/covid19_DTM/preprocessing/QALY_model/long_COVID/prevalence_first_fit_no_AD.png b/results/covid19_DTM/preprocessing/QALY_model/long_COVID/prevalence_first_fit_no_AD.png deleted file mode 100644 index f11df7b4c..000000000 Binary files a/results/covid19_DTM/preprocessing/QALY_model/long_COVID/prevalence_first_fit_no_AD.png and /dev/null differ diff --git a/src/covid19_DTM/data/model_parameters.py b/src/covid19_DTM/data/model_parameters.py index f91b3956d..ec844bc11 100644 --- a/src/covid19_DTM/data/model_parameters.py +++ b/src/covid19_DTM/data/model_parameters.py @@ -411,7 +411,7 @@ def errorfcn(multiplier, n, n_desired): if not agg: # Set the average values for beta, seasonality, contact effectivities and mentality according to 'BASE' calibration dictionary samples_path = '../../../data/covid19_DTM/interim/model_parameters/calibrations/national/' - base_dict_name = 'national_REF_SAMPLES_2023-06-09.json' + base_dict_name = 'national_REF_SAMPLES_2023-06-28.json' base_samples_dict = load_samples_dict(os.path.join(abs_dir, samples_path+base_dict_name), age_stratification_size=age_stratification_size) pars_dict.update({ 'eff_home': 1, diff --git a/src/covid19_DTM/data/utils.py b/src/covid19_DTM/data/utils.py index c966e8bdf..602d825a2 100644 --- a/src/covid19_DTM/data/utils.py +++ b/src/covid19_DTM/data/utils.py @@ -39,10 +39,10 @@ def construct_initN(age_classes=None, agg=None): if age_classes is not None: age_struct['age_class'] = pd.cut(age_struct.age, bins=age_classes) - age_piramid = age_struct.groupby(['NIS','age_class']).sum().reset_index() + age_piramid = age_struct.groupby(['NIS','age_class'], observed=False).sum().reset_index() initN = age_piramid.pivot(index='NIS', columns='age_class', values='number') else: - age_piramid = age_struct.groupby(['NIS','age']).sum().reset_index() + age_piramid = age_struct.groupby(['NIS','age'], observed=False).sum().reset_index() initN = age_piramid.pivot(index='NIS', columns='age', values='number') initN = initN.fillna(0) diff --git a/src/covid19_DTM/models/QALY.py b/src/covid19_DTM/models/QALY.py index ddd05ea88..9b360c71d 100644 --- a/src/covid19_DTM/models/QALY.py +++ b/src/covid19_DTM/models/QALY.py @@ -56,7 +56,7 @@ def __init__(self, comorbidity_parameters=None): # Define absolute path abs_dir = os.path.dirname(__file__) # Import life table (q_x) - self.life_table = pd.read_csv(os.path.join(abs_dir, '../../../data/covid19_DTM/interim/QALY_model/Life_table_Belgium_2019.csv'),sep=',',index_col=0) + self.life_table = pd.read_csv(os.path.join(abs_dir, '../../../data/covid19_DTM/interim/QALY_model/life_table_model/Life_table_Belgium_2019.csv'),sep=',',index_col=0) # Define mu_x explictly to enhance readability of the code self.mu_x = self.life_table['mu_x'] # Define overall Belgian QoL scores (source: ...) @@ -460,9 +460,9 @@ def lost_QALYs(out,AD_non_hospitalised=False): out_enlarged = xr.concat([out]*int(np.ceil(100/sim_draws)), dim='draws') if AD_non_hospitalised: - hospitalisation_groups = ['Non-hospitalised','Cohort','ICU'] + hospitalisation_groups = ['Non-hospitalised','Hospitalised (no IC)','Hospitalised (IC)'] else: - hospitalisation_groups = ['Non-hospitalised (no AD)','Cohort','ICU'] + hospitalisation_groups = ['Non-hospitalised (no AD)','Hospitalised (no IC)','Hospitalised (IC)'] # Load average QALY losses abs_dir = os.path.dirname(__file__) @@ -474,7 +474,7 @@ def lost_QALYs(out,AD_non_hospitalised=False): multi_index = pd.MultiIndex.from_product([hospitalisation_groups,ages,np.arange(QALY_draws)],names=['hospitalisation','age','draw']) QALY_long_COVID_per_age = pd.Series(index = multi_index, dtype=float) for idx,(hospitalisation,age,draw) in enumerate(QALY_long_COVID_per_age.index): - QALY_long_COVID_per_age[idx] = np.random.normal(average_QALY_losses['mean'][hospitalisation,age], + QALY_long_COVID_per_age.iloc[idx] = np.random.normal(average_QALY_losses['mean'][hospitalisation,age], average_QALY_losses['sd'][hospitalisation,age]) # bin data