Learn R Programming

astsa (version 2.2)

lap.xts: LA Pollution-Mortality Study: Sampled Daily

Description

Original data from a study of the effects of pollution and weather on mortality, LA, 1970-1979. These are 3652 daily observations for the 10 year period. The data set is an xts object indexed by Date.

Arguments

Format

The format is: An xts object on 1970-01-01 / 1979-12-31 containing:
Data: double [3652, 11]
Columns: Tmort, Rmort, Cmort, Temp, Rhumid, CO, SO2, NO2, HC, Ozone, Part
Index: Date [3652] (TZ: "UTC")

Details

columns are time serieswith names
(1) Total MortalityTmort
(2) Respiratory MortalityRmort
(3) Cardiovascular MortalityCmort
(4) TemperatureTemp
(5) Relative HumidityRhumid
(6) Carbon MonoxideCO
(7) Sulfur DioxideSO2
(8) Nitrogen DioxideNO2
(9) HydrocarbonsHC
(10) OzoneOzone
(11) ParticulatesPart

References

You can find demonstrations of astsa capabilities at FUN WITH ASTSA.

The most recent version of the package can be found at https://github.com/nickpoison/astsa/.

In addition, the News and ChangeLog files are at https://github.com/nickpoison/astsa/blob/master/NEWS.md.

The webpages for the texts and some help on using R for time series analysis can be found at https://nickpoison.github.io/.

See Also

lap

Examples

Run this code
if (FALSE) {

library(xts)   # assumes package has been installed

plot(lap.xts$Cmort, col=4)   

lapw = apply.weekly(lap.xts, FUN=colMeans)  # get weekly averages

plot(lapw[,c('Cmort', 'Temp', 'Part')], col=astsa.col(2:4, .7), main=NA) 
addLegend(col=2:4, lty=1, lwd=2, ncol=3, bty="white")

sarima(lapw$Cmort, 0,1,1, no.constant=TRUE)  # fit ARIMA(0,1,1) to weekly Cmort
}

Run the code above in your browser using DataLab