load_vcftools_stats: Load a VCFtools results file.

In nwvankuren/gwplotting: Utilities for plotting statistics from genome-wide analyses

Load a VCFtools results file.

Description

This function will read a (possibly gzipped) VCFtools results file from windowed Fst analysis (.windowed.weir.fst), nucleotide diversity analysis (.windowed.pi), or Tajima's D (.Tajima.D) and return a tibble containing just the scaffold, position, and statistic.

Usage

load_vcftools_stats(
  file,
  stat = "mean_fst",
  min_sites = 0,
  position = "midpoint"
)

Arguments

file	VCFtools results file. Can be gzipped.
stat	Name of the column to use as statistic. This will be mean_fst, weighted_fst, pi, tajimas_d, ...
min_sites	Minimum number of sites required to keep a window
position	Coordinate to use as position: start, end, or midpoint

Value

A three-column tibble containing scaffold, bin position, and statistic

Examples

a <- system.file("extdata", "test.vcftools_pi.txt.gz",
                 package = "gwplotting")
b <- load_vcftools_stats( a, stat = 'pi' )
b

nwvankuren/gwplotting documentation built on July 5, 2024, 4:24 p.m.