To understand if the size of the systems affects the accuracy of the methods, we performed separate statistical analyses on systems with 1-to-4, and 5-6 non-hydrogen atoms. We collected the related indicators in Tables S16-S18 in the SI. From this analysis, we observe small variations in the trends for each method when moving from smaller to larger molecules. Interestingly, we observe how Q-CCSD, E-CCSD, Q-CC2, and I-ADC2 methods reduce their performances. For example the MAPE of Q-CCSD passes from 6.5% to 18.4%, of Q-CC2 from 40.4% to 45.1%, and of I-ADC2 from 46.1% to 60.8% when going from small to large molecules. In contrast, the performances of TD-DFT methods tend to improve, with 0 50 100 150 for instance the MAPE of CAM-B3LYP going from 38.7% to 32.8%. An exception to this trend is the LC-BLYP47 functional, which shows a MAPE of 36.8% for small molecules and 45.6% for the large ones.