Thanks.
In case it might help, I'm putting here a rude but hopefully correct LOO
sample code for STFDF objects:

-------------------------------------
N_POINTS <- nrow(mySTFDF[,<temporal-selection>])
krigeST.predictions <- array(dim=N_POINTS)

for (point in 1:N_POINTS) {
    # Remove one station
    mySTFDF.lessone <- mySTFD[c(sequence(point-1), seq(min(point
+1,N_POINTS),N_POINTS)),]
    mySTFDF.one <- SpatialPointsDataFrame(
                       SpatialPoints(matrix(coordinates(mySTFDF)[point
,],ncol=2),
                                     proj4string=CRS(PROJ.STRING)),
                       as.data.frame(mySTFDF[point,pred.day]))
    stked.predictions[point] <- krigeST(
                                    formula=<formula>,
                                    data=mySTFDF.lessone,

newdata=STFDF(as(mySTFDF.one,"SpatialPoints"),
                                                  stPM.month.proj@time
[<temporal-selection],
                                                  mySTFDF.one@data,
                                                  endTime=...),
                                    modelList=list(space=...,time=...),
                                    computeVar=FALSE)$pred
}
krigeST.residuals <- krigeST.predictions - mySTFDF
[,<temporal-selection>,<variable-selection>]@data
krigeST.rmse <- round(sqrt(mean(krigeST.residuals^2, na.rm=TRUE)),3)
-------------------------------------

Piero

On 13 January 2013 18:20, Edzer Pebesma <edzer.pebesma@uni-muenster.de>wrote:

>
>
> On 01/13/2013 01:53 PM, Piero Campalani wrote:
> > Dear list,
> >
> > I was wondering if there is a way to automatically compute x-validation
> > with the new spacetime objects (as is in gstat::gstat.cv).
> > Looking through the available functions it does not look so, thus I bet I
> > just need to manually run it. ?
> >
> I guess so, yes - this is what we've done so far.
>
> > Thanks,
> > Piero
> >
> >       [[alternative HTML version deleted]]
> >
> > _______________________________________________
> > R-sig-Geo mailing list
> > R-sig-Geo@r-project.org
> > https://stat.ethz.ch/mailman/listinfo/r-sig-geo
> >
>
> --
> Edzer Pebesma
> Institute for Geoinformatics (ifgi), University of Münster
> Weseler Straße 253, 48151 Münster, Germany. Phone: +49 251
> 8333081, Fax: +49 251 8339763  http://ifgi.uni-muenster.de
> http://www.52north.org/geostatistics      e.pebesma@wwu.de
>
> _______________________________________________
> R-sig-Geo mailing list
> R-sig-Geo@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>

	[[alternative HTML version deleted]]

