1

我对包的fit.variogram 功能有一些问题。gstat我的代码、数据、错误消息和traceback()消息可以在下面找到。任何解决此问题的建议将不胜感激。

library(gstat)
mydf<-as.data.frame(read.table("Sample_data.xlsx",header = TRUE,sep = 
",",na.strings = "EMPTY"))
coordinates(mydf) = ~x+y
c.vgm.exp<-variogram(z~1,data=mydf,cutoff=60000, width = 60000/15)
c.vgm.fit<-fit.variogram(c.vgm.exp,vgm(nugget=1000000,psill=6000000,model 
="Exp",range = 40000),fit.method = 7)

数据文件可以从: https ://gtvault-my.sharepoint.com/personal/srathore6_gatech_edu/_layouts/15/guestaccess.aspx?docid=0b5983d45c9544203a5081cd501c64de4&authkey=AYhC8_ZVvGFG18JHHWMgUa4&expiration=2017-1452-021T10%

错误信息:

c.vgm.fit<-fit.variogram(c.vgm.exp,vgm(nugget=1000000,psill=6000000,model 
="Exp",range = 40000),fit.method = 7)
Error in switch(model, exponential = fit.exponential(v.object, c0 = nugget,: 
EXPR must be a length 1 vector

traceback()
1: fit.variogram(c.vgm.exp, vgm(nugget = 1e+06, psill = 6e+06, model = 
"Exp", range = 40000), fit.method = 7)
4

0 回答 0