dista: Distance between vectors and a matrix

Description Usage Arguments Details Value Author(s) See Also Examples

View source: R/dista.R

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 calculaitons, 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.

Author(s)

Michail Tsagris

R implementation and documentation: Michail Tsagris <[email protected]> and Manos Papadakis <[email protected]>.

See Also

mahala, Dist, total.dist, total.dista

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( iris[-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

Rfast documentation built on Nov. 18, 2017, 4 a.m.