swh:1:snp:ffdd0a7d2c8ea15ad41d45b3b178f668bd942287
Tip revision: bea2788f0c0fbbd150a73565cd2689e0970629e5 authored by Derek Young on 01 May 2007, 04:51:41 UTC
version 0.2.0
version 0.2.0
Tip revision: bea2788
depth.r
######################################################################
# the following depth function can compute the depth for #
# multi points at one time #
######################################################################
#dyn.load("~fxuan/depth/spherical/sphericaldepth.so")
depth <- function(pts,x,Cx=var(x)){
x=x%*%t(chol(solve(Cx)))
pts=pts%*%t(chol(solve(Cx)))
y <- .Fortran("mudepth",
as.integer(nrow(x)),
as.integer(nrow(pts)),
as.integer(ncol(x)),
as.single(pts),
as.single(x),
as.integer(1:nrow(pts)),
as.single(1:nrow(pts)),
PACKAGE="mixtools")
count <- y[[6]]
n <- nrow(x)
# depth <- sqrt(8)*(count-n*(n-1)/4)/sqrt((n*(n-1)))#this is to standardize the depth
# depth <- (factorial(n)/(2*factorial(n-2)))^(-1)*count
depth <- choose(n,2)^(-1)*count
depth
}