Вариограмма сюжета для Силла, Самородка, Спектра
Я работаю над данными Temporal-spatio. Данные в структуре STFDF.
Я нахожусь на этапе применения вариограммы и пытаюсь определить Nugget, Sill и Range из графика вариограммы, который должен выглядеть следующим образом 
,
что-то вроде этого: 
Однако, когда я строю свою вариограмму, она показывает:
plot(vario, main="Flow")
А также plot(vario,map=FALSE, main="Spatio-temporal correlation") 
А также plot(vario,wireframe=T, main="Spatio-temporal correlation") 
Как получить сюжет, который может определить диапазон, подоконник и самородок!! Я не мог определить свои параметры из этих графиков. Я использовал функцию variogramST и попробовал слишком функцию variogram, но оба не могут создать график, который мне нужен.
vario<-variogramST(Flow~1,data=data,tunit="hours",assumeRegular=F,na.omit=T)
Есть ли способ определить Range, Nugget и Sill?
####### separable model, least squares fit
separableModel <- vgmST("separable",space=vgm(0.9,"Exp",100,0.1),time=vgm(0.9,"Exp",1000,0.1),sill=40)
separable_fit <- fit.StVariogram(model=separableModel,object=vario)
plot(vario,separable_fit,all=T,map=F)
attr(separable_fit,"MSE") # calculate the Mean Absolute Error
### product sum model, least squares fit
ProductSum <- vgmST("productSum",space =vgm(psill= 15000,"Exp",range= 0.07, nugget= 0),time=vgm(psill= 7500,"Exp",range= 0.07, nugget=0),k=0.000649771173341919)
ProductSum_fit <- fit.StVariogram(model=ProductSum,object=vario)
ProductSum_fit
attr(ProductSum_fit,"MSE")
plot(vario,ProductSum_fit,all=T,map=F)
plot(vario,ProductSum_fit,all=T)
plot(vario,ProductSum_fit, all=T, wireframe=T)
#### product sum model, manual fit
ProductSum_man <- vgmST("productSum",space=vgm(7500,"Exp",2e5,0,add.to=vgm(10,"Exp",9e3,0)), time=vgm(7500,"Exp",800,9.5,add.to=vgm(0,"Exp",850,0)),k=0.035)
ProductSum_man_fit <- fit.StVariogram(model=ProductSum_man,object=vario)
attr(ProductSum_man_fit,"MSE")
plot(vario,ProductSum_man_fit,all=T,map=F)
plot(vario,ProductSum_man_fit,all=T)