Eu estou querendo implementar uma função de distância euclideana ponderada pelo erro em R (similar, mas não exatamente a mesma que Kumar & Patel. 2005. Rutcor Research Report. RRR 12-2005). Eu já sei que quero usar o pr_DB dentro do proxy para isso. Aqui é onde eu me deparo com problemas. A primeira coisa que sempre faço quando olho para um novo pacote é recortar e colar os exemplos em R para ver se eles fazem o que eu espero que eles façam. Este não é o que aconteceu neste caso. Comecei com o exemplo simples incluído, a saber:

mydist <- function(x,y) x * y
pr_DB$set_entry(FUN = mydist, names = c("test", "mydist"))

Ok, até agora tudo bem. Eu acho que, se funcionar como uma função normal do tipo "dist", eu adicionarei uma informação muito simples.

>toydat
        x           y
a -0.12817993 -1.03238513
b  1.56200731  0.93826937
c -1.24051847 -1.31005852
d -1.12892553 -1.57133401
e -1.10098308  0.06577006

Apenas um brinquedinho, claro. SOOOOOO, eu experimento o seguinte:

toydist <- dist(toydat,method="mydist")

Eu recebo a seguinte mensagem:

Error in do.call(".External", c(list(CFUN, x, y, pairwise, if (!is.function(method)) get(method) else method),  : 
not a scalar return value

Em um palpite, eu tentei:

toydist <- dist(toydat,method="Euclidean")
toydist <- dist(toydat,method="Manhatten")

e outros. Todos eles funcionam como esperado. Estou presumindo que há algo especial que deve ser feito com a fórmula básica para que ele funcione de maneira adequada para calcular as matrizes de distância. O que eu quero calcular é sqrt ((xi - xj)2) + (yi - yj)2) + ... (ni - nj)2) ÷ sqrt ((σxi2 + σxj2) + (σyi2 + σyj2) + ... (σni2 + σnj2), em todas as combinações de pares para o meu conjunto de dados.

Eu percebo que eu poderia fazer isso por planilha para um único conjunto de dados, mas eu quero fazer alguns bootstrapping de cluster, que eu sei fazer em R. Não é prático com uma planilha.