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]]

