svylme (version 1.5-1)

pisa: Data from the PISA international school survey

Description

Data from the PISA survey of schools, obtained from Stata, who obtained it from Rabe-Hesketh & Skrondal.

Usage

data("pisa")

Format

A data frame with 2069 observations on the following 11 variables.

female

1 for female

isei

socioeconomic index

w_fstuwt

student sampling weight (total)

wnrschbw

school sampling weight

high_school

1 if highest level of parents' education is high school

college

1 if highest level of parents' education is college/uni

one_for

1 if one parent is foreign-born

both_for

1 if both parents are foreign-born

test_lang

1 if the test language is spoken at home

pass_read

1 if the student passed a reading proficiency test

id_school

school (sampling unit) identifier

References

Rabe-Hesketh, S., and A. Skrondal. 2006. Multilevel modelling of complex survey data.Journal of the Royal Statistical Society, Series A. 169: 805-827

Examples

data(pisa)

## This model doesn't make a lot of sense, but it's the one in the
## Stata documentation because the outcome variable is numeric.

pisa$w_condstuwt<-with(pisa, w_fstuwt/wnrschbw)
pisa$id_student<-1:nrow(pisa)

dpisa<-survey::svydesign(id=~id_school+id_student, weight=~wnrschbw+w_condstuwt, data=pisa)


svy2lme(isei~(1|id_school)+female+high_school+college+one_for+both_for+test_lang,
	design=dpisa)