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