[R] raster time series statistics

[Alexander.Herr at csiro.au]

Hi List,

The following code returns an "Error in as.POSIXlt.character(x, tz, ...) :   character string is not in a standard unambiguous format"

require(raster)
require(rts)
require(stringi)
r <- raster(ncol=100, nrow=100)
values(r) <- runif(ncell(r))
list(ID=seq(1:24),month=rep(str_pad(1:12, pad = 0,width = 2 , "left"),2),year=sort(rep(2016:2017,12)))->dt
stack(r)->s
r->rs
for(i in 1:23){
rs[]<-r[]*i
  addLayer(s,rs)->s
print(nlayers(s))
}
timelst<-paste0(unlist(dt['year']),'-',unlist(dt['month']))
rts(s,time=as.yearmon(timelst))->rsts
str(rsts at time)
apply.monthly(rsts,mean)

I was expecting that the statistics accept the yearmonth format of the timeslot, as it allows subsetting.
Any suggestions?
Thanks
Herry

	[[alternative HTML version deleted]]