subtract marginal (spatial and temporal) means from observations
subtract marginal (spatial and temporal) means from observations
stInteraction(x, ...)
object of class STFDF with each attribute replaced by its residual, computed by $y_ij=x_ij-m_.jmi./m$ with $m$ the grand mean, $m_.j$ the temporal mean, $m_i.$ the spatial mean and $m$ the grand mean.
if (require(gstat)) {
library(sp)
data(wind)
wind.loc$y = as.numeric(char2dms(as.character(wind.loc[["Latitude"]])))
wind.loc$x = as.numeric(char2dms(as.character(wind.loc[["Longitude"]])))
coordinates(wind.loc) = ~x+y
proj4string(wind.loc) = "+proj=longlat +datum=WGS84"
# match station order to names in wide table:
stations = 4:15
wind.loc = wind.loc[match(names(wind[stations]), wind.loc$Code),]
row.names(wind.loc) = wind.loc$Station
wind$time = ISOdate(wind$year+1900, wind$month, wind$day, 0)
space = list(values = names(wind)[stations])
wind.st = stConstruct(wind[stations], space, wind$time, SpatialObj = wind.loc)
wind.sti = stInteraction(wind.st)
# temporal means for any station should be zero:
c(mean(wind.sti[3,]),
# spatial mean for each time step should be zero:
mean(wind.sti[,5][[1]]))
}Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.