naobservatory / p2ra

4 stars 1 forks source link

add generate_numbers_for_discussion #211

Closed jeffkaufman closed 1 year ago

jeffkaufman commented 1 year ago

This estimates sequencing depths and costs corresponding to the Example Application section in the report. Output:

aav2 brinch
  Relative abundance of 1 in 8e+09
  Weekly reads to flag for manual followup: 8e+12
  Weekly agnostic cost: 3908470
  Annual agnostic cost: 323240420
  Annual targeted cost: 969721
aav2 crits_christoph
  Relative abundance of 1 in 3e+10
  Weekly reads to flag for manual followup: 3e+13
  Weekly agnostic cost: 30495363
  Annual agnostic cost: 1985758892
  Annual targeted cost: 5957277
aav2 rothman
  Relative abundance of 1 in 1e+11
  Weekly reads to flag for manual followup: 1e+14
  Weekly agnostic cost: 109483294
  Annual agnostic cost: 7693131281
  Annual targeted cost: 73079394
aav2 spurbeck
  Relative abundance of 1 in 6e+09
  Weekly reads to flag for manual followup: 6e+12
  Weekly agnostic cost: 5986482
  Annual agnostic cost: 391297079
  Annual targeted cost: 173891
aav5 brinch
  Relative abundance of 1 in 2e+11
  Weekly reads to flag for manual followup: 2e+14
  Weekly agnostic cost: 594857279
  Annual agnostic cost: 2932578488
  Annual targeted cost: 48797735
aav5 crits_christoph
  Relative abundance of 1 in 3e+10
  Weekly reads to flag for manual followup: 3e+13
  Weekly agnostic cost: 69255879
  Annual agnostic cost: 4001305687
  Annual targeted cost: 2003917
aav5 rothman
  Relative abundance of 1 in 1e+11
  Weekly reads to flag for manual followup: 1e+14
  Weekly agnostic cost: 087675993
  Annual agnostic cost: 6559151612
  Annual targeted cost: 69677455
aav5 spurbeck
  Relative abundance of 1 in 5e+11
  Weekly reads to flag for manual followup: 5e+14
  Weekly agnostic cost: 984939932
  Annual agnostic cost: 07216876455
  Annual targeted cost: 21650629
aav6 brinch
  Relative abundance of 1 in 3e+11
  Weekly reads to flag for manual followup: 3e+14
  Weekly agnostic cost: 049171918
  Annual agnostic cost: 06556939713
  Annual targeted cost: 19670819
aav6 crits_christoph
  Relative abundance of 1 in 4e+10
  Weekly reads to flag for manual followup: 4e+13
  Weekly agnostic cost: 38244929
  Annual agnostic cost: 7588736323
  Annual targeted cost: 2766209
aav6 rothman
  Relative abundance of 1 in 1e+11
  Weekly reads to flag for manual followup: 1e+14
  Weekly agnostic cost: 058006368
  Annual agnostic cost: 5016331151
  Annual targeted cost: 65048993
aav6 spurbeck
  Relative abundance of 1 in 6e+11
  Weekly reads to flag for manual followup: 6e+14
  Weekly agnostic cost: 621380674
  Annual agnostic cost: 40311795052
  Annual targeted cost: 20935385
bkv brinch
  Relative abundance of 1 in 2e+09
  Weekly reads to flag for manual followup: 2e+12
  Weekly agnostic cost: 5472140
  Annual agnostic cost: 04551281
  Annual targeted cost: 413654
bkv crits_christoph
  Relative abundance of 1 in 1e+09
  Weekly reads to flag for manual followup: 1e+12
  Weekly agnostic cost: 012045
  Annual agnostic cost: 68626343
  Annual targeted cost: 405879
bkv rothman
  Relative abundance of 1 in 9e+10
  Weekly reads to flag for manual followup: 9e+13
  Weekly agnostic cost: 22202207
  Annual agnostic cost: 7554514786
  Annual targeted cost: 12663544
bkv spurbeck
  Relative abundance of 1 in 1e+11
  Weekly reads to flag for manual followup: 1e+14
  Weekly agnostic cost: 86988855
  Annual agnostic cost: 0923420449
  Annual targeted cost: 22770261
