possible bug in diag()

Why is this a bug ? help(diag) says

     If 'x' is a vector (or 1D array) of length two or more, then
     'diag(x)' returns a diagonal matrix whose diagonal is 'x'.

     If 'x' is a vector of length one then 'diag(x)' returns an
     identity matrix of order the nearest integer to 'x'.  The
     dimension of the returned matrix can be specified by 'nrow' and
     'ncol' (the default is square).

There's no way of knowing what rnorm(1) returned in your particular call, but 
my guess is that the integer "nearest" to it was 0.

It concerns trival diagonal matrices:

diag(1)

     [,1]
[1,]    1

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.

Regards,

John Marsland


**********************************************************************
This is a commercial communication from Commerzbank AG.\ \ T...{{dropped}}

______________________________________________
R-devel@stat.math.ethz.ch mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-devel