Hello Geo-Rs!
I'm using the latest version of gstat and spacetime. I'm trying to fit a
variogramST, but I get the following:
> city1b.st <-
> STFDF(sp=tempb,time=liltime,data=data.frame(index=city1b.df$index))
>
> str(city1b.st)
Formal class 'STFDF' [package "spacetime"] with 4 slots
..@ data :'data.frame': 253 obs. of 1 variable:
.. ..$ index: num [1:253] 95.3 92.5 83.2 89.8 92 ...
..@ sp :Formal class 'SpatialPoints' [package "sp"] with 3 slots
.. .. ..@ coords : num [1:23, 1:2] 1525044 1333239 1716737 2064910
1751258 ...
.. .. .. ..- attr(*, "dimnames")=List of 2
.. .. .. .. ..$ : NULL
.. .. .. .. ..$ : chr [1:2] "lon" "lat"
.. .. ..@ bbox : num [1:2, 1:2] 887073 6868534 2064910 7901351
.. .. .. ..- 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=aea +lat_1=27.5 +lat_2=35 +lat_0=18
+lon_0=-100 +x_0=1500000 +y_0=6000000 +ellps=GRS80 +datum=NAD83 +units=m
+no_defs +to"| __truncated__
..@ time :An xts object on 2001-01-01/2011-01-01 containing:
Data: int [1:11, 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:11], format: "2002-01-01" "2003-01-01" ...
> vv <- variogramST(index~1,city1b.st)
|
| | 0%
|
|==== | 6%
|
|========= | 12%
|
|============= | 19%
|
|================== | 25%
|
|====================== | 31%
|
|========================== | 38%
|
|=============================== | 44%
|
|=================================== | 50%
|
|======================================= | 56%
|
|============================================ | 62%
|
|================================================ |
69%Error in `[.xts`(x@time, j) : subscript out of bounds
In addition: Warning message:
In variogramST(index ~ 1, city1b.st) :
strictly irregular time steps were assumed to be regular
>
Does this look familiar, please?
Thanks for any help.
Sincerely,
Erin
[[alternative HTML version deleted]]
_______________________________________________
R-sig-Geo mailing list
R-sig-Geo@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
