library(sf)
library(spdep)
library(mapview)
head(Aragonez)
summary(Aragonez)
AragonezSF
Aragonez3763Vor
plot(Aragonez3763Vor[,"yield"], key.pos=4)
Aragonez3763 <- st_transform(AragonezSF, crs=3763)
AragonezPoints <- as_Spatial(Aragonez3763)
str(AragonezPoints)
spplot(AragonezPoints ,zcol="yield", key.space="right")
AragonezPoints$yieldct <- AragonezPoints$yield - mean(AragonezPoints$yield)
Arag.lm <- lm(yield ~ rowm + colm , data=AragonezPoints)
AragonezPoints$yieldldt <- AragonezPoints$yield - fitted(Arag.lm)
head(AragonezPoints@data)
spplot(AragonezPoints, layout=c(2,1), zcol=c("yieldct", "yieldldt"))
bubble(AragonezPoints, zcol="yieldct")
bubble(AragonezPoints, zcol="yieldldt")
Aragonez3763Vor$yieldct <- Aragonez3763Vor$yield-mean(Aragonez3763Vor$yield)
Aragonez3763Vor$yieldldt <- residuals(Arag.lm)
plot(Aragonez3763Vor[,c("yieldldt")], key.pos=2)
Aragonez3763Grid$yieldct <- Aragonez3763Grid$yield-mean(Aragonez3763Grid$yield)
Aragonez3763Grid$yieldldt <- residuals(Arag.lm)
plot(Aragonez3763Grid[,c("yieldldt")], key.pos=4)
vignette("nb")
cell2nb(3,3)
dnearneigh(AragonezPoints, d1=0, d2=3)
plot(dnearneigh(Aragonez3763, d1=0, d2=3), coords=st_geometry(Aragonez3763), col="red", pch=16)
dnearneigh(AragonezPoints, d1=0, d2=5)
plot(dnearneigh(AragonezPoints, d1=0, d2=5),coord=coordinates(AragonezPoints), col="red")
knn2nb(knearneigh(Aragonez3763, k=4))
plot(knn2nb(knearneigh(Aragonez3763, k=4)), st_geometry(Aragonez3763), col="red", pch=16)
poly2nb(Aragonez3763Grid)
plot(poly2nb(Aragonez3763Grid), coords=st_geometry(Aragonez3763), col="red", pch=16)
nb2mat(cell2nb(3,3))
nb2listw(cell2nb(3,3))
str(nb2listw(cell2nb(3,3)))
nb2listw(poly2nb(Aragonez3763Grid), style="C")
Wd3 <- nb2listw(dnearneigh(AragonezPoints, d1=0, d2=3))
moran.test(AragonezPoints$yield, listw=Wd3)
moran.test(AragonezPoints$yield, listw=Wd3, randomisation=FALSE)
geary.test(AragonezPoints$yield, listw=Wd3)
moran.mc(AragonezPoints$yield, listw=Wd3, nsim=10000)
lm.morantest(Arag.lm, listw=Wd3)
nb.k4 <- knn2nb(knearneigh(AragonezPoints, k=4))
nblag(nb.k4, maxlag=3)
sp.correlogram(nb.k4, var=AragonezPoints$yieldldt, method="I", order=3)
plot(sp.correlogram(nb.k4, var=AragonezPoints$yieldldt, method="I", order=10))
variogram(yield ~ 1, data=AragonezPoints)
variogram(yield ~ 1, data=AragonezPoints, cutoff=75)
plot(variogram(yield ~ 1, data=AragonezPoints, cutoff=75))
variogram(yield ~ colm + rowm, data=AragonezPoints)
library(geoR)
variog(coords=coordinates(AragonezPoints), data=AragonezPoints$yieldldt)
plot(variog(coords=coordinates(AragonezPoints), data=AragonezPoints$yieldldt))
show.vgms()
AragVarioLin <- variogram(yield ~ colm + rowm, data=AragonezPoints, locations=coordinates(AragonezPoints))
m.fit <- fit.variogram(AragVarioLin,  model=vgm(psill=0.3,"Exp", range=7, nugget=0.8))
m.fit
plot(AragVarioLin, m.fit)
m2.fit <- fit.variogram(AragVarioLin, model=vgm(psill=0.3,"Sph", range=7, nugget=0.8))
m2.fit
plot(AragVarioLin, m2.fit)
AragVariog <- variog(coords=coordinates(AragonezPoints), data=AragonezPoints\$yieldldt, max.dist=80)
plot(AragVariog)
lines.variomodel(cov.model="exp", cov.pars=c(0.2, 20), nugget=0.9)
variogram(yield ~ colm + rowm , data=AragonezPoints, alpha=c(0,90))
head(meteo)
round(cor(meteo[,3:7]),d=2)
meteo.sf <- st_as_sf(meteo, coords=c("lon","lat"), crs=4326)
meteo.sp <- as_Spatial(meteo.sf)
plot(meteo.sf, pch=16)
library(mapview)
mapView(meteo.sf, zcol="stl1")
library(dismo)
meteo.voronoi <- voronoi(meteo.sp)
mapView(meteo.voronoi, zcol="tp")
gobj <- gstat(NULL, "t2m", t2m ~ coords.x1 + coords.x2, meteo.sp)
gobj <- gstat(gobj, "stl1" , stl1 ~ coords.x1 + coords.x2, meteo.sp)
gobj <- gstat(gobj, "stl2" , stl2 ~ coords.x1 + coords.x2, meteo.sp)
gobj <- gstat(gobj, "sund" , sund ~ coords.x1 + coords.x2, meteo.sp)
gobj <- gstat(gobj, "tp" , tp ~ coords.x1 + coords.x2, meteo.sp)
gobj
vario.meteo <- variogram(gobj)
plot(vario.meteo)
vmeteo.fit <- fit.lmc(vario.meteo, gobj, vgm(psill=1, "Sph", range=800, nugget=1))
vmeteo.fit
plot(vario.meteo, vmeteo.fit)