按时间设置空时：：SDFDF

Question

spacetime::SDFDF 如何按时间划分(具有全时空网格的时空数据)？

索瓦尔，我试过：

library("maps")
library("maptools")
library("spacetime")
library("plm")

states.m <- map("state", plot = FALSE, fill = TRUE)
IDs <- sapply(strsplit(states.m$names, ":"), function(x) x[1])
states <- map2SpatialPolygons(states.m, IDs = IDs)
yrs <- 1970:1986
time <- as.POSIXct(paste(yrs, "-01-01", sep = ""), tz = "GMT")
data("Produc")

Produc.st <- STFDF(states[-8], time, Produc[order(Produc[2], Produc[1]),])
Produc.st@time[c(1,5,17)]
Produc.st[Produc.st@time[c(1,5,17)]]

但这给了我error: ncol(i) == 2 is not TRUE。

有什么想法吗？

Answer 1

请尝试

Produc.st[,index(Produc.st@time[c(1,5,17)])]

也就是说，时间选择是在,之后完成的，不要像xts对象那样选择Produc.st@time[c(1,5,17)])，而是使用时间(POSIXct)向量进行选择。