パッケージ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-12-0301aZ519 から取得できます 。
エラーメッセージ：
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)