data(pork) # pork tasting data, 10 individual paired comparison matrices
btl1 <- eba(apply(pork, 1:2, sum)) # fit a BTL model
b <- boot(pork, 200) # resample 200 times
plot(btl1$estimate, b$stat[,'mean'], log='xy')
abline(0, 1, lty=2)
Run the code above in your browser using DataLab