Update PL_analysis.R