In our simulations, we use this (in Matlab/Octave):
function l=linearityindex(matrix)
m=(matrix>matrix')+.5*(matrix==matrix');
rowsums=sum(m,2)-.5;
n=length(rowsums);
l=12/(n^3-n)*sum((rowsums-(n-1)/2).^2);
You may want to check out Landau's 1951 paper, "On dominance relations and the structure of animal societies: I. Effect of inherent characteristics".
From http://www.springerlink.com/content/366920782vwr0648/
Chat with our AI personalities