Pi_N.par: Calculate Pi_N for loci with a high impact or moderate impact...

View source: R/picalc.R

Pi_N.parR Documentation

Calculate Pi_N for loci with a high impact or moderate impact mutation using parallel computation

Description

Calculate Pi_N for loci with a high impact or moderate impact mutation using parallel computation

Usage

Pi_N.par(TH, exon.sel, exon.mut, nSNP = NA, datatype = "full", cores = 15)

Arguments

TH

Site-wise theta calculated using ANGSD

exon.sel

Selected gene models

exon.mut

Coordinates of the non-synonymous positions.

nSNP

If given, divides the total Pi with the given number of SNPs.

datatype

"full" assumes that all positions are called Pi is a mean of observed values, "SNP" assumes only variant positions and the Pi is divided by the total number of non-synonymous positions in the selected gene models.

cores

Number of cores for parallel computation

Value

a list with Pi, chromosome-wise cumulative sums and numbers of SNPs.


jsalojar/PiNSiR documentation built on Nov. 1, 2023, 1:47 p.m.