How to select significant cell type eQTL
zhangpicb opened this issue · 6 comments
Dear @npklein
Thanks for your great tool!
I want to use Decon-qtl to find cell-type eQTL with cibersort results,But I find no significant p-value after B-H correction. Most p-value in results are 1,like README
systemsgenetics/Decon2/README.md
Lines 151 to 153 in a7fc9e2
How to select significant cell type eQTL?
Thanks in advanced!
Dear @zhangpicb
How many samples did you use? The reason that they are all not significant in the README is because it is a toy example with a very small data set. To identify interactions at least a couple hundred samples are necesarry.
I find the method in Decon-qtl is similar to GTEx cell type interaction eQTL .Becuase our lab own GTEx v8 data,without BLUPRINT data.So I replicate https://www.science.org/doi/10.1126/science.aaz8528 Fig5c , But the pvalue is 1.
Our method is quite a bit different from GTEx method. In our paper we show that we trade higher specificity for lower sensitivity, so that could be why it does no replicate.
Alternatively, can also be due to strict pvalue correction? For how many SNP-Gene pairs did you run it? We run it only for top SNPs, so if you ran for all snp gene combinations that could also be the issue
Dear @npklein
Thanks for your quick reply!
I use all snp-gene combinations,around 1500 snp-gene pairs per gene.The SNPs are from gene tss±1Mb.
How to select top SNPs for each gene?
Thanks in advanced!
Okay that is not too many tests, in that case I think more samples are necessary