final script and Rmarkdown