exp4_code.R