Samyak

Hello! If I am not mistaken the MultiComparison allows tukey test of several factors easily `import statsmodels.stats.multicomp as mc` `tukey=mc.MultiComparison(df['value'],df['Genotype']+df['years'])` `tukey_test=tuk.tukeyhsd()` I believe this will provide you with the solution