1#
2#This code returns the clusters found by KUlldorff's test
3#
4
5get.knclusters<-function(d, knresults)
6{
7
8	ncl<-nrow(knresults)
9	clusters<-as.list( rep(NA, ncl) )
10
11	for(i in 1:ncl)
12	{
13		x<-knresults$x[i]
14		y<-knresults$y[i]
15		s<-knresults$size[i]
16
17		dst<-(d$x-x)^2+(d$y-y)^2
18		clusters[[i]]<-(order(dst))[1:s]
19	}
20
21	return(clusters)
22}
23
24