R/get_eacf.R

In spmodel: Spatial Statistical Modeling and Prediction

get_eacf <- function(residual_vector2, dist_vector, bins, cutoff, formula) {

  # compute semivariogram classes
  dist_classes <- cut(dist_vector, breaks = seq(0, cutoff, length.out = bins + 1))

  # compute squared differences within each class
  acov <- tapply(residual_vector2, dist_classes, function(x) mean(x))

  # compute pairs within each class
  np <- tapply(residual_vector2, dist_classes, length)

  # set as zero if necessary
  np <- ifelse(is.na(np), 0, np)

  # compute average distance within each class
  dist <- tapply(dist_vector, dist_classes, mean)

  # return output
  eacf_out <- tibble::tibble(
    bins = factor(levels(dist_classes), levels = levels(dist_classes)),
    dist = as.numeric(dist),
    acov = as.numeric(acov),
    np = as.numeric(np)
  )

  # set row names to NULL
  # row.names(eacf_out) <- NULL

  eacf_out
}


get_eacf_cloud <- function(residual_vector2, dist_vector) {

  eacf_out <- tibble::tibble(dist = dist_vector, acov = residual_vector2)

  # set row names to NULL
  # row.names(eacf_out) <- NULL

  eacf_out
}