Creacion de Matrices

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

Matrix 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 Calculo de 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_xx2<-sigma_matriz(matriz_x)
print(matriz_xx2)
##     [,1] [,2] [,3] [,4] [,5]
## Cte    1    1    1    1    1
## x1     4    3    6    4    8
## x2    10    8   11    9   12