plot.Knet {SANTA} | R Documentation |
Knet
function
Plot the observed Knet curve against the quantiles of the permuted Knet curves and the observed AUK against the permuted AUKs.
## S3 method for class 'Knet' plot(x, sequential=FALSE, ...)
x |
Results from the Knet function. |
sequential |
Logical, if |
... |
Additional arguments to be passed to |
If the high-weight vertices are clustered, then the observed Knet curve and AUK will be high relative to the permuted Knet curves and AUKs. The greater the strength of clustering, the greater the difference between the observed and permuted statistics. If the strength of clustering is low, then the observed and permuted curves and AUKs will likely overlap.
The first plot displays the the observed curve in red and the quantiles of the permuted curves in yellow. The quantile boundaries are displayed as grey lines. These boundaries are specified in the Knet function. The second plot displays the observed AUK as a red line and the distribution of permuted AUKs in grey.
Alex J. Cornish a.cornish12@imperial.ac.uk and Florian Markowetz
Cornish, A.J. and Markowetz, F. (2014) SANTA: Quantifying the Functional Content of Molecular Networks.. PLOS Computational Biology. 10:9, e1003808.
# plot results with hit clustering g.clustered <- barabasi.game(100, directed=FALSE) g.clustered <- SpreadHits(g.clustered, h=10, lambda=10) res.clustered <- Knet(g.clustered, nperm=10, vertex.attr="hits") res.clustered$pval plot(res.clustered) # plot results without hit clustering g.unclustered <- barabasi.game(100, directed=FALSE) g.unclustered <- SpreadHits(g.unclustered, h=10, lambda=0) res.unclustered <- Knet(g.unclustered, nperm=10, vertex.attr="hits") res.unclustered$pval plot(res.unclustered)