To do so use the following commands. > j<-c(4.1, 3.2, 4.0, 5.1, 2.3, 2.9, 3.2, 3.9, 4.0, 4.2, 4.9, 4.8) > j2<-randtestoneway(j,c(4,4,4),5000) You should see something similar to the following printed out on your screen (the actual test statistic value should be exactly the same, but the pvalue will dier slightly because you are using randomly sampled arrangements not the same ones that I got). Run this same command several times to see how variable the reported pvalues are. How many possible treatment arrangements are there for this experiment?