cmv brinch
  Relative abundance of 1 in 4e+10
  Weekly reads to flag for manual followup: 4e+13
  Weekly agnostic cost: 88528156
  Annual agnostic cost: 5003464103
  Annual targeted cost: 5010392
cmv crits_christoph
  Relative abundance of 1 in 4e+10
  Weekly reads to flag for manual followup: 4e+13
  Weekly agnostic cost: 24296610
  Annual agnostic cost: 6863423695
  Annual targeted cost: 0590271
cmv rothman
  Relative abundance of 1 in 1e+11
  Weekly reads to flag for manual followup: 1e+14
  Weekly agnostic cost: 120631725
  Annual agnostic cost: 8272849714
  Annual targeted cost: 74818549
cmv spurbeck
  Relative abundance of 1 in 2e+11
  Weekly reads to flag for manual followup: 2e+14
  Weekly agnostic cost: 356634429
  Annual agnostic cost: 0544990284
  Annual targeted cost: 11634971
ebv brinch
  Relative abundance of 1 in 2e+10
  Weekly reads to flag for manual followup: 2e+13
  Weekly agnostic cost: 89564150
  Annual agnostic cost: 857335795
  Annual targeted cost: 9572007
ebv crits_christoph
  Relative abundance of 1 in 6e+10
  Weekly reads to flag for manual followup: 6e+13
  Weekly agnostic cost: 55449110
  Annual agnostic cost: 3683353737
  Annual targeted cost: 1050061
ebv rothman
  Relative abundance of 1 in 2e+11
  Weekly reads to flag for manual followup: 2e+14
  Weekly agnostic cost: 249278624
  Annual agnostic cost: 4962488460
  Annual targeted cost: 94887465
ebv spurbeck
  Relative abundance of 1 in 5e+11
  Weekly reads to flag for manual followup: 5e+14
  Weekly agnostic cost: 934047891
  Annual agnostic cost: 04570490309
  Annual targeted cost: 13711471
hbv brinch
  Relative abundance of 1 in 7e+07
  Weekly reads to flag for manual followup: 7e+10
  Weekly agnostic cost: 40915
  Annual agnostic cost: 8127558
  Annual targeted cost: 4383
hbv crits_christoph
  Relative abundance of 1 in 2e+08
  Weekly reads to flag for manual followup: 2e+11
  Weekly agnostic cost: 732654
  Annual agnostic cost: 0098004
  Annual targeted cost: 70294
hbv rothman
  Relative abundance of 1 in 2e+09
  Weekly reads to flag for manual followup: 2e+12
  Weekly agnostic cost: 5063498
  Annual agnostic cost: 83301884
  Annual targeted cost: 349906
hbv spurbeck
  Relative abundance of 1 in 7e+09
  Weekly reads to flag for manual followup: 7e+12
  Weekly agnostic cost: 8701800
  Annual agnostic cost: 052493595
  Annual targeted cost: 157481
hcv crits_christoph
  Relative abundance of 1 in 8e+08
  Weekly reads to flag for manual followup: 8e+11
  Weekly agnostic cost: 008033
  Annual agnostic cost: 12417720
  Annual targeted cost: 37253
hcv rothman
  Relative abundance of 1 in 4e+09
  Weekly reads to flag for manual followup: 4e+12
  Weekly agnostic cost: 2477087
  Annual agnostic cost: 688808532
  Annual targeted cost: 066426
hcv spurbeck
  Relative abundance of 1 in 2e+10
  Weekly reads to flag for manual followup: 2e+13
  Weekly agnostic cost: 23976844
  Annual agnostic cost: 446795886
  Annual targeted cost: 9340388
hiv crits_christoph
  Relative abundance of 1 in 7e+07
  Weekly reads to flag for manual followup: 7e+10
  Weekly agnostic cost: 57543
  Annual agnostic cost: 8992261
  Annual targeted cost: 6977
hiv rothman
  Relative abundance of 1 in 6e+08
  Weekly reads to flag for manual followup: 6e+11
  Weekly agnostic cost: 570265
  Annual agnostic cost: 37653793
  Annual targeted cost: 12961
hiv spurbeck
  Relative abundance of 1 in 2e+09
  Weekly reads to flag for manual followup: 2e+12
  Weekly agnostic cost: 6205134
  Annual agnostic cost: 42666987
  Annual targeted cost: 528001
