data(pu)
 kNN <- similarity(pu@data[2:ncol(pu)], k = 4, frequency = FALSE, ID = pu@data$UNIT_ID)  
## Not run:   
#  kNN <- similarity(pu@data[2:ncol(pu)], k = 4, frequency = TRUE, ID = pu@data$UNIT_ID)  
#  p <- kNN$freq   
#  clr <- c("#3288BD", "#99D594", "#E6F598", "#FEE08B", 
#           "#FC8D59", "#D53E4F")   
#  p <- ifelse(p <= 0, clr[1], 
#         ifelse(p > 0 & p < 10, clr[2],
#           ifelse(p >= 10 & p < 20, clr[3],
#  	       ifelse(p >= 20 & p < 50, clr[4],
#  	         ifelse(p >= 50 & p < 100, clr[5],
#  	           ifelse(p >= 100, clr[6], NA))))))
#  plot(pu, col=p, border=NA)
#    legend("topleft", legend=c("None","<10","10-20",
#           "20-50","50-100",">100"),
#           fill=clr, cex=0.6, bty="n") 
#    box()
# ## End(Not run)
 
Run the code above in your browser using DataLab