##Numeral 1
options(scipen = 99999999)
load("C:/Users/Molin/Downloads/modelo_ventas.RData")
#Construyendo matriz A
Matriz_X<-model.matrix(modelo_ventas)
Matriz_XX<-t(Matriz_X)%*%Matriz_X
Matriz_A<-solve(Matriz_XX)%*%t(Matriz_X)
#Construyendo matriz P
Matriz_P<-(Matriz_X)%*%(Matriz_A)
#Constuyendo matriz M
n<-nrow(Matriz_X)
Matriz_M<-diag(n)-Matriz_P##Numeral 2
library(magrittr)## Warning: package 'magrittr' was built under R version 4.2.3Residuos_modelo_ventas<-modelo_ventas$residuals
Datos_modelo<-modelo_ventas$model
Residuos_matrices<-Matriz_M%*%Datos_modelo$ventas
cbind(Residuos_matrices,Residuos_modelo_ventas, Residuos_modelo_ventas-Residuos_matrices)%>% round(digits = 2)%>% as.data.frame()->comparacion
names(comparacion)<-c("Por matrices", "En el modelo", "Diferencia")
comparacion## Por matrices En el modelo Diferencia
## 1 -15.93 -15.93 0
## 2 19.33 19.33 0
## 3 38.02 38.02 0
## 4 -15.43 -15.43 0
## 5 5.16 5.16 0
## 6 80.22 80.22 0
## 7 -16.35 -16.35 0
## 8 -22.89 -22.89 0
## 9 -34.40 -34.40 0
## 10 46.09 46.09 0
## 11 -40.56 -40.56 0
## 12 9.25 9.25 0
## 13 5.82 5.82 0
## 14 -19.64 -19.64 0
## 15 -2.74 -2.74 0
## 16 -20.58 -20.58 0
## 17 -4.89 -4.89 0
## 18 -0.83 -0.83 0
## 19 -36.61 -36.61 0
## 20 -8.11 -8.11 0
## 21 11.00 11.00 0
## 22 56.50 56.50 0
## 23 -7.55 -7.55 0
## 24 31.71 31.71 0
## 25 -40.48 -40.48 0
## 26 86.80 86.80 0
## 27 -1.44 -1.44 0
## 28 35.92 35.92 0
## 29 24.51 24.51 0
## 30 -38.94 -38.94 0
## 31 25.44 25.44 0
## 32 -17.24 -17.24 0
## 33 -44.11 -44.11 0
## 34 47.62 47.62 0
## 35 -38.06 -38.06 0
## 36 104.56 104.56 0
## 37 -22.01 -22.01 0
## 38 6.26 6.26 0
## 39 -20.85 -20.85 0
## 40 -8.14 -8.14 0
## 41 3.52 3.52 0
## 42 7.13 7.13 0
## 43 39.33 39.33 0
## 44 29.22 29.22 0
## 45 -13.40 -13.40 0
## 46 8.21 8.21 0
## 47 -42.24 -42.24 0
## 48 -15.16 -15.16 0
## 49 36.15 36.15 0
## 50 -42.33 -42.33 0
## 51 29.60 29.60 0
## 52 -27.94 -27.94 0
## 53 -25.78 -25.78 0
## 54 -16.15 -16.15 0
## 55 20.52 20.52 0
## 56 -34.10 -34.10 0
## 57 32.84 32.84 0
## 58 -8.84 -8.84 0
## 59 -25.90 -25.90 0
## 60 4.38 4.38 0
## 61 -56.13 -56.13 0
## 62 -5.97 -5.97 0
## 63 34.06 34.06 0
## 64 -22.00 -22.00 0
## 65 -22.92 -22.92 0
## 66 -33.91 -33.91 0
## 67 -25.98 -25.98 0
## 68 -24.10 -24.10 0
## 69 15.18 15.18 0
## 70 -12.96 -12.96 0
## 71 -2.78 -2.78 0
## 72 -39.34 -39.34 0
## 73 9.10 9.10 0
## 74 -23.58 -23.58 0
## 75 16.53 16.53 0
## 76 41.92 41.92 0
## 77 -61.13 -61.13 0
## 78 -11.80 -11.80 0
## 79 42.30 42.30 0
## 80 -28.28 -28.28 0
## 81 -18.10 -18.10 0
## 82 55.90 55.90 0
## 83 -32.84 -32.84 0
## 84 1.99 1.99 0
## 85 -7.85 -7.85 0
## 86 2.06 2.06 0
## 87 -18.80 -18.80 0
## 88 -12.42 -12.42 0
## 89 -34.44 -34.44 0
## 90 -0.71 -0.71 0
## 91 -19.86 -19.86 0
## 92 -71.04 -71.04 0
## 93 5.47 5.47 0
## 94 -1.83 -1.83 0
## 95 -22.68 -22.68 0
## 96 -9.60 -9.60 0
## 97 30.46 30.46 0
## 98 2.80 2.80 0
## 99 -1.39 -1.39 0
## 100 -6.92 -6.92 0
## 101 47.13 47.13 0
## 102 13.67 13.67 0
## 103 74.24 74.24 0
## 104 8.29 8.29 0
## 105 5.20 5.20 0
## 106 -23.16 -23.16 0
## 107 -39.16 -39.16 0
## 108 -36.55 -36.55 0
## 109 -51.39 -51.39 0
## 110 14.93 14.93 0
## 111 31.77 31.77 0
## 112 0.85 0.85 0
## 113 3.55 3.55 0
## 114 20.05 20.05 0
## 115 1.25 1.25 0
## 116 -6.04 -6.04 0
## 117 -5.24 -5.24 0
## 118 -59.49 -59.49 0
## 119 -9.08 -9.08 0
## 120 -15.33 -15.33 0
## 121 -22.74 -22.74 0
## 122 -6.98 -6.98 0
## 123 52.60 52.60 0
## 124 -7.38 -7.38 0
## 125 6.99 6.99 0
## 126 -36.29 -36.29 0
## 127 51.55 51.55 0
## 128 -45.71 -45.71 0
## 129 -38.99 -38.99 0
## 130 -49.74 -49.74 0
## 131 137.77 137.77 0
## 132 68.15 68.15 0
## 133 29.49 29.49 0
## 134 4.04 4.04 0
## 135 2.44 2.44 0
## 136 24.96 24.96 0
## 137 14.83 14.83 0
## 138 18.65 18.65 0
## 139 -13.45 -13.45 0
## 140 -17.13 -17.13 0
## 141 -36.00 -36.00 0
## 142 -7.79 -7.79 0
## 143 -5.46 -5.46 0
## 144 -38.02 -38.02 0
## 145 -33.90 -33.90 0
## 146 -7.02 -7.02 0
## 147 57.02 57.02 0
## 148 -24.42 -24.42 0
## 149 5.76 5.76 0
## 150 -21.33 -21.33 0
## 151 60.96 60.96 0
## 152 -37.40 -37.40 0
## 153 3.70 3.70 0
## 154 -11.53 -11.53 0
## 155 5.57 5.57 0
## 156 20.34 20.34 0
## 157 -6.39 -6.39 0
## 158 0.10 0.10 0
## 159 35.21 35.21 0
## 160 -12.86 -12.86 0
## 161 -0.28 -0.28 0
## 162 -5.49 -5.49 0
## 163 7.07 7.07 0
## 164 -11.48 -11.48 0
## 165 -17.07 -17.07 0
## 166 43.48 43.48 0
## 167 30.49 30.49 0
## 168 32.76 32.76 0
## 169 9.33 9.33 0
## 170 79.25 79.25 0
## 171 -32.15 -32.15 0
## 172 -2.32 -2.32 0
## 173 -19.07 -19.07 0
## 174 11.45 11.45 0
## 175 56.81 56.81 0
## 176 -20.90 -20.90 0
## 177 10.72 10.72 0
## 178 10.72 10.72 0
## 179 99.16 99.16 0
## 180 0.68 0.68 0
## 181 8.07 8.07 0
## 182 43.24 43.24 0
## 183 -53.93 -53.93 0
## 184 -14.41 -14.41 0
## 185 33.06 33.06 0
## 186 -26.80 -26.80 0
## 187 -11.78 -11.78 0
## 188 2.02 2.02 0
## 189 75.55 75.55 0
## 190 -31.29 -31.29 0
## 191 11.46 11.46 0
## 192 -34.64 -34.64 0
## 193 -47.21 -47.21 0
## 194 -20.94 -20.94 0
## 195 -16.29 -16.29 0
## 196 -54.55 -54.55 0
## 197 -31.98 -31.98 0
## 198 8.34 8.34 0
## 199 -9.58 -9.58 0
## 200 49.59 49.59 0##Numeral 3
eigen(x=Matriz_XX)->descomposicion
auto_valores<-descomposicion$values
print(auto_valores)## [1] 311421698.6388 70252.5341 40973.4590 3714.3627 12.7735print(auto_valores>0)## [1] TRUE TRUE TRUE TRUE TRUE