In eldafani/intsvy: International Assessment Data Manager

knitr::opts_chunk$set(echo = TRUE)

library("tidyverse")
library("haven")

library("intsvy")
options(width=120)

Read data

dir <- "/home/eldani/eldani/International LSA/SEA-PLM/Data/SEA-PLM_Regional_Students-sav"
student <- read_spss(file.path(dir, "SEA-PLM_Regional_Students.sav"))
dir <- "/home/eldani/eldani/International LSA/SEA-PLM/Data/SEA-PLM_Regional_Schools-sav"
school <- read_spss(file.path(dir, "SEA-PLM_Regional_Schools.sav"))

Merge data

sea <- left_join(student, school, by = c("CNT", "SchID"))

Create grouping variables

sea$Sex <- factor(sea$Gender, 1:2, labels = c("Female", "Male"))

sea$Language <- factor(sea$S_LANG, 0:1, labels = c("No", "Yes"))

sea$Location <- factor(sea$SC09Q01, 1:5, labels = c("Rural", rep("Urban", 4)))

Proficiency reading levels by CNT

intsvy.ben.pv(pvnames= paste0("PV", 1:5, "_R"), cutoff = c(304, 317),
              by= c("CNT"), data=sea, config = sea_conf)

Proficiency reading levels by CNT and sex

intsvy.ben.pv(pvnames= paste0("PV", 1:5, "_R"), cutoff =   c(304, 317), 
              by= c("CNT", "Sex"), data=sea, config = sea_conf)

Mean reading achievement by groups

intsvy.mean.pv(pvnames=paste0("PV", 1:5, "_R"), by= c("CNT"), data=sea, config = sea_conf)

intsvy.mean.pv(pvnames=paste0("PV", 1:5, "_R"), by= c("CNT", "Location"), data=sea, config = sea_conf)

intsvy.mean.pv(pvnames=paste0("PV", 1:5, "_R"), by= c("CNT", "Sex"), 
               data=sea[!is.na(sea$Sex), ], config = sea_conf)

intsvy.mean.pv(pvnames= paste0("PV", 1:5, "_R"), by= c("CNT", "Language"), data=sea, config = sea_conf)

Mean writing and math achievement by CNT

intsvy.mean.pv(pvnames=paste0("PV", 1:5, "_W"), by= c("CNT"), data=sea, config = sea_conf)
intsvy.mean.pv(pvnames=paste0("PV", 1:5, "_M"), by= c("CNT"), data=sea, config = sea_conf)

Regression of reading achievement on SES

intsvy.reg.pv(pvnames= paste0("PV", 1:5, "_R"), x= "SES", by= c("CNT"), data=sea, config = sea_conf)

Regression of math achievement on SES and sex

intsvy.reg.pv(pvnames= paste0("PV", 1:5, "_M"), x= c("SES", "Sex"), by= c("CNT"), data=sea, config = sea_conf)

eldafani/intsvy documentation built on March 18, 2024, 7:54 p.m.