Maybe I'm wrong, but aren't you calculating just the same as cor(data
frame,method ="spearman"), with some further parameters being monitored?
Could you please provide a commented, minimal, self-contained, reproducible
code, so that we can see what is actually going on and what is the way you
want it?
Dieter Best wrote:
>
> Hi there,
>
> I have a multiple for loop over a list of data frames
>
> for ( i in 1:(N-1) ) {
> for ( j in (i+1):N ) {
> for ( p in 1:M ) {
> v_i[p] = alist[[p]][i,"v"]
> v_j[p] = alist[[p]][j,"v"]
> }
> rho_s = cor(v_i, v_j, method = "spearman")
> rho_p = cor(v_i, v_j, method = "pearson" )
> iv = c( iv, min(i, j) )
> jv = c( jv, max(i, j) )
> rho_sv = c( rho_sv, rho_s)
> rho_pv = c( rho_pv, rho_p)
> }
> }
>
> N is of the order of 400, M about 800.
>
> This takes me an entire day basically. Is there anything I could do to
> speed things up or is cor really that slow?
>
> -- D
>
>
>
> ---------------------------------
>
>
> [[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.
>
>
--
View this message in context:
http://www.nabble.com/How-to-speed-up-multiple-for-loop-over-list-of-data-frames-tf4638585.html#a13249895
Sent from the R help mailing list archive at Nabble.com.
______________________________________________
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.
