# Loading of the verbal data
data(verbal)
# Computing proportions of success
props <- cbind(colMeans(verbal[verbal[,26]==0,1:24]),
colMeans(verbal[verbal[,26]==1,1:24]) )
# Perpendicular distances
trItemDiff(props)
# Removing item 6 from the set of anchor items
trItemDiff(props, anchor=c(1:5,7:24))Run the code above in your browser using DataLab