PCR positvity rates

trying URL 'https://cran.ma.imperial.ac.uk/bin/macosx/big-sur-arm64/contrib/4.1/tsibble_1.1.1.tgz'
Content type 'application/x-gzip' length 1047751 bytes (1023 KB)
==================================================
downloaded 1023 KB

The downloaded binary packages are in
    /var/folders/bk/jrqs03tx5mq9s28mhml5xzhm0000gn/T//RtmpehX7ad/downloaded_packages

tsibble installed
also installing the dependency ‘fabletools’

trying URL 'https://cran.ma.imperial.ac.uk/bin/macosx/big-sur-arm64/contrib/4.1/fabletools_0.3.2.tgz'
Content type 'application/x-gzip' length 661125 bytes (645 KB)
==================================================
downloaded 645 KB

trying URL 'https://cran.ma.imperial.ac.uk/bin/macosx/big-sur-arm64/contrib/4.1/feasts_0.2.2.tgz'
Content type 'application/x-gzip' length 1255459 bytes (1.2 MB)
==================================================
downloaded 1.2 MB

The downloaded binary packages are in
    /var/folders/bk/jrqs03tx5mq9s28mhml5xzhm0000gn/T//RtmpehX7ad/downloaded_packages

feasts installed
trying URL 'https://cran.ma.imperial.ac.uk/bin/macosx/big-sur-arm64/contrib/4.1/fable_0.3.1.tgz'
Content type 'application/x-gzip' length 1048077 bytes (1023 KB)
==================================================
downloaded 1023 KB

The downloaded binary packages are in
    /var/folders/bk/jrqs03tx5mq9s28mhml5xzhm0000gn/T//RtmpehX7ad/downloaded_packages

fable installed

Get positivity data

Series: uniqueCasePositivityBySpecimenDateRollingSum 
Model: ETS(M,A,N) 
  Smoothing parameters:
    alpha = 0.9998995 
    beta  = 0.1506035 

  Initial states:
     l[0]     b[0]
 1.703391 2.926506

  sigma^2:  0.024

     AIC     AICc      BIC 
3674.047 3674.139 3696.485 

Age specific cases

Healthcare data by region

`geom_smooth()` using formula 'y ~ s(x, bs = "cs")'
Warning: Removed 14 rows containing non-finite values (stat_smooth).
Warning: Removed 14 row(s) containing missing values (geom_path).

Warning: Removed 1 row(s) containing missing values (geom_path).

