Skip to content
Back to formatted view

Raw Message

Message-ID: <de2296f30504290833bf9f251@mail.gmail.com>
Date: 2005-04-29T15:33:36Z
From: Gunnar Hellmund
Subject: Subarrays

Define an array

> v<-1:256
> dim(v)<-rep(4,4)

Subarrays can be obtained as follows:

> v[3,2,,2]
[1]  71  87 103 119
> v[3,,,2]
     [,1] [,2] [,3] [,4]
[1,]   67   83   99  115
[2,]   71   87  103  119
[3,]   75   91  107  123
[4,]   79   95  111  127

In the general case this procedure is very tedious. 

Given an array 
A, dim(A)=(dim_1,dim_2,...,dim_d) 
and two vectors
v1=(n_i1,...n_ik), v2=(int_1,...,int_k) ('marginals' and relevant
'interval numbers')
is there a smart way to obtain 
A[,...,int_1,....,int_2,....,....,int_k,....]
?

Best wishes
Gunnar Hellmund