FW: [R] Optimising code

Crispin Miller CMiller at PICR.man.ac.uk
Tue Oct 7 13:21:19 CEST 2003


Yup, my mistake for not being very clear! Alas the computer is already very fast - looks like it's C for me...

> -----Original Message-----
> From: Patrick Burns [mailto:pburns at pburns.seanet.com]
> Sent: 07 October 2003 12:00
> Cc: Crispin Miller
> Subject: Re: FW: [R] Optimising code
> 
> 
> I had misunderstood also.  I agree with Uwe either move
> to C or buy a faster computer -- or at least another one.
> 
> Pat
> 
> Uwe Ligges wrote:
> 
> > Crispin Miller wrote:
> >
> >>>> I have a function that applies a wilcoxon test to 12 
> sets of about 
> >>>> a quarter of a million pairs
> >>>
> >>
> >>
> >>> ... and let me guess: everything is significiant to an almost 
> >>> arbitrary value of \alpha?
> >>
> >>
> >>
> >> :-) For each of quarter of a million sets, I do a wilcoxon between 
> >> two pairs each containing twenty numbers...
> >> I do this 12 times...
> >
> >
> > Ah. Sorry for misunderstanding.
> >
> >
> >>>> (and takes about 3 hours). I've replaced the inner loop I 
> >>>
> >>>
> >>> had originally with a function call via mapply, and also 
> considered 
> >>> different approximations  of the wilcoxon, rather than 
> that which is 
> >>> implemented in wilcox.test, but that makes little 
> difference 9and if 
> >>> anything slows things down :-).
> >>>
> >>> Are you using wilcox.test(), or a self written function?
> >>>
> >>
> >>
> >> I started off with wilcox.test, and tried to speed it up 
> by hacking 
> >> out as much from the original function as possible, since a fair 
> >> proportion of the code being executed was logic deciding 
> between the 
> >> different flavours of the test (I didn't expect that to make much 
> >> difference, but tried it anyway), there was also a call to pnorm() 
> >> which I replaced with a numeric approximation. This is not inlined 
> >> into the original function.
> >
> >
> >
> > In that case you cannot do very much except for writing it in a 
> > language like C / Fortran, if you really need the speed.
> >
> > Uwe Ligges
> >
> >> Crispin
> >>  
> >> --------------------------------------------------------
> >>
> >>  
> >> This email is confidential and intended solely for the use 
> >> o...{{dropped}}
> >>
> >> ______________________________________________
> >> R-help at stat.math.ethz.ch mailing list
> >> https://www.stat.math.ethz.ch/mailman/listinfo/r-help
> >
> >
> > ______________________________________________
> > R-help at stat.math.ethz.ch mailing list
> > https://www.stat.math.ethz.ch/mailman/listinfo/r-help
> >
> >
> 
> 
> 
>
 
--------------------------------------------------------

 
This email is confidential and intended solely for the use o...{{dropped}}




More information about the R-help mailing list