Martin Maechler has pointed out to me that I omitted to include my values for x when circulating my query. They are a simple rank list from 1 to 26: x <- 1:26 Thanks, Andrew