Karolis Koncevičius
2019-Jun-18 21:37 UTC
[Rd] Small bug in fligner.test - constant values can produce significant results (patch attached)
In specific cases fligner.test() can produce a small p-value even when both
groups have constant variance.
Here is an illustration:
fligner.test(c(1,1,2,2),
c("a","a","b","b"))
# p-value = NA
But:
fligner.test(c(1,1,1,2,2,2),
c("a","a","a","b","b","b"))
# p-value < 2.2e-16
This can potentially get dangerous if people perform lots of parallel tests of
this type (i.e. when doing a test for each gene in genomic studies).
Submitted a proposed patch that should solve the issue by producing an error
"data is essentially constant"
- which is the same error message found in t-test under similar conditions.
P.S. First time writing to this list. Read all the guides of posting, but sorry
in advance if I still missed any rules.
---
Karolis Koncevi?ius
-------------- next part --------------
An embedded and charset-unspecified text was scrubbed...
Name: patch.diff
URL:
<https://stat.ethz.ch/pipermail/r-devel/attachments/20190619/0642547d/attachment.ksh>
Apparently Analagous Threads
- [Suggested patch] to fligner.test - constant values can produce significant results
- Should '@" now be listed in tools:::.get_internal_S3_generics() ?
- Should '@" now be listed in tools:::.get_internal_S3_generics() ?
- Should '@" now be listed in tools:::.get_internal_S3_generics() ?
- Should '@" now be listed in tools:::.get_internal_S3_generics() ?
