[R-sig-Geo] kriging output

[Dave Depew]

Thanks,
This worked.

I'm still confused why the if else statement didn't work...
If one wanted to do conditional arithmetic would a for statement bee needed?


e.g.
 meuse.grid[["class"]] = for(i in 1:length(meuse.grid[["dist"]])){
if (meuse.grid[["dist"]]<0.5) 
{meuse.grid[["class"]]=10*meuse.grid[["dist]]} else
{meuse.grid[["class"]]=100* meuse.grid[["dist"]]}
 }




Edzer Pebesma wrote:
> Dave Depew wrote:
>> Hi all,
>> I've got a question regarding kriging outputs. I have an interpolated 
>> dataset which due to the nugget effect contains some negative values 
>> as the predictions. I would like to truncate these @ "0", rather than 
>> having them as a negative prediction.
>> I've tried something similar with the meuse data set...
>>
>> data(meuse.grid) coordinates(meuse.grid) = c("x", "y") 
>> gridded(meuse.grid) <- TRUE
>> meuse.grid[["idist"]] = 1 - meuse.grid[["dist"]]
>> Somewhat arbitrary, but essentially, I'd like to set all "dist" data 
>> less than 0.5 to "0", and have the rest remain as they are...
>>
>>
>> meuse.grid[["class"]] = for(i in 1:length(meuse.grid[["dist"]])){
>> if (meuse.grid[["dist"]]<0.5) {meuse.grid[["class"]]==0} else
> I guess you meant an = rather than a == here
>> {meuse.grid[["class"]]=meuse.grid[["dist"]]}
>> }
>>
>> All I've managed to get is the new attribute "class" all filled with 0.
>>
> Try:
>
> meuse.grid$class = meuse.grid$dist
> meuse.grid$class[meuse.grid$dist < .5] = 0
> plot(class~dist,meuse.grid)
> -- 
> Edzer
>