[R] < 0 x 0 matrix >

[Barry Rowlingson]

On Fri, Sep 4, 2009 at 11:25 AM, Markku
Karhunen<markku.karhunen at helsinki.fi> wrote:
> Hi,
>
> Does anybody know, what is going on here?
>
>> diag(sqrt(1))
>
>     [,1]
> [1,]    1
>>
>> diag(sqrt(0.3333))
>
> <0 x 0 matrix>
>>
>> sqrt(1)
>
> [1] 1
>>
>> sqrt(0.3333)
>
> [1] 0.5773214
>

 Read the help for diag yet?

    'diag' has four distinct usages:
...
        3.  'x' is a scalar (length-one vector) and the only argument
           it a square identity matrix of size given by the scalar.
...

So diag(0.1) becomes diag(0) which is a 0-size matrix. Try diag(2.4)

Barry