R_script_for_VPT