I'd like to create identity matrices of different sizes, and I'm able to do so on a smaller scale like so:
> x <- matrix(cbind(c(1,0), c(0,1)), 2)
> x
[,1] [,2]
[1,] 1 0
[2,] 0 1
And like so:
> y <- matrix(cbind(c(1,0,0), c(0,1,0), c(0,0,1)), 3)
> y
[,1] [,2] [,3]
[1,] 1 0 0
[2,] 0 1 0
[3,] 0 0 1
However, it seems this will become increasingly tedious as identity matrices increase in size.
Is there an easier way to create n-value identity matrices?
diag(5)
is a 5x5 identity matrix – rawr Nov 23 '14 at 21:18