load_abbababa: Load a results file from Simon Martin's ABBABABAwindows.py...

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

Load a results file from Simon Martin's ABBABABAwindows.py script.

Description

This function was developed strictly following his GitHub description of the output file. The columns in this output file are: scaffold, start, end, sites, sitesUsed, ABBA, BABA, D, fd, fdM. You can choose to plot any of the statistics, and filter windows by the number of sites used.

Usage

load_abbababa(file, stat = "D", position = "midpoint", min_sites = 0)

Arguments

file	Input file name
stat	Statistic(s) to keep. This can be a vector of column names.
position	Use the start, end, or midpoint as the position?
min_sites	Minimum number of sitesUsed required for window inclusion.

Value

A tibble.

Examples

a <- system.file("extdata", "test.abbababa.csv.gz",
                 package = "gwplotting")
b <- load_abbababa( a, stat = 'D' )
b

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