DEar Alessandra, In the function variog(), you complete with uvec=(??), this allow you define the lag distance and max.dist= ??, define the max distance for your experimental variogram or domain, An Example: uvec=seq(0,37040,by=926), max.dist=37040
regards, Carolina Lang Carolina Lang Abarzúa Investigador Instituto de Fomento Pesquero Av. Blanco #839, Valparaiso, Chile. Fono: (32) 2151572 www.ifop.cl ----- Mensaje original ----- De: "Alessandra Carioli" <alessandracari...@gmail.com> Para: r-sig-geo@r-project.org Enviados: Viernes, 18 de Julio 2014 4:00:18 Asunto: [R-sig-Geo] semivariogram geoR Dear bloggers, I am using geoR package to do semivariograms (variog). I would like to specify 50km distances with a maximum of 2000km distance and I am at loss on how to do it. The Coordinates of my shape file are projected longitude and latitude A 13.515094 37.44574 B 8.664217 44.82956 C 13.155760 43.50457 D 7.393502 45.73274 E 11.856243 43.53362 F 13.555405 42.88887 G 8.187188 44.87850 [�] Here is my code: lat <- coordinates(shape.shp)[,1] lon <- coordinates(shape.shp)[,2] data <- cbind(lat,lon,fert) dist <- dist(data[,1:2]) summary(dist) breaks=seq(0,12,l=21) v1 <- variog(coords = fert.l[,1:2], data = fert[,5], breaks = breaks,na.rm = TRUE) v1.summary <- cbind(c(1:20), v1$v, v1$n) colnames(v1.summary) <- c("lag", "semi-variance", "# of pairs") plot(v1, type = "b", main = "Variogram") Any help on the matter is greatly appreciated. Ale [[alternative HTML version deleted]] _______________________________________________ R-sig-Geo mailing list R-sig-Geo@r-project.org https://stat.ethz.ch/mailman/listinfo/r-sig-geo ------------------------------------------------------------------------------------------------ Certificaci�n ISO 9001/2008: Sistema de Datos Bio-Pesqueros ( Arica, Iquique, Coquimbo, Valpara�so, San Antonio, Talcahuano y Calbuco, pesquer�as industriales y artesanales)
_______________________________________________ R-sig-Geo mailing list R-sig-Geo@r-project.org https://stat.ethz.ch/mailman/listinfo/r-sig-geo