[R] How to create a function returning an array ?

[Duncan Murdoch]

On 27/06/2013 11:38 PM, David Winsemius wrote:
> On Jun 27, 2013, at 8:04 PM, Kaptue Tchuente, Armel wrote:
>
> > Hi there,
> >
> > I would like to know how to change the line "img=as.single(rnorm(m)))" such that instead of being a vector of length m as it is now, img is an array of dimension c=(n,m,o) for instance
> >
> > ---------------------------------
> > read_ts<-function(n,m,o,img) {
> >   out<-.Fortran("read_ts",
> >                as.integer(n),
> >                as.integer(m),
> >                as.integer(o),
> >                img=as.single(rnorm(n)))
> >   return(out$img)
> > ------------------------------------------
> >
>
> Well, assuming that  the 'out$img' object has a R-length of n*m*o , wouldn't if be simpler to just change the return call to:

In fact, out$img has a length of n, same as on input.  .Fortran won't 
change the length of its arguments.

Duncan Murdoch

>
> return( array( out$img, dim=c(n,m,o) )
>
> I don't think you wnat start naming your dimension vectors "c".
>