On Fri, Sep 4, 2009 at 11:25 AM, Markku
Karhunen<markku.karhu...@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
______________________________________________
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.
