可以如下实现方式

你好，感谢您的回复，我试了依然不行。我用的包是
library(raster)
library(rgdal)
library(dismo)
library(gstat)
我想做克里金插值，可是插值单次运行就可以，循环起来就出错，运行代码(出错部分）：library(gstat)
v<- variogram(log(rain) ~ 1, data =dsp)
出错代码：
Error in if (rat > 1e+06 || rat < 1e-06) print("a possible solution MIGHT be to scale semivariances and/or distances") :   需要TRUE/FALSE值的地方不可以用缺少值

也就是建半方差函数的过程中，log(rain)这里我改成我的列名log(X1)~1，因为我有很多数据要算，依次X1,X2...X3652这么多列。这个式子表示就成问题。
我看了帮助，log(X1)~1的部分它是variogram.formula，可能跟一般formula的定义不一样~！！
请用过这个的帮我看看解决方法~谢谢！