Message-ID: <786850361.1555668.1405690870104.JavaMail.zimbra@ifop.cl>
Date: 2014-07-18T13:41:10Z
From: carolina.lang at ifop.cl
Subject: semivariogram geoR
In-Reply-To: <A04E92E9-75C8-44A3-B9D0-E559783C4622@gmail.com>
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)