2. model code in R.zip