Description Usage Arguments Value Examples
The PsigA package provides a set of functions to perform Principal Component analysis based on genetic pathways. Moreover PSigA incorporates an algorithm to measure the spread of the data in the PCA space and rank the gene-signatures based on their ability to split the data in distinct entities.
PsigA blah blah
1 2 3 |
data |
data frame of matrix with gene expression values where rows represent genes and columns represent samples. |
signatures |
list where every entry is a character vector of
the gene names that correspond to a gene-pathway. If the gene format is
different from the one in |
threshold |
low density cutoff. Used in |
n |
Number of grid points in each direction. Can be scalar or a length-2 integer vector. See kde2d. |
magnitude |
When TRUE the score is multiplied by the cluster density. Default: FALSE. |
scale |
a logical value indicating whether the variables should be scaled to have unit variance before the analysis takes place. See prcomp for more details. |
A data frame with 2 columns; 1) score
, the PsigA
score of a given signature and 2) size
, the number of genes in the
signature
that were found in the data
. If
parsed = FALSE
, the function returns a list, where the first entry
contains the data frame described above and the second entry the parsed
data
.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 | require(Biobase)
require(breastCancerVDX)
data(vdx)
data(RAPIN)
VDX <- parseData(data = exprs(vdx), geneIds = fData(vdx)$Gene.symbol)
scores <- PsigA(VDX, RAPIN, threshold = 0.003)
head(scores)
#plot top signature
sigPlot(VDX, RAPIN[[rownames(scores)[1]]])
sigBiplot(VDX, RAPIN[[rownames(scores)[1]]], factor(pData(vdx)$grade),
main = rownames(scores)[1])
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.