R/levels.incidence.R

Defines functions levels.incidence

Documented in levels.incidence

levels.incidence <-
function(x) {
  n <- nrow(x)
  lev <- rep(NA, n)
  names(lev) <- rownames(x)
  i <- 1
  sub <- (1:n)[is.na(lev)]
  while(length(sub)>0) {
    tmp_z <- x[sub, sub]
    class(tmp_z) <- class(x)
    quali <- which(depths(tmp_z)==1)
    lev[sub[quali]] <- i
    i <- i + 1
    sub <- (1:n)[is.na(lev)]
  }
  lev
}

Try the parsec package in your browser

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

parsec documentation built on Aug. 19, 2023, 5:07 p.m.