Hi I have analysed a dataset in Graphpad prism with a 2 Way Anova.

1) In the experimental design i select "Each column represents a different time point, so matched values are spread accross the row."

2) Then in Multiple Comparisons tab we select "With each column, compare rows".

3) And then we correct for multiple comparisons using Bonferroni.

Eventually i get a p value for each row comparing different treatments.

I can not seem to be able to replicate same results or even same pattern programmatically.

So far I have tried pipelines:

http://jpktd.blogspot.com/2013/03/multiple-comparison-and-tukey-hsd-or_25.html

https://pythonfordatascience.org/anova-python/