# [Rd] Error in funcion distancia() in package dprep(v1.0) (PR#9745)

From: <kyousan_at_sina.com>
Date: Wed, 20 Jun 2007 04:26:13 +0200 (CEST)

There is a bug in function distancia() in package dprep. See the description below.

> distancia

1 function (x, y)
2 {

```3      if (class(y) == "matrix") {
4          distancia = drop(sqrt(colSums((x - t(y))^2)))
5          distancia = t(distancia)
6      }
7      else distancia = sqrt(sum((x - y)^2))
8      distancia
```

9 }

The t(y) in line 4 should be modified to y. Before the change, see example: > (x=rnorm(4))
 0.02670846 1.24512505 -0.22145290 1.99760521
> (y=matrix(rnorm(12),4,3))

[,1] [,2] [,3]
[1,] 0.4181184 1.8980350 -1.7398861
[2,] -0.1906483 -0.1777747 0.6666913
[3,] 0.2239468 -0.2631944 -0.3586917
[4,] -1.0316230 -0.7121346 0.3587542

> distancia(x,y)

[,1] [,2] [,3] [,4]
[1,] 1.698568 2.272631 2.336264 2.847838

After the change:
a<-function (x, y)
{

```    if (class(y) == "matrix") {
distancia = drop(sqrt(colSums((x - y)^2)))
distancia = t(distancia)
```

}
else distancia = sqrt(sum((x - y)^2))     distancia
}
> a(x,y)

[,1] [,2] [,3]
[1,] 3.404299 3.587609 2.481957

The latter result is resonable according to the description of the function return value: "distancia numeric value representing the Euclidean distance between x and y, or a row matrix representing the Euclidean distance between x and each column of y."

Thanks for attention.

R-devel_at_r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-devel Received on Wed 20 Jun 2007 - 13:30:27 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Wed 20 Jun 2007 - 16:35:00 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-devel. Please read the posting guide before posting to the list.