Full_Name: John Hendrickx
Version: 1.6.0
OS: Windows 98
Submission from: (NULL) (137.224.174.216)
CrossTable in the "gregmisc" package fails when the fisher.exact test produces
an error (I suspect this is because the number of cases is too large). This can
be fixed using "FTt <- try(fisher.test(t, alternative = "two.sided"))" or by
making the test optional.
bugtab <- array(c(
24, 31, 2, 50,
4, 101, 16, 130,
0, 29, 3, 52,
3, 39, 11, 210))
dim(bugtab)<-c(4,4)
bugtab<-t(bugtab)
rownames(bugtab)<-c("a","b","c","d")
colnames(bugtab)<-c("A","B","C","D")
bugtab
thisworks<-trunc(bugtab/10)
thisworks
library(gregmisc)
# fisher.exact test works fine on small tables
CrossTable(thisworks)
# I suspect the bugtab table is too large
CrossTable(bugtab)
bug in CrossTable (package:gregmisc) (PR#2480)
2 messages · John_Hendrickx@yahoo.com, Marc Schwartz
-----Original Message-----
From: r-devel-admin@stat.math.ethz.ch
[mailto:r-devel-admin@stat.math.ethz.ch] On Behalf Of
John_Hendrickx@yahoo.com
Sent: Tuesday, January 21, 2003 6:51 AM
To: r-devel@stat.math.ethz.ch
Cc: R-bugs@biostat.ku.dk
Subject: [Rd] bug in CrossTable (package:gregmisc) (PR#2480)
Full_Name: John Hendrickx
Version: 1.6.0
OS: Windows 98
Submission from: (NULL) (137.224.174.216)
CrossTable in the "gregmisc" package fails when the
fisher.exact test produces an error (I suspect this is
because the number of cases is too large). This can be fixed
using "FTt <- try(fisher.test(t, alternative = "two.sided"))"
or by making the test optional.
bugtab <- array(c(
24, 31, 2, 50,
4, 101, 16, 130,
0, 29, 3, 52,
3, 39, 11, 210))
dim(bugtab)<-c(4,4)
bugtab<-t(bugtab)
rownames(bugtab)<-c("a","b","c","d")
colnames(bugtab)<-c("A","B","C","D")
bugtab
thisworks<-trunc(bugtab/10)
thisworks
library(gregmisc)
# fisher.exact test works fine on small tables
CrossTable(thisworks)
# I suspect the bugtab table is too large
CrossTable(bugtab)
John, Thank you for bringing this to my attention. I was out of town all day on a business trip and will look at this issue further with fresh eyes in the morning. I'll get an update to Greg as soon as I have a chance to incorporate a fix. I also have some pending updates for barplot2(), so I will coordinate the updates with Greg to conserve his time. Best regards, Marc Schwartz