LS0tCnRpdGxlOiAiUiBOb3RlYm9vayIKb3V0cHV0OiBodG1sX25vdGVib29rCi0tLQoKIyMgUENSIHBvc2l0dml0eSByYXRlcwoKYGBge3J9CmtuaXRyOjpvcHRzX2NodW5rJHNldChjYWNoZSA9IFRSVUUsIGVjaG8gPSBGQUxTRSwgbWVzc2FnZSA9IEZBTFNFLCB3YXJuaW5nID0gRkFMU0UpCgpsaWJyYXJ5KHBhY21hbikKcF9sb2FkKHRpZHl2ZXJzZSwgdHNpYmJsZSwgZmVhc3RzLCBmYWJsZSwgZmFibGV0b29scywgc2FzcywgdXJjYSkKCgpgYGAKCiMjIyBHZXQgcG9zaXRpdml0eSBkYXRhCgpgYGB7cn0KCmRmMSA8LSByZWFkX2NzdigiaHR0cHM6Ly9hcGkuY29yb25hdmlydXMuZGF0YS5nb3YudWsvdjIvZGF0YT9hcmVhVHlwZT1sdGxhJm1ldHJpYz11bmlxdWVDYXNlUG9zaXRpdml0eUJ5U3BlY2ltZW5EYXRlUm9sbGluZ1N1bSZmb3JtYXQ9Y3N2Iiwgc2hvd19jb2xfdHlwZXMgPSBGQUxTRSkKCmRmMQoKIyMgY29lcmNlIHRvIHRzaWJibGUKCgpkZjEgPC0gZGYxICU+JQogIGFzX3RzaWJibGUoa2V5ID0gYyhhcmVhTmFtZSwgYXJlYUNvZGUpLCBpbmRleCA9IGRhdGUpCgpkZjEgPC0gZGYxICU+JQogIHRzaWJibGU6OmZpbGxfZ2FwcyguZnVsbCA9IEZBTFNFKSAlPiUKICBkcm9wX25hKCkKCmRmMQoKcCA8LSBkZjEgJT4lCiAgZmlsdGVyKHN0cl9kZXRlY3QoYXJlYU5hbWUsICJOb3J0aCIpKSAlPiUKICBhdXRvcGxvdCgpICsKICB0aGVtZShsZWdlbmQucG9zaXRpb24gPSAiYm90dG9tIikKCmxpYnJhcnkocGxvdGx5KQoKZ2dwbG90bHkocCkKCmBgYAoKCgoKCmBgYHtyfQoKCm1vZCA8LSBUU0xNKHVuaXF1ZUNhc2VQb3NpdGl2aXR5QnlTcGVjaW1lbkRhdGVSb2xsaW5nU3VtIH4gdHJlbmQoKSkKCmZpdCA8LSBkZjEgJT4lCiAgIyBkcm9wX25hKCkgJT4lCiAgZmFibGV0b29sczo6bW9kZWwoZXRzID0gRVRTKHVuaXF1ZUNhc2VQb3NpdGl2aXR5QnlTcGVjaW1lbkRhdGVSb2xsaW5nU3VtIH4gdHJlbmQoIkEiKSksIAogICAgICAgICAgICAgICAgICAgIGFyaW1hID0gQVJJTUEodW5pcXVlQ2FzZVBvc2l0aXZpdHlCeVNwZWNpbWVuRGF0ZVJvbGxpbmdTdW0pKQoKZml0ICU+JQogIGF1Z21lbnQoKSAlPiUKICBmaWx0ZXIoIWlzLm5hKC5maXR0ZWQpLCAKICAgICAgICAgYXJlYU5hbWUgPT0gIkxld2lzaGFtIikgJT4lCiAgZ2dwbG90KCkgKwogIGdlb21fcG9pbnQoYWVzKGRhdGUsIHVuaXF1ZUNhc2VQb3NpdGl2aXR5QnlTcGVjaW1lbkRhdGVSb2xsaW5nU3VtKSwgcGNoID0gMSkgKwogIGdlb21fbGluZShhZXMoZGF0ZSwgeSA9IC5maXR0ZWQpLCBjb2xvdXIgPSAicmVkIiwgbHR5ID0gImRvdHRlZCIpICsKICBmYWNldF93cmFwKH4gLm1vZGVsKQogIApmaXQgJT4lCiAgZmlsdGVyKGFyZWFOYW1lID09ICJCbGFja3Bvb2wiKSAlPiUKICBzZWxlY3QoZXRzKSAlPiUKICByZXBvcnQoKQoKJT4lCiAgYWNjdXJhY3koKQogIAoKYGBgCgoKCmBgYHtyfQoKZGYxICU+JQogIGZlYXR1cmVzKGxpc3QobWVhbiA9IG1lYW4pKQoKCmRmMiA8LSBkZjEgJT4lCiAgZ3JvdXBfYnlfa2V5KGMoImFyZWFOYW1lIiwgImFyZWFDb2RlIikpICU+JQogIGluZGV4X2J5KHllYXJfd2VlayA9IH4geWVhcndlZWsoLikpICU+JQogIHN1bW1hcmlzZSgKICAgIAogICAgbWVhbl9wb3MgPSByb3VuZChtZWFuKHVuaXF1ZUNhc2VQb3NpdGl2aXR5QnlTcGVjaW1lbkRhdGVSb2xsaW5nU3VtLCBuYS5ybSA9IFRSVUUpLCAyKSkKICAgIAogICkKICAKZGYyICU+JQogIHNlcGFyYXRlKHllYXJfd2VlaywgYygieWVhciIsICJ3ZWVrIiksICJcXHMiKSAlPiUKICBmaWx0ZXIoc3RyX2RldGVjdChhcmVhQ29kZSwgIkUwNiIpKSAlPiUKICBnZ3Bsb3QoYWVzKHdlZWssIG1lYW5fcG9zLCBjb2xvdXIgPSB5ZWFyLCAgZ3JvdXAgPSB5ZWFyKSkgKwogIGdlb21fbGluZSgpICsKICBzY2FsZV94X2Rpc2NyZXRlKCkgKwogIGZhY2V0X3dyYXAofmFyZWFOYW1lKQoKCmRmMiAlPiUKICBmaWx0ZXIoYXJlYU5hbWUgJWluJSBjKCJCYWJlcmdoIiwgIkFkdXIiLCAiU3V0dG9uIiwgIlRvd2VyIEhhbWxldHMiKSkgJT4lCiAgYXV0b3Bsb3QoKQoKZGYyICU+JQogIGZlYXR1cmVzKG1lYW5fcG9zLCBmZWF0X3N0bCkgJT4lCiAgZmlsdGVyKGxpbmVhcml0eSA9PSBtaW4obGluZWFyaXR5KSkKICAjZmlsdGVyKGFyZWFOYW1lICVpbiUgYygiQmFiZXJnaCIsICJBZHVyIiwgIlN1dHRvbiIsICJUb3dlciBIYW1sZXRzIikpICU+JQogIGdncGxvdChhZXMoeCA9IGN1cnZhdHVyZSwgeSA9IHNwaWtpbmVzcykpICsKICBnZW9tX3BvaW50KCkgCgoKYGBgCgpgYGB7cn0KCmRmMSAlPiUKICBmaWx0ZXIoYXJlYU5hbWUgJWluJSBjKCJCYWJlcmdoIikpICU+JQogIEFDRih1bmlxdWVDYXNlUG9zaXRpdml0eUJ5U3BlY2ltZW5EYXRlUm9sbGluZ1N1bSwgbGFnX21heCA9IDQ4KSAlPiUKICBhdXRvcGxvdCgpCgpgYGAKCmBgYHtyfQpkY21wIDwtIGRmMiAlPiUKICBmaWx0ZXIoYXJlYU5hbWUgJWluJSBjKCJMYW1iZXRoIikpIAoKZGNtcCAlPiUgCiAgbXV0YXRlKAogICAgYDctTUFgID0gc2xpZGVyOjpzbGlkZV9kYmwobWVhbl9wb3MsIG1lYW4sCiAgICAgICAgICAgICAgICAuYmVmb3JlID0gMywgLmFmdGVyID0gMywgLmNvbXBsZXRlID0gVFJVRSkKICApICU+JQogIGF1dG9wbG90KCkgKwogIGdlb21fbGluZShhZXMoeSA9IGA3LU1BYCksIGNvbG91ciA9ICJyZWQiKQoKCmRmMSAlPiUKICBtb2RlbCgKICAgIGNsYXNzaWNhbF9kZWNvbXBvc2l0aW9uKHVuaXF1ZUNhc2VQb3NpdGl2aXR5QnlTcGVjaW1lbkRhdGVSb2xsaW5nU3VtLCB0eXBlID0gImFkZGl0aXZlIikKICApICU+JQogIGNvbXBvbmVudHMoKSAlPiUKICBhdXRvcGxvdCgpIAoKYGBgCgojIyBBZ2Ugc3BlY2lmaWMgY2FzZXMKCmBgYHtyfQoKbGlicmFyeShkYXRhLnRhYmxlKQoKYXNfY2FzZXMgPC0gZnJlYWQoImh0dHBzOi8vYXBpLmNvcm9uYXZpcnVzLmRhdGEuZ292LnVrL3YyL2RhdGE/YXJlYVR5cGU9bHRsYSZtZXRyaWM9bmV3Q2FzZXNCeVNwZWNpbWVuRGF0ZUFnZURlbW9ncmFwaGljcyZmb3JtYXQ9Y3N2IikKCmFzX2Nhc2VzX3QgPC0gYXNfY2FzZXMgJT4lCiAgYXNfdHNpYmJsZShrZXkgPSBjKGFyZWFOYW1lLCBhcmVhQ29kZSwgYWdlKSwgaW5kZXggPSBkYXRlKQoKYXNfY2FzZXNfdAoKIyMgcmF0ZXMgaW4gc2Nob29sIGNoaWxkcmVuCgphc19jYXNlc190ICU+JQogIGZpbHRlcihhZ2UgJWluJSBjKCIxMF8xNCIsICIxNV8xOSIsICIyMF8yNCIpLCAKICAgICAgICAgZGF0ZSA+PSAiMjAyMS0xMS0wMSIsCiAgc3RyX2RldGVjdChhcmVhQ29kZSwgIkUwOSIpKSAlPiUKICBnZ3Bsb3QoYWVzKGRhdGUsIHJvbGxpbmdSYXRlLCBjb2xvdXIgPSBhZ2UpKSArCiAgZ2VvbV9saW5lKGFlcyhncm91cCA9IGFnZSkpICsKICBnZW9tX3ZsaW5lKHhpbnRlcmNlcHQgPSBhcy5EYXRlKCIyMDIxLTEyLTIwIikpICsKICBmYWNldF93cmFwKH5hcmVhTmFtZSkKCiMjIHJhdGVzIGluIG9sZGVyIGFkdWx0cwoKYXNfY2FzZXNfdCAlPiUKICBmaWx0ZXIoYWdlICVpbiUgYygiNjArIiksIAogICAgICAgICAjZGF0ZSA+PSAiMjAyMS0xMS0wMSIsCiAgc3RyX2RldGVjdChhcmVhQ29kZSwgIkUwOSIpKSAlPiUKICBnZ3Bsb3QoYWVzKGRhdGUsIHJvbGxpbmdSYXRlLCBjb2xvdXIgPSBhZ2UpKSArCiAgZ2VvbV9saW5lKGFlcyhncm91cCA9IGFnZSkpICsKICBmYWNldF93cmFwKH5hcmVhTmFtZSkgKwogIGdlb21fdmxpbmUoeGludGVyY2VwdCA9IGFzLkRhdGUoIjIwMjEtMTItMjAiKSkKCgpgYGAKCgoKYGBge3IgZm9yZWNhc3R9CgpkZjIgJT4lCiAgCiAgICBmaWx0ZXIoYXJlYU5hbWUgJWluJSBjKCJCYWJlcmdoIiwgIkFkdXIiLCAiU3V0dG9uIiwgIlRvd2VyIEhhbWxldHMiKSkgJT4lCiAgCiAgbW9kZWwoCiAgICBldHMgPSBFVFMoYm94X2NveCh1bmlxdWVDYXNlUG9zaXRpdml0eUJ5U3BlY2ltZW5EYXRlUm9sbGluZ1N1bSwgMC4zKSksCiAgICBhcmltYSA9IEFSSU1BKGxvZyh1bmlxdWVDYXNlUG9zaXRpdml0eUJ5U3BlY2ltZW5EYXRlUm9sbGluZ1N1bSkpLAogICAgc25haXZlID0gU05BSVZFKHVuaXF1ZUNhc2VQb3NpdGl2aXR5QnlTcGVjaW1lbkRhdGVSb2xsaW5nU3VtKQogICkgJT4lCiAgZm9yZWNhc3QoaCA9ICIyIHllYXJzIikgCgolPiUgCiAgYXV0b3Bsb3QoZmlsdGVyKGF1c19yZXRhaWwsIHllYXIoTW9udGgpID4gMjAxMCksIGxldmVsID0gTlVMTCkKCgoKCmBgYAoKCiMjIEhlYWx0aGNhcmUgZGF0YSBieSByZWdpb24KCgpgYGB7cn0KCmhvc3AgPC0gZnJlYWQoImh0dHBzOi8vYXBpLmNvcm9uYXZpcnVzLmRhdGEuZ292LnVrL3YyL2RhdGE/YXJlYVR5cGU9bmhzUmVnaW9uJm1ldHJpYz1ob3NwaXRhbENhc2VzJm1ldHJpYz1uZXdBZG1pc3Npb25zJm1ldHJpYz1jb3ZpZE9jY3VwaWVkTVZCZWRzJmZvcm1hdD1jc3YiKQoKaG9zcAoKaG9zcCAlPiUKICBmaWx0ZXIoZGF0ZSA+PSAiMjAyMS0xMS0wMSIpICU+JQogIGdncGxvdChhZXMoZGF0ZSwgbG9nKG5ld0FkbWlzc2lvbnMpLCBncm91cCA9IGFyZWFOYW1lKSApICsKICBnZW9tX2xpbmUoKSArCiAgZ2VvbV9zbW9vdGgobWV0aG9kID0gImdhbSIpICsKICBnZW9tX2xpbmUoYWVzKHkgPSBsb2coaG9zcGl0YWxDYXNlcykpKSArCiAgZmFjZXRfd3JhcCh+YXJlYU5hbWUpCgpgYGAKYGBge3J9CgoKCmFkX2FnZSA8LSBmcmVhZCgiaHR0cHM6Ly9hcGkuY29yb25hdmlydXMuZGF0YS5nb3YudWsvdjIvZGF0YT9hcmVhVHlwZT1uaHNSZWdpb24mbWV0cmljPWN1bUFkbWlzc2lvbnNCeUFnZSZmb3JtYXQ9Y3N2IikKCgphZF9hZ2UgJT4lCiAgYXNfdHNpYmJsZShrZXkgPSBjKGFyZWFDb2RlLCBhcmVhTmFtZSwgYWdlKSwgaW5kZXggPSBkYXRlKSAlPiUKICBncm91cF9ieShhZ2UsIGFyZWFOYW1lLCBhcmVhQ29kZSkgJT4lCiAgbXV0YXRlKGxlYWRfcmF0ZSA9IC1yYXRlICsgbGVhZChyYXRlKSkgJT4lCiAgZ2dwbG90KGFlcyhkYXRlLCBsZWFkX3JhdGUpKSArCiAgZ2VvbV9saW5lKCkgKwogICNnZW9tX3Ntb290aChtZXRob2QgPSAibG9lc3MiKSArCiAgZmFjZXRfZ3JpZChhZ2UgfiBhcmVhTmFtZSwgc2NhbGVzID0gImZyZWVfeSIpCgpgYGAKCg==