Displaying 1 result from an estimated 1 matches for "pisqr".
Did you mean:
pisar
2009 Mar 27
1
General help for a function I'm attempting to write
...combinations
of columns, perhaps returning a matrix of distances. Thanks for any help or
suggestions.
Cheers
Colin
Function:
JostD <- function(DF, X, Y) {
Ni1 <- DF[,X]
Ni2 <- DF[,Y]
N1 <- sum(Ni1)
N2 <- sum(Ni2)
pi1 <-Ni1/N1
pi2 <-Ni2/N2
pisqr <- ((pi1+pi2)/2)^2
H1 <- 1 - sum(pi1^2)
H2 <- 1 - sum(pi2^2)
Ha <- 0.5*(H1+H2)
Da <- 1/(1-Ha)
Dg <- 1/sum(pisqr)
Db <- Dg/Da
D <- -2*((1/Db) - 1)
D
}
Sample data:
e<-c(0,0,0,4,27)
r<-c(0,1,0,7,16)
t<-c(1,0,0,16,44)
y&l...