Create joint.R - Script to analyze 40k patients