How about something like:
wilcox.bonf <- function(a,b) {
n <- length(a)
mapply( function(x,y) wilcox.test(x,y)$p.value*n,a,b)
}
On 11/09/06, Raj, Towfique <[EMAIL PROTECTED]> wrote:
> Dear all,
>
> I am trying to run Wilcoxon Rank-Sum Test with Bonferroni's
> correction. I have two lists: l0, l1:
>
> mapply(function(x,y)wilcox.test(x,y)$p.value, l0, l1)
>
> How do I run Bonferroni's correction on mapply? Any help is much apperciated.
> Thanks,
>
> -Raj
>
> ______________________________________________
> R-help@stat.math.ethz.ch mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.
>
--
=================================
David Barron
Said Business School
University of Oxford
Park End Street
Oxford OX1 1HP
______________________________________________
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
