SpatioTemporal (version 1.1.2)

predictNaive: Naive Temporal Predictions

Description

Computes naive predictions that are based on a few sites. These predictions can then be used, e.g. in summary.predCVSTmodel, to evaluate how much better the spatial-temporal model performs compared to simple (temporal) predictions. The function requires one of location and type to be specified, if both are given location will be used over type. If type is given locations such that as.character(STmodel$locations$type) type will be used.

Usage

predictNaive(STmodel, locations = NULL, type = NULL)

Arguments

STmodel

STmodel object for which to compute simple predictions.

locations

Locations on which to base the naive predictions.

type

The type of sites to base the predictions on, uses the (optional) field STmodel$locations$type.

Value

• A list with items:
• pred
  A (number of observations) - by - (6) data.frame containing the four naive predictions described under details, along with dates and IDs.
• locations
  The locations used for the naive predictions.

Details

Given a set of locations the function computes 4 sets of naive prediction for the observations in STmodel: [object Object],[object Object],[object Object],[object Object]

See Also

Other cross-validation functions: computeLTA, createCV, dropObservations, estimateCV.STmodel, predictCV.STmodel

Other STmodel functions: createCV, createDataMatrix, createSTmodel, dropObservations, loglikeST, loglikeSTdim, loglikeSTnaive, processLocation, processLUR, processST, updateCovf

Examples

##load data
data(mesa.model)

##naive predictions based on either AQS,
pred.aqs <- predictNaive(mesa.model, type="AQS")
##...FIXED sites,
pred.fixed <- predictNaive(mesa.model, type="FIXED")
##...or only one sites,
pred.1site <- predictNaive(mesa.model, locations="60372005")

##plot the predictions - The two cases that are constant in space
par(mfcol=c(3,1), mar=c(4.5,4.5,1,.5))
##observations as a function of date
plot(mesa.model, "loc.obs", type=as.factor(mesa.model$locations$ID),
     legend.loc=NULL, pch=19, cex=.25)

##Add the predictions based on the smooth fitted to all sites
lines(pred.aqs$pred$date, pred.aqs$pred$smooth.fixed, col=1, lwd=2)
lines(pred.fixed$pred$date, pred.fixed$pred$smooth.fixed, col=2, lwd=2)
lines(pred.1site$pred$date, pred.1site$pred$smooth.fixed, col=3, lwd=2)


plot(mesa.model, "loc.obs", type=as.factor(mesa.model$locations$ID),
     legend.loc=NULL, pch=19, cex=.25)
##Add predictions based on the temporal average
lines(pred.aqs$pred$date, pred.aqs$pred$avg.fixed, col=1, lwd=2)
lines(pred.fixed$pred$date, pred.fixed$pred$avg.fixed, col=2, lwd=2)

##plot the predictions - One of the cases that vary in space, i.e. the smooth
##fit to the closest site.
##first extract as a data matrix
D <- createDataMatrix(obs=pred.aqs$pred$smooth.closest.fixed,
                      date=pred.aqs$pred$date, ID=pred.aqs$pred$ID)

##observations as a function of date
##(only five sites for clarity)
mesa.model <- dropObservations(mesa.model, !(mesa.model$obs$idx %in% c(1:5)))
plot(mesa.model, "loc.obs", type=as.factor(mesa.model$locations$ID),
     legend.loc=NULL, pch=19, cex=.25)
##Add the predictions based on the smooth
##fitted to the closest site
for(i in 1:5){
  lines(as.Date(rownames(D)), D[,mesa.model$locations$ID[i]], col=i, lwd=2)
}