simulation
sample size = 20000
num of genes = 5000
10 annotations, with logRR = c(4,3,2,1.5,1,-8,-6,-5,-3,-1)
prior=0.05
100 rounds of simulation
statistics of estimated parameters
select annotations 1:5
1 | summary(rr.dff) |
ptv>0.995 spidex mpc2 RADAR_RBP_top5%
Min. :3.786 Min. :2.808 Min. :1.764 Min. :1.316
1st Qu.:3.912 1st Qu.:2.899 1st Qu.:1.949 1st Qu.:1.381
Median :3.941 Median :2.956 Median :1.993 Median :1.454
Mean :3.972 Mean :2.954 Mean :1.982 Mean :1.455
3rd Qu.:4.048 3rd Qu.:3.023 3rd Qu.:2.040 3rd Qu.:1.524
Max. :4.204 Max. :3.069 Max. :2.109 Max. :1.628
RBP-VarDB prior
Min. :0.2228 Min. :0.04103
1st Qu.:0.6487 1st Qu.:0.04701
Median :0.7710 Median :0.04989
Mean :0.7356 Mean :0.05097
3rd Qu.:0.8348 3rd Qu.:0.05501
Max. :1.0120 Max. :0.06288 FDR
1 | readRDS("/storage11_7T/fuy/TADA-A/cell_WES/DNM/simulation/10annota_sim/org.fdr.check.res.rds") |
| cutoff.lst | fdr.average |
|---|---|
| <dbl> | <dbl> |
| 0.02 | 0.01922601 |
| 0.05 | 0.05140520 |
| 0.08 | 0.08159085 |
| 0.10 | 0.10391246 |
| 0.15 | 0.16152949 |
| 0.20 | 0.21270958 |
| 0.30 | 0.31172867 |
| 0.50 | 0.51087254 |