[R-sig-Geo] semivariogram geoR

[carolina.lang at ifop.cl]

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" <alessandracarioli at gmail.com>
Para: r-sig-geo at 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 at 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)