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 BR, Markku Karhunen researcher University of Helsinki ______________________________________________ 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.