Description Usage Arguments Value Author(s) References Examples
View source: R/networkBasedSVM.R
mapping must be a data.frame with at least two
columns. The column names have to be
c('probesetID','graphID'). Where 'probesetID' is
the probeset ID present in the expression matrix (i.e.
colnames(x)) and 'graphID' is any ID that
represents the nodes in the diffusionKernel (i.e.
colnames(diffusionKernel) or
rownames(diffusionKernel)). The purpose of the
this mapping is that a gene or protein in the network
might be represented by more than one probe set on the
chip. Therefore, the algorithm must know which
genes/protein in the network belongs to which probeset on
the chip.
1 2 3 |
exps |
a p x n matrix of expression measurements with p samples and n genes. |
y |
a factor of length p comprising the class labels. |
DEBUG |
should debugging information be plotted. |
n.inner |
number of fold for the inner cross-validation. |
scale |
a character vector defining if the data
should be centered and/or scaled. Possible values are
center and/or scale. Defaults to
|
sd.cutoff |
a cutoff on the standard deviation (sd) of genes. Only genes with sd > sd.cutoff stay in the analysis. |
lambdas |
a set of values for lambda regularization parameter of the L_∞-Norm. Which, if properly chosen, eliminates factors that are completely irrelevant to the response, what in turn leads to a factor-wise (subnetwork-wise) feature selection. The 'best' lambda is found by an inner-cross validation. |
adjacencyList |
a adjacency list representing the
network structure. The list can be generated from a
adjacency matrix by using the function
|
a networkBasedSVM object containing
features |
the selected features |
lambda.performance |
overview how different values of lambda performed in the inner cross validation |
fit |
the fitted network based SVM model |
Marc Johannes JohannesMarc@gmail.com
Zhu Y. et al. (2009). Network-based support vector machine for classification of microarray samples. BMC Bioinformatics
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 | ## Not run:
library(Biobase)
data(sample.ExpressionSet)
x <- t(exprs(sample.ExpressionSet))
y <- factor(pData(sample.ExpressionSet)$sex)
# create the mapping
library('hgu95av2.db')
mapped.probes <- mappedkeys(hgu95av2REFSEQ)
refseq <- as.list(hgu95av2REFSEQ[mapped.probes])
times <- sapply(refseq, length)
mapping <- data.frame(probesetID=rep(names(refseq), times=times), graphID=unlist(refseq),
row.names=NULL, stringsAsFactors=FALSE)
mapping <- unique(mapping)
library(pathClass)
data(adjacency.matrix)
matched <- matchMatrices(x=x, adjacency=adjacency.matrix, mapping=mapping)
ad.list <- as.adjacencyList(matched$adjacency)
res.nBSVM <- crossval(matched$x, y, theta.fit=fit.networkBasedSVM, folds=3, repeats=1, DEBUG=TRUE,
parallel=FALSE, adjacencyList=ad.list, lambdas=10^(-1:2), sd.cutoff=50)
## End(Not run)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.