Entering edit mode
Jose Duarte
▴
100
@jose-duarte-555
Last seen 10.2 years ago
I'm getting an error while using webbioc in the multtest step.
All works fine when doing a two-sample t-test (equal or unequal
variances) and choosing "abs" or "upper" as Rejection region.
However when choosing "lower" I get the following error back:
Error in pt(teststat, df1, df2, lower.tail = FALSE) :
Object "df1" not found
After looking a bit into the R code that webbioc is using it seems
that
the error is in the mt.wrapper function of the multtest package. It
should only use one "df" and not "df1" and "df2" when doing the
two-sample t-test. But it looks like for "lower" is using both "df1"
and
"df2" while is alright for "upper" and "abs".
Is this a bug? Or am I missing something?
Versions are:
Linux RedHat 8
R 1.8.0
webbioc 0.2
multtest 1.3.1
Thanks in advance
Jose
--
Jose Duarte
MRC Functional Genetics Unit
Department of Human Anatomy and Genetics
University of Oxford
South Parks Road
Oxford OX1 3QX
Phone: +44 (0)1865 282650
e-mail: jose.duarte@anat.ox.ac.uk