HEMDAG (version 2.7.4)

obozinski.heuristic.methods: Obozinski heuristic methods

Description

Implementation of the Obozinski's heuristic methods Max, And, Or (Obozinski et al., Genome Biology, 2008, 10.1186/gb-2008-9-s1-s6).

Usage

obozinski.max(S, g, root = "00")
obozinski.and(S, g, root = "00")
obozinski.or(S, g, root = "00")

Arguments

S

a named flat scores matrix with examples on rows and classes on columns.

g

a graph of class graphNEL. It represents the hierarchy of the classes.

root

name of the class that it is the top-level of the hierarchy (def:00).

Value

A matrix with the scores of the classes corrected according to the chosen Obozinski's heuristic algorithm.

Details

Obozinski's heuristic methods:

1. Max: reports the largest logistic regression (LR) value of self and all descendants: \(p_i = max_{j \in descendants(i)} \hat{p_j}\);

2. And: reports the product of LR values of all ancestors and self. This is equivalent to computing the probability that all ancestral terms are "on" assuming that, conditional on the data, all predictions are independent: \(p_i = \prod_{j \in ancestors(i)} \hat{p_j}\);

3. Or: computes the probability that at least one of the descendant terms is "on" assuming again that, conditional on the data, all predictions are independent: \(1 - p_i = \prod_{j \in descendants(i)} (1 - \hat{p_j})\);

Examples

# NOT RUN {
data(graph);
data(scores);
root  <- root.node(g);
S.max <- obozinski.max(S,g,root);
S.and <- obozinski.and(S,g,root);
S.or  <- obozinski.or(S,g,root);
# }