230 likes | 240 Views
This lecture covers cross-classified models, multiple membership models, and examples such as Danish chickens and educational data from ALSPAC. The lecture also includes an AI example and slides by Jon Rasbash.
E N D
Lecture 21 Cross-classified and Multiple membership models
Lecture Contents • Cross classified models • AI example • Multiple membership models • Danish chickens example • More complex structures • ALSPAC educational example Thanks to Jon Rasbash for slides!
Hospital H1 H2 H3 H4 Patient P1 P2 P3 P4 P5 P6 P7 P8 P9 P10 P11 P12 Nbhd N1 N2 N3 Cross-classification For example, hospitals by neighbourhoods. Hospitals will draw patients from many different neighbourhoods and the inhabitants of a neighbourhood will go to many hospitals. No pure hierarchy can be found and patients are said to be contained within a cross-classification of hospitals by neighbourhoods :
Other examples of cross-classifications • pupils within primary schools by secondary schools. • patients within GPs by hospitals. • interviewees within interviewers by surveys. • repeated measures within raters by individual. (e.g. patients by nurses)
Notation With hierarchical models we use a subscript notation that has one subscript per level and nesting is implied reading from the left. For example, subscript pattern ijk denotes the i’th level 1 unit within the j’th level 2 unit within the k’th level 3 unit. If models become cross-classified we use the term classification instead of level. With notation that has one subscript per classification, that captures the relationship between classifications, notation can become very cumbersome. We propose an alternative notation introduced in Browne et al. (2001) that only has a single subscript no matter how many classifications are in the model.
i nbhd(i) hosp(i) 1 1 1 2 2 1 3 1 1 4 2 2 5 1 2 6 2 2 Hospital H1 H2 H3 H4 7 2 3 Patient P1 P2 P3 P4 P5 P6 P7 P8 P9 P10 P11 P12 8 3 3 Nbhd N1 N2 N3 9 3 4 10 2 4 11 3 4 12 3 4 Single subscript notation We write the model as Where classification 2 is neighbourhood and classification 3 is hospital. Classification 1 always corresponds to the classification at which the response measurements are made, in this case patients. For patients 1 and 11 equation (1) becomes:
Neighbourhood Hospital Patient Classification diagrams In the single subscript notation we lose information about the relationship(crossed or nested) between classifications. A useful way of conveying this information is with the classification diagram. Which has one node per classification and nodes linked by arrows have a nested relationship and unlinked nodes have a crossed relationship. Hospital Neighbourhood Patient Cross-classified structure where patients from a hospital come from many neighbourhoods and people from a neighbourhood attend several hospitals. Nested structure where hospitals are contained within neighbourhoods
Donor Donation Woman Cycle Example : Artificial insemination by donor 1901 women 279 donors 1328 donations 12100 ovulatory cycles response is whether conception occurs in a given cycle In terms of a unit diagram: Or a classification diagram:
Parameter Description Estimate(se) intercept -4.04(2.30) azoospermia * 0.22(0.11) semen quality 0.19(0.03) womens age>35 -0.30(0.14) sperm count 0.20(0.07) sperm motility 0.02(0.06) insemination to early -0.72(0.19) insemination to late -0.27(0.10) women variance 1.02(0.21) donation variance 0.644(0.21) donor variance 0.338(0.07) Model for artificial insemination data We can write the model as Results:
Multiple membership models • When level 1 units are members of more than one higher level unit we describe a model for such data as a multiple membership model. • For example, • Pupils change schools/classes and each school/class has an effect on pupil outcomes. • Patients are seen by more than one nurse during the course of their treatment.
n1(j=1) n2(j=2) n3(j=3) p1(i=1) 0.5 0 0.5 p2(i=2) 1 0 0 p3(i=3) 0 0.5 0.5 p4(i=4) 0.5 0.5 0 Here patient 1 was seen by nurse 1 and 3 but not nurse 2 and so on. If we substitute the values of w(2)i,j, i and j. from the table into (1) we get the series of equations : Notation Note that nurse(i) now indexes the set of nurses that treat patient i and w(2)i,jis a weighting factor relating patient i to nurse j. For example, with four patients and three nurses, we may have the following weights:
nurse Here patients are multiple members of nurses, nurses are nested within hospitals and GP practice is crossed with both nurse and hospital. patient hospital nurse GP practice patient Classification diagrams for multiple membership relationships Double arrows indicate a multiple membership relationship between classifications. We can mix multiple membership, crossed and hierarchical structures in a single model.
Farm House Parent flock Child flock Example involving nesting, crossing and multiple membership – Danish chickens Production hierarchy 10,127 child flocks 725 houses 304 farms Breeding hierarchy 10,127 child flocks 200 parent flocks As a unit diagram: As a classification diagram:
Parameter Description Estimate(se) intercept -2.322(0.213) 1996 -1.239(0.162) 1997 -1.165(0.187) Results: hatchery 2 -1.733(0.255) hatchery 3 -0.211(0.252) hatchery 4 -1.062(0.388) parent flock variance 0.895(0.179) house variance 0.208(0.108) farm variance 0.927(0.197) Model and results
ALSPAC data All the children born in the Avon area in 1990 followed up longitudinally. Many measurements made including educational attainment measures. Children span 3 school year cohorts(say 1994,1995,1996). Suppose we wish to model development of numeracy over the schooling period. We may have the following attainment measures on a child : m1 m2 m3 m4 m5 m6 m7 m8 primary school secondary school
Structure for primary schools Primary school P School Cohort Area Pupil P. Teacher M. Occasion • Measurement occasions within pupils. • At each occasion there may be a different teacher. • Pupils are nested within primary school cohorts. • All this structure is nested within primary school. • Pupils are nested within residential areas.
A mixture of nested and crossed relationships Primary school P School Cohort Area Pupil P. Teacher M. occasions Nodes directly connected by a single arrow are nested, otherwise nodes are cross-classified. For example, measurement occasions are nested within pupils. However, cohort are cross-classified with primary teachers, that is teachers teach more than one cohort and a cohort is taught by more than one teacher.
Multiple membership Primary school P School Cohort Area Pupil P. Teacher m1 m2 m3 m4 M. occasions t1 t2 t3 t4 It is reasonable to suppose the attainment of a child in a particualr year is influenced not only by the current teacher, but also by teachers in previous years. That is measurements occasions are “multiple members” of teachers. We represent this in the classification diagram by using a double arrow.
What happens if pupils move area? Primary school Area P School Cohort P. Teacher Primary school Pupil P School Cohort P. Teacher M. occasions Area Pupil M. occasions Classification diagram without pupils moving residential areas. If pupils move area, then pupils are no longer nested within areas. Pupils and areas are cross-classified. Also it is reasonable to suppose that pupils measured attainments are effected by the areas they have previously lived in. So measurement occasions are multiple members of areas. Classification diagram where pupils move between residential areas. BUT…
If pupils move area they will also move schools Primary school Primary school P School Cohort P. Teacher P School Cohort P. Teacher Area Pupil Area Pupil M. occasions M. occasions Classification diagram where pupils move between areas but not schools. If pupils move schools they are no longer nested within primary school or primary school cohort. Also we can expect, for the mobile pupils, both their previous and current cohort and school to effect measured attainments. Classification diagram where pupils move between schools and areas.
If pupils move area they will also move schools cnt’d Primary school P School Cohort P. Teacher Area Pupil M. occasions And secondary schools… We could also extend the above model to take account of Secondary school, secondary school cohort and secondary school teachers.
Other predictor variables Remember we are partitioning the variability in attainment over time between primary school, residential area, pupil, p. school cohort, teacher and occasion. We also have predictor variables for these classifications, eg pupil social class, teacher training, school budget and so on. We can introduce these predictor variables to see to what extent they explain the partitioned variability.
Information for the practicals • We have two MLwiN practicals taken from chapters of Browne (2003). We firstly look at a cross-classified model for education data (primary schools and secondary schools. We secondly look at a multiple membership model for a (simulated) earnings dataset.