R2MLwiN (version 0.8-6)

lips1: Lips data

Description

Observed counts of male lip cancer for the 56 regions of Scotland over the period 1975-1980.

Format

A data frame with 56 observations on the following 41 variables:

area

Region ID.

cons

Constant (=1).

obs

Observed cases of lip cancer.

exp

Expected count.

perc_aff

Percentage of the region who work in agriculture, fishing and forestry.

offs

Log of the expected count.

pcons

Constant (=1).

denom

Constant (=1).

neigh1

First neighbours.

neigh2

Second neighbours.

neigh3

Third neighbours.

neigh4

Fourth neighbours.

neigh5

Fifth neighbours.

neigh6

Sixth neighbours.

neigh7

Seventh neighbours.

neigh8

Eighth neighbours.

neigh9

Ninth neighbours.

neigh10

Tenth neightbours.

neigh11

Eleventh neightbours.

weight1

First neighbours' weights.

weight2

Second neighbours' weights.

weight3

Third neighbours' weights.

weight4

Fourth neighbours' weights.

weight5

Fifth neighbours' weights.

weight6

Sixth neighbours' weights.

weight7

Seventh neighbours' weights.

weight8

Eighth neighbours' weights.

weight9

Ninth neighbours' weights.

weight10

Tenth neightbours' weights.

weight11

Eleventh neightbours' weights.

wcar1

First neighbours' CAR weights.

wcar2

Second neighbours' CAR weights.

wcar3

Third neighbours' CAR weights.

wcar4

Fourth neighbours' CAR weights.

wcar5

Fifth neighbours' CAR weights.

wcar6

Sixth neighbours' CAR weights.

wcar7

Seventh neighbours' CAR weights.

wcar8

Eighth neighbours' CAR weights.

wcar9

Ninth neighbours' CAR weights.

wcar10

Tenth neightbours' CAR weights.

wcar11

Eleventh neightbours' CAR weights.

Details

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.

Examples

# 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))

# }