gnomad.v2.1.1.all_lofs

RR=2.2615612

[1.6800949, 2.8430276]

PTVs in paper

RR(average by genes) = 1.949428

MPC_gt2

RR=2.0834977

[1.8200768,1.5566559] 

1
log(22.15)  ### MisB [2,~]

3.09783749649114

1
log(4.18) ### MisA [1,2)

1.43031124653666

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
annota = c("gnomad.v2.1.1.all_lofs_snv (PTVs)",
           "MPC_gt2 (MisB)",
           "lofs_pLI_nlt_995",
           "lofs_pLI_05_995",
           "MPC_1-2 (MisA)")

RR_our = c(2.2615612,
           2.0834977,
           3.1078241,
           -1.0000000,
           -1.000000)

RR_paper = c(1.949428,
             3.097837,
             3.931826,
             1.94591, 
             1.43031)

df = data.frame(RR_our,annota,RR_paper)
df
A data.frame: 5 × 3
RR_ourannotaRR_paper
<dbl><fct><dbl>
2.261561gnomad.v2.1.1.all_lofs_snv (PTVs)1.949428
2.083498MPC_gt2 (MisB) 3.097837
3.107824lofs_pLI_nlt_995 3.931826
-1.000000lofs_pLI_05_995 1.945910
-1.000000MPC_1-2 (MisA) 1.430310