[Rd] possible bug in diag()

Simon Urbanek simon.urbanek at math.uni-augsburg.de
Wed Sep 17 21:01:14 MEST 2003

On Wednesday, September 17, 2003, at 06:51 PM, Marsland, John wrote:

> 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:

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 


More information about the R-devel mailing list