File: Harding.R

package info (click to toggle)
r-cran-ksamples 1.2-10-1
  • links: PTS, VCS
  • area: main
  • in suites: forky, sid, trixie
  • size: 456 kB
  • sloc: ansic: 1,321; makefile: 2
file content (21 lines) | stat: -rw-r--r-- 404 bytes parent folder | download | duplicates (2) 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
 
Harding <-
function (nn) 
{
nn <- sort(as.integer(nn))
nvec <- rev(cumsum(rev(nn)))
k <- length(nn)
L1 <- sum(nn[1:(k-1)]*nvec[2:k])+1
freq <- double(L1)

out <- .C("Harding0", k=as.integer(k), L1=as.integer(L1),
		nn=as.integer(nn), nvec=as.integer(nvec),
		freq=as.double(freq), PACKAGE = "kSamples")

freq <- out$freq

if(is.nan(sum(freq)) ||
        abs(sum(freq)-1) > .0000001) freq <- NaN

freq
}