diag                  package:base                  R Documentation

_M_a_t_r_i_x _D_i_a_g_o_n_a_l_s

_D_e_s_c_r_i_p_t_i_o_n:

     Extract or replace the diagonal of a matrix, or construct a
     diagonal matrix.

_U_s_a_g_e:

     diag(x, nrow, ncol)
     diag(x) <- value

_V_a_l_u_e:

     If `x' is a matrix then `diag(x)' returns the diagonal of `x'. 
     The resulting vector will have `names' if the matrix `x' has
     matching column and row names.

     If `x' is a vector (or a 1-d array) then `diag(x)' returns a
     diagonal matrix whose diagonal is `x'.  If `x' is an integer then
     `diag(x)' returns an identity matrix of order `x'. The dimension
     of the returned matrix can be specified by `nrow' and `ncol' (the
     default is square).

     The assignment form sets the diagonal of the matrix `x' to the
     given value(s).

_S_e_e _A_l_s_o:

     `matrix'.

_E_x_a_m_p_l_e_s:

     dim(diag(3))
     diag(10,3,4) # guess what?
     all(diag(1:3) == {m <- matrix(0,3,3); diag(m) <- 1:3; m})

     diag(var(M <- cbind(X=1:5, Y=rnorm(5))))#-> vector with names "X" and "Y"
     rownames(M) <- c(colnames(M),rep("",3));
     M; diag(M) #  named as well

