[R-sig-Geo] Spatio-temporal variogram error

Jonesmus Mutua jonesmus.mutua at gmail.com
Sat Feb 16 20:08:26 CET 2013


Dear Ben,
Thanks a lot, it actually solves the problem. What I was wondering is whether the time lag is only allowed upto lag 15. The data am using has 21 time points but the lag on the graphs is upto 15, no matter what I change.

Is there a way of  determining what to use as width and cutoff in the formula to compute the spatio-temporal variogram instead of trying different values severally?

Regards,
Jonesmus 

-----Original Message-----
From: Benedikt Gräler [mailto:ben.graeler at uni-muenster.de] 
Sent: 15 February 2013 18:44
To: Jonesmus Wambua
Cc: Edzer Pebesma
Subject: Re: [R-sig-Geo] Spatio-temporal variogram error

Dear Jonesmus,

I just submitted the updated version to r-forge. You will find a package build (.tar.gz) and windows binaries (.zip) attached to check whether this solves your problem.

Best wishes,

 Ben

On Freitag, 15. Februar 2013 13:31:18, Jonesmus Wambua wrote:
>
> Hi ben,
> what is the solution to the problem am having?
>
> On Fri, Feb 15, 2013 at 3:28 PM, Jonesmus Wambua 
> <jonesmus.mutua at gmail.com <mailto:jonesmus.mutua at gmail.com>> wrote:
>
>     no, he has not.
>
>
>     On Fri, Feb 15, 2013 at 3:28 PM, Edzer Pebesma
>     <edzer.pebesma at uni-muenster.de
>     <mailto:edzer.pebesma at uni-muenster.de>> wrote:
>
>         It seems Ben solved your problem, has he told you?
>
>
>
>         -------- Original Message --------
>         Subject:        Re: [R-sig-Geo] Spatio-temporal variogram error
>         Date:   Fri, 15 Feb 2013 15:26:53 +0300
>         From:   Jonesmus Wambua <jonesmus.mutua at gmail.com
>         <mailto:jonesmus.mutua at gmail.com>>
>         To:     Edzer Pebesma <edzer.pebesma at uni-muenster.de
>         <mailto:edzer.pebesma at uni-muenster.de>__>
>
>
>
>         Hi,
>         here is the data, coordinates and code. the spatial variogram
>         part of
>         the code works but the spatio-temporal part doesn't work.
>
>         Thanks for the help.
>         regards.
>
>
>         On Fri, Feb 15, 2013 at 2:16 PM, Edzer Pebesma
>         <edzer.pebesma at uni-muenster.de
>         <mailto:edzer.pebesma at uni-muenster.de>
>         <mailto:edzer.pebesma at uni-__muenster.de
>         <mailto:edzer.pebesma at uni-muenster.de>>>
>
>         wrote:
>
>             I can have a look if you send me the data, as .RData file,
>         off-list.
>
>
>             On 02/15/2013 09:26 AM, Jonesmus Wambua wrote:
>
>                 hi, am doing a spatio-temporal variogram and am
>         getting this
>                 error below
>                 that i have no response variable in my data set even
>         though it
>                 is there. my
>                 data am calling gdata and the response variable is
>         pfpr. i have
>                 tried the
>                 code below with no access. the summary of my data is
>         as follows:
>
>                     summary(gdata)
>
>                 Object of class STFDF
>                    with Dimensions (s, t, attr): (438, 21, 1)
>                 [[Spatial:]]
>                 Object of class SpatialPoints
>                 Coordinates:
>                            min      max
>                 Long -16.718 -13.8500
>                 Lat   13.176  13.7915
>                 Is projected: NA
>                 proj4string : [NA]
>                 Number of points: 438
>                 [[Temporal:]]
>                        Index              timeIndex
>                    Min.   :1990-01-01   Min.   : 1
>                    1st Qu.:1994-12-31   1st Qu.: 6
>                    Median :1999-12-30   Median :11
>                    Mean   :1999-12-30   Mean   :11
>                    3rd Qu.:2004-12-28   3rd Qu.:16
>                    Max.   :2009-12-27   Max.   :21
>                 [[Data attributes:]]
>                      Min. 1st Qu.  Median    Mean 3rd Qu.    Max.    NA's
>                       0.0     1.0     4.0    10.1    13.0   330.0    8640
>
>
>                     s=variogram(gdata at data~1, gdata)
>
>
>                     |
>
>                     |
>                                                  |   0%Error in
>         .subset2(x, i,
>                 exact = exact)
>                 : no such index at level 1
>
>
>                     s=variogram(pfpr~1, gdata)
>
>
>                     |
>
>                     |
>                                                  |   0%Error in
>                 gstat.formula(d$formula,
>                 d$data) :
>                     no response variable present in formula
>                 In addition: Warning messages:
>                 1: In min(x) : no non-missing arguments to min;
>         returning Inf
>                 2: In max(x) : no non-missing arguments to max;
>         returning -Inf
>                 3: In min(x) : no non-missing arguments to min;
>         returning Inf
>                 4: In max(x) : no non-missing arguments to max;
>         returning -Inf
>
>
>
>                 can anyone help?
>                 J
>
>                          [[alternative HTML version deleted]]
>
>                 ___________________________________________________
>                 R-sig-Geo mailing list
>         R-sig-Geo at r-project.org <mailto:R-sig-Geo at r-project.org>
>         <mailto:R-sig-Geo at r-project.__org
>         <mailto:R-sig-Geo at r-project.org>>
>         https://stat.ethz.ch/mailman/____listinfo/r-sig-geo
>         <https://stat.ethz.ch/mailman/__listinfo/r-sig-geo>
>
>                 <https://stat.ethz.ch/mailman/__listinfo/r-sig-geo
>         <https://stat.ethz.ch/mailman/listinfo/r-sig-geo>>
>
>
>             --
>             Edzer Pebesma
>             Institute for Geoinformatics (ifgi), University of Münster
>             Weseler Straße 253, 48151 Münster, Germany. Phone: +49 251
>             8333081, Fax: +49 251 8339763 <tel:%2B49%20251%208339763>
>         <tel:%2B49%20251%208339763>
>         http://ifgi.uni-muenster.de
>         http://www.52north.org/____geostatistics
>         <http://www.52north.org/__geostatistics>
>             <http://www.52north.org/__geostatistics
>         <http://www.52north.org/geostatistics>> e.pebesma at wwu.de
>         <mailto:e.pebesma at wwu.de>
>             <mailto:e.pebesma at wwu.de <mailto:e.pebesma at wwu.de>>
>
>
>
>
>
>

--
Benedikt Gräler

ifgi - Institute for Geoinformatics
University of Muenster

http://ifgi.uni-muenster.de/graeler

Phone: +49 251 83-33082
Mail: ben.graeler at uni-muenster.de



More information about the R-sig-Geo mailing list