MyFisher24

A very compact routine to compute Fisher's exact test on a 2x4 matrix

Description

Fisher's exact test of 2x4 contingency tables permits calculation of precise probabilities in situation where, as a consequence of small cell frequencies, the much more rapid normal approximation and chi-square calculations are liable to be inaccurate. The Fisher's exact test involves the computations of several factorials to obtain the probability of the observed and each of the more extreme tables. Factorials growth quickly, so it's necessary use logarithms of factorials. This computations is very easy in Matlab because x!=gamma(x+1) and log(x!)=gammaln(x+1). This function is now fully vectorized to speed up the computation. Actually, the function also computes the mid-P correction to make the test less conservative.

You can visit my homepage http://home.tele2.it/cardillo
My profile on XING http://www.xing.com/go/invita/13675097
My profile on LinkedIN http://it.linkedin.com/in/giuseppecardillo

• MATLAB release: MATLAB 7.5 (R2007b)

[Jon]

I get error messages for some data.
For instance:
y =
    27 13 9 15
    12 10 3 16
>> myfisher24(y)
??? Attempted to access idxBstart(15); index out of bounds because
numel(idxBstart)=14.
Error in ==> myfisher24 at 116
            Tables(idxBstart(J+1):idxBstop(J+1),2)=J;

I hope this can be fixed?
Jon

[Giuseppe Cardillo]

Jon was right. I have fixed the bug and uploaded the new version

[Jon]

Giuseppe thank you for your quick response.
Your routine takes all my data now.
Jon