backlog

[Yetta Jager <zij@ornl.gov>]

Gstat users,

Thanks to Edzer for the suggestions on automating variogram
estimation.

We're having very good success using Gstat to generate alternative
realizations that have similar structure to an original map.  Now,
we're trying to understand how to control the error rates and what
the lower bounds on generated map errors are. 

This question resurfaced during a meeting yesterday.  During conditional
simulation, we subsample a grid or mask to provide conditioning data.
Gstat won't continue when there is overlap between the input data sites 
and the output mask sites.  We have been tricking Gstat into giving us
predictions by scooting them over slightly.

Is there any reason, in principle, why we shouldn't be able to
get predictions at zero distance?  We've also been discussing
whether to expect "honoring of the data" at these locations when
there is a nugget in simulations -- I suppose the kriging variance is 
still zero?

Any ideas?

Thanks,

Yetta
------------------------------------------------------
Yetta Jager                                             
Environmental Sciences Division   
Oak Ridge National Laboratory         
P.O. Box 2008, MS 6036             
Oak Ridge, TN  37831-6036          
OFFICE: 423/574-8143
FAX:    423/576-8543 
Work email: jagerhi@ornl.gov 
Home email: hjager@utk.edu
WEBpage: http://www.esd.ornl.gov/~zij/
-----------------------------------------------------

Descartes goes into a bar where the bar tender asks him if he would like a
beer. Descartes says "I think not" and then disappears.