R help - Nov 2007 - fisher.test, chisq.test

Hi,

I want to analyse a contigency table (3 x 12) with a fisher.test
beacause there are cells that are less than 5.

        ?mmen Anken Baf Belchen H?chi Hof Porti R?m Schmutz Sch?n Sissa Tann
class14     7    26 150       2    46  68   126  66       3    31     7   61
class24     7     6  55       5    49  71    93  90       1    18    16   79
class34     1     1   4       3    19   8    29  61       4     4     5   37

So in made a matrix and used
>fisher.test(classtable4).
But then this error occurs:

Fehler in fisher.test(classtable4) : FEXACT error 7.
LDSTP is too small for this problem.
Try increasing the size of the workspace.

Well I then increased the workspace:
>fisher.test(classtable4, workspace=1e8)
But then the same error occurs.

Do someone knows an alternative to the fisher.test or the chisq.test?

Thank you very much!

All the best,
Denis

PS: I use R 2.6.0-patched on a PC with Windows XP.

On Fri, 9 Nov 2007, Denis Aydin wrote:
> Hi,
>
> I want to analyse a contigency table (3 x 12) with a fisher.test
> beacause there are cells that are less than 5.
...expected cells. Anyways, your data set yields a statistic > 160
and I wouldn't worry that this could be the result of random variation
under the null hypothesis. Also looking at the data, e.g., via a
spineplot() of the (transposed) table would clearly indicate a deviation
from the null hypothesis.
>         ?mmen Anken Baf Belchen H?chi Hof Porti R?m Schmutz Sch?n Sissa
Tann
> class14     7    26 150       2    46  68   126  66       3    31     7  
61
> class24     7     6  55       5    49  71    93  90       1    18    16  
79
> class34     1     1   4       3    19   8    29  61       4     4     5  
37
>
> So in made a matrix and used
>
> >fisher.test(classtable4).
>
> But then this error occurs:
>
> Fehler in fisher.test(classtable4) : FEXACT error 7.
> LDSTP is too small for this problem.
> Try increasing the size of the workspace.
>
> Well I then increased the workspace:
>
> >fisher.test(classtable4, workspace=1e8)
>
> But then the same error occurs.
You don't really want/need the exact distribution...
> Do someone knows an alternative to the fisher.test or the chisq.test?
Both functions have an argument simulate.p.value which you could set to
TRUE, possibly increasing B (the number of replications). Even with a
million replications (under independence) I did not obtain a single one
that was larger than the observed statistic. So I think it is pretty safe
to reject independence...

hth,
Z