rankSizeSlope:

In JusteRaimbault/r-spatialmorphology: What the package does (short line)

Usage

rankSizeSlope(m = NULL)

Arguments

m

Examples

##---- Should be DIRECTLY executable !! ----
##-- ==>  Define data, use random,
##--	or do  help(data=index)  for the standard data sets.

## The function is currently defined as
function (m = NULL) 
{
    if (is.null(m)) {
        return(list(rankSizeAlpha = NA, rankSizeRSquared = NA))
    }
    r_pop = raster(m)
    size = cellStats(r_pop, function(x, ...) {
        na.omit(log(x))
    })
    size = size[size > 0]
    size = sort(size, decreasing = TRUE)
    rank = log(1:length(size))
    if (length(size) > 0) {
        reg = lm(size ~ rank, data.frame(rank, size))
        return(list(rankSizeAlpha = reg$coefficients[2], rankSizeRSquared = summary(reg)$r.squared))
    }
    else {
        return(list(rankSizeAlpha = NA, rankSizeRSquared = NA))
    }
  }

JusteRaimbault/r-spatialmorphology documentation built on May 22, 2019, 5:07 p.m.