error using krige function

[brilliana]

Hi, i am new to rstudio, and i tried the krige function. But i always got text warning "There were 50 or more warnings (use warnings() to see the first 50)" after i run my krige function and keep getting my var1.pred and var1.var to NA. Can someone pls help me? Thanks

[nirgrahamuk]

This is not a standard function in R. Therefore please declare the package.
Have you looked at the warnings ?
There is likely a problem with the data you provide to this function, but, we can't see that to judge it.

Please read this guide to understand how to share representative data.

FAQ: How to do a minimal reproducible example ( reprex ) for beginners

[brilliana]

hi, thanks for answering.
im sorry if my english and the format for sharing data is bad
i used this package :

library(ggplot2)
library(gstat)
library(sp)

I tried to do the kriging interpolation using the example data below. X and Y are the location coordinates. while Depth is the data that I have used to make kriging interpolation.

No	X	Y	Depth
1	436162	9225207	72
2	421918	9228627	54
3	421796	9228805	66
4	425039	9226890	55
5	438778	9224380	60
6	428920	9219980	15
7	435402	9215807	72
8	444398	9218928	62
9	424909	9225120	55
10	430204	9225798	17

and the code is

koordinat <- read.csv("coba_koordinat432.csv")
coordinates(koordinat) = ~X+Y
plot(koordinat)
x.range <- as.integer(range(koordinat@coords[,1]))
y.range <- as.integer(range(koordinat@coords[,2]))
grd <- expand.grid(x=seq(from=x.range[1], to=x.range[2], by=200),
y=seq(from=y.range[1], to=y.range[2], by=200))
coordinates(grd) <- ~ x+y
gridded(grd) <- TRUE
plot(grd, cex=1.5)
points(koordinat, pch=1, col='red', cex=1)
variogcloud<-variogram(Depth~1, locations=koordinat, data=koordinat, cloud=TRUE)
plot(variogcloud)
semivariog<-variogram(Depth~1, locations=koordinat, data=koordinat)
plot(semivariog)
semivariog
model.variog<-vgm(psill=600, model="Exp", nugget=80, range=8675)
fit.variog<-fit.variogram(semivariog, model.variog)
plot(semivariog, fit.variog)
krig<-krige(formula=Depth ~ 1, locations=koordinat, newdata=grd, model=model.variog)
krig.output=as.data.frame(krig)
names(krig.output)[1:3]<-c("long","lat","var1.pred")

and after i run the code above, i got warning massage like this,

Warning messages:
1: In flat_str(content, breaks) : Coercing content to character
2: c(421779, 421979, 422179, 422379, 422579, 422779, 422979,....,etc