| bdiag {Matrix} | R Documentation |
Build a block diagonal matrix given several building block matrices.
bdiag(...) .bdiag(lst)
... |
individual matrices or a |
lst |
non-empty |
For non-trivial argument list, bdiag() calls .bdiag().
The latter maybe useful to programmers
A sparse matrix obtained by combining the arguments into a block diagonal matrix.
The value of bdiag() inheris from class
CsparseMatrix, whereas
.bdiag() returns a TsparseMatrix.
Martin Maechler, built on a version posted by Berton Gunter to
R-help; earlier versions have been posted by other authors, notably
Scott Chasalow to S-news. Doug Bates's faster implementation builds
on TsparseMatrix objects.
Diagonal for constructing matrices of
class diagonalMatrix, or kronecker
which also works for "Matrix" inheriting matrices.
bandSparse constructs a banded sparse matrix from
its non-zero sub-/super - diagonals.
Note that other CRAN R packages have own versions of bdiag()
which return traditional matrices.
bdiag(matrix(1:4, 2), diag(3))
## combine "Matrix" class and traditional matrices:
bdiag(Diagonal(2), matrix(1:3, 3,4), diag(3:2))
mlist <- list(1, 2:3, diag(x=5:3), 27, cbind(1,3:6), 100:101)
bdiag(mlist)
stopifnot(identical(bdiag(mlist),
bdiag(lapply(mlist, as.matrix))))
ml <- c(as(matrix((1:24)%% 11 == 0, 6,4),"nMatrix"),
rep(list(Diagonal(2, x=TRUE)), 3))
mln <- c(ml, Diagonal(x = 1:3))
stopifnot(is(bdiag(ml), "lsparseMatrix"),
is(bdiag(mln),"dsparseMatrix") )