R/AF2threshold.R

In parsec: Partial Orders in Socio-Economics

AF2threshold <-
function(mpi, prof, zeta = NULL) {
#     if (class(mpi)!="ophi")
      if (!is(mpi, "ophi"))
        stop("mpi is not an object of class 'ophi'")
#     if (class(prof)!="wprof")
      if (!is(prof, "wprof"))
        stop("prof is not an object of class 'wprof'")
    if (!is.null(zeta)) {
        # if (class(zeta)!="incidence")
        if (!is(zeta, "incidence"))
          stop("zeta is not an object of class 'incidence'")
    } else zeta <- getzeta(prof)
    ds <- apply(prof$profiles, 1, mpi$rho_k)
    gen.downset(zeta, ds)
}