similar to: Block-Diagonal Matrix and Multivariate Skew Normal

Dear R users, Suppose i have an A square matrix rxr. I want to obtain a block matrix B (pxr,pxr) where the p diagonal blocks are A and the others values are 0. I would like to do something like : diag(A,ncol=pr, nrow=pr) How can i do it ?? Thanks in advance, Olivier. -- -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- Olivier MARTIN PhD student phone: (33) 04

Dear R-users, Does anybody know how to construct a block-diagonal matrix (with the blocks being different matrixs, concerning the dimension and the values), without use loops ? Thanks all, Caio --------------------------------- [[alternative HTML version deleted]]

Hello List I have just written a little function that takes two matrices as arguments and returns a large matrix that is composed of the two input matrices in upper-left position and lower-right position with a padding value everywhere else. (function definition and toy example below). I need nonsquare matrices and rowname() and colname() inherited appropriately. Two questions: (1) Is there a

> -----Original Message----- > From: r-help-bounces at stat.math.ethz.ch > [mailto:r-help-bounces at stat.math.ethz.ch]On Behalf Of Caio Lucidius > Naberezny Azevedo > Sent: 01 September 2005 12:09 > To: Help mailing list - R > Subject: [R] Multivariate Skew Normal distribution > > > Hi all, > > Could anyone tell me if there is any package (or function)

Hello, I am trying to do some simulation using the rmetasim package and I've run to this problem. --beginning of error msg-- Error in "[<-"(`*tmp*`, slice[l, ], slice[l, ], value = c(0.200000002980232, : number of items to replace is not a multiple of replacement length --end of error msg-- Here is the script I used. --script starts here-- ## load 'rmetasim'

Hi R users: How can I built a direct sum function of matrices in R? I mean A(mxn), B(pxq), C(rxs),... X<-ds(A,B,C,...) X = [ A, 0, 0 0, B, 0 0, 0, C] ((m+p+r+...) x (n+q+s+...)) Thank you for your help. Kenneth -- Using Opera's revolutionary e-mail client: http://www.opera.com/m2/

Is there a way to add covariates to the system equation in a time-varying approach: Y[t] = F'[t]theta[t] + v[t], v[t] ~ N[0,V] #observation equation theta[t] = theta[t-1] + psi*Z[t] + w[t], w[t] ~ N[0,W] #system equation While F[t] is a matrix of regressors to capture the short term effect on the response series Y, Z[t] measures the long-term effect of either (1) two policies by a step

Hi I have two arbitrarily dimensioned arrays, "a" and "b", with length(dim(a))==length(dim(b)). I want to form a sort of "corner-to-corner" version of abind(), or a multidimensional version of blockdiag(). In the case of matrices, the function is easy to write and if a=matrix(1,3,4) and b=matrix(2,2,2), then adiag(a,b) would return: [,1] [,2] [,3] [,4] [,5]

Dear R experts, I have a matrix (from some sort of classification) like this: object group [1,] 1 1 [2,] 2 2 [3,] 3 1 [4,] 4 1 [5,] 5 3 And I need something like this: [,1] [,2] [,3] [,4] [,5] [1,] 1 0 1 1 0 [2,] 0 1 0 0 0 [3,] 1 0 1 1 0 [4,] 1 0 1 1 0 [5,] 0 0 0 0 1 where all

Hi all - I am NEW to R and NEW to any type of programming. I am making heatmaps using the heatmap.2 function within gplots package. At present, when the heatmap is plotted it uses the row identifiers as 1,2,3,4...etc. However, I much rather use my own labels. I was told my another well-versed R programmer to use the follow script: x<-as.matrix(test1[,-1]) ## skip column 1 rownames(x)<-

Hi, I'm new about R. My problem is: I have 2 lists of proteins and I would compare them. I import the 2 lists like 2 different matrixs and I would that the first entry in the matrix 1 match with all entries of the matrix 2, then the second entry of the matrix 1 matchs with all entries of the matrix 2 and repeat this process with all entries of the matrix 1. Anybodies know can I do

Dear Dr.Goslee and anyone may intrested in matrix manipulate, I am using your ecodist to do mantel and partial mantel test, I have locality data and shape variation data, and the two distance matrixs are given as belowings. When I run the analysis, it is always report that the matrix is not square, but I didn't know what's wrong with my data. Would you please help me on this. I am quite

3.5 was already released, and we only merge bug fixes to release branches. Caio, do you really need it on 3.5 specifically or are you just asking if folks are interested to get it into the trunk? - Hans On Mon, Oct 27, 2014 at 3:33 AM, Yaron Keren <yaron.keren at gmail.com> wrote: > +Hans > > Caio, even if this patch does not go into 3.5, maybe you'd like to post it > on

Hello there! I'm intending to port LLVM old JIT API for PPC64le once I wiil need this support in one of my applications. Would a port of this old piece of code be merged on LLVM 3.5, even though it is on the verge of obsoleteness? Thank you! -------------- next part -------------- An HTML attachment was scrubbed... URL:

Hi all, I am using the function "plotCI" with the following command: plotCI(m.residuos.p.2 [1:41],li=m.residuos.p.3 [1:41],ui=m.residuos.p.4 [1:41],lty=1,ylab="") This generates exactly what I want except for the fact that I wanna drawn a line linking the points (m.residuos). How could I do that? Thanks a lot in advance, Caio [[alternative HTML version deleted]]

Dear R-users, I would like to know what are the commands to install (from a local zip file) a package and then to load it. Thaks all, Bests, Caio "Perco a consciencia, mas n併o importa, encontro a maior serenidade na alucina佺ao. 伾 curioso como n併o sei dizer quem sou. Quer dizer, sei-o bem, mas n併o posso dizer. Sobretudo tenho medo de dizer, porque no momento em que tento falar, n併o

Hello everyone, Does anyone know what would cause the skewness() function (from e1071), as well as skew() from psych, to return a value of NaN? I have a vector of positively-skewed data (https://docs.google.com/file/d/0B6-m45Jvl3ZmYzlHRVRHRURzbVk/edit?usp=sharing) which these functions return a value for like normal: > skewness( data ) # returns 1.400405 but when I instead give those

I'm exploring pedigree() in lme4 and found some strange results trying to multiply a dtCMatrix by a diag() matrix using %*%. Studying the Matrix documentation, I stumbled upon the following error in the example on the ddiMatrix-class {Matrix} help page. I wonder if my problem could be related to this. > matrix(cbind(1, 2:4),3,2) %*% diag(c(10,1)) # OK [,1] [,2] [1,] 10 2 [2,]

Hi, is there any direct R function to write an diagonal matrix in an opposite way? for example I want to get like:   > diag(rnorm(5))[,5:1]            [,1]       [,2]       [,3]      [,4]       [,5] [1,]  0.0000000  0.0000000  0.0000000  0.000000 -0.1504687 [2,]  0.0000000  0.0000000  0.0000000 -2.139669  0.0000000 [3,]  0.0000000  0.0000000 -0.2102133  0.000000  0.0000000 [4,]  0.0000000

Hi, I have the following problem with ggplot2: When I specify black contours for bars in a barplot, ggplot automatically shows diagonal lines in the legend boxes. Is there a way, to remove these diagonal lines? Here's a simple example: library(ggplot2) df<-data.frame(x=gl(6,1),y=rnorm(6,100,50)) ggplot(df,aes(y=y,x=x,fill=x))+geom_bar(aes(group=x),colour='black',size=1)+