[STATSGRASS] gstat probability maps

[Pierluigi De Rosa]

Dear All,

I need to prepare a probability map of exceeding a value...
I can make the simulation, for example I tried to make it with data meuse:

v <- variogram(log(zinc)~1,loc= meuse) 
m <- vgm(.59, "Sph", 874, .04)
sim <- krige(formula = log(zinc)~1, locations = meuse, model = m, newdata = 
meuse.grid, nmax = 15, nsim = 5)
names(sim)
[1] "sim1" "sim2" "sim3" "sim4" "sim5"

Now How can I create a probabilility map of exceding a value (for example 5)??

Thanks
Pierluigi

-- 

Pierluigi De Rosa

Department of Enviromental & Civil engineering 
Faculty of Engineering 
Perugia (Italy)