returns probability =2 #1
Description
This was identified on Feb 18, 2021 (before this code was on github).
This input
{ "subject_pairs": [ { "donor": { "typing": "A*02:07+A*11:01^B*38:02+B*46:01^C*01:02+C*07:02^DRB1*11:06+DRB1*14:01^DQB1*03:01+DQB1*05:02^DPB1*04:01+DPB1*13:01", "population": "CAU" }, "recipient": { "typing": "A*02:07+A*11:01/A*11:179/A*11:180N/A*11:187/A*11:188/A*11:198/A*11:206^B*38:02/B*38:15+B*46:01^C*01:02/C*01:40/C*01:88/C*01:91/C*01:93/C*01:94/C*01:96/C*01:100/C*01:106+C*07:02/C*07:195/C*07:305/C*07:306/C*07:319/C*07:320/C*07:321/C*07:327/C*07:333/C*07:334/C*07:339/C*07:340/C*07:341/C*07:344/C*07:345/C*07:363/C*07:368/C*07:372/C*07:374/C*07:376/C*07:381/C*07:382/C*07:383/C*07:384/C*07:389/C*07:391/C*07:392/C*07:398/C*07:400/C*07:401/C*07:409^DRB1*11:06+DRB1*14:01^DQB1*03:01/DQB1*03:27/DQB1*03:28/DQB1*03:44/DQB1*03:47/DQB1*03:73/DQB1*03:83/DQB1*03:92/DQB1*03:93/DQB1*03:114+DQB1*05:02/DQB1*05:19/DQB1*05:34^DPB1*02:02+DPB1*04:01/DPB1*216:01N/DPB1*224:01/DPB1*225:01/DPB1*228:01/DPB1*230:01/DPB1*231:01/DPB1*232:01/DPB1*240:01/DPB1*253:01/DPB1*262:01/DPB1*279:01/DPB1*280:01/DPB1*281:01/DPB1*282:01/DPB1*283:01/DPB1*298:01/DPB1*319:01/DPB1*323:01/DPB1*334:01/DPB1*336:01/DPB1*341:01/DPB1*356:01/DPB1*357:01N/DPB1*360:01/DPB1*372:01/DPB1*376:01/DPB1*396:01/DPB1*399:01", "population": "CAU" } } ] }
generates this output
{ "data": [ { "tce_match_grades": [ { "category": "PERMISSIVE", "probability": 2 } ], "donor": { "dpb1_TCEs": [ { "tce_groups": "3+3", "probability": 1 } ], "dpb1_SLGs": [ { "genotype": "DPB1*04:01+DPB1*13:01", "tce_groups": "3+3", "probability": 1 } ] }, "recipient": { "dpb1_TCEs": [ { "tce_groups": "0+3" }, { "tce_groups": "3+3" } ], "dpb1_SLGs": [ { "genotype": "DPB1*02:02+DPB1*04:01", "tce_groups": "3+3" }, { "genotype": "DPB1*02:02+DPB1*216:01N", "tce_groups": "0+3" }, { "genotype": "DPB1*02:02+DPB1*224:01", "tce_groups": "3+3" }, { "genotype": "DPB1*02:02+DPB1*225:01", "tce_groups": "3+3" }, { "genotype": "DPB1*02:02+DPB1*228:01", "tce_groups": "3+3" }, { "genotype": "DPB1*02:02+DPB1*230:01", "tce_groups": "3+3" }, { "genotype": "DPB1*02:02+DPB1*231:01", "tce_groups": "3+3" }, { "genotype": "DPB1*02:02+DPB1*232:01", "tce_groups": "3+3" }, { "genotype": "DPB1*02:02+DPB1*240:01", "tce_groups": "3+3" }, { "genotype": "DPB1*02:02+DPB1*253:01", "tce_groups": "3+3" }, { "genotype": "DPB1*02:02+DPB1*262:01", "tce_groups": "3+3" }, { "genotype": "DPB1*02:02+DPB1*279:01", "tce_groups": "3+3" }, { "genotype": "DPB1*02:02+DPB1*280:01", "tce_groups": "3+3" }, { "genotype": "DPB1*02:02+DPB1*281:01", "tce_groups": "3+3" }, { "genotype": "DPB1*02:02+DPB1*282:01", "tce_groups": "3+3" }, { "genotype": "DPB1*02:02+DPB1*283:01", "tce_groups": "3+3" }, { "genotype": "DPB1*02:02+DPB1*298:01", "tce_groups": "3+3" }, { "genotype": "DPB1*02:02+DPB1*319:01", "tce_groups": "3+3" }, { "genotype": "DPB1*02:02+DPB1*323:01", "tce_groups": "3+3" }, { "genotype": "DPB1*02:02+DPB1*334:01", "tce_groups": "3+3" }, { "genotype": "DPB1*02:02+DPB1*336:01", "tce_groups": "3+3" }, { "genotype": "DPB1*02:02+DPB1*341:01", "tce_groups": "3+3" }, { "genotype": "DPB1*02:02+DPB1*356:01", "tce_groups": "3+3" }, { "genotype": "DPB1*02:02+DPB1*357:01N", "tce_groups": "0+3" }, { "genotype": "DPB1*02:02+DPB1*360:01", "tce_groups": "3+3" }, { "genotype": "DPB1*02:02+DPB1*372:01", "tce_groups": "3+3" }, { "genotype": "DPB1*02:02+DPB1*376:01", "tce_groups": "3+3" }, { "genotype": "DPB1*02:02+DPB1*396:01", "tce_groups": "3+3" }, { "genotype": "DPB1*02:02+DPB1*399:01", "tce_groups": "3+3" } ] } } ] }
The output of a PERMISSIVE probability of 2 is not mathematically valid.
This may be due to the presence of a NULL allele in the genotype and the results somehow adding the probability for the 0+3 and 3+3 cases together but thats just speculation.
Removing the null DBP1 alleles does lead to the service returning 1.0