sharp (version 1.4.6)

PFER: Per Family Error Rate

Description

Computes the Per Family Error Rate upper-bound of a stability selection model using the methods proposed by Meinshausen and Bühlmann (2010) or Shah and Samworth (2013). In stability selection, the PFER corresponds to the expected number of stably selected features that are not relevant to the outcome (i.e. False Positives).

Usage

PFER(q, pi, N, K, PFER_method = "MB")

Value

The estimated upper-bound in PFER.

Arguments

q

average number of features selected by the underlying algorithm.

pi

threshold in selection proportions.

N

total number of features.

K

number of resampling iterations.

PFER_method

method used to compute the upper-bound of the expected number of False Positives (or Per Family Error Rate, PFER). If PFER_method="MB", the method proposed by Meinshausen and Bühlmann (2010) is used. If PFER_method="SS", the method proposed by Shah and Samworth (2013) under the assumption of unimodality is used.

References

stabilityselectionMBsharp

stabilityselectionSSsharp

See Also

Other stability metric functions: ConsensusScore(), FDP(), StabilityMetrics(), StabilityScore()

Examples

# Computing PFER for 10/50 selected features and threshold of 0.8
pfer_mb <- PFER(q = 10, pi = 0.8, N = 50, K = 100, PFER_method = "MB")
pfer_ss <- PFER(q = 10, pi = 0.8, N = 50, K = 100, PFER_method = "SS")