Description Usage Arguments Value Author(s) Examples
View source: R/screen_bam_multi.R
An efficient function to identify sites of interestes from a list of longitudinal sequence files with parallel computing######
1 | screen_bam_multi(file_list, ref = NULL, region, priorPA = 0.001, core = 1)
|
file_list |
A list of bam files. |
ref |
A referece sequence. |
region |
An optional region from the BAM file (as a |
priorPA |
Prior assocation to possible models. |
core |
Number of cores for parallel computation. |
PPA association for each site.
Fei Xiang (xf3087@gmail.com)
1 2 3 4 5 6 7 8 9 10 11 12 13 | p4305.d3=system.file("extdata/4305-d3.bam",package = "popseq")
p4305.d4=system.file("extdata/4305-d4.bam",package = "popseq")
p4305.d5=system.file("extdata/4305-d5.bam",package = "popseq")
ref=system.file("extdata/A_England_195_2009.fa",package = "popseq")
pigbam=list(p4305.d3,p4305.d4,p4305.d5)
res1=screen_bam_multi(file_list=pigbam,ref=ref,region= GRanges("GQ166656_1",IRanges(1,2280)),priorPA=.001,core=8)
res2=screen_bam_multi(file_list=pigbam,ref=ref,region = GRanges("GQ166655_2",IRanges(1,2274)),priorPA=.001,core=8)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.