Performs tree inference under "RAxML"
for aligned fasta sequences in
a given folder (default is "2.Alignments"
).
Usage
tree.constraint(
taxonomy_folder = "1.CuratedSequences",
targetColumns = c("kingdom", "phylum", "class", "order", "family", "genus",
"species_names"),
Topology = "((ingroup), outgroup);",
outgroup = NULL
)
Arguments
- taxonomy_folder
Name of the folder where the 1.Taxonomy file is stored.
- targetColumns
Where to find
"RAxML"
or how to run it from the console? (string).- Topology
A string summarizing the desired topological constraint in newick format.
- outgroup
Optional and only required when the topology argument is not "((ingroup), outgroup);".
Examples
if (FALSE) {
sq.retrieve.direct(
clades = c("Felis", "Vulpes", "Phoca"),
species = "Manis_pentadactyla",
genes = c("ADORA3", "CYTB")
)
sq.curate(
filterTaxonomicCriteria = "Felis|Vulpes|Phoca|Manis",
kingdom = "animals", folder = "0.Sequences"
)
tree.constraint(
taxonomy_folder = "1.CuratedSequences",
targetColumns = c("kingdom", "phylum", "class", "order", "family",
"genus", "species_names"),
Topology = "((ingroup), outgroup);",
outgroup = "Manis_pentadactyla"
)
tree.constraint(
taxonomy_folder = "1.CuratedSequences",
targetColumns = c("kingdom", "phylum", "class", "order", "family",
"genus", "species_names"),
Topology = "((Felis), (Phoca));"
)
}