Create a Phylogenetic Tree Object
genPhytree.Rd
Create a phylo tree object using a manipulated DNA string set object.
Examples
# create and plot a phylo tree
specdf_Anth <- querySpecData("Antheraea polyphemus")[1:10,]
DNABin_Anth <- genDNABin(specdf_Anth)
DNAStringset_Anth <- genDNAStringSet(DNABin_Anth)
#> Warning: cannot remove file 'new_DNAString_file.fasta', reason 'Permission denied'
DNAStringSet_Anth_manipulated <- ManipStringSet(DNAStringset_Anth)
#> ========================================================================================================================================================================================================
#>
#> Time difference of 0.05 secs
#> Determining distance matrix based on shared 9-mers:
#> ================================================================================
#>
#> Time difference of 0 secs
#>
#> Clustering into groups by similarity:
#> ================================================================================
#>
#> Time difference of 0 secs
#>
#> Aligning Sequences:
#> ================================================================================
#>
#> Time difference of 0.05 secs
#>
#> Iteration 1 of 2:
#>
#> Determining distance matrix based on alignment:
#> ================================================================================
#>
#> Time difference of 0 secs
#>
#> Reclustering into groups by similarity:
#> ================================================================================
#>
#> Time difference of 0 secs
#>
#> Realigning Sequences:
#> ================================================================================
#>
#> Time difference of 0.05 secs
#>
#> Alignment converged - skipping remaining iteration.
#>
Phytree_Anth <- genPhytree(DNAStringSet_Anth_manipulated)
Phytree_Anth$edge.length <- Phytree_Anth$edge.length + (mean(Phytree_Anth$edge.length) * 0.1)
plot(Phytree_Anth)