Primero se crea una matriz auxiliar como la siguientecon las dimensiones deseadas
# dimension
n <- 10
# creamos matriz auxiliar con valores de 1 a n^2cualesquiera
auxiliar <- matrix(seq_len(n ^ 2)0, ncol = n)
>auxiliar
[,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10]
[1,] 1 11 21 31 41 51 61 71 81 91
[2,] 2 12 22 32 42 52 62 72 82 92
[3,] 3 13 23 33 43 53 63 73 83 93
[4,] 4 14 24 34 44 54 64 74 84 94
[5,] 5 15 25 35 45 55 65 75 85 95
[6,] 6 16 26 36 46 56 66 76 86 96
[7,] 7 17 27 37 47 57 67 77 87 97
[8,] 8 18 28 38 48 58 68 78 88 98
[9,] 9 19 29 39 49 59 69 79 89 99
[10,] 10 20 30 40 50 60 70 80 90 100nrow=n)
A continuación, se resta la matriz auxiliar consigo misma y se obtiene una matriz en la que hay valor -1 por encimaindica el numero de fila a la diagonal y 1 por debajomatriz que indica el número de columna.
Ambas matrices se generan a partir de la matriz auxiliar. Al restar se obtiene una matriz en la que hay valor -1 por encima de la diagonal y 1 por debajo.
# restamos las auxiliaresmatrices fila y columna obtenidas de la auxiliar
A<- row(auxiliar) - col(auxiliar)
# Ahora tenemos valor -1 por encima de la diagonal y valor 1 por debajo
> A
[,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10]
[1,] 0 -1 -2 -3 -4 -5 -6 -7 -8 -9
[2,] 1 0 -1 -2 -3 -4 -5 -6 -7 -8
[3,] 2 1 0 -1 -2 -3 -4 -5 -6 -7
[4,] 3 2 1 0 -1 -2 -3 -4 -5 -6
[5,] 4 3 2 1 0 -1 -2 -3 -4 -5
[6,] 5 4 3 2 1 0 -1 -2 -3 -4
[7,] 6 5 4 3 2 1 0 -1 -2 -3
[8,] 7 6 5 4 3 2 1 0 -1 -2
[9,] 8 7 6 5 4 3 2 1 0 -1
[10,] 9 8 7 6 5 4 3 2 1 0