File 3: Statistical analysis R scripts