R/tfl2spc.R

In zipfR: Statistical Models for Word Frequency Distributions

tfl2spc <- function (tfl)
{
  if (! inherits(tfl, "tfl")) stop("argument must be object of class 'tfl'")
  if (attr(tfl, "incomplete")) stop("incomplete type frequency lists are not supported")

  f <- tfl$f
  if (!is.integer(f)) {
    if (! all(tfl$f == floor(tfl$f))) stop("type frequencies in 'tfl' must be integer values")
    f <- as.integer(f)
  }

  idx <- tfl$f > 0
  f <- if (!all(idx)) tfl$f[idx] else tfl$f
  x <- rle(sort(f))
  spc(Vm=x$lengths, m=x$values)
}