Dear all, I think earlier i didn't clear myself. Now i am sending my data file and the code. My data set is attached with the mail where first and second column is X and Y coordinates and the column is the value of samples. I hope this time i am able to communicate. Thanking you Pujan Code: > mol.gr <- read.table(“D:/mol.dat”) > mol.gr[,"x"] <-mol.gr[,1] > mol.gr[,"y"] <-mol.gr[,2] > mol.gr[,"Green Band Reflectance"] <-mol.gr[,3] >map<-as.geodata(mol.gr) # make a geoR object to use various ready-made R >unctions > coal.var <- variog(map, max.dist =283.6) # omnidirectional semivariogram up > to half max distance > plot(coal.var) _________________________________________________________________ Helping your favorite cause is as easy as instant messaging. You IM, we give.
mol.rar
Description: Binary data
_______________________________________________ R-sig-Geo mailing list R-sig-Geo@stat.math.ethz.ch https://stat.ethz.ch/mailman/listinfo/r-sig-geo