exp1_code.R