Re: error in gstat universal kriging example

["Edzer J. Pebesma" <e.pebesma@geog.uu.nl>]

femke wrote:

> Hello,
>  
> I'm trying to play with the universal kriging example, and using the 
> code as specified in the demos folder, but i get the follwing error 
> message:
>  
> > x <- krige(log(zinc) ~ -1 + sqrt(dist)+ part.a + part.b,
> + ~ x + y, meuse, meuse.grid,
> + model = vgm(.149, "Sph", 700, .0674))
> Error in eval(expr, envir, enclos) : Object "part.a" not found
>  
>  
> Trying to fix it with the following solved that problem, but gave me 
> another:
>  
>  
> > x <- krige(log(zinc) ~ -1 + sqrt(meuse.grid$dist)+ meuse.grid$part.a 
> + meuse.grid$part.b,
> + ~ x + y, meuse, meuse.grid,
> + model = vgm(.149, "Sph", 700, .0674))
> Error in model.frame(formula, rownames, variables, varnames, extras, 
> extranames,  :
>         variable lengths differ
>  
> I'm not sure where the problem is as the lengths of the variables are 
> the same:
>  
> > length(meuse.grid$part.a)
> [1] 3103
> > length(meuse.grid$part.b)
> [1] 3103
> > length(meuse.grid$dist)
> [1] 3103
>
> Any ideas about what might be going wrong?
>  
> Cheers,
>  
> femke

Femke, did you try a full

demo(examples)
?
--
Edzer