معرفی شرکت ها

تبلیغات ما

مشتریان به طور فزاینده ای آنلاین هستند. تبلیغات می تواند به آنها کمک کند تا کسب و کار شما را پیدا کنند.

مشاهده بیشتر

تبلیغات ما

مشتریان به طور فزاینده ای آنلاین هستند. تبلیغات می تواند به آنها کمک کند تا کسب و کار شما را پیدا کنند.

مشاهده بیشتر

تبلیغات ما

مشتریان به طور فزاینده ای آنلاین هستند. تبلیغات می تواند به آنها کمک کند تا کسب و کار شما را پیدا کنند.

مشاهده بیشتر

تبلیغات ما

مشتریان به طور فزاینده ای آنلاین هستند. تبلیغات می تواند به آنها کمک کند تا کسب و کار شما را پیدا کنند.

مشاهده بیشتر

تبلیغات ما

مشتریان به طور فزاینده ای آنلاین هستند. تبلیغات می تواند به آنها کمک کند تا کسب و کار شما را پیدا کنند.

مشاهده بیشتر

توضیحات

Discrete mixture model based ploidy inference tool

ویژگی	مقدار
سیستم عامل	OS Independent
نام فایل	PloidPy-1.2.1
نام	PloidPy
نسخه کتابخانه	1.2.1
نگهدارنده	[]
ایمیل نگهدارنده	[]
نویسنده	Oluwatosin Olayinka
ایمیل نویسنده	oaolayin@live.unc.edu
آدرس صفحه اصلی	https://github.com/floutt/PloidPy
آدرس اینترنتی	https://pypi.org/project/PloidPy/
مجوز	-

# PloidPy ## Introduction PloidPy is a program written in Python designed to infer ploidy from next-generation reads aligned to a haploid reference genome. The program makes use of a stastical model representing the distribution of specific nucleotide counts and selects the most probable ploidy on the basis of a minimum AIC. ## Installation Installation of PloidPy is relatively simple and can be done easily using pip. The most recent stable version can be installed using the following command: ``` pip install PloidPy ``` And the developmental one can be installed directly from the git repository: ``` pip install git+git://github.com/floutt/PloidPy ``` ### Dependencies In order to run PloidPy, the following dependencies are required: - `Python 3.6+` - `NumPy` - `SciPy` - `Statsmodels` - `matplotlib` - `seaborn` ## Getting Started Once installed, PloidPy can be run on an indexed BAM file. For this tutorial we will download some BAM files created using simulated data from the [ploidyNGS](https://github.com/diriano/ploidyNGS) repository. ``` wget https://github.com/diriano/ploidyNGS/raw/master/test_data/simulatedDiploidGenome/Ploidy2.bowtie2.sorted.bam wget https://github.com/diriano/ploidyNGS/raw/master/test_data/simulatedTriploidGenome/Ploidy3.bowtie2.sorted.bam wget https://github.com/diriano/ploidyNGS/raw/master/test_data/simulatedTetraploidGenome/Ploidy4.bowtie2.sorted.bam ``` After downloading these files, we have to index all of these files using samtools. ``` samtools index Ploidy2.bowtie2.sorted.bam samtools index Ploidy3.bowtie2.sorted.bam samtools index Ploidy4.bowtie2.sorted.bam ``` From here we can extract the minor allele count (MAC) and total read coverage from each of the BAM files using `PloidPy process_bam`. In the following example the Phred-like mapping threshold is set to 15 ``` PloidPy process_bam --bam Ploidy2.bowtie2.sorted.bam --out diploid.count --quality 15 PloidPy process_bam --bam Ploidy3.bowtie2.sorted.bam --out triploid.count --quality 15 PloidPy process_bam --bam Ploidy4.bowtie2.sorted.bam --out tetraploid.count --quality 15 ``` These commands will produce two different files per command - a gzip-compressed archive containing the count data (`*count.gz`) and a file containing important metadata needed for ploidy evaluation (`*count.info`). After this has been done, one of two steps can be taken. Either inferring ploidy using the unfiltered data with an error component incorporated into the model or filtering the data and running the model with no error component. One can filter out the data from each individual using the following commands: ``` PloidPy filter --count_file diploid.count.gz --out diploid.count.filtered PloidPy filter --count_file triploid.count.gz --out triploid.count.filtered PloidPy filter --count_file tetraploid.count.gz --out tetraploid.count.filtered ``` This will create gzip-compressed files with count data filtering out data resulting from sequencing error and, with this, we can evaluate the ploidy of each individual using the `PloidPy assess` subcommand. PloidPy automatically detects whether or not the count file is filtered (by the presence or absence of the `*info` file) and adjusts the model accordingly. In this example, ploidies of 2n to 8n will be evaluated, although this can be done with any set of ploidies greater than or equal to 2n. This can be done with the following commands: ``` PloidPy assess --count_file diploid.count.filtered.gz --out diploid.filtered.tsv --ploidies 2 3 4 5 6 7 8 PloidPy assess --count_file triploid.count.filtered.gz --out triploid.filtered.tsv --ploidies 2 3 4 5 6 7 8 PloidPy assess --count_file tetraploid.count.filtered.gz --out tetraploid.filtered.tsv --ploidies 2 3 4 5 6 7 8 ``` The evaluation on unfiltered can similarly be done with the following commands: ``` PloidPy assess --count_file diploid.count.gz --out diploid.tsv --ploidies 2 3 4 5 6 7 8 PloidPy assess --count_file triploid.count.gz --out triploid.tsv --ploidies 2 3 4 5 6 7 8 PloidPy assess --count_file tetraploid.count.gz --out tetraploid.tsv --ploidies 2 3 4 5 6 7 8 ``` As you can see, in both cases all of the predictions were correct! Additional information can be found in the \*tsv files. Each column represents the following: |Column |Meaning| |-------------------|-------| |Ploidy |Ploidy model| |Log_Likelihood |Log likelihood of ploidy model| |AIC |AIC value of ploidy model| |Het_Weights |Weight parameter of each heterozygous state component| |Uniform Weight |Weight parameter of the uniform component of the model| |Binomial_Err_Weight|Weight parameter of the binomial error component (only present in unfiltered data)| Additionally, the filtered data can be used to produce helpful figures to visualize the joint distribution of TRC and MAC values. This can be done with the following commands: ``` PloidPy jdist --count_file diploid.count.filtered --out diploid PloidPy jdist --count_file triploid.count.filtered --out triploid PloidPy jdist --count_file tetraploid.count.filtered --out tetraploid ``` A visual overview of this process can be found [here](https://github.com/floutt/PloidPy/blob/master/figures/PloidPy_visual_guide.pdf) as well as some example [joint](https://github.com/floutt/PloidPy/blob/master/figures/diploid_joint_dist.pdf) [distribution](https://github.com/floutt/PloidPy/blob/master/figures/triploid_joint_dist.pdf) [figures](https://github.com/floutt/PloidPy/blob/master/figures/tetraploid_joint_dist.pdf).

زبان مورد نیاز

مقدار	نام
>=3.6	Python

نحوه نصب

نصب پکیج whl PloidPy-1.2.1:

pip install PloidPy-1.2.1.whl

نصب پکیج tar.gz PloidPy-1.2.1:

pip install PloidPy-1.2.1.tar.gz