Defines functions `dist.eval`

`dist.eval` <-
function(x, dist){
    x <- as.matrix(x)
    tots <- rowSums(x)
    if(any(tots==0)) {
        cat("Warning: the community matrix contains some sites with only zero abundances\n")
        cat("You may want to use functions removezerospecies or dist.zeroes from Biodiversity.R\n")
        tests <- c("manhattan", "euclidean", "canberra", "clark", "bray", "kulczynski", "jaccard", "gower", "altGower", "morisita", "horn", "mountford", "raup" , "binomial", 
            "chao", "cao", "mahalanobis", "hellinger")
        op <- options()
        if (any(dist==tests)) {
            dist1 <- vegdist(x, method=dist)
            if((dist %in% c("w", "-1", "c", "wb", "r", "I", "e", "t", "me", "j", "sor", "m", "-2", "co", "cc", "g", "-3", "l", "19", "hk", "rlb", "sim", "gl", "z")) == F) {
                stop("Provide acceptable method for betadiver")
            dist1 <- betadiver(x, method=dist)
        dist2 <- no.shared(x)
        list1 <- (dist2==0)
        list2 <- (dist2==1)
        max <- max(dist1[list1])
        min <- min(dist1[list2])
        if(min<max) {
            cat("Warning: min distance for sites with no shared species(",min,") < max dist for other sites(", max, ")\n")
            cat("Choose other distance measure or stepacross\n")

Try the BiodiversityR package in your browser

Any scripts or data that you put into this service are public.

BiodiversityR documentation built on Jan. 6, 2023, 5:18 p.m.