mc19167, Jairo Rodrigo Mendez Carrillo
Matrix Y
#Generación de la Matriz y
Matriz_y<-matrix(data = c(30,20,36,24,40),
nrow = 5,ncol = 1,byrow = TRUE)
colnames(Matriz_y)<-c("y")
print(Matriz_y)
## y
## [1,] 30
## [2,] 20
## [3,] 36
## [4,] 24
## [5,] 40
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 se X’X (Sigma matriz)
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
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