r - spacetime::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.

何か案は？