Description Usage Format Source Examples

The data arose from a single hospital (at Campbelltown) as part of a larger (ongoing) study into the relationship between atmospheric pollution and the number of asthma cases presenting themselves to various emergency departments in local hospitals in the South West region of Sydney, Australia.

1 |

A data frame containing the following columns:

[, 1] | Count | Daily counts of asthma at Campbelltown Hospital. |

[, 2] | Intercept | A vector of ones, providing the intercept in the model. |

[, 3] | Sunday | Takes value one for Sundays, otherwise zero. |

[, 4] | Monday | Takes value one for Mondays, otherwise zero. |

[, 5] | CosAnnual | cos((2*pi*t)/365), annual cosine term. |

[, 6] | SinAnnual | sin((2*pi*t)/365), annual sine term. |

[, 7] | H7 | Scaled lagged and smoothed humidity variable. |

[, 8] | NO2max | Maximum daily nitrogen dioxide. |

[, 9:16] | T1.1990 - T2.1993 | Smooth shapes to capture school terms in each year. |

Davis, Richard A and Dunsmuir, William TM and Streett, Sarah B (2003)
Observation-driven models for Poisson counts. *Biometrika*,
**90**, 777–790.

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 | ```
### Example with asthma data
data(Asthma)
y <- Asthma[,1]
X <- as.matrix(Asthma[,2:16])
## Model in Davis, Dunsmuir and Streett (2003)
## MA(7) specification - see Davis, Dunsmuir and Streett (2003)
## Pearson Residuals, Fisher Scoring
glarmamod <- glarma(y, X, thetaLags = 7, type = "Poi", method = "FS",
residuals = "Pearson", maxit = 100, grad = 1e-6)
glarmamod
summary(glarmamod)
likTests(glarmamod)
plot.glarma(glarmamod)
## Not run:
## Example is specified as \dontrun because it takes too long
## for package inclusion on CRAN
## Pearson Residuals, Newton Raphson, Negative Binomial
## Initial value of the shape parameter take to be zero
glarmamod <- glarma(y, X, thetaLags = 7, type = "NegBin", method = "NR",
residuals = "Pearson", alphaInit = 0,
maxit = 100, grad = 1e-6)
glarmamod
summary(glarmamod)
likTests(glarmamod)
plot.glarma(glarmamod)
## End(Not run)
``` |

glarma documentation built on Feb. 9, 2018, 6:08 a.m.

Embedding an R snippet on your website

Add the following code to your website.

For more information on customizing the embed code, read Embedding Snippets.