The goal is to create a function that can perform the Critchlow-Fligner simultaneous confidence intervals for pairwise comparisons (as in section 6.9 of the textbook).

Looking for a function that can

• Accept data as data.frame or as a matrix where each treatment corresponds to a column.
• Other functional arguments should be "alpha" and whether it is "two sided" or "one sided".
• Create all pairwise comparisons and output the results in a nice coherent way.

Feel free to look at other functions from NSM3 or other packages, to get ideas.

(Optional) Other considerations are:

• For a large data set (decide on what "large" is) should the function use the large sample approximation