On Sun, May 19, 2013 at 10:33 AM, mary <mary.d...@libero.it> wrote: > > Sorry in the previuos message I've inverted: I would like generate ... > "same eigenvectors and different eigenvalues"
Try this: m0 <- 5 + matrix(c(.05, -.05, -.05, .05), 2) m0.values <- eigen(m0)$values with(eigen(matrix(rnorm(4), 2)), vectors %*% diag(m0.values) %*% solve(vectors)) ______________________________________________ R-help@r-project.org 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.