It concerns trival diagonal matrices: [...]diag(rnorm(1))<0 x 0 matrix>
diag(rnorm(1),nrow=1)[,1] [1,] 0.4843697
There's an obvious work around... but I thought it was worth notifying the
list.
I fail to see any reason for you posting this here, since this is a perfectly documented behavior - ?diag says:
Note: Using 'diag(x)' can have unexpected effects if 'x' is a vector that could be of length one. Use 'diag(x, nrow = length(x))' for consistent behaviour.
Hence you can safely assume that the "list" needs no notification of this.
Cheers, Simon
______________________________________________ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-devel