R/rad.R

In piklprado/sads: Maximum Likelihood Models for Species Abundance Distributions

rad <- function(x){
    if(is(x, "fitsad"))
        y <- new("rad", data.frame(rank=1:length(x@data$x), abund=sort(x@data$x, decreasing=T)))
    else if(is(x,"fitrad"))
        y <- x@rad.tab
    else if(is(x,"numeric")){
        z <- x[x>0]
        y <- new("rad", data.frame(rank=1:length(z), abund=sort(z, decreasing = TRUE)))
    }
    else if(is(x, "histogram")){
        ab <- rep(x$mids, x$counts)
        y <- new("rad", data.frame(rank=1:length(ab), abund=sort(ab, decreasing = TRUE)))
    }
    else if(is(x, "fitsadC"))
        y <- rad(x@hist)
    return(y)
}