Nothing
bdiag <-
function (...) {
mlist <- list(...)
## handle case in which list of matrices is given
if (length(mlist) == 1)
mlist <- unlist(mlist, recursive = FALSE)
csdim <- rbind(c(0, 0), apply(sapply(mlist, dim), 1, cumsum))
ret <- array(0, dim = csdim[length(mlist) + 1, ])
add1 <- matrix(rep(1:0, 2), ncol = 2)
for (i in seq_along(mlist)) {
indx <- apply(csdim[i:(i+1), ] + add1, 2, function(x) x[1]:x[2])
if(is.null(dim(indx))) {
## non-square matrix
ret[indx[[1]], indx[[2]]] <- mlist[[i]]
} else {
## square matrix
ret[indx[, 1], indx[, 2]] <- mlist[[i]]
}
}
colnames(ret) <- unlist(lapply(mlist, colnames))
ret
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.