Binned data of number of abnormal rotation rates and number of flares is illustrated in Fig. 3a. The red continuous line is a linear least squares fit of the form nF_bin = (−1.489 ± 0.561) + (1.292 ± 0.151)nARR_bin (where nF_bin and the nARR_bin are the number of flares and the number of abnormal rotation rates respectively). Binned data of number of non-abnormal rotation rates and number of flares is illustrated in Fig. 3b. The red continuous line is a linear least squares fit of the form nF_bin = (1.136 ± 0.264) + (0.051 ± 0.048)nNARR_bin (where nF_bin and nNARR_bin are the number of flares and the number of binned non-abnormal rotation rates respectively). The correlation coefficient and its significance for the least square fits are 0.95 and 99.74% respectively for Fig. 3a and, 0.25 and 93.2% for Fig. 3b, respectively.