sample_from_angsd_probs: Sample a genotype (0,1,2) from probabilities for all the snps...

Description Usage Arguments

View source: R/RcppExports.R

Description

Sample a genotype (0,1,2) from probabilities for all the snps and individuals in x.

Usage

1
sample_from_angsd_probs(x, samples, snps, thresh)

Arguments

x

a numeric vector that holds the posterior probs from an ANGSD geno.gz file (from doGenos 32)

samples

character vector of sample names

snps

character vector of SNP names

thresh

any snp in an individual must have one posterior greater than T in order to have a non-missing value simulated for it. For example, if T is 0.6 and the posteriors are 0.33, 0.33, 0.33, this locus will just get a missing (-1).


eriqande/genoscapeRtools documentation built on Dec. 27, 2021, 8:01 a.m.