Daniel:Notebook/Haplotyping/SISSORPipeline/2014-12-5
Segmentation Analysis Pipeline (Started 12/2/2014)[edit]
On 12/3, I identified that the cutoff of 20 reads used in the pipeline is actually pretty arbitrary, and is based off the expected number of reads per bin based on the amount of sequencing we do. As a result, I analyzed the read count distribution from Chamber 1 of PGP1-22 (arbitrarily chosen). I loaded the readcount data into MATLAB, and attempted to fit the distribution using the fitdist command. The code is here The two most closely fitting distributions are the Negative Binomial and the Gamma, which makes sense based on the underlying assumptions of the probability distributions. For reference, I have included the CDFs of the two distributions below.
- Negative Binomial Fit.jpg
Negative Binomial distribution fit
- Gamma Fit Chamber 1.jpg
Gamma distribution fit
Conclusion[edit]
After talking with Dr. Zhang, we have decided to keep the 20 for now, and that I may look in to the cutoff later, if it becomes a bigger problem.