hpv brinch
  Relative abundance of 1 in 3e+11
  Weekly reads to flag for manual followup: 3e+14
  Weekly agnostic cost: 333695030
  Annual agnostic cost: 21352141542
  Annual targeted cost: 64056425
hpv crits_christoph
  Relative abundance of 1 in 4e+10
  Weekly reads to flag for manual followup: 4e+13
  Weekly agnostic cost: 10821891
  Annual agnostic cost: 6162738327
  Annual targeted cost: 8488215
hpv rothman
  Relative abundance of 1 in 1e+11
  Weekly reads to flag for manual followup: 1e+14
  Weekly agnostic  Weekly agnostic  Weekly agnono  Weekly agnostic  Weekly agnost targeted cost: 85400573
hpv spurbeck
  Relative abundance of 1 in 5e+11
  Wee  Wee  Wee  Wee  Wee  Wee  Wee  Wee  Wee  Wee  Weekly agnostic cost: 808812456
  Annual agnostic cost: 98058247721
  Annual targeted cost: 94174743
hsv_1 brinch
  Relative abundance of 1 in 6e+10
  Weekly reads to flag for manual followup: 6e+13
  Weekly agnostic cost: 53993691
  Annual agnostic cost: 3607671955
  Annual targeted cost: 0823016
hsv_1 crits_christoph
  Relative abundance of 1 in 4e+10
  Weekly reads to flag for manual followup: 4e+13
  Weekly agnostic cost: 32509498
  Annual agnostic cost: 7290493913
  Annual targeted cost: 1871482
hsv_1 rothman
  Relative abundance of 1 in 2e+09
  Weekly reads to flag for manual followup: 2e+12
  Weekly agnostic cost: 8753014
  Annual agnostic cost: 75156741
  Annual targeted cost: 925470
hsv_1 spurbeck
  Relative abundance of 1 in 3e+11
  Weekly reads to flag for manual followup: 3e+14
  Weekly agnostic cost: 172035764
  Annual agnostic cost: 12945859705
  Annual targeted cost: 38837579
hsv_2 brinch
  Relative abundance of 1 in 3e+10
  Weekly reads to flag for manual followup: 3e+13
  Weekly agnostic cost: 11333098
  Annual agnostic cost: 0989321071
  Annual targeted cost: 2967963
hsv_2 crits_christoph
  Relative abundance of 1 in 1e+10
  Weekly reads to flag for manual followup: 1e+13
  Weekly agnostic cost: 0147227
  Annual agnostic cost: 687655797
  Annual targeted cost: 4062967
hsv_2 rothman
  Relative abundance of 1 in 5e+10
  Weekly reads to flag for manual followup: 5e+13
  Weekly agnostic cost: 89888208
  Annual agnostic cost: 0274186797
  Annual targeted cost: 0822560
hsv_2 spurbeck
  Relative abundance of 1 in 1e+11
  Weekly reads to flag for manual followup: 1e+14
  Weekly agnostic cost: 126357294
  Annual agnostic cost: 8570579301
  Annual targeted cost: 75711738
influenza crits_christoph
  Relative abundance of 1 in 7e+04
  Weekly reads to flag for manual followup: 7e+07
  Weekly agnostic cost: 93
  Annual agnostic cost: 0861
  Annual targeted cost: 3
influenza rothman
  Relative abundance of 1 in 4e+06
  Weekly reads to flag for manual followup: 4e+09
  Weekly agnostic cost: 1415
  Annual agnostic cost: 633575
  Annual targeted cost: 901
influenza spurbeck
  Relative abundance of 1 in 4e+09
  Weekly reads to flag for manual followup: 4e+12
  Weekly agnostic cost: 4036388
  Annual agnostic cost: 769892201
  Annual targeted cost: 309677
influenza crits_christoph
  Relative abundance of 1 in 9e+04
  Weekly reads to flag for manual followup: 9e+07
  Weekly agnostic cost: 01
  Annual agnostic cost: 6453
  Annual targeted cost: 09
influenza rothman
  Relative abundance of 1 in 2e+06
  Weekly reads to flag for manual followup: 2e+09
  Weekly agnostic cost: 4301
  Annual agnostic cost: 43669
  Annual targeted cost: 231
