Description Usage Arguments Details Value note Author(s) References See Also Examples
This function computes the k pairs of genes that achieved the maximum difference between the sensitivity and the specificity (in absolute value) between two specific groups based on the comparison of the expressions of the two genes present in the pairs. The function ktspcalc() is the general function and uses the function kts.pair once the dataset has been prepared for this function. The function ktspcalc() is also able to deal with eSets.
1 |
dat |
A matrix of m lines (the gene expressions) and n columns (the observations). |
grp |
A vector of 0 and 1 for the groups of the observations. |
k |
The number of pairs of genes that the function will select. |
display |
Allows the user to avoid the function kts.pair() to print warning message (mainly used in the function crossvalidation). |
length |
This paramters allows the used to control the length of the list used in the C code. |
med |
If the mean of the median between the two groups for each gene should be substracted to the dataset or not. |
This function only works with matrices and vector of group containing only 0 and 1. For a more general use (eSets and labels for the groups) see the function ktspcalc(). This classifier can only be used for classifications with two groups. The k-TSP was introduced in Tan et al. (2005) and is an extension of the TSP, which was presented in Geman et al. (2004).
A ktsp object with the following elements:
index |
A k by 2 matrix composed of genes where the ith row stands for the ith best pair of genes with the restriction that a gene can appear in only one pair. The pairs are seleceted with respect to the score Delta and Gamma (in case of ties), see Tan et al. (2005) for more details about the k-TSP. |
ktspcore |
A vector of size k containing the scores Delta achieved by each selected pair of genes. The score Delta is based on the sensitivity and the specificity of a pair, see Geman et al. (2004) for more details. |
grp |
The group for each observation in a binary form |
ktspdat |
The row i and the row i+k represents the expressions of the genes present in the ith pair. |
k |
The number of pairs of genes. |
labels |
The name of the two groups that were present in the original variable grp. |
rankscore |
The score Gamma achieved by each pair of genes, for more details on this score see Geman et al. (2004). |
accuracy |
A vector of the estimated percentage of correct prediction for the k-TSP with k=1,3,5,7,9. |
accuracy_k |
The estimated percentage of correct prediction of the k-TSP with the selected k. |
sensitivity |
A vector of the estimated sensitivity for the k-TSP with k=1,3,5,7,9. |
specificity |
A vector of the estimated specificity for the k-TSP with k=1,3,5,7,9. |
med |
If the mean of the medians within each group has been substracted to the dataset return the values of the mean of the median, return FALSE otherwise |
The length sets to the list used in the C code (defined by the paramter length) has to be at least as big as k.
Julien Damond julien.damond@gmail.com
D. Geman, C. d'Avignon, D. Naiman and R. Winslow, "Classifying gene expression profiles from pairwise mRNA comparisons," Statist. Appl. in Genetics and Molecular Biology, 3, 2004.
A.C. Tan, D.Q. Naiman, L. Xu, R.L. Winslow, D. Geman, "Simple decision rules for classifying human cancers from gene expression profiles," Bioinformatics, 21: 3896-3904, 2005.
J. Damond, supervised by S. Morgenthaler and S. Hosseinian, "Presentation and study of robustness for several methods to classify individuals based on their gene expressions", Master thesis, Swiss Federal Institute of Technology Lausanne (Switzerland), 2011.
J. Damond, S. Morgenthaler, S. Hosseinian, "The robustness of the TSP and the k-TSP and the computation of ROC curves", paper is submitted in Bioinformatics, December 2011.
Jeffrey T. Leek <jtleek@jhu.edu> (). tspair: Top Scoring Pairs for Microarray Classification. R package version 1.10.0.
ktspcalc
, ktspplot
,predict.ktsp
, summary.ktsp
1 2 3 4 5 6 7 8 9 10 11 12 13 |
Loading required package: Biobase
Loading required package: BiocGenerics
Loading required package: parallel
Attaching package: 'BiocGenerics'
The following objects are masked from 'package:parallel':
clusterApply, clusterApplyLB, clusterCall, clusterEvalQ,
clusterExport, clusterMap, parApply, parCapply, parLapply,
parLapplyLB, parRapply, parSapply, parSapplyLB
The following objects are masked from 'package:stats':
IQR, mad, sd, var, xtabs
The following objects are masked from 'package:base':
Filter, Find, Map, Position, Reduce, anyDuplicated, append,
as.data.frame, cbind, colMeans, colSums, colnames, do.call,
duplicated, eval, evalq, get, grep, grepl, intersect, is.unsorted,
lapply, lengths, mapply, match, mget, order, paste, pmax, pmax.int,
pmin, pmin.int, rank, rbind, rowMeans, rowSums, rownames, sapply,
setdiff, sort, table, tapply, union, unique, unsplit, which,
which.max, which.min
Welcome to Bioconductor
Vignettes contain introductory material; view with
'browseVignettes()'. To cite Bioconductor, see
'citation("Biobase")', and for packages 'citation("pkgname")'.
Warning message:
In kts.pair(dat, grp, 3) : NAs introduced by coercion
[,1] [,2]
[1,] 1 2
[2,] 3 4
[3,] 5 6
[1] 0 0 0
The k-TSP cannot be computed for k = 9
The value of k has been reduced to k = 5
Warning message:
In kts.pair(dat, grp, 9) : NAs introduced by coercion
k-TSP object with: 5 TSPs
Pair: TSP Score Indices
TSP 1 : 0 1 2
TSP 2 : 0 3 4
TSP 3 : 0 5 6
TSP 4 : 0 999 1000
TSP 5 : 0 7 8
The k-TSP cannot be computed for k = 9
The value of k has been reduced to k = 5
Warning message:
In kts.pair(dat, grp, 9, length = 40) : NAs introduced by coercion
k-TSP object with: 5 TSPs
Pair: TSP Score Indices
TSP 1 : 0 1 2
TSP 2 : 0 3 4
TSP 3 : 0 5 6
TSP 4 : 0 999 1000
TSP 5 : 0 7 8
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.