[R] elegant matrix creation

Robin Hankin rksh at soc.soton.ac.uk
Wed Jul 28 16:49:19 CEST 2004


Spencer

thank you!

just what I wanted.   I have a feeling that jj3 is qualitatively 
different from the other two.
Nevertheless, I'm sure it'll crack sooner or later!

best wishes

rksh




At 07:34 am -0700 28/07/04, Spencer Graves wrote:
>      It's not obvious what pattern you want, but some variants of 
>the following would work for jj1 and jj2:
>      1+outer(1:9, 1:9, "+")%%3
>
>In particular the following are equal to your jj1 and jj2: 
>jj1. <- 1+outer(0:8, rep(0:2, e=3), "+")%%3
>
>jj2. <- 1+outer(0:8, c(1,2,1,3,1,3,2,3,2)-1, "+")%%3    
>
>      I couldn't figure out jj3, but this system may not work so 
>easily for that.  hope this helps.  spencer graves
>
>
>Robin Hankin wrote:
>
>>Hello everybody.
>>
>>I am trying to reproduce a particular matrix in an elegant way.  If I
>>have
>>
>>jj1 <-
>>structure(c(1,2,3,1,2,3,1,2,3,1,2,3,1,2,3,1,2,
>>3,1,2,3,1,2,3,1,2,3,2,3,1,2,3,1,2,3,1,2,3,
>>1,2,3,1,2,3,1,2,3,1,2,3,1,2,3,1,3,1,2,3,1,
>>2,3,1,2,3,1,2,3,1,2,3,1,2,3,1,2,3,1,2,3,1,
>>2),.Dim = as.integer(c(9,9)))
>>
>>[ image(jj1) is good here ] then I can get this with
>>
>>kronecker(matrix(1,3,1),kronecker(1+outer(0:2,0:2,"+")%%3,matrix(1,1,3)))
>>
>>I want to reproduce the following matrices in an equivalent way:
>>
>>jj2 <- matrix(c(1,2,3,1,2,3,1,2,3,2,3,1,2,3,1,2,3,1,
>>1,2,3,1,2,3,1,2,3,3,1,2,3,1,2,3,1,2,1,2,3,1,2,
>>3,1,2,3,3,1,2,3,1,2,3,1,2,2,3,1,2,3,1,2,3,1,3,
>>1,2,3,1,2,3,1,2,2,3,1,2,3,1,2,3,1),9,9)
>>
>>jj3 <- structure(c(1,2,3,2,3,1,3,1,2,1,2,1,2,3,2,3,1,
>>3,1,3,1,2,1,2,3,2,3,2,3,1,3,1,2,1,2,3,2,3,
>>2,3,1,3,1,2,1,2,1,2,3,2,3,1,3,1,3,1,2,1,2,
>>3,2,3,1,3,1,3,1,2,1,2,3,2,3,2,3,1,3,1,2,1, 2),.Dim =
>>as.integer(c(9,9)))
>>
>>[ note that jj1-jj3 each have precisely 3 occurrences of A, B, and C
>>along each row, column and (broken) diagonal ].
>>
>>Can anyone give me a nice elegant way of creating jj2 and jj3 please?
>>
>>


-- 
Robin Hankin
Uncertainty Analyst
Southampton Oceanography Centre
SO14 3ZH
tel +44(0)23-8059-7743
initialDOTsurname at soc.soton.ac.uk (edit in obvious way; spam precaution)




More information about the R-help mailing list