# NOT RUN {
## Data as a square table, as in Agresti p448
citations
##
## Convert frequencies to success/failure data:
##
citations.sf <- countsToBinomial(citations)
names(citations.sf)[1:2] <- c("journal1", "journal2")
## Standard Bradley-Terry model fitted to these data
citeModel <- BTm(cbind(win1, win2), journal1, journal2,
data = citations.sf)
# }
Run the code above in your browser using DataLab