[R] How to generating diagnal blocks ?
Robin Hankin
r.hankin at noc.soton.ac.uk
Thu Sep 21 09:28:18 CEST 2006
Hello
You need adiag(), which
is now part of the magic package:
> library(magic)
> adiag(matrix(1:6,2,3),matrix(1:9,3,3))
[,1] [,2] [,3] [,4] [,5] [,6]
[1,] 1 3 5 0 0 0
[2,] 2 4 6 0 0 0
[3,] 0 0 0 1 4 7
[4,] 0 0 0 2 5 8
[5,] 0 0 0 3 6 9
>
the function also has extra functionality over the one Vito
points to in the archives.
HTH
rksh
On 21 Sep 2006, at 08:14, vito muggeo wrote:
> If I remember well, there should be a package including the function
> bdiag() making the job..but I am not able to remember its name..
>
> However a quick search via RSiteSearch("bdiag") yields
>
> http://finzi.psych.upenn.edu/R/Rhelp02a/archive/40393.html
>
> Hope this helps you,
>
> vito
>
>
> Tong Wang wrote:
>> Hi,
>> I am trying to creat a matrix with diagnal blocks, say, I
>> have a matrix X of any dimension (nxm) ,and would like to have:
>>
>> X
>> X
>> X
>>
>> .....
>> otherwise space filled with 0's. Is there a handy way to do this ?
>>
>> Thanks a lot in advance.
>>
>> best
>>
>> ______________________________________________
>> R-help at stat.math.ethz.ch 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.
>>
>>
>
> --
> ====================================
> Vito M.R. Muggeo
> Dip.to Sc Statist e Matem `Vianelli'
> Università di Palermo
> viale delle Scienze, edificio 13
> 90128 Palermo - ITALY
> tel: 091 6626240
> fax: 091 485726/485612
>
> ______________________________________________
> R-help at stat.math.ethz.ch 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.
--
Robin Hankin
Uncertainty Analyst
National Oceanography Centre, Southampton
European Way, Southampton SO14 3ZH, UK
tel 023-8059-7743
More information about the R-help
mailing list