201019-tadaA-wes-interaction-res

separate estimation

rr

A data.frame: 15 × 4
log_RR	lower_bound	upper_bound	annota
<dbl>	<dbl>	<dbl>	<chr>
2.5987195	2.1512099	3.046229	annovar_lof & ccr_gt90
1.7289160	0.5694150	2.888416	ago_adult_brain.BA4.hg19 & ago_adult_brain.Cingulate.gyrus.hg19
2.0650330	0.2667945	3.863271	ago_adult_brain.BA4.hg19 & ccr_gt90
2.6081985	1.5657429	3.650654	elavl_Adult_brain.BA9.hg19 & ccr_gt90
-1.0000000	-12.1245250	10.124525	ago_adult_brain.Cingulate.gyrus.hg19 & ccr_gt90
1.3859830	-1.9550871	4.727053	elavl_Adult_brain.all_human_samples.hg19 & ccr_gt90
-1.0000000	-9.7084343	7.708434	elavl_Adult_brain.BA9_Alzheimer.hg19 & ccr_gt90
2.3654829	1.2929190	3.438047	RBP.ccr90
0.5516327	-8.2777934	9.381059	ribosnitch.ccr90
1.8249369	-2.2379312	5.887805	CLIPdb.ccr90
2.4386057	1.2020223	3.675189	alt3.spliced.exon.3SS.hg19 in top5% with ccr_gt90
3.3635423	2.6357400	4.091344	alt3.spliced.exon.5SS.hg19 in top5% with ccr_gt90
2.4734789	1.0323978	3.914560	alt5.spliced.exon.3SS.hg19 in top5% with ccr_gt90
1.8634173	-0.4201538	4.146988	alt5.spliced.exon.5SS.hg19 in top5% with ccr_gt90

joint estimation

test1

1 2	z = readRDS(paste0("/storage11_7T/fuy/TADA-A/annotation/results/",Sys.Date(),"_wes_interaction_sig_joint_RR.rds")) z

A data.frame: 4 × 4
logRR	lower_bound	upper_bound	annota
<dbl>	<dbl>	<dbl>	<chr>
1.2712498	0.7477706	1.794729	annovar_lof & ccr_gt90
1.8316680	1.0435827	2.619753	elavl_Adult_brain.BA9.hg19 & ccr_gt90
0.9908704	0.7841269	1.197614	pathogenic_rare_missense_MVP_gt75
1.2161730	0.9857013	1.446645	splidex_lower10pct

test2

z2

A data.frame: 12 × 4
logRR	lower_bound	upper_bound	annota
<dbl>	<dbl>	<dbl>	<chr>
-0.46660527	NaN	NaN	RBP.ccr90
0.09961991	-1.3043522	1.5035920	alt3.spliced.exon.3SS.hg19.top5%.ccr90
1.64550461	0.5611467	2.7298625	alt3.spliced.exon.5SS.hg19.top5%_altA.ccr90
0.01852650	-3.0703185	3.1073715	alt5.spliced.exon.3SS.hg19.top5%_altA.ccr90
1.10205708	0.5595139	1.6446002	annovar_lof & ccr_gt90
-0.93862380	NaN	NaN	ago_adult_brain.BA4.hg19 & ccr_gt90
1.26739765	0.3275906	2.2072047	elavl_Adult_brain.BA9.hg19 & ccr_gt90
-1.81554186	NaN	NaN	pLI_nlt_995_gnomad.v2.1.1.all_lofs_snv
0.70539583	0.4553003	0.9554914	MVP_rare_missense_pathogen_gt75
0.08359544	-0.2266574	0.3938483	primateAI_exome_mutation_pathogen_rank_gt_80
1.22933849	1.0103475	1.4483294	spidex_lower10pct
0.96205800	0.6231595	1.3009565	MPC_gt2

1	readRDS(paste0("/storage11_7T/fuy/TADA-A/annotation/results/",Sys.Date(),"_wes_interaction_test2_g_BF2.rds"))

A data.table: 23 × 9
genename	logBF_coding	logBF_additional	logBF_all	BF_coding	BF_additional	BF_all	FDR_prior	FDR_all
<chr>	<dbl>	<dbl>	<dbl>	<dbl>	<dbl>	<dbl>	<dbl>	<dbl>
SCN2A	28.240698	0	28.240698	1.839837e+12	1	1.839837e+12	0.94	8.515300e-12
CHD8	14.733305	0	14.733305	2.503762e+06	1	2.503762e+06	0.94	3.128610e-06
SLC6A1	14.474444	0	14.474444	1.932731e+06	1	1.932731e+06	0.94	4.787710e-06
FOXP1	10.314673	0	10.314673	3.017209e+04	1	3.017209e+04	0.94	1.333343e-04
GRIN2B	8.942270	0	8.942270	7.648539e+03	1	7.648539e+03	0.94	5.154943e-04
SYNGAP1	8.117830	0	8.117830	3.353736e+03	1	3.353736e+03	0.94	1.204526e-03
WBSCR17	6.639490	0	6.639490	7.647052e+02	1	7.647052e+02	0.94	3.900437e-03
DEAF1	6.352399	0	6.352399	5.738675e+02	1	5.738675e+02	0.94	6.734714e-03
LRRC43	5.987375	0	5.987375	3.983675e+02	1	3.983675e+02	0.94	1.019075e-02
KCNQ3	5.545837	0	5.545837	2.561688e+02	1	2.561688e+02	0.94	1.493497e-02
PRKAR1B	5.538171	0	5.538171	2.542126e+02	1	2.542126e+02	0.94	1.885457e-02
TBR1	5.473951	0	5.473951	2.384004e+02	1	2.384004e+02	0.94	2.242199e-02
KDM5B	5.107110	0	5.107110	1.651923e+02	1	1.651923e+02	0.94	2.736058e-02
CHD2	4.772796	0	4.772796	1.182495e+02	1	1.182495e+02	0.94	3.376258e-02
SETD5	4.707779	0	4.707779	1.108058e+02	1	1.108058e+02	0.94	3.977002e-02
IQSEC2	4.465698	0	4.465698	8.698168e+01	1	8.698168e+01	0.94	4.682343e-02
SPOCK3	4.331342	0	4.331342	7.604624e+01	1	7.604624e+01	0.94	5.411752e-02
TFAP4	4.257804	0	4.257804	7.065467e+01	1	7.065467e+01	0.94	6.119390e-02
PAX8	4.236171	0	4.236171	6.914262e+01	1	6.914262e+01	0.94	6.769571e-02
RHOBTB2	4.223501	0	4.223501	6.827210e+01	1	6.827210e+01	0.94	7.364312e-02
LRP10	4.065766	0	4.065766	5.830955e+01	1	5.830955e+01	0.94	8.022106e-02
DPYSL2	3.947032	0	3.947032	5.178145e+01	1	5.178145e+01	0.94	8.713270e-02
STXBP1	3.833372	0	3.833372	4.621810e+01	1	4.621810e+01	0.94	9.435122e-02