Perhaps this might help clarify: sample A: 10 15 20 sample B: 12 15 22
Median of sample A = 15; median of sample B = 15. Sample medians are =. But: B-A differences are 2,0,2 with a median of 2. So the median difference does not equal the difference of the medians. Clarity in what you wish to test (and why!) is essential to determine how. Bert Gunter Genentech Nonclinical Statistics -----Original Message----- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of cheba meier Sent: Friday, May 07, 2010 12:46 PM To: Joris Meys Cc: R-help@r-project.org; Thomas Lumley Subject: Re: [R] median of two groups Hi all, Thank you for your reply. if done properly! What does this mean? The R-code I have is using the R-function sample without replacement. Am I doing this properly? median of the differences is zero! Does this mean if I run 1000 permutation and for each permutation I compute the median difference and as a result I have 1000 differences. Is the the H0: median(1000 differences) =0? If yes, which conclusion one would have from this H0? Best wishes, Cheba 2010/5/7 Joris Meys <jorism...@gmail.com> > depends on how you interprete "absolute median difference". Is that the > absolute difference of the medians, or the median of the absolute > differences. Probably the latter one, so you would be right. If it's the > former one, then it is testing whether the difference of the medians is > zero. > > Cheers > Joris > > > On Fri, May 7, 2010 at 6:52 PM, Thomas Lumley <tlum...@u.washington.edu>wrote: > >> On Fri, 7 May 2010, cheba meier wrote: >> >> Dear Thomas, >>> >>> I have been running simulations in order me to understand this problem! I >>> have found something online where the absolute median difference is >>> computed >>> and permutations are ran to compute a p-value. Is such a test (if I can >>> call >>> it a test) tests the null hypothesis that median group 1 = median group >>> 2? >>> >> >> No, that is testing whether the median of the differences is zero. This >> is not the same as testing whether the difference of the medians is zero. >> >> -thomas >> >> >> >> Thank you in advance for your help. >>> >>> Regards, >>> Cheba >>> >>> 2010/4/6 Thomas Lumley <tlum...@u.washington.edu> >>> >>> >>>> >>>> None of them. >>>> >>>> - mood.test() looks promising until you read the help page and see that >>>> it >>>> does not do Mood's test for equality of quantiles, it does Mood's test >>>> for >>>> equality of scale parameters. >>>> - wilcox.test() is not a test for equal medians >>>> - ks.test() is not a test for equal medians. >>>> >>>> >>>> Mood's test for the median involves dichotomizing the data at the pooled >>>> median and then doing Fisher's exact test to see if the binary variable >>>> has >>>> the same mean in the two samples. >>>> >>>> median.test<-function(x,y){ >>>> z<-c(x,y) >>>> g <- rep(1:2, c(length(x),length(y))) >>>> m<-median(z) >>>> fisher.test(z<m,g)$p.value >>>> } >>>> >>>> Like most exact tests, it is quite conservative at small sample sizes. >>>> >>>> -thomas >>>> >>>> >>>> On Tue, 6 Apr 2010, cheba meier wrote: >>>> >>>> Dear all, >>>> >>>>> >>>>> What is the right test to test whether the median of two groups are >>>>> statistically significant? Is it the wilcox.test, mood.test or the >>>>> ks.test? >>>>> In the text book I have got there is explanation for the Wilcoxon (Mann >>>>> Whitney) test which tests ob the two variable are from the same >>>>> population >>>>> and also ks.test! >>>>> >>>>> Regards, >>>>> Cheba >>>>> >>>>> [[alternative HTML version deleted]] >>>>> >>>>> ______________________________________________ >>>>> R-help@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. >>>>> >>>>> >>>>> Thomas Lumley Assoc. Professor, Biostatistics >>>> tlum...@u.washington.edu University of Washington, Seattle >>>> >>>> >>>> >>> [[alternative HTML version deleted]] >>> >>> ______________________________________________ >>> R-help@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. >>> >>> >> Thomas Lumley Assoc. Professor, Biostatistics >> tlum...@u.washington.edu University of Washington, Seattle >> >> ______________________________________________ >> R-help@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. >> > > > > -- > Joris Meys > Statistical Consultant > > Ghent University > Faculty of Bioscience Engineering > Department of Applied mathematics, biometrics and process control > > Coupure Links 653 > B-9000 Gent > > tel : +32 9 264 59 87 > joris.m...@ugent.be > ------------------------------- > Disclaimer : http://helpdesk.ugent.be/e-maildisclaimer.php > [[alternative HTML version deleted]] ______________________________________________ R-help@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. ______________________________________________ R-help@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.
