All codes have been tested on local Windows machine runnign R version 4.0.2.

To rerun the codes 1-4, you have to modify the R working directory by the function setwd()

The code 5 has been tested on this local fold. 