Matriz Y
#Geración de la Matriz Y.
Matriz_Y <- matrix(data = c(20,30,36,24,40),
nrow =5,
ncol = 1,
byrow = TRUE)
colnames(Matriz_Y) <-c("Y")
print(Matriz_Y)
## Y
## [1,] 20
## [2,] 30
## [3,] 36
## [4,] 24
## [5,] 40
Matriz X
#Generación de la Matriz X.
Matriz_X <-cbind(rep(1,5),
matrix(data=c(4,10,3,8,6,11,4,9,8,12),
nrow = 5,
ncol = 2,
byrow = TRUE))
colnames(Matriz_X) <-c("Cte","X1","X2")
print(Matriz_X)
## Cte X1 X2
## [1,] 1 4 10
## [2,] 1 3 8
## [3,] 1 6 11
## [4,] 1 4 9
## [5,] 1 8 12
Producto de Matrices
Cálculo de X´X (Sigma Matriz)
#Generación de la Matriz Transpuesta.
#De manera manual.
Matriz_XX <-t(Matriz_X)%*%Matriz_X
print(Matriz_XX)
## Cte X1 X2
## Cte 5 25 50
## X1 25 141 262
## X2 50 262 510
#Generación de la Matriz Transpuesta.
#Apartir de la creación de una función.
Sigma_matriz <-function(matriz){t(matriz)%*%matriz}
Matriz_XX2 <-Sigma_matriz(Matriz_X)
print(Matriz_XX2)
## Cte X1 X2
## Cte 5 25 50
## X1 25 141 262
## X2 50 262 510