piquelab / QuASAR

Quantitative Allele Specific Analysis of Reads. Joint genotyping & ASE inference for RNA-seq data
MIT License
25 stars 10 forks source link

Any way to use different indeviduals as replicates #15

Open IanSudbery opened 5 years ago

IanSudbery commented 5 years ago

Issue #14 talks about the situation where there are no replicates. I have a similar situation - for each individual there is only one replicate, but I have many individuals. Not all individuals are going to be hetrozygous at the same positions, but where two individuals are hetrozygous at the same position, presumably you'd be more confident of the significance of any imbalance at that location if it were present in both individuals, especially if we were talking about 20 or 30 individuals rather than 2.

Is there anyway the QuASAR algorithm can be altered to allow this sort of replication?