y11 <- c(21, 39, 32, 29, 27, 17, 27, 21, 28, 17, 12, 27)
y12 <- c(20, 36, 33, 27, 28, 14, 30, 20, 27, 15, 11, 22)
y21 <- c(21, 36, 30, 27, 28, 15, 27, 18, 29, 16, 11, 22)
y22 <- c(18, 34, 29, 28, 28, 17, 27, 21, 26, 16, 14, 23)
ci.2x2.stdmean.ws(.05, y11, y12, y21, y22)
# Should return:
# Estimate adj Estimate SE LL UL
# AB: 0.1725 0.1645 0.13655 -0.0951 0.4401
# A: 0.1092 0.1042 0.05753 -0.0035 0.2220
# B: 0.0747 0.0713 0.05921 -0.0413 0.1908
# A at b1: 0.1955 0.1864 0.08461 0.0297 0.3613
# A at b2: 0.0230 0.0219 0.09372 -0.1607 0.2067
# B at a1: 0.1610 0.1535 0.09457 -0.0244 0.3463
# B at a2: -0.0115 -0.0110 0.08596 -0.1800 0.1570
Run the code above in your browser using DataLab