dista: Distance between vectors and a matrix In Rfast: A Collection of Efficient and Extremely Fast R Functions

Description

Distance between vectors and a matrix.

Usage

 1 dista(xnew, x, type = "euclidean", k = 0, index = FALSE, trans = TRUE, square = FALSE)

Arguments

 xnew A matrix with some data or a vector. x A matrix with the data, where rows denotes observations (vectors) and the columns contain the variables. type This can be either "euclidean" or "manhattan". k Should the k smaller distances or their indices be returned? If k > 0 this will happen. index In case k is greater than 0, you have the option to get the indices of the k smallest distances. trans Do you want the returned matrix to be transposed? TRUE or FALSE. square If you choose "euclidean" as the method, then you can have the optino to return the squared Euclidean distances by setting this argument to TRUE.

Details

The target of this function is to calculate the distances between xnew and x without having to calculate the whole distance matrix of xnew and x. The latter does extra calculations, which can be avoided.

Value

A matrix with the distances of each xnew from each vector of x. The number of rows of the xnew and and the number of columns of xnew are the dimensions of this matrix.

Michail Tsagris.

Examples

 1 2 3 4 5 6 7 8 9 10 11 12 13 xnew <- as.matrix( iris[1:10, 1:4] ) x <- as.matrix( -c(1:10), 1:4] ) a <- dista(xnew, x) b <- as.matrix( dist( rbind(xnew, x) ) ) b <- b[ 1:10, -c(1:10) ] sum( abs(a - b) ) ## see the time x <- matrix( rnorm(1000 * 4), ncol = 4 ) system.time( dista(xnew, x) ) system.time( as.matrix( dist( rbind(xnew, x) ) ) ) x<-b<-a<-xnew<-NULL