Create JM_analysis.R