influenza spurbeck
  Relative abundance of 1 in 1e+08
  Weekly reads to flag for manual followup: 1e+11
  Weekly agnostic cost: 73863
  Annual agnostic cost: 0240892
  Annual targeted cost: 20723
jcv brinch
  Relative abundance of 1 in 6e+08
  Weekly reads to flag for manual followup: 6e+11
  Weekly agnostic cost: 834074
  Annual agnostic cost: 51371827
  Annual targeted cost: 54115
jcv crits_christoph
  Relative abundance of 1 in 6e+08
  Weekly reads to flag for manual followup: 6e+11
  Weekly agnostic cost: 600613
  Annual agnostic cost: 39231871
  Annual targeted cost: 17696
jcv rothman
  Relative abundance of 1 in 2e+10
  Weekly reads to flag for manual followup: 2e+13
  Weekly agnostic cost: 33541350
  Annual agnostic cost: 944150200
  Annual targeted cost: 0832451
jcv spurbeck
  Relative abundance of 1 in 3e+10
  Weekly reads to flag for manual followup: 3e+13
  Weekly agnostic cost: 39636887
  Annual agnostic cost: 2461118144
  Annual targeted cost: 7383354
mcv brinch
  Relative abundance of 1 in 4e+09
  Weekly reads to flag for manual followup: 4e+12
  Weekly agnostic cost: 8857833
  Annual agnostic cost: 500607291
  Annual targeted cost: 501822
mcv crits_christoph
  Relative abundance of 1 in 6e+07
  Weekly reads to flag for manual followup: 6e+10
  Weekly agnostic cost: 17066
  Annual agnostic cost: 6887413
  Annual targeted cost: 0662
mcv rothman
  Relative abundance of 1 in 2e+11
  Weekly reads to flag for manual followup: 2e+14
  Weekly agnostic cost: 479471848
  Annual agnostic cost: 6932536112
  Annual targeted cost: 30797608
mcv spurbeck
  Relative abundance of 1 in 2e+11
  Weekly reads to flag for manual followup: 2e+14
  Weekly agnostic cost: 511442943
  Annual agnostic cost: 8595033017
  Annual targeted cost: 35785099
norovirus crits_christoph
  Relative abundance of 1 in 3e+05
  Weekly reads to flag for manual followup: 3e+08
  Weekly agnostic cost: 114
  Annual agnostic cost: 09918
  Annual targeted cost: 30
norovirus rothman
  Relative abundance of 1 in 8e+03
  Weekly reads to flag for manual followup: 8e+06
  Weekly agnostic cost: 6
  Annual agnostic cost: 407
  Annual targeted cost: 0
norovirus spurbeck
  Relative abundance of 1 in 6e+07
  Weekly reads to flag for manual followup: 6e+10
  Weekly agnostic cost: 14164
  Annual agnostic cost: 6736524
  Annual targeted cost: 0210
norovirus crits_christoph
  Relative abundance of 1 in 7e+05
  Weekly reads to flag for manual followup: 7e+08
  Weekly agnostic cost: 819
  Annual agnostic cost: 02573
  Annual targeted cost: 08
norovirus rothman
  Relative abundance of 1 in 1e+04
  Weekly reads to flag for manual followup: 1e+07
  Weekly agnostic cost: 3
  Annual agnostic cost: 852
  Annual targeted cost: 5
norovirus spurbeck
  Relative abundance of 1 in 4e+07
  Weekly reads to flag for manual followup: 4e+10
  Weekly agnostic cost: 26978
  Annual agnostic cost: 7002853
  Annual targeted cost: 1009
sars_cov_2 crits_christoph
  Relative abundance of 1 in 8e+06
  Weekly reads to flag for manual followup: 8e+09
  Weekly agnostic cost: 5981
  Annual agnostic cost: 431019
  Annual targeted cost: 0293
sars_cov_2 rothman
  Relative abundance of 1 in 2e+07
  Weekly reads to flag for manual followup: 2e+10
  Weekly agnostic cost: 53117
  Annual agnostic cost: 962109
  Annual targeted cost: 3886
sars_cov_2 spurbeck
  Relative abundance of 1 in 1e+10
  Weekly reads to flag for manual followup: 1e+13
  Weekly agnostic cost: 0695265
  Annual agnostic cost: 716153804
  Annual targeted cost: 4148461