ethseq.Analysis {EthSEQ}

Ancestry analysis from whole-exome and targeted sequencing data

Description

This function performs ancestry analysis of a set of samples ad reports the results.

Usage

ethseq.Analysis(
  target.vcf = NA,
  target.gds = NA,
  bam.list = NA,
  out.dir = tempdir(),
  model.gds = NA,
  model.available = NA,
  model.assembly = "hg38",
  model.pop = "All",
  model.folder = tempdir(),
  run.genotype = FALSE,
  aseq.path = tempdir(),
  mbq = 20,
  mrq = 20,
  mdc = 10,
  cores = 1,
  verbose = TRUE,
  composite.model.call.rate = 1,
  refinement.analysis = NA,
  space = "2D",
  bam.chr.encoding = FALSE
)

Arguments

target.vcf

Path to the sample's genotypes in VCF format

target.gds

Path to the sample's genotypes in GDS format

bam.list

Path to a file containing a list of BAM files paths

out.dir

Path to the folder where the output of the analysis is saved

model.gds

Path to a GDS file specifying the reference model

model.available

String specifying the pre-computed reference model to use

model.assembly

String value indicating the assembly version to download for the pre-build models

model.pop

String value indicating the population to download for the pre-build models

model.folder

Path to the folder where reference models are already present or downloaded when needed

run.genotype

Logical values indicating whether the ASEQ genotype should be run

aseq.path

Path to the folder where ASEQ binary is available or is downloaded when needed

mbq

Minmum base quality used in the pileup by ASEQ

mrq

Minimum read quality used in the piluep by ASEQ

mdc

Minimum read count acceptable for genotype inference by ASEQ

cores

Number of parallel cores used for the analysis

verbose

Print detailed information

composite.model.call.rate

SNP call rate used to run Principal Component Analysis (PCA)

refinement.analysis

Matrix specifying a tree of ancestry sets

space

Dimensions of PCA space used to infer ancestry (2D or 3D)

bam.chr.encoding

Logical value indicating whether input BAM files have chromosomes encoded with "chr" prefix

Value

Logical value indicating the success of the analysis