[Gstat-info] variogram

[Hoogland, Tom]

Thanks, Edzer. This works fine, but how does the doubling of locations
occur?

Tom

-----Original Message-----
From: Edzer Pebesma [mailto:edzer.pebesma at uni-muenster.de] 
Sent: donderdag 20 maart 2008 16:17
To: Hoogland, Tom
Cc: gstat-info at geo.uu.nl
Subject: Re: [Gstat-info] variogram

Try

variogram(GHG_NAP~1,pnt[!is.na(pnt$GHG_NAP), ])
--
Edzer

Hoogland, Tom wrote:
>
> I use gstat in R for the first time and in the first steps I encounter

> this error:
>
> > library(gstat)
> Loading required package: sp
> > pnt<-read.table("pnt-norg.csv", header=TRUE, sep=",")
> > coordinates(pnt) = ~X + Y
> > varioghg <- variogram(GHG_NAP ~ 1, pnt)
> Error: dimensions do not match: locations 426 and data 159
> > summary(pnt)
> Object of class SpatialPointsDataFrame
> Coordinates:
>      min    max
> X 228793 229679
> Y 566722 567816
> Is projected: NA
> proj4string : [NA]
> Number of points: 213
>
> The number of locations is indeed 213 as given by summary(pnt) but the

> variogram finctions seems to find 426 locations. Can anyone tell me 
> whats wrong...
>
> Tom Hoogland
> Alterra, Centrum Bodem
> Tel: 0317-484251
> E-mail: Tom.Hoogland at wur.nl
>
> ----------------------------------------------------------------------
> --
>
> _______________________________________________
> Gstat-info mailing list
> Gstat-info at geo.uu.nl
> http://mailman.geo.uu.nl/mailman/listinfo/gstat-info
>