y11 <- c(14, 15, 11, 7, 16, 12, 15, 16, 10, 9)
y12 <- c(18, 24, 14, 18, 22, 21, 16, 17, 14, 13)
y21 <- c(16, 11, 10, 17, 13, 18, 12, 16, 6, 15)
y22 <- c(18, 17, 11, 9, 9, 13, 18, 15, 14, 11)
ci.2x2.mean.bs(.05, y11, y12, y21, y22)
# Should return:
# Estimate SE t df p LL UL
# AB: -5.10 2.224860 -2.2923 35.48 0.02793 -9.6145264 -0.5854736
# A: 1.65 1.112430 1.4832 35.48 0.14684 -0.6072632 3.9072632
# B: -2.65 1.112430 -2.3822 35.48 0.02270 -4.9072632 -0.3927368
# A at b1: -0.90 1.545244 -0.5824 17.56 0.56768 -4.1522367 2.3522367
# A at b2: 4.20 1.600694 2.6239 17.94 0.01725 0.8362274 7.5637726
# B at a1: -5.20 1.536952 -3.3833 17.61 0.00339 -8.4341379 -1.9658621
# B at a2: -0.10 1.608657 -0.0622 17.92 0.95112 -3.4807927 3.2807927
Run the code above in your browser using DataLab