Práctica de matrices en R y Publicación en Rpubs

Creación de matrices

Matriz Y

#Generació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

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")
#cte=constante
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

Resultado de matrices

Cálculo 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

creación de la 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