R help - Feb 2009 - wilcoxon test with bonferroni correction

Hi!
I need to run a wilcoxon (Mann-whitly, in fact) test with bonferroni  
correction, as I am running 10 consecutive wilcoxon test not  
independent, and I know that bonferroni will partially correct for  
this problem, but I have no idea how to do it with R, I have been  
looking in the archive but couldn't understand how to do it.

The format I am using at the moment is

r4_o <-

  [1] 1.05 2.60 1.57 3.07 1.20 1.00 2.11 1.10 0.10

r4_m <-

[1] 0 0 0 0 0 0 0 0 0

wilcoxon.test (r4_o, r3_m)

Does any body know how to make the bonferroni correction when I  
compare them with the wilcoxon test?

Thank you very much.

Lucia






  Lucia Prieto Godino
PhD student.
Department of Zoology,
Downing street
University of Cambridge.
UK

On 2/1/2009 8:32 PM, Laura Lucia Prieto Godino wrote:
> Hi!
> I need to run a wilcoxon (Mann-whitly, in fact) test with bonferroni
> correction, as I am running 10 consecutive wilcoxon test not
> independent, and I know that bonferroni will partially correct for this
> problem, but I have no idea how to do it with R, I have been looking in
> the archive but couldn't understand how to do it.
> 
> The format I am using at the moment is
> 
> r4_o <-
> 
>  [1] 1.05 2.60 1.57 3.07 1.20 1.00 2.11 1.10 0.10
> 
> r4_m <-
> 
> [1] 0 0 0 0 0 0 0 0 0
> 
> wilcoxon.test (r4_o, r3_m)
> 
> Does any body know how to make the bonferroni correction when I compare
> them with the wilcoxon test?
# Ten p-values
> X <- seq(.001, .10, len=10)
> X
 [1] 0.001 0.012 0.023 0.034 0.045 0.056 0.067 0.078 0.089 0.100

# Same ten p-values adjusted by the Bonferroni method
> p.adjust(X, method="bonferroni")
 [1] 0.01 0.12 0.23 0.34 0.45 0.56 0.67 0.78 0.89 1.00

?p.adjust
> Thank you very much.
> 
> Lucia 
> 
>  Lucia Prieto Godino
> PhD student.
> Department of Zoology,
> Downing street
> University of Cambridge.
> UK
> 
> ______________________________________________
> R-help at r-project.org 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. 
-- 
Chuck Cleland, Ph.D.
NDRI, Inc. (www.ndri.org)
71 West 23rd Street, 8th floor
New York, NY 10010
tel: (212) 845-4495 (Tu, Th)
tel: (732) 512-0171 (M, W, F)
fax: (917) 438-0894