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]]

