sampleplot {globaltest} | R Documentation |
Produces a plot to show the influence of individual
samples on the test result produced by globaltest
.
sampleplot(gt, geneset, samplesubset, scale = TRUE, drawlabels = TRUE, labelsize = 0.6, plot = TRUE, addlegend = TRUE, ...)
gt |
The output of a call to globaltest . |
geneset |
The name or number of the geneset to be plotted (only necessary if multiple genesets were tested). |
samplesubset |
A vector of names or numbers of samples to be plotted. Default: plot all samples |
scale |
Logical: should the bars be scaled to unit standard deviation? |
drawlabels |
Logical value to control drawing of the samplenames on the x-axis of the plot. |
labelsize |
Relative size of the labels on the x-axis. If it is NULL ,
the current value for par("cex.axis") is used |
plot |
If FALSE : does not plot, but only returns a
gt.barplot object. |
addlegend |
If FALSE : does not add a legend to the plot
or to the gt.barplot object. |
... |
Any extra arguments will be forwarded to the plotting function. |
The sampleplot shows a bar and a reference line for each sample. The bar shows the influence of each gene on the test statistic. Samples with a positive influence carry evidence against the null hypothesis (in favour of a significant p-value), because they are are similar in expression profile to samples with a similar clinical outcome. Samples with a negative influence bar supply evidence in favour of the null hypothesis and of a non-significant p-value: they are relatively similar in expression profile to samples with a different clinical outcome.
The influence varies around zero if the tested geneset is not associated with the outcome. Marks on the bars show the standarddeviation of the influence under the null hypothesis for those samples which are more than one standard deviation away from zero.
The color of the bar indicates the sign of the residual of Y. In a logistic model the coloring this distinguishes the original groups.
The bottom margin is adjusted to allow enough space for the longest samplename to draw under the axis.
An object of type gt.barplot
.
Jelle Goeman: j.j.goeman@lumc.nl; Jan Oosting
globaltest
, geneplot
,
regressionplot
, checkerboard
.
# Breast cancer data (ExpressionSet) from the Netherlands Cancer # Institute with annotation: data(vandeVijver) data(annotation.vandeVijver) gt <- globaltest(vandeVijver, "StGallen", annotation.vandeVijver) if (interactive()){ sampleplot(gt[1]) } sp <- sampleplot(gt[1], plot = FALSE) if (interactive()){ plot(sort(sp)) }