From: Roastgsa: a comparison of rotation-based scores for gene set enrichment analysis
Competitive score | |
---|---|
Meanrank | \(T_{meanrank} = \frac{1}{m_o}\sum \limits _{i\in S}\frac{(q+1)/2 - \rho _i}{q}\) |
Ksmax | \(T_{ksmax} = \hbox {I}(A>|a|)A + \hbox {I}(A\le |a|)a,\) |
\(A = \max \limits _{l \in S} ks(l|S), a = \min \limits _{l \in S} ks(l|S),\) | |
\(ks(l|S) = \frac{\sum \limits _{i\in S} |\gamma _i|^k I(\rho _i \le l)}{\sum \limits _{i\in S} |\gamma _i|^k} - \frac{\sum \limits _{i\not \in S} I(\rho _i \le l)}{q - m_o}\) | |
Ksmean | \(T_{ksmean} = \max \limits _{l \in S} ks(l|S) + \min \limits _{l \in S} ks(l|S),\) |
\(ks(l|S) = \frac{\sum \limits _{i\in S} |\eta _i|^k I(\rho _i \le l)}{\sum \limits _{i\in S} |\eta _i|^k} - \frac{\sum \limits _{i\not \in S} I(\rho _i \le l)}{q - m_o}\) | |
Notation | \([\delta _i] \equiv \hbox {modt-statistics, } i\in \Omega = [1,\ldots ,q]\) |
\([\rho _i] \equiv \hbox {rank for } [\delta _i] \hbox { in decreasing order}\) | |
\(S \equiv \hbox {Testing gene set, } S\subset \Omega , \,\,\, m_0 = |S|,\,\, C = \Omega \setminus S\) | |
\(\gamma _i = \delta _i - \bar{\delta },\,\,\, \eta _i = (2\rho _i+q+1)/2\) |