Search Archives

Yes, I would like every generated matrix to be drawn from a uniform distribution. Martin Maechler's solution was interesting but when I compute the product of the obtained matrix with any N-vector Y, the resulting vector is most...

Thanks, it's perfect for my needs. ----- Message d'origine ---- De : Martin Maechler <maechler at stat.math.ethz.ch> ? : Florent Bresson <f_bresson at yahoo.fr> Cc : Richard M. Heiberger <rmh at temple.edu>; r-help at stat.math.ethz...

I just have a try with sapply. The problem is that my function pbeta2 has two parameters z and p (wich is a vector of two parameters). If I use sapply , R returns a message incicating that parameter p is...

Thanks, it is exactly what I was looking for. Florent Bresson --- Gavin Simpson <gavin.simpson at ucl.ac.uk> a ?crit?: > On Wed, 2006-07-19 at 16:12 +0200, Florent Bresson > wrote: > > Hi, I have to use a loop...

Thanks, I think it's a shrewd solution, but the problem is that it cannot generate every N*N bistochastic matrix and every cell tends to 1/N as B tends to infinity Florent Bresson ----- Message d'origine ---- De : Dimitris...

Thanks, I tried someting like this, but computation takes times for large matrices btransf <- function(y,X=length(y)^4) { N <- length(y) bm <- matrix(rep(1/N,N^2),N,N) for(j in 1:X){ coord <- sample(1...

It works really faster with mapply. I just had to change my fonction from pbeta2(z,p) with p=c(p1,p2) to pbeta2(z,p1,p2). Thanks for the tip. --- Thomas Lumley <tlumley at u.washington.edu> a ??crit...

The problem was with the use of the integrate command for the definition of the cdf of the generalized beta of the second kind. I solved the problem with a transformation of the function qbeta. I then used an optimize...

Thanks Brian for your help. I am sorry, I should have red the help page more carefully. I am working on Kubuntu 8.04 with R 2.8.0. The problem with the file is that I do not know...