pm10: Realistic data which illustrate the usage of the package Stem
In Stem: Spatio-temporal models in R

Description Usage Format Author(s) References Examples

This simple data set is a list of three objects and refers to 22 spatial locations and 366 time points.

1	data(pm10)

A list of with three objects with the following components:

coords: the coordinates of the 22 spatial locations.
covariates: it is a 8052 by 3 matrix referring to the following covariates: intercept, emissions (g/s) and altitude (km). The first 366 rows refer to the first spatial location, the rows from 367 to 732 refer to the second spatial location and so on.
z: it is a 366 by 22 observation matrix referring to PM10 concentration measurements (log scale).

Michela Cameletti michela.cameletti@unibg.it

Fasso', A., Cameletti, M. (2007) A general spatio-temporal model for environmental data. Tech.rep. n.27 Graspa - The Italian Group of Environmental Statistics - http://www.graspa.org.

data(pm10)
names(pm10)

#plot the coordinates
dim(pm10$coords)
plot(pm10$coords[,1],pm10$coords[,2],xlab=colnames(pm10$coords)[1],
	ylab=colnames(pm10$coords)[2])

#plot the data
dim(pm10$z)

#summary by station
apply(pm10$z,2,summary)

#plot the time series for station n.22
plot(pm10$z[,22],t="l",xlab="Days",ylab="PM10 concentrations (log)")

#plot the station altitude
plot(pm10$covariates[,3],ylab=colnames(pm10$covariates)[3],xaxt="n",xlab="")
positions = seq(1,8052,366)+366/2
axis(1, at=positions, labels=rownames(pm10$coords), las=2)

Loading required package: mvtnorm
Loading required package: MASS
[1] "coords"     "covariates" "z"         
[1] 22  2
[1] 366  22
           Stat1    Stat2    Stat3    Stat4    Stat5    Stat6    Stat7    Stat8
Min.    1.945910 1.609438 2.197225 1.386294 1.386294 1.945910 1.609438 2.302585
1st Qu. 3.146134 3.433987 3.135494 3.258097 2.833213 3.332205 3.222742 3.178054
Median  3.433987 3.761200 3.555348 3.663562 3.496508 3.737670 3.624252 3.526361
Mean    3.485402 3.752521 3.607629 3.632213 3.332381 3.736664 3.575426 3.532639
3rd Qu. 3.923922 4.128504 4.158883 4.077537 3.806662 4.158883 3.970292 3.926875
Max.    5.017280 5.176150 5.099866 5.327876 4.997212 5.247024 5.187386 5.049856
           Stat9   Stat10   Stat11   Stat12   Stat13   Stat14   Stat15   Stat16
Min.    1.386294 1.098612 1.609438 2.079442 1.386294 1.386294 1.945910 1.609438
1st Qu. 3.265650 3.091042 2.890372 3.434967 3.295837 3.091042 2.833213 2.944439
Median  3.684844 3.465736 3.367296 3.806662 3.610918 3.465736 3.258097 3.295837
Mean    3.639113 3.476453 3.305686 3.801897 3.658579 3.418776 3.216801 3.326619
3rd Qu. 4.002746 3.931826 3.731645 4.230457 4.002746 3.850148 3.555348 3.711628
Max.    5.198497 4.919981 5.056246 5.176150 5.187386 5.141664 4.727388 4.744932
          Stat17   Stat18   Stat19   Stat20   Stat21   Stat22
Min.    1.386294 2.197225 2.079442 2.197225 2.197225 1.609438
1st Qu. 2.890372 3.806662 3.433987 3.465736 3.375771 2.772589
Median  3.218876 4.195758 3.737670 3.912023 3.828641 3.198465
Mean    3.219606 4.127616 3.743610 3.900893 3.826566 3.180278
3rd Qu. 3.610918 4.509436 4.110874 4.369448 4.262680 3.610918
Max.    4.969813 5.389072 5.225747 5.252273 5.192957 4.983607