## Wilcoxon for 2x2x4. "Not many things easier than that"? [Nonparametrics]

Dear Smart People!
Could you please advise me how to extend Wilcoxon test for 4-period replicate study properly?
1. Do we need to use means of two periods of T's and R's? If yes, why do we need arithmetic means, why don't we need neither medians nor harmonic means nor so on?
2. How does Phoenix perform this calculation in Crossover object? According to guideline it is used only for 2x2 design, but it also works for 2x2x4 somehow.
Below is my suggested code (first two periods coincide with the example from the post). For this data Phoenix gave p=0.6863 - does it make any sense?
 data<-c("Subject Sequence P1 P2 P3 P4 1 TRTR 2.0 2.0 3.0 3.0 2 RTRT 1.5 3.0 1.5 3.0 3 TRTR 2.0 2.0 2.0 3.0 4 RTRT 2.0 2.0 1.5 2.0 5 TRTR 3.0 2.0 2.0 1.5 6 RTRT 3.0 2.0 2.0 3.0 7 TRTR 1.5 2.0 1.5 2.0 8 RTRT 3.0 2.0 2.0 2.0 9 TRTR 2.0 3.0 1.5 3.0 10 RTRT 2.0 1.5 2.0 2.0 11 TRTR 1.5 2.0 2.0 2.0 12 RTRT 2.0 3.0 3.0 2.0 13 TRTR 1.5 3.0 3.0 2.0 14 RTRT 3.0 3.0 1.5 2.0") tc <- textConnection(data) PKt  <- read.table(tc, header=TRUE, strip.white=TRUE, as.is=TRUE) close(tc)   PKt$pdiff <- PKt$P1+PKt$P3-PKt$P2-PKt$P4 library(coin) PKt$Sequence <- as.factor(PKt\$Sequence) t1 <- wilcox_test(pdiff ~ Sequence, data=PKt, conf.int=TRUE, conf.level=0.9, distribution="exact") print(t1) #p=0.4353 t2 <- wilcox_test(pdiff ~ Sequence, data=PKt, conf.int=TRUE, conf.level=0.9, distribution="asymptotic") print(t2) #p=0.3983 

"Being in minority, even a minority of one, did not make you mad"  Ing. Helmut Schütz 