2021-05-24-nonfs-rate-fs-gama
1234567891011121314151617181920212223242526272829303132library(data.table)all = fread("/storage11_7T/fuy/TADA-A/cell_WES/denovo_all.csv")[,c ...
2021-05-20-fix-pi-VS-estim-pi-effect-size-num-enrich
effect size
add a new annotation mixture of spidex_low3% and merged spliceai, without PTVs, and remove the original ones.
the difference of effect s ...
2021-05-14-tadaA-report
估参算法改进
将sapply全部改为Rcpp进行循环，整体的估参速度提升很多，如5组annotation联合估参时，改进前需要7h，而改进后只需要1.2h
多个参数估参时，optim method 由默认的Nelder-Mead改为BFGS，避免陷入局部最优
prior由原来固定的0.05， ...
2021-05-13-indel+snv-novel-CHD-GO
12library(data.table)fread("/storage11_7T/fuy/TADA-A/CHD/indel_rate/indel+snv_novel_GO.lst")[,-7]
A data.table: 13 × 14
MyListGene S ...
2021-05-13_md_indel-snv_ASD_riskG_GO_DisGeNET_17
123library(data.table)f = fread("/storage11_7T/fuy/TADA-A/cell_WES/DNM/indel_rate/0513/17_indel-snv_ASD_risk_gene_GO.lst")[,c(-2,-6)]f
A ...
2021-05-11-rr-all-annota-ASD-CHD
ASD1readRDS("/storage11_7T/fuy/TADA-A/cell_WES/DNM/results/0511-all-rr.rds")
A data.table: 23 × 4
separate_logRRupper_boundlower_bou ...
2021-05-13-indel-model-ASD
indel rate123456789101112131415library(data.table)options(digits=22)sp = 6430 ### sample sizedf = read.table("/storage11_7T/fuy/TADA-A/cell_WES/D ...
2021-05-10-tadaA-report
2021-05-10-tadaA-report
/*!
*
* Twitter Bootstrap
*
*/
/*!
* Bootstrap v3.3.7 (http://getbootstrap.com)
* Copyright 2011-2016 Twitter ...
2021-04-19-tadaA-report
2021-04-19-tadaA-report
/*!
*
* Twitter Bootstrap
*
*/
/*!
* Bootstrap v3.3.7 (http://getbootstrap.com)
* Copyright 2011-2016 Twitter ...
2021-04-11-optim-update-posterior-pi
2021-04-11-optim-update-posterior-pi
/*!
*
* Twitter Bootstrap
*
*/
/*!
* Bootstrap v3.3.7 (http://getbootstrap.com)
* Copyright 2011 ...