Read.Directory.PhylipFiles: Read.Directory.PhylipFiles.R Function to read the alignments...

In radamsRHA/Misc.Statistics.Genetics.Models:

Description

The function returns a list of alignments obtained from a directory that contains phylip formatted alignments

Usage

Read.Directory.PhylipFiles(phy.dir, popID)

Arguments

phy.dir	directory containing phylip formats
popID	vector containing all of the samples ID that are extracted from the alignments

Examples

library(AlnStats)
library(stringr)

EXAMPLE
wanted <- c(">CV0031_0", ">CV0031_1", ">CV0310_0",">CV0310_1")
aln.data <- read.alleles2(alleles.file = '~/Desktop/EXAMPLE.alleles', popID = wanted)
S.data <- Calc.S(alleles.data = aln.data)
hist(S.data$theta.hat.vec)
hist(S.data$S.vec)

radamsRHA/Misc.Statistics.Genetics.Models documentation built on May 5, 2019, 6:56 p.m.