Dear all, I want to apply a weighted median on a huge dataset, and I remember a function from the package R.basic that could do this using an internal sorting algorithm qsort. This speeded things up quite a bit. Alas, I can't find that package anywhere anymore. There is a weighted.median function in the package limma too, but I didn't use that before. Anybody who knows what happened to R.basic? Cheers Joris -- 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.Meys@Ugent.be ------------------------------- Disclaimer : http://helpdesk.ugent.be/e-maildisclaimer.php [[alternative HTML version deleted]]
Hi, good memory. weightedMedian() is now available in the aroma.light package (it was moved there from R.basic in Feb 2006). /Henrik (author of both packages) On Tue, Mar 30, 2010 at 6:30 PM, Joris Meys <jorismeys at gmail.com> wrote:> Dear all, > > I want to apply a weighted median on a huge dataset, and I remember a > function from the package R.basic that could do this using an internal > sorting algorithm qsort. This speeded things up quite a bit. Alas, I can't > find that package anywhere anymore. There is a weighted.median function in > the package limma too, but I didn't use that before. > > Anybody who knows what happened to ?R.basic? > > Cheers > Joris > > -- > 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.Meys at Ugent.be > ------------------------------- > Disclaimer : http://helpdesk.ugent.be/e-maildisclaimer.php > > ? ? ? ?[[alternative HTML version deleted]] > > ______________________________________________ > 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. >
While perhaps not the solution you were looking for, you might consider estimating weighted medians with linear quantile regression (just specify an intercept for single sample analysis, tau=0.50, and weights = your weights) in the quantreg package. Quantile regression does not require sorting to estimate medians (minimizes and objective function) and thus might require less computing time on a large data set. Brian Brian S. Cade, PhD U. S. Geological Survey Fort Collins Science Center 2150 Centre Ave., Bldg. C Fort Collins, CO 80526-8818 email: brian_cade@usgs.gov tel: 970 226-9326 From: Joris Meys <jorismeys@gmail.com> To: R mailing list <r-help@r-project.org> Date: 03/30/2010 10:39 AM Subject: [R] weighted.median function from package R.basic Sent by: r-help-bounces@r-project.org Dear all, I want to apply a weighted median on a huge dataset, and I remember a function from the package R.basic that could do this using an internal sorting algorithm qsort. This speeded things up quite a bit. Alas, I can't find that package anywhere anymore. There is a weighted.median function in the package limma too, but I didn't use that before. Anybody who knows what happened to R.basic? Cheers Joris -- 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.Meys@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. [[alternative HTML version deleted]]