Hi all Can you help me label the points Si' s accordingly with subscripts, also to write a distance on between location S0 (0.5,0.5 ) and all order points Si lines if i have this kind of code: library(geoR) coords<-matrix(c(0,0,1,1,0,1,0,1), nrow=4, ncol=2) data<-c(3,5,6,4) coordata<-data.frame(coords=coords,data=data) geo<-as.geodata(coordata) plot(geo$coords, xlim = c(0, 1.2), ylim = c(0,1.2), xlab = "X-axis", ylab = "Y-axis") locations<- matrix(c(0.5, 0.5), ncol = 2) text(locations, as.character(0), col = "red") polygon(x = c(0, 1, 1, 0), y = c(0, 0, 1, 1),lty = 2) szero<- krige.conv(geo, locations=locations, krige=krige.control(cov.pars=c(1,1))) [[alternative HTML version deleted]]