如何在 R 中为 autoKrige {automap} 循环编写公式
How to write a formula in a loop for autoKrige {automap} in R
我有一个大矩阵,列名如下:
colid=vector(length = 60)
for(i in 1"60) {
colid[i]=paste0("V",i)
}
当我在automap
中使用autoKrige
函数时,公式必须写成这样的格式:V1~1。当我尝试使用循环执行此操作时,出现错误:
library(automap)
value=list()
for(i in 1:60) {
value[[i]]=autoKrige(colid[i]~1,Mydata,new_data = newgrid)
}
Error: too many spatial dimensions: 3068 In addition: Warning message:
NAs introduced by coercion
然后我尝试解决这个问题,我测试了公式:
> colid[10]~1
colid[10] ~ 1
所以,问题是因为公式的格式不正确。我想知道如何解决这个问题?非常感谢。
使用函数as.formula
将字符转换为公式:
colid=vector(length = 60)
value=list()
for(i in 1:60) {
colid[i] <- paste0("V",i, " ~ 1")
value[[i]]=autoKrige(as.formula(colid[i]),Mydata,new_data = newgrid)
}