[R] Looping multiple dimensions

[jim holtman]

I would suggest that you use a 'list' for the data.  Something like
this should work:

objs <- grep("^mrunoff_", ls(), values=TRUE)  # get the names of the objects
result <- lapply(objs, function(.obj){
    .obj <- get(.obj)  # get the value
    lapply(1:12, function(. iter)  .obj[,, .iter])
})



On Fri, Jan 22, 2010 at 11:32 AM, Steve Murray <smurray444 at hotmail.com> wrote:
>
> Dear all,
>
> I have 30 arrays, each with dimensions 720,360,12. The naming format for each of these 30 objects is: mrunoff_5221, mrunoff_5222... mrunoff_5250.
>
>
> For example:
>> str(mrunoff_5221)
>   num [1:720, 1:360, 1:12] NA NA NA NA NA NA NA NA NA NA ...  (the initial NA's are nothing to worry about)
>
>
> I am looking for a way by which I can extract each of the third dimension of these grids (1:12) in turn, along with the first and second dimensions, to create new objects in the following style:
>
>
> #2071
> mrunoff_207101 <- mrunoff_5221[,,1]
> mrunoff_207102 <- mrunoff_5221[,,2]
> mrunoff_207103 <- mrunoff_5221[,,3]
> mrunoff_207104 <- mrunoff_5221[,,4]
> mrunoff_207105 <- mrunoff_5221[,,5]   ...(etc. - up to [,,12])
>
> #2072
> mrunoff_207201 <- mrunoff_5222[,,1]
> mrunoff_207202 <- mrunoff_5222[,,2]
> mrunoff_207203 <- mrunoff_5222[,,3]
> mrunoff_207204 <- mrunoff_5222[,,4]
> mrunoff_207205 <- mrunoff_5222[,,5]  ...(etc. - up to [,,12]) and mrunoff_ continues to 2100 and 5250 respectively.
>
>
> Clearly, this is a cumbersome and non-sustainable way to proceed! There will be 360 new objects in total, and I imagine that there must be a more effective way of achieving this, either via a loop or, possibly, one of the 'apply' functions. Yet my attempts to date have so far resulted in... well, a complete mess!
>
> If anyone has any suggestions as to a more efficient means of achieving this, then I'd be very grateful to hear them.
>
> Many thanks,
>
> Steve
>
>
> _________________________________________________________________
> Tell us your greatest, weirdest and funniest Hotmail stories
>
> ______________________________________________
> 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.
>



-- 
Jim Holtman
Cincinnati, OH
+1 513 646 9390

What is the problem that you are trying to solve?