Laura Lucia Prieto Godino
2009-Feb-02 01:32 UTC
[R] 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
Seemingly Similar Threads
- Chi-squared test adjusted for multiple comparisons? Harbe's test?
- Wilcoxon Rank-Sum Test with Bonferroni's correction
- Bonferroni simultaneous confidence intervals for multiple regression
- calculating p-values from t-values for a Bonferroni adjustment
- Bonferroni correction for multiple correlation tests