Finds genes that have similar patterns of expression.
ExpressionSet or a
matrix of gene expressions, and the
indices of the genes of interest,
genefinder returns a
list of the
numResults closest genes.
The user can specify one of the standard distance measures listed
The number of values to return can be specified. The return value is a
list with two components:
genes (measured through the desired distance method) to the genes
of interest (where X is the number of desired results returned) and
Number of results to display, starting from the least distance to the greatest.
One of "none", "range", or "zscore". Scaling is carried out separately on each row.
A vector of weights applied across the columns of
One of "euclidean", "maximum", "manhattan", "canberra", "correlation", "binary".
scale option is "range", then the input matrix is scaled using
genescale(). If it is "zscore", then the input matrix is scaled using
scale builtin with no arguments.
The method option specifies the metric used for gene comparisons. The
metric is applied, row by row, for each gene specified in
The "correlation" option for the distance method will return a value equal to 1-correlation(x).
dist for a more detailed description of the distances.
The returned value is a
list containing an entry for each gene
list entry contains an array of
distances for that gene of interest.
J. Gentry and M. Kajen
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21
set.seed(12345) #create some fake expression profiles m1 <- matrix (1:12, 4, 3) v1 <- 1 nr <- 2 #find the 2 rows of m1 that are closest to row 1 genefinder (m1, v1, nr, method="euc") v2 <- c(1,3) genefinder (m1, v2, nr) genefinder (m1, v2, nr, scale="range") genefinder (m1, v2, nr, method="manhattan") m2 <- matrix (rnorm(100), 10, 10) v3 <- c(2, 5, 6, 8) nr2 <- 6 genefinder (m2, v3, nr2, scale="zscore")
Want to suggest features or report bugs for rdrr.io? Use the GitHub issue tracker.