Learn R Programming

ktspair (version 1.0)

cv2: Crossvalidation with chosen pairs of genes

Description

This function computes the performance of the model where the pairs of genes has been selected by the user.

Usage

cv2(indice, dat, grp, cross = 5, display = FALSE, length = 40, seed = NULL, med = FALSE, healthy = NULL)

Arguments

indice
The pairs that should be used to construct the ktsp object
dat
Can either be (a) a matrix of m lines (the gene expressions) and n columns (the observations) or (b) an eSet object.
grp
Can either be (a) a character (or numeric) vector indicating the group of each observations or (b) an integer indicating the column of pData(dat) that represents the group of the observations.
cross
The number of fold that should be used in the crossvalidation.
display
Allows the user to avoid the function ktspcalc() to print waring message over the loop.
length
This paramters allows the used to control the length of the list used in the C code.
seed
Allow the user to set a seed.
med
If the mean of the median between the two groups for each gene should be substracted to the dataset or not.
healthy
This variable is used to determine which group will be considerer as the healthy group (specificity). Need to give the label of the group.

Value

k
The selected value for k
accuracy
The estimated percentage of correct prediction reached by the k-TSP with the chosen pairs of genes.
sensitivity
The estimated sensitivity reached by the k-TSP with the chosen pairs of genes.
specificity
The estimated specificity reached by the k-TSP with the chosen pairs of genes.

Details

This function is used to compute the performance reached by the model where the pairs of genes have been chosen by the user.

References

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 (). tspair: Top Scoring Pairs for Microarray Classification. R package version 1.10.0.

See Also

ktspcalc2, cv

Examples

Run this code
  ## Not run: 
#   ## Load data
#   data(ktspdata) 
#   cv <- cv(dat, grp,cross =10)
#   ktsp <- ktspcalc(dat, grp, cv$k)
#   ktsp
#   cv
#  ## End(Not run)

Run the code above in your browser using DataLab