Dear R-users,
I would like to know whether you know some trick for skipping some of
the steps in the example below (especially the last step in a way that
would make easier to be written succinctly in a text).
I could try to explain in words the whole process, but I'm sure the code
below would be clearer.
Thanks in advance for your help,
Giancarlo
## data in matrices
D <- matrix(1:15, 3, 5)
T <- matrix(0, 3, 3)
T[c(2,4,6,8)] <- 1
## how to place the diag matrices of each row
M0 <- matrix(0, nrow(T), sum(T))
wr <- which(T==1, arr.ind=TRUE)[,2]
wc <- 1:ncol(M0)
M0[cbind(wr,wc)] <- 1
## number of columns
m <- ncol(D)
## non-zero positions
M <- kronecker(M0, diag(m))
## which rows to take
pos <- which(T==1, arr.ind=TRUE)[,1]
## filling up with data
M[M!=0] <- t(D[wr,])
[[alternative HTML version deleted]]
______________________________________________
[email protected] 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.
