y11 <- c(18.3, 19.5, 20.1, 17.4, 20.5, 16.1)
y12 <- c(19.2, 16.4, 16.5, 14.0, 16.9, 18.3)
y21 <- c(19.1, 18.4, 19.8, 20.0, 17.2, 16.8)
y22 <- c(16.5, 10.2, 12.7, 9.9, 13.5, 15.0)
ci.2x2.median.mixed(.05, y11, y12, y21, y22)
# Should return:
# Estimate SE LL UL
# AB: -3.450 1.6317863 -6.6482423 -0.2517577
# A: 1.875 0.8158931 0.2758788 3.4741212
# B: 3.925 1.4262367 1.1296274 6.7203726
# A at b1: 0.150 1.4243192 -2.6416144 2.9416144
# A at b2: 3.600 0.7962670 2.0393454 5.1606546
# B at a1: 2.200 1.5812792 -0.8992503 5.2992503
# B at a2: 5.650 1.7027101 2.3127496 8.9872504
Run the code above in your browser using DataLab