filtering the component in simtub?

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

filtering the component in simtub?

Postby MaximeBeauchamp » Thu Dec 07, 2017 4:00 pm

Hi,

I would like to carry out a lot of conditional simulations of one of the basic components of my variogram. I saw in the kriging function the parameter cov.extract. Am I right to say that it enables you to estimate spatially the required component of the process? If so, is the same possible with the simtub function because I don't see this option in the help of the function.

Thanks in advance,

Maxime
MaximeBeauchamp
 
Posts: 9
Joined: Mon May 26, 2014 10:12 am

Re: filtering the component in simtub?

Postby MaximeBeauchamp » Fri Dec 08, 2017 11:07 am

Hi again,

To clarify my question, what I want is to produce a conditional simulation of my process without the nugget effect. If it is not possible via simtub but the option is available in kriging (basically I just want to remove the Nugget effect that is always the first component of the variogram in my setup), I was thinking to do the conditional simulation in a more indirect way by doing something like:

Code: Select all
# 0) model without the nugget effect (i.e removing the first basic structure)
data.model2=new('model',ndim=2,nvar=1,basics=data.model@basics[seq(2,length(data.model@basics))])
# 1) non conditional simulations on x_alpha's and x_0's
grid.db=simtub(,grid.db,data.model2,nbsimu=1,nbtuba=1000)
data.db=simtub(,data.db,data.model2,nbsimu=1,nbtuba=1000)
# 2) kriging of the true data on the grid x_0's
data.db=db.locate(data.db,"Gaussian.thick",'z')
grid.db=kriging(data.db,grid.db,data.model,data.u.neigh,cov.extract=c(2:length(data.model@basics)),radix="KS.Part")
# 3) kriging of the non conditional simulated data on the grid x_0's
data.db=db.locate(data.db,"Simu.V1.S1",'z')
grid.db=kriging(data.db,grid.db,data.model,data.u.neigh,cov.extract=c(2:length(data.model@basics)),radix="KS.Part")
# 4) final value of the conditional simulations on x_0's
grid.db=db.add(grid.db,Simu.Gaussian.thick.S1=db.extract(grid.db,"KS.Part.Gaussian.thick.estim")+(db.extract(grid.db,"Simu.V1.S1")-db.extract(grid.db,"KS.Part.Simu.V1.S1.estim")))


Do you think it is correct? Even if it is the case, I am still open to suggestions for having a direct solution with simtub.

Thanks for your help,

Maxime
MaximeBeauchamp
 
Posts: 9
Joined: Mon May 26, 2014 10:12 am


Return to Other Troubleshooting

Who is online

Users browsing this forum: No registered users and 12 guests

cron