add R script to show how to prepare functional phenotype