'variogramMoLibrary not found
usgs
usgs<- read.csv(file.choose(),header=T)
lat.range <- as.integer(range(usgs[,2]))
lat.range
long.range <- as.integer(range(usgs[,3]))
long.range
lat=seq(from=lat.range[1], to=lat.range[2], by=2000)
long=seq(from=long.range[1], to=long.range[2], by=2000)
plot(grd)
points(usgs)
g <- gstat(id="data", formula = data~1, locations = ~lat+long, data = usgs)
q <- variogram(g, cressie=TRUE, cutoff=210000,width=21000)
plot(q)
VariogramMolibrary(gstat)
grd <- expand.grid(lat=seq(from=lat.range[1], to=lat.range[2], by=2000))
del<- vgm(psill=0.4,model="Sph",nugget=0.05,range=150000,kappa=0.5)
plot(q,model=VariogramModel,main="Semi Variogram")
VariogramModel
pr_ok <- krige(id="data", formula = data~1, locations=~lat+long, model=VariogramModel, data=usgs, newdata=grd)
options(max.print = 1000000)
pr_ok
cbind(pr_ok$lat[1:10000],pr_ok$long[1:10000],pr_ok$data.pred[1:10000], pr_ok$data.var[1:10000])
my<-cbind(pr_ok$lat[1:10000],pr_ok$long[1:10000],pr_ok$data.pred[1:10000], pr_ok$data.var[1:10000])
write.table(my,file="PREDICTION-VAR gstat1.csv",row.names=F,sep= ",")
qqq <- matrix(pr_ok$data.pred,
length(seq(from=lat.range[1], to=lat.range[2], by=2000)),
length(seq(from=long.range[1], to=long.range[2], by=2000)) )
image(seq(from=lat.range[1], to=lat.range[2], by=2000),
seq(from=long.range[1],to=long.range[2], by=2000), qqq,
xlab="West to East", ylab="South to North", main="Predicted values")
points(usgs)
contour(seq(from=lat.range[1], to=lat.range[2], by=2000),
seq(from=long.range[1],to=long.range[2], by=2000), qqq,
add=TRUE, col="black", labcex=1)
kr_cv = krige.cv(formula = data~1, locations=~lat+long, model=VariogramModel, data=usgs)
kr_cv[1:15,]
plot(kr_cv)
write.table(kr_cv,file="CROSS VALIDATION gstat1.csv",row.names=F,sep= ",")
this is the code and when i run it it shows variogramMoLibrary error help me with that
it R language and i also have the data file with me ## Heading ## here is data file link https://drive.google.com/file/d/1qkjgxL8AE0IwCnMqsUjirSs9RIoyUOB-/view?usp=sharing
Sources
This article follows the attribution requirements of Stack Overflow and is licensed under CC BY-SA 3.0.
Source: Stack Overflow
| Solution | Source |
|---|
