[R-sig-Geo] question about variogramST

[Edzer Pebesma]

Erin, see

https://stat.ethz.ch/pipermail/r-sig-geo/2013-February/017577.html

(where the same error message appeared, but in Italian!)

On 02/23/2013 04:57 AM, Hodgess, Erin wrote:
> Hello everyone:
>
> I have an STFDF object and I would like to produce a variogramST.  I'm using the version of spacetime from 2/19.
>
> However, I'm getting a weird error message.  Here is the output:
>
>> str(city4.st)
> Formal class 'STFDF' [package "spacetime"] with 4 slots
>    ..@ data   :'data.frame':     230 obs. of  1 variable:
>    .. ..$ index: num [1:230] 100 100 100 100 100 100 100 100 100 100 ...
>    ..@ sp     :Formal class 'SpatialPoints' [package "sp"] with 3 slots
>    .. .. ..@ coords     : num [1:23, 1:2] -99.7 -101.8 -97.7 -94.1 -97.5 ...
>    .. .. .. ..- attr(*, "dimnames")=List of 2
>    .. .. .. .. ..$ : NULL
>    .. .. .. .. ..$ : chr [1:2] "lon" "lat"
>    .. .. ..@ bbox       : num [1:2, 1:2] -106.5 25.9 -94 35.2
>    .. .. .. ..- attr(*, "dimnames")=List of 2
>    .. .. .. .. ..$ : chr [1:2] "lon" "lat"
>    .. .. .. .. ..$ : chr [1:2] "min" "max"
>    .. .. ..@ proj4string:Formal class 'CRS' [package "sp"] with 1 slots
>    .. .. .. .. ..@ projargs: chr " +proj=lcc +lat_1=27.5 +lat_2=35 +lat_0=18 +lon_0=-100 +x_0=1500000 +y_0=5000000 +ellps=GRS80 +datum=NAD83 +units=m +no_defs +t"| __truncated__
>    ..@ time   :An ‘xts’ object on 2001-01-01/2010-01-01 containing:
>    Data: int [1:10, 1] 1 2 3 4 5 6 7 8 9 10
>   - attr(*, "dimnames")=List of 2
>    ..$ : NULL
>    ..$ : chr "timeIndex"
>    Indexed by objects of class: [Date] TZ: UTC
>    xts Attributes:
>   NULL
>    ..@ endTime: POSIXct[1:10], format: "2002-01-01" "2003-01-01" "2004-01-01" "2005-01-01" ...
>> variogramST(index~1,city4.st)
> Error in if (nrow(x) < nrow(y)) matrix(t(apply(x, 1, spDiN1, y = y, ll = longlat)),  :
>    argument is of length zero
>>
>
> Any suggestions would be much appreciated.
>
> Thanks,
> Sincerely,
> Erin
>
>
> 	[[alternative HTML version deleted]]
>
>
>
> _______________________________________________
> R-sig-Geo mailing list
> R-sig-Geo at r-project.org
> 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  http://ifgi.uni-muenster.de
http://www.52north.org/geostatistics      e.pebesma at wwu.de