Learn R Programming

parsec (version 1.1.1)

AF2threshold: Poset threshold equivalent to the AF counting approach

Description

The function computes the profile poset threshold equivalent to the AF counting approach described in argument mpi.

Usage

AF2threshold(mpi, prof, zeta = NULL)

Arguments

mpi
an object of class ophi, see AF for details.
prof
an object of class wprof.
zeta
an object of class incidence.

See Also

AF

Examples

Run this code
vl <- c(2, 3, 3, 2)
prof <- var2prof(varlen = vl, labtype = "progressive")

res <- AF(prof, z = c(1, 2, 1, 1), k = 1)

thr <- AF2threshold(res, prof)

plot(prof, col = 1 + thr, lwd = 1 + res$c,
     main = "Comparison between OPHI and parsec",
     sub = "bold: deprived profiles identified by OPHI, red: parsec threshold")

eval <- evaluation(prof, thr, maxint = 10^4, nit = 10^7)

ord <- order(eval$idn_f, res$c)
plot(eval$idn_f[ord], col = "red", lwd=2, type = "l", xlab="",
     ylab = "", axes = FALSE, frame.plot = TRUE,
     main = "Comparison between OPHI and parsec",
     sub = "red: identification function, black: OPHI deprived profiles")
points(res$c[ord], type="l", lwd=2)
axis(2)

Run the code above in your browser using DataLab