Skip to content

Arrays of variable dimensionality

7 messages · Jim Lemon, Brian Ripley, William Dunlap +3 more

Original
WRAY NICHOLAS
# 
Hello folks

Supposing I have a multidimensional array in an R prog, say a 4D array.

I want the coordinate quantities to be read off from a vector.  The values
in the vector (vec) are generated by a function.

This is easy if the number of dimensions is fixed for both the array and
the number of elements in the  vector (say 4):

X<-array[vec[1],vec[2],vec[3],vec[4]]

But if the number of dimensions of the array is not fixed and can change
during the course of the prog I am stuck as to how to do this, as I don?t
know a way of feeding a vector or list of unspecified beforehand length
into the coordinates for an array

I can?t find anything useful about this on the net

Does anyone have any ideas?
Thanks, Nick
Jim Lemon
# 
Hi Nick,
Does length(dim(X)) help by getting the current dimensions of the
array before you pass the vectors of indices?

Jim


On Sat, May 30, 2015 at 8:29 PM, WRAY NICHOLAS
<nicholas.wray at ntlworld.com> wrote:
Brian Ripley
#
On 30/05/2015 11:29, WRAY NICHOLAS wrote:
It is not clear what you want to do.
Is 'array' an array and not the function of that name?

I am guess you are missing do.call, the way to generate calls with 
variable-length argument lists.  My guess of your intention could be written

arr <- array(1:256, rep(4,4))
vec <- 1:4
arr[1, 2, 3, 4] # same as
do.call(`[`, c(list(arr), as.list(vec)))
PLEASE do.
William Dunlap
# 
I don't know for sure what result you want, but it may be that
using a length(dim(array))-column matrix as your subscript will
do what you want.  E.g.,
  > a <- array(101:124, dim=4:2)
  > subscriptMatrix <- cbind(c(1,3), c(2,2), c(2,1))
  > subscriptMatrix
       [,1] [,2] [,3]
  [1,]    1    2    2
  [2,]    3    2    1
  > a[subscriptMatrix]
  [1] 117 107
  > a[1, 2, 2]
  [1] 117
  > a[3, 2, 1]
  [1] 107


Bill Dunlap
TIBCO Software
wdunlap tibco.com

On Sat, May 30, 2015 at 3:29 AM, WRAY NICHOLAS <nicholas.wray at ntlworld.com>
wrote:

  
  


  


      

    

    

      
      

        


  

        

      Jeff Newmiller
    

    

      
      #
    

  


  

      
Don't.

Arrays in R don't have variable dimensions [1]. To the extent that you pretend otherwise, you will degrade performance and complicate your program. I would argue that similar effects arise in languages that do let you pretend that arrays have variable dimensions.

The main reason variable-sized arrays seem useful is that some algorithms or data sources yield variable amounts of data. There are various ways to handle these cases, and the use of lists to hold intermediate units of data followed by transfer of that data into an array (strategically handled after you know how big the array has to be) is one of the more common ones.

Another technique is to use one of the ragged array or sparse matrix data structures, but they generally work best when the data actually are sparse.

Being more specific about your immediate problem can elicit more specific help.

BTW please post plain text on this list... HTML is not supported by the list reflector and leads to misunderstandings.

[1] More precisely, the total number of elements in the underlying vector is fixed, though you can redefine how the dimensions use the elements in that vector. For example, study the "aperm" function.
---------------------------------------------------------------------------
Jeff Newmiller                        The     .....       .....  Go Live...
DCN:<jdnewmil at dcn.davis.ca.us>        Basics: ##.#.       ##.#.  Live Go...
                                      Live:   OO#.. Dead: OO#..  Playing
Research Engineer (Solar/Batteries            O.O#.       #.O#.  with
/Software/Embedded Controllers)               .OO#.       .OO#.  rocks...1k
--------------------------------------------------------------------------- 
Sent from my phone. Please excuse my brevity.

        
On May 30, 2015 3:29:09 AM PDT, WRAY NICHOLAS <nicholas.wray at ntlworld.com> wrote:

            

      
      
      
    

        

      
      
    

  
  


  


      

    

    

      
      

        


  

        

      Henrik Bengtsson
    

    

      
      #
    

  


  

      
Basically dropping by saying what's already been said, but also that
you probably want to do whatever you're trying to achieve using
vectorized operations.  A small summary with some tweaks:

## Generate an array with a "random" number of dimensions
dim <- c(4, 5, 6, 7)
dimnames <- lapply(1:4, FUN=function(i) sprintf("%s%d", letters[i], 1:dim[i])
)
x <- array(seq_len(prod(dim)), dim=dim, dimnames=dimnames)

# A single cell

            

      
      
      
    

        
[1] 405

            

      
      
      
    

        
# A "block" subset of the array

            

      
      
      
    

        
c3  c4
a1 405 425
a2 406 426

            

      
      
      
    

        
# A "block" subset of the array (preserving all 4 dimensions)

            

      
      
      
    

        
int [1:2, 1, 1:2, 1] 405 406 425 426
 - attr(*, "dimnames")=List of 4
  ..$ : chr [1:2] "a1" "a2"
  ..$ : chr "b2"
  ..$ : chr [1:2] "c3" "c4"
  ..$ : chr "d4"

            

      
      
      
    

        
# Why alist() and not list()?  Because if you also need to do ...

            

      
      
      
    

        
int [1:2, 1:5, 1:2, 1:7] 1 2 5 6 9 10 13 14 17 18 ...
 - attr(*, "dimnames")=List of 4
  ..$ : chr [1:2] "a1" "a2"
  ..$ : chr [1:5] "b1" "b2" "b3" "b4" ...
  ..$ : chr [1:2] "c1" "c2"
  ..$ : chr [1:7] "d1" "d2" "d3" "d4" ...

            

      
      
      
    

        
# Any subset of cells - not a "block" (using what is called "matrix
indexing" in R)

            

      
      
      
    

        
[1] 385 247

            

      
      
      
    

        
For a (possibly) more readable subsetting by "blocks" than using
do.call(), there's also:

            

      
      
      
    

        
It also has some other bells and whistles, but it cannot do more than
what already shown above.

/Henrik
(author of R.utils)

On Sat, May 30, 2015 at 10:34 AM, Jeff Newmiller

        
<jdnewmil at dcn.davis.ca.us> wrote:

            

      
      
      
    

        

      
      
    

            

      
      
      
    

        

      
      
    

  
  


  


      

    

    

      
      

        


  

        

      WRAY NICHOLAS
    

    

      
      #
    

  


  

      
Thanks for all the ideas   I will run through them next week when I'm back
at work... Nick
---------- Forwarded message ----------
From: WRAY NICHOLAS <nicholas.wray at ntlworld.com>
Date: 30 May 2015 at 11:29
Subject: Arrays of variable dimensionality
To: r-help <r-help at r-project.org>


Hello folks

Supposing I have a multidimensional array in an R prog, say a 4D array.

I want the coordinate quantities to be read off from a vector.  The values
in the vector (vec) are generated by a function.

This is easy if the number of dimensions is fixed for both the array and
the number of elements in the  vector (say 4):

X<-array[vec[1],vec[2],vec[3],vec[4]]

But if the number of dimensions of the array is not fixed and can change
during the course of the prog I am stuck as to how to do this, as I don?t
know a way of feeding a vector or list of unspecified beforehand length
into the coordinates for an array

I can?t find anything useful about this on the net

Does anyone have any ideas?
Thanks, Nick