returns number of DNA sites used
rcpp_pairwiseDeletionDNA(dnavector, ncores = 1L, symmetric = 1L)
StringVector [mandatory]
number of cores [default: 1]
symmetric score matrix [default: 1]
list
## load example sequence data
data("woodmouse", package="ape")
w <- woodmouse |> dnabin2dnastring() |> as.character()
rcpp_pairwiseDeletionDNA(dnavector=w, ncores=1)
#>
Computing: [========================================] 100% (done)
#> $sitesUsed
#> No305 No304 No306 No0906S No0908S No0909S No0910S No0912S No0913S
#> No305 962 959 960 958 958 958 958 958 958
#> No304 959 962 960 959 959 959 959 959 959
#> No306 960 960 963 959 959 959 959 959 959
#> No0906S 958 959 959 961 961 961 961 961 961
#> No0908S 958 959 959 961 961 961 961 961 961
#> No0909S 958 959 959 961 961 961 961 961 961
#> No0910S 958 959 959 961 961 961 961 961 961
#> No0912S 958 959 959 961 961 961 961 961 961
#> No0913S 958 959 959 961 961 961 961 961 961
#> No1103S 958 959 959 961 961 961 961 961 961
#> No1007S 958 959 959 961 961 961 961 961 961
#> No1114S 914 913 913 915 915 915 915 915 915
#> No1202S 958 959 959 961 961 961 961 961 961
#> No1206S 958 959 959 961 961 961 961 961 961
#> No1208S 956 956 956 958 958 958 958 958 958
#> No1103S No1007S No1114S No1202S No1206S No1208S
#> No305 958 958 914 958 958 956
#> No304 959 959 913 959 959 956
#> No306 959 959 913 959 959 956
#> No0906S 961 961 915 961 961 958
#> No0908S 961 961 915 961 961 958
#> No0909S 961 961 915 961 961 958
#> No0910S 961 961 915 961 961 958
#> No0912S 961 961 915 961 961 958
#> No0913S 961 961 915 961 961 958
#> No1103S 961 961 915 961 961 958
#> No1007S 961 961 915 961 961 958
#> No1114S 915 915 915 915 915 915
#> No1202S 961 961 915 961 961 958
#> No1206S 961 961 915 961 961 958
#> No1208S 958 958 915 958 958 958
#>