lips1 | R Documentation |
Observed counts of male lip cancer for the 56 regions of Scotland over the period 1975-1980.
lips1
A data frame with 56 observations on the following 41 variables:
Region ID.
Constant (=1).
Observed cases of lip cancer.
Expected count.
Percentage of the region who work in agriculture, fishing and forestry.
Log of the expected count.
Constant (=1).
Constant (=1).
First neighbours.
Second neighbours.
Third neighbours.
Fourth neighbours.
Fifth neighbours.
Sixth neighbours.
Seventh neighbours.
Eighth neighbours.
Ninth neighbours.
Tenth neightbours.
Eleventh neightbours.
First neighbours' weights.
Second neighbours' weights.
Third neighbours' weights.
Fourth neighbours' weights.
Fifth neighbours' weights.
Sixth neighbours' weights.
Seventh neighbours' weights.
Eighth neighbours' weights.
Ninth neighbours' weights.
Tenth neightbours' weights.
Eleventh neightbours' weights.
First neighbours' CAR weights.
Second neighbours' CAR weights.
Third neighbours' CAR weights.
Fourth neighbours' CAR weights.
Fifth neighbours' CAR weights.
Sixth neighbours' CAR weights.
Seventh neighbours' CAR weights.
Eighth neighbours' CAR weights.
Ninth neighbours' CAR weights.
Tenth neightbours' CAR weights.
Eleventh neightbours' CAR weights.
The lips1
dataset is one of the sample datasets provided with the
multilevel-modelling software package MLwiN (Rasbash et al., 2009), and was
analysed in Clayton & Kaldor (1987); see also Browne (2012) for more
details.
Browne, W. J. (2012) MCMC Estimation in MLwiN Version 2.26. University of Bristol: Centre for Multilevel Modelling.
Clayton, D., Kaldor, J. (1987) Empirical Bayes estimates of age-standardized relative risks for use in disease mapping. Biometrics 43: 671-681.
Rasbash, J., Charlton, C., Browne, W.J., Healy, M. and Cameron, B. (2009) MLwiN Version 2.1. Centre for Multilevel Modelling, University of Bristol.
## Not run:
data(lips1, package = "R2MLwiN")
(mymodel <- runMLwiN(log(obs) ~ 1 + perc_aff + offset(offs) + (0 | neigh1) + (1 | area),
D = "Poisson", estoptions = list(EstM = 1), data = lips1))
## End(Not run)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.