sampleplot {globaltest}R Documentation

Sample Plot for Global Test

Description

Produces a plot to show the influence of individual samples on the test result produced by globaltest.

Usage

sampleplot(gt, geneset, samplesubset, scale = TRUE, drawlabels = TRUE,
  labelsize = 0.6, plot = TRUE, addlegend = TRUE, ...)

Arguments

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.

Details

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.

Value

An object of type gt.barplot.

Author(s)

Jelle Goeman: j.j.goeman@lumc.nl; Jan Oosting

See Also

globaltest, geneplot, regressionplot, checkerboard.

Examples

    # Breast cancer data (exprSet) 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))
    }

[Package globaltest version 4.2.0 Index]