Log Likelihood ratio em Octave

Calculando o log likelihood ratio de uma tabela de contingência (na matriz M) usando Octave:

function [LL, p] = logLikelihood(M)
  LL = zeros(size(M,1),1);
  EF = zeros(size(M,1), size(M,2));

  EF = [sum(M)']  * [sum(M,2)'] ./ sum(sum(M));
  LL = (2 * sum(M .* log(M./EF' + (M == 0)),2));
   p = 1 - chi2cdf(LL, (size(M,1)-1)*(size(M,2)-1));
end


Você pode achar interessante também essa calculadora online:
http://ucrel.lancs.ac.uk/llwizard.html