[SOLVED] Specify lower bounds for parameter model.auto

If you have problems using RGeostats package or if you simply need some help on an unlisted subject.

[SOLVED] Specify lower bounds for parameter model.auto

Postby markm0705 » Wed Jun 26, 2019 1:56 am

I would like to set the nugget effect for a two structure model to a little above zero using model.auto.
I'm afraid the help is a little complicated I'm now sure how to fill out the 'lower' parameter for the nugget

Code: Select all
  vg.fit  <- model.auto(tmp,c("Nugget effect","Gaussian"),niter = 1000,
                        wmode = 0,draw = FALSE,flag.norm.sill = TRUE,
                        flag.noreduce = TRUE,lower = "?")
Last edited by markm0705 on Wed Jun 26, 2019 2:06 pm, edited 1 time in total.
markm0705
 
Posts: 22
Joined: Sun Mar 17, 2019 9:39 am

Re: Specify lower bounds for parameter model.auto

Postby Fabien Ors » Wed Jun 26, 2019 10:39 am

Dear

According the model.param.define function help and example here: http://rgeostats.free.fr/html/model.param.define.html
To make sure that the sill of your nugget effect structure remains a little above zero (say 0.001), you should use the following :
Code: Select all
  vg.fit <- model.auto(tmp,c("Nugget effect","Gaussian"),niter = 1000,
                       wmode = 0,draw = FALSE,flag.norm.sill = TRUE,
                       flag.noreduce = TRUE,lower = "M1V=0.001")


Hope this helps.
Best regards
Fabien Ors
Administrateur du site
 
Posts: 226
Joined: Thu Sep 20, 2012 1:07 pm

Re: Specify lower bounds for parameter model.auto

Postby markm0705 » Wed Jun 26, 2019 2:05 pm

Thanks
markm0705
 
Posts: 22
Joined: Sun Mar 17, 2019 9:39 am


Return to Other Troubleshooting

Who is online

Users browsing this forum: No registered users and 5 guests

cron