pisa {svylme} | R Documentation |
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
Source
Data downloaded from https://www.stata-press.com/data/r15/pisa2000.dta
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)