Modelación y Simulación PEP 1.2 PDF
Modelación y Simulación PEP 1.2 PDF
Modelación y Simulación PEP 1.2 PDF
MODELACIÓN Y SIMULACIÓN
1.- Los clientes se encargan de evaluar los diseños preliminares de varios productos. En el
pasado, el 95 % de los productos que con mayor éxito en el mercado recibieron buenas
evaluaciones, el 60 % de los productos con éxito moderado recibieron buenas evaluaciones, y el
80% de los productos con baja aceptación recibieron malas evaluaciones. Además, el 40 % de
los productos ha tenido mucho éxito, el 35 % un éxito moderado, y el 25 % una baja aceptación.
Determine:
a.- ¿Cuál es la probabilidad de que un producto obtenga una buena evaluación?
b.- Si un nuevo diseño obtiene una buena evaluación, ¿cuál es la probabilidad de que se convierta en
un producto de gran éxito?
c.- Si un producto no tiene una buena evaluación ¿Cuál es la probabilidad de que se convierta en un
producto de gran éxito?
a) Solución: Sean:
𝐴 = “Producto con mayor éxito” 𝐸 = “Producto con Buenas Evaluaciones”
𝐵 = “Producto con éxito moderado” 𝐸 𝐶 = “Producto con Malas Evaluaciones”
𝐶 = “Producto con baja aceptación”
0.40 0.05 EC
0.60 E
0.35
B
A C
0.40 E
0.25 0.20 E
C
C
0.80 E
Por lo que las probabilidades que entrega el problema están dadas por:
𝑃 (𝐴) = 0.40 𝑃 (𝐵) = 0.35 𝑃 (𝐶 ) = 0.25
𝑃(𝐸/𝐴) = 0.95 𝑃(𝐸/𝐵) = 0.60 𝑃(𝐸 𝐶 /𝐶) = 0.80
c) Solución:
𝑃(𝐸 𝐶 /𝐴) ∙ 𝑃 (𝐴) 𝑃(𝐸 𝐶 /𝐴) ∙ 𝑃 (𝐴) 0.05 ∙ 0.40
𝑃(𝐴/𝐸 𝐶 ) = = = = 0.056
𝑃 (𝐸 𝐶 ) 1 − 𝑃 (𝐸 ) 1 − 0.64
2.- Tres líneas de producción de ánodos, A, B, C producen el 40%, 35%, 25%, respectivamente,
del total de la producción. Los porcentajes de ánodos defectuosos de estas líneas son del 2%,
3%, 4%. Determinar:
a) La probabilidad de que sea defectuoso un ánodo al azar.
b) Al tomar al azar un ánodo, y resulta ser defectuoso; calcule la probabilidad de haber sido
producido en la línea B.
c) ¿Qué línea tiene la mayor probabilidad de haber producido el ánodo defectuoso?
a) Solución: Sean:
𝐴 = “Línea de producción de ánodos A” 𝐷 = “Ánodos defectuosos”
𝐵 = “Línea de producción de ánodos B” 𝐷𝐶 = “Ánodos sin defectos”
𝐶 = “Línea de producción de ánodos C”
0.40 0.05 DC
0.60 D
0.35
B
A C
0.40 D
0.25 0.20 D
C
C
0.80 D
Por lo que las probabilidades que entrega el problema están dadas por:
𝑃 (𝐴) = 0.40 𝑃 (𝐵) = 0.35 𝑃 (𝐶 ) = 0.25
𝑃(𝐷/𝐴) = 0.02 𝑃(𝐷/𝐵) = 0.03 𝑃(𝐷/𝐶) = 0.04
c) Solución:
𝑃(𝐷/𝐴) ∙ 𝑃 (𝐴) 0.02 ∙ 0,40
𝑃(𝐴/𝐷) = = = 0.280
𝑃 (𝐷) 0.0285
Respuesta: La línea con mayor probabilidad de haber producido el ánodo defectuoso corresponde a la
B.
3.- Tras un test de cultura general se observa que las puntuaciones obtenidas siguen una
distribución N(65; 182). Se desea clasificar a los examinados en tres grupos (de baja cultura
general, de cultura general aceptable y de excelente cultura general) de modo que haya en el
primero un 20 % de la población, un 65 % en el segundo y un 15 % en el tercero. ¿Cuáles han de
ser las puntuaciones que marcan el paso de un grupo al otro?
Respuesta: El puntaje que marca el paso de baja cultura general a cultura general aceptable es de 48.2
puntos, y el paso de cultura general aceptable a excelente cultura general son 71.93 puntos.
Es decir que las piezas que cumplen con las especificaciones corresponden al 95.44%. Luego por
propiedad de complemento es posible determinar la cantidad de piezas que no cumplen las
especificaciones, quedando de la siguiente forma:
%𝑃(𝐶 𝑐 ) = 100% − %𝑃 (𝐶 ) = 100% − 95.44% = 4.56%
Respuesta: El 4.56% de las piezas se descartarán por no cumplir con las especificaciones establecidas
por el comprador.
5.- Se utilizan medidores para rechazar todos los componentes donde cierta dimensión no está
dentro de la especificación 1,50 +/-d .Se sabe que esta medición se distribuye de forma normal
con media 1,50 y desviación estándar 0,2. Determinar los valores de x que cubran el 95 % de las
mediciones
Se sabe que la distribución normal es simétrica, por lo que el 95% de las mediciones se encontrarán
entre el 2.5% y 97.5% de los datos, lo que es posible apreciarlo a continuación:
𝑑 𝑑
95%
2.5% 97.5%
𝜇 = 1.5
Luego se procede a determinar el valor de 𝑑:
(1.5 + 𝑑 ) − 1.5
𝑃(𝑥 < 1.5 + 𝑑) = 𝑃 (𝑧 < ) = 0.975
0.2
(1.5 + 𝑑) − 1.5
= 1.96 → 𝑑 = 0.392
0.2
Respuesta: Los componentes que cumplen con las especificaciones se encuentran entre 1.108 y 1.892.
6.- Los nuevos contenedores de emulsión para explosivo deben contener 72 lt. La fábrica ajusta
la máquina llenadora para que los recipientes contengan 72.6 lt (este valor se denomina set point).
Si la distribución del peso de llenado es aproximadamente normal con una desviación estándar
de 0.4 lt.
6) Solución: Sea: 𝑥 = “Peso de llenado de emulsión para explosivos”; 𝑥 ~ 𝑁(𝜇 = 7.2 ; 𝜎 2 = 0.42 )
72 − 72.6
𝑃(𝑥 < 72) = 𝑃 (𝑧 < ) = 𝑧(−1.5) = 0.0668
0.4
Respuesta: De un lote de 100 000 contenedores, 6890 unidades contendrán menos de 72 litros.
b) Solución: Lo primero será determinar la proporción entre los 300 tarros y los 100 00 contenedores, lo
que por principio de probabilidad corresponde a:
300
𝑝= = 0.003
100 000
72 − 𝜇
= −2.75 → 𝜇 = 73.1
0.4
Respuesta: El valor de media que debería usarse como set point para que no más de 300 tarros en
100000 contengan menos de 72 lt es 73.1 lt.
7.- En una faena minera se realiza una tronadura al día durante el mes de noviembre, la media de
tiros por tronadura es de 23 tiros, con una desviación estándar de 5, suponiendo se distribuyen
de manera normal, Calcule el número de tronaduras al mes en los que se espera quemar entre 21
y 27 tiros.
21 − 23 27 − 23
𝑃(21 < 𝑥 < 27) = 𝑃 ( <𝑧< ) = 𝑃 (−0.40 < 𝑧 < 0.80) = 𝑧(0.80) − 𝑧(−0.40)
5 5
= 0.7881 − 0.3446 = 0.4435
Respuesta: El número de tronaduras al mes de noviembre que se esperan quemar entre 21 y 27 tiros,
corresponde a 13 tronaduras.
X Y
0 9.1
1 7.3
2 3.2
3 4.6
4 4.8
5 2.9
6 5.7
7 7.1
8 8.8
9 10.2
Solución: Lo primero que se debe realizar es calcular las sumatorias necesarias para establecer el
sistema de ecuaciones expresado a continuación, el cual se determina previamente mediante método
de los mínimos cuadrados:
1) 𝑛 ∙ 𝑏𝑜 + ∑ 𝑥 ∙ 𝑏1 + ∑ 𝑥 2 ∙ 𝑏2 = ∑ 𝑦
2) ∑ 𝑥 ∙ 𝑏𝑜 + ∑ 𝑥 2 ∙ 𝑏1 + ∑ 𝑥 3 ∙ 𝑏2 = ∑ 𝑥𝑦
3) ∑ 𝑥 2 ∙ 𝑏𝑜 + ∑ 𝑥 3 ∙ 𝑏1 + ∑ 𝑥 4 ∙ 𝑏2 = ∑ 𝑥 2 𝑦
En seguida se procede a reemplazar los valores obtenidos en el sistema de ecuaciones antes descrito:
1) 10 ∙ 𝑏𝑜 + 45 ∙ 𝑏1 + 285 ∙ 𝑏2 = 63.7
2) 45 ∙ 𝑏𝑜 + 285 ∙ 𝑏1 + 2025 ∙ 𝑏2 = 307.3
3) 285 ∙ 𝑏𝑜 + 2025 ∙ 𝑏1 + 15333 ∙ 𝑏2 = 2153.3
8.- Se tiene la siguiente tabla de frecuencias que representa el tonelaje total extraído de una mina
expresada en miles de toneladas
a) Dibuje el Histograma
b) Determine si se trata de una normal aplicando el test de Chi- cuadrado
a) Solución: Lo primero que debe realizarse es determinar la marca de clase, la que corresponde a la
semisuma del límite inferior con el límite superior de cada intervalo
Intervalo de Clase 𝑛𝑖 𝑥𝑖
170 – 180 3 175
180 – 190 7 185
190 – 200 13 195
200 – 210 8 205
210 – 220 5 215
220 – 230 4 225
En seguida se grafica dicha información. Cabe destacar que todos los intervalos poseen el mismo ancho,
por lo que no es necesario recalcular las alturas de cada intervalo.
Histograma
14
12
10
Frecuencia
8
6
Series1
4
2
0
175 185 195 205 215 225
Marca de Clase
b) Solución: El primer paso a realizar es determinar las marcas de clases para poder ingresar los datos
a la calculadora
𝐿𝐼 − 𝜇 𝐿𝑆 − 𝜇
𝑍1 = 𝑦 𝑍2 =
𝜎 𝜎
Donde LI es el límite inferior de cada intervalo, LS es el límite superior de cada intervalo, 𝜇 es la media
muestral, 𝜎 es la desviación estándar muestral.
Donde 𝑁(𝑍1 ) corresponde a la normalización del límite inferior, y 𝑁(𝑍2 ) a la normalización del límite
superior.
𝐸𝑗 = Á𝑟𝑒𝑎 2 − Á𝑟𝑒𝑎 1
𝑛𝑖
∙ 100 𝑂𝑗 =
𝑛
Donde 𝑂𝑗 es la frecuencia observada que presentará la muestra, 𝑛𝑖 la frecuencia absoluta de la muestra,
y 𝑛 el tamaño de la muestra
𝑛 2
2 (𝐸𝑗 − 𝑂𝑗 )
𝑜𝑏𝑠 =∑
𝐸𝑗
𝑖=1
2
Y 𝑜𝑏𝑠 es estadístico de prueba del test Chicuadrado. Posteriormente se procede a determinar los
valores antes señalados, lo que se muestra en la siguiente tabla:
2
(𝐸𝑗 − 𝑂𝑗 )
𝐿𝐼 𝐿𝑆 𝑍1 𝑍2 Á𝑟𝑒𝑎 1 Á𝑟𝑒𝑎 2 𝐸𝑗 𝑛𝑖 𝑂𝑗
𝐸𝑗
170 180 -2.10 -1.38 1.79 8.38 6.59 3 7.5 0.1248
180 190 -1.38 -0.66 8.38 25.46 17.08 7 17.5 0.0102
190 200 -0.66 0.05 25.46 51.99 26.53 13 32.5 1.3428
200 210 0.05 0.77 51.99 77.94 25.94 8 20 1.3607
210 220 0.77 1.49 77.94 93.19 15.25 5 12.5 0.4971
220 230 1.49 2.21 93.19 98.64 5.46 4 10 3.7846
SUMA 7.1201
Donde los grados de libertad están dados por = 𝑘 − − 1, con 𝑘 es el número de intervalos y el
número de estimadores. Para este caso 𝑘 = 6 y el número de estimadores corresponde a dos, ya que
fue estimada la desviación estándar y media.
2 2
Luego como 0,95;3 = 14.067 y 0,05;7 = 2.167, se cumple:
2 2 2
0,95;7 = 7.815 > 0,05;7 = 0.352 > 𝑜𝑏𝑠 = 7.1201
Por lo tanto, los datos presentan un ajuste bueno a la distribución normal con media 𝜇 y desviación
estándar de 𝜎.
9.- En un yacimiento de fierro se muestrearon los pozos de tronadura registrando las siguientes
leyes:
Preguntas
a) Demuestre por el test de Chi cuadrado que es una normal
b) Determine el área un 10% y 90 %
a) Solución: El primer paso a realizar es determinar las marcas de clases para poder ingresar los datos
a la calculadora
𝐿𝐼 − 𝜇 𝐿𝑆 − 𝜇
𝑍1 = 𝑦 𝑍2 =
𝜎 𝜎
Donde LI es el límite inferior de cada intervalo, LS es el límite superior de cada intervalo, 𝜇 es la media
muestral, 𝜎 es la desviación estándar muestral.
Donde 𝑁(𝑍1 ) corresponde a la normalización del límite inferior, y 𝑁(𝑍2 ) a la normalización del límite
superior.
𝐸𝑗 = Á𝑟𝑒𝑎 2 − Á𝑟𝑒𝑎 1
2
(𝐸𝑗 − 𝑂𝑗 )
LI LS 𝑍1 𝑍2 Á𝑟𝑒𝑎 1 Á𝑟𝑒𝑎 2 𝐸𝑗 𝑛𝑖 𝑂𝑗
𝐸𝑗
27.5 31.5 -3.43 -2.75 0.03 0.30 0.27 1 0.1 0.1070
31.5 35.5 -2.75 -2.07 0.30 1.92 1.62 14 1.4 0.0299
35.5 39.5 -2.07 -1.39 1.92 8.23 6.31 56 5.6 0.0799
39.5 43.5 -1.39 -0.72 8.23 23.58 15.35 172 17.2 0.2230
43.5 47.5 -0.72 -0.04 23.58 48.40 24.82 245 24.5 0.0041
47.5 51.5 -0.04 0.64 48.40 73.89 25.49 263 26.3 0.0257
51.5 55.5 0.64 1.32 73.89 90.66 16.77 156 15.6 0.0816
55.5 59.5 1.32 2.00 90.66 97.72 7.06 67 6.7 0.0184
59.5 63.5 2.00 2.68 97.72 99.63 1.91 23 2.3 0.0796
63.5 67.5 2.68 3.36 99.63 99.96 0.33 3 0.3 0.0027
SUMA 0.6520
Donde los grados de libertad están dados por = 𝑘 − − 1, con 𝑘 es el número de intervalos y el
número de estimadores. Para este caso 𝑘 = 10 y el número de estimadores corresponde a dos, ya que
fue estimada la desviación estándar y media.
2 2
Luego como 0,95;7 = 14.067 y 0,05;7 = 2.167, se cumple:
2 2 2
0,95;7 = 14.067 > 0,05;7 = 2.167 > 𝑜𝑏𝑠 = 0.6520
Por lo tanto, los datos presentan un ajuste muy bueno a la distribución normal con media 𝜇 y desviación
estándar de 𝜎.
b) Solución:
Área z 𝑥 =𝑧∙𝜎+𝜇
10% 0.1 48.3009
20% 0.2 48.8898
30% 0.3 49.4787
40% 0.4 50.0676
50% 0.5 50.6565
60% 0.6 51.2454
70% 0.7 51.8343
80% 0.8 52.4232
90% 0.9 53.0121
Nota: Corresponde a la respuesta planteada por el profesor
Solución: Lo primero que se debe realizar es determinar el factor de utilización del sistema (𝜌), el que
se encuentra dado por:
𝜆 4
𝜌= = = 0.4
𝜇 10
Esto significa que en promedio el 40% del sistema permanece ocupado. En base a esto es posible
determinar la probabilidad que ninguna unidad en el sistema este ocupado (𝑃0 ), lo cual se determina a
partir del complemento, como ve a continuación:
𝜆 4
𝑃0 = 1 − = 1 − = 0.6
𝜇 10
Esto significa que en promedio el 60% del sistema permanece desocupado.
Luego se determinar el número promedio de unidades en el sistema en estudio (𝐿𝑆 ), cuya fórmula es:
𝜆 4
𝐿𝑆 = = = 0.67
𝜇 − 𝜆 10 − 4
En seguida se determina el tiempo promedio que una unidad permanezca en el sistema (𝑊𝑆 ), en horas,
siendo este:
1 1
𝑊𝑆 = = = 0.167
𝜇 − 𝜆 10 − 4
Finalmente el tiempo promedio que una unidad espera en la cola (𝑊𝑞 ), en horas, es:
𝜆
𝑊𝑞 = = 𝜌 ∙ 𝑊𝑆 = 0.4 ∙ 0.167 = 0.0668
𝜇 (𝜇 − 𝜆)
Los valores determinados significan que el factor 50% del sistema se encuentra ocupa, por ende el 50%
restante se encontrará desocupado. La cantidad de autos promedio en la línea de espera es de 0.25
autos. En promedio la cantidad de autos en el sistema es de 0.75 autos. El tiempo de espera promedio
en la línea de espera antes de ser atendido un auto es de 3.75 minutos. En promedio un auto espera en
el sistema 11.25 minutos antes de ser atendido
12.-Una máquina fotocopiadora es utilizada por 3 secretarias de una oficina para obtener las
copias que su sección requiere; como la magnitud del trabajo difiere de acuerdo al número de
copias que cada quien traiga, se hizo un análisis el cual dejó concluir que la máquina tiende a un
proceso de Poisson con un promedio de 8 trabajos por hora. Los requerimientos de utilización
son también aleatorios de acuerdo a un proceso Poissoniano con una tasa media de 5 trabajos
por hora. Calcular las características de utilización de la fotocopiadora.
12) Solución: Lo primero será identificar que se trata de un sistema multicanal o M/M/S, el cual
presenta la siguiente fórmula para la determinación de la probabilidad que el sistema este ocioso es:
1 1
𝑃0 = 𝑛 = 0 1 = 0.1496
∑𝑘𝑛=0 [
𝑘!
(
𝜆
) ] 3! 5 3! 5 2! 5 2 3! 5 3
(𝑘 − 𝑛)! 𝜇 (
(3 − 0)! 8 ) + (
(3 − 1)! 8 ) + (
(3 − 2)! 8 ) + (
(3 − 3)! 8 )
Por lo que en promedio el tiempo gastado en la cola para llegar a obtener un servicio en la fotocopiadora
es de 9.36 minutos.Finalmente 𝑊𝑆 corresponde a:
𝑊𝑆 = 𝑊 − 𝑊𝑞 = 0.32 − 0.156 = 0.164 ℎ𝑟 = 9.84 𝑚𝑖𝑛
13.- Suponga que un cajero bancario puede atender a los clientes a una velocidad promedio de
diez clientes por hora (µ = 10). Además, suponga que los clientes llegan a la ventanilla del cajero
a una tasa promedio de 7 por hora (𝝀 = 7). Se considera que el patrón de arribo sigue la
Distribución de Poisson y el patrón de servicio sigue la Distribución exponencial.
Determine el % de clientes que perderá, debido a que los clientes se alejan del cajero siempre
que existan 3 o más clientes antes que ellos en el sistema.
13) Solución: Primeramente se determinará el factor de utilización del sistema, el cual se encuentra dado
por:
𝜆 7
𝜌= = = 0.7 = 70%
𝜇 10
Por lo que la probabilidad que ninguna unidad este en el sistema será:
𝜆 7
𝑃0 = 1 − = 1 − = 0.3 = 30%
𝜇 10
Luego se sabe que la probabilidad que existan tres o más clientes en el sistema se encuentra dado por:
𝑃(𝑥 ≥ 3) = 1 − [𝑃0 + 𝑃1 + 𝑃2 + 𝑃3 ] = 1 − [𝑃0 + (𝑃0 ∙ 𝜌1 ) + (𝑃0 ∙ 𝜌2 ) + (𝑃0 ∙ 𝜌3 )]
= 1 − [0.3 + (0.3 ∙ 0.71 ) + (0.3 ∙ 0.72 ) + (0.3 ∙ 0.73 )] = 1 − [0.21 + 0.147 + 0.1029] = 0.2401
Otra forma de determinar este valor es a partir de la fórmula de probabilidad que más de “k” clientes
estén en el sistema, la que presenta la siguiente forma:
𝜆 𝑘+1 7 3+1
𝑃𝑛>𝑘 = ( ) → 𝑃𝑛>3 = ( ) = 0.2401
𝜇 10
Finalmente el porcentaje de clientes que perderá, corresponde al 24.01%
14.- Un promedio de 10 automóviles por hora llegan a un cajero con un solo servidor que
proporciona servicio sin que uno descienda del automóvil. Suponga que el tiempo de servicio
promedio por cada cliente es 4 minutos, y que tanto los tiempos entre llegadas y los tiempos de
servicios son exponenciales. Conteste las preguntas siguientes:
a) Solución: El ejercicio dice que en promedio 10 automóviles llegan a un cajero por hora, por lo que al
transformarlo a minutos queda dado:
𝑐𝑙𝑖𝑒𝑛𝑡𝑒𝑠 1 ℎ𝑜𝑟𝑎
𝜆 = 10 ( ) = 1/6 𝑐𝑙𝑖𝑒𝑛𝑡𝑒𝑠/𝑚𝑖𝑛
ℎ𝑜𝑟𝑎 60 𝑚𝑖𝑛
Luego como el tiempo promedio por cada cliente es 4 minutos, por lo que 𝜇 toma este valor.
Para determinar lo pedido se debe determinar el factor de utilización, el que se encuentra dado por:
𝜆 1/6 2
𝜌= = = = 0.6667 = 66.67%
𝜇 1/4 3
Es decir la probabilidad que el cajero este de ocioso corresponderá a:
𝜆 2
𝑃0 = 1 − = 1 − = 0.3333 = 33.33%
𝜇 3
b) Solución: El número promedio de automóviles que están en la cola del cajero es:
𝜆 1/6 4
𝐿𝑞 = = = = 1.33 𝑎𝑢𝑡𝑜𝑠
𝜇 (𝜇 − 𝜆) 1/4 ∙ (1/4 − 1/6) 3
c) Solución: El tiempo promedio que un cliente pasa en el estacionamiento del banco (incluyendo el
tiempo de servicio) es:
1 1 1
𝑊𝑆 = = = = 12 𝑚𝑖𝑛
𝜇 − 𝜆 1/4 − 1/6 1/12
d) Solución: Si el cajero estuviera ocupado atendería en promedio 15 clientes por hora, ya que:
1 𝑐𝑙𝑖𝑒𝑛𝑡𝑒 60 𝑚𝑖𝑛
𝜇= ( ) = 15 𝑐𝑙𝑖𝑒𝑛𝑡𝑒𝑠/ℎ𝑟
4 𝑚𝑖𝑛 1 ℎ𝑟
Luego, la probabilidad que este ocupado es 0.6667, por lo que los clientes que se atenderán en promedio
por hora serán:
𝜌 ∙ 𝜇 = 0.6667 ∙ 15 = 10 𝑐𝑙𝑖𝑒𝑛𝑡𝑒𝑠
ANOVA 2k
15. Aplicando el diseño 2K resuelva el siguiente problema:
Un ingeniero mecánico estudia la rugosidad superficial de una pieza producida en una operación
de corte de metal. El interés recae en tres factores: la rapidez con la que se hace el corte (A), la
profundidad de éste (B), y el ángulo de la herramienta (C). A los tres factores se les asigna dos
niveles, y se realizan dos réplicas del diseño factorial. Los datos se presentan en la siguiente
tabla.
Profundidad del Corte B
0.025 pulgadas 0.040 pulgadas ángulo de la
Ángulo de la herramienta C herramienta C
Rapidez de corte (A) 15" 25" 15" 25"
12 13 24 15
20 pulg/min
13 17 17 18
15 18 20 20
30 pulg/min
18 22 23 22
Solución: Debido a que tiene tres factores de interés: A (Rapidez de Corte), B (Profundidad de Corte) y
C (Ángulo de la Herramienta), los que presentan dos niveles cada uno. Las ocho corridas o tratamientos
de diseño 23 se pueden representar geométricamente como un cubo similar al que se muestra en la
siguiente figura:
Cada arista del cubo corresponde a una corrida o combinación de tratamientos diferentes, las que se
representan por (1), a, b, ab, c, ac, bc y abc.
El efecto principal A puede estimarse promediando las cuatro combinaciones de tratamiento de la cara
derecha del cubo, donde el nivel A es alto, y después restando de esta cantidad el promedio de las
cuatro combinaciones de tratamientos que están en la cara izquierda del cubo A tiene el nivel bajo, por
lo que se tiene la siguiente notación:
𝑎 + 𝑎𝑏 + 𝑎𝑐 + 𝑎𝑏𝑐 (1) + 𝑏 + 𝑐 + 𝑏𝑐
𝐴 = 𝑦̅𝐴+ − 𝑦̅𝐴− = −
4𝑛 4𝑛
1
𝐴= [𝑎 + 𝑎𝑏 + 𝑎𝑐 + 𝑎𝑏𝑐 − (1) − 𝑏 − 𝑐 − 𝑏𝑐]
4𝑛
De manera similar, el efecto de B se puede determinar como la diferencia en promedios entre las cuatro
combinaciones de tratamientos de la cara posterior del cubo y las cuatro combinaciones de la cara
anterior, por lo que se tiene:
𝑏 + 𝑎𝑏 + 𝑏𝑐 + 𝑎𝑏𝑐 (1) + 𝑎 + 𝑐 + 𝑎𝑐
𝐵 = 𝑦̅𝐵+ − 𝑦̅𝐵− = −
4𝑛 4𝑛
1
𝐵= [𝑏 + 𝑎𝑏 + 𝑏𝑐 + 𝑎𝑏𝑐 − (1) − 𝑎 − 𝑐 − 𝑎𝑐]
4𝑛
𝑐 + 𝑎𝑐 + 𝑏𝑐 + 𝑎𝑏𝑐 (1) + 𝑎 + 𝑏 + 𝑎𝑏
𝐶 = 𝑦̅𝐶+ − 𝑦̅𝐶− = −
4𝑛 4𝑛
1
𝐶= [𝑐 + 𝑎𝑐 + 𝑏𝑐 + 𝑎𝑏𝑐 − (1) − 𝑎 − 𝑏 − 𝑎𝑏]
4𝑛
Posteriormente, los efectos de interacción también se pueden obtener con facilidad, donde la interacción
entre A y B se puede obtener como la diferencia entre los promedios de los efectos de A en los dos
niveles de B, es decir:
(𝑎𝑏𝑐 − 𝑏𝑐) + (𝑎𝑏 − 𝑏)
𝐵 𝑎𝑙𝑡𝑜 (+) → 𝐸𝑓𝑒𝑐𝑡𝑜 𝐴𝑃𝑎𝑟𝑎 𝐵+ =
2𝑛
(𝑎𝑐 − 𝑐) + (𝑎 − (1))
𝐵 𝑎𝑙𝑡𝑜 (−) → 𝐸𝑓𝑒𝑐𝑡𝑜 𝐴𝑃𝑎𝑟𝑎 𝐵− =
2𝑛
En la siguiente figura se muestra de forma gráfica los términos involucrados en el cálculo de los efectos
principales y de interacción en un diseño factorial 23
La figura anteriormente expuesta, muestra los vértices del cubo (tratamientos) que deben ir con signo
positivo o negativo para la determinación de cada efecto, los cuales concuerdan con los obtenidos con
anterioridad.
A pesar que puede representarse geométricamente en un cubo, el cálculo de los contrastes y de los
efectos se pueden obtener de forma más sencilla a partir de la tabla de signos mostrada a continuación:
Tabla 1: Signos para el diseño 23
Tratamiento A B C AB AC BC ABC
1 (1) - - - + + + -
2 a + - - - - + +
3 b - + - - + - +
4 ab + + - + - - -
5 c - - + + + + -
6 ac + - + - - + +
7 bc - + + - + - +
8 abc + + + + - - -
Cabe destacar que la columna A cambia de signo por cada tratamiento, luego la columna B cambia de
signo por cada dos tratamientos y finalmente la columna C cambia de signo por cada cuatro tratamientos.
En seguida se pueden obtener los contrastes de interacción, multiplicando las columnas adecuadas de
los efectos principales.
La suma de cuadrados y efectos para la construcción del ANOVA se pueden obtener a partir de las
fórmulas:
𝐶𝑜𝑛𝑡𝑟𝑎𝑠𝑡𝑒
𝐸𝑓𝑒𝑐𝑡𝑜 =
4𝑛
(𝐶𝑜𝑛𝑠𝑡𝑟𝑎𝑠𝑡𝑒)2
𝑆𝑆 =
8𝑛
Donde la fórmula general corresponde a:
𝐶𝑜𝑛𝑡𝑟𝑎𝑠𝑡𝑒
𝐸𝑓𝑒𝑐𝑡𝑜 =
𝑛 2𝑘−1
(𝐶𝑜𝑛𝑠𝑡𝑟𝑎𝑠𝑡𝑒)2
𝑆𝑆 =
𝑛 2𝑘
Siguiendo con el ejercicio, se dispone a determinar las corridas, las que se encontrarán dadas por la
suma de los valores entregados. A modo de ejemplo se muestra el cálculo de la corrida (1)
(1) = 12 + 13 = 25
12 13 24 15
20 pulg/min
13 17 17 18
(1) = 25 c = 30 b = 41 bc = 33
15 18 20 20
30 pulg/min
18 22 23 22
a = 33 ac = 40 ab = 43 abc = 42
• ∑ ∑ 𝑦𝑖𝑗 2 = 122 + 132 + 132 + 172 + 242 + 172 + 152 + 182 + 152 + 182 + 182 + 222 + 202 + 232
𝑖=1 𝑗=1
+ 202 + 222 = 5355
𝑦•• 2 12 + 13 + 13 + 17 + 24 + 17 + 15 + 18 + 15 + 18 + 18 + 22 + 20 + 23 + 20 + 22 2872
• = =
𝑁 16 16
Por lo tanto
𝑞 𝑛
𝑦•• 2 2872
𝑆𝑆𝑇 = ∑ ∑ 𝑦𝑖𝑗 2 − = 5355 − = 206.9375
𝑁 16
𝑖=1 𝑗=1
Posteriormente se calcula la suma de cuadrados del error, la que se encuentra dada por:
𝑆𝑆𝐸 = 𝑆𝑆𝑇 − 𝑆𝑆𝑇𝑟𝑎𝑡𝑎𝑚𝑖𝑒𝑛𝑡𝑜𝑠 = 𝑆𝑆𝑇 − (𝑆𝑆𝑎 + 𝑆𝑆𝑏 + 𝑆𝑆𝑐 + 𝑆𝑆𝑎𝑏 + 𝑆𝑆𝑎𝑐 + 𝑆𝑆𝑏𝑐 + 𝑆𝑆𝑎𝑏𝑐 )
= 206.9375 − (52.5625 + 60.0625 + 0.5625 + 3.0625 + 5.0625 + 27.5625 + 1.5625) = 56.5000
El valor de “𝑎” es el número de subclasificaciones que presenta la rapidez de corte, las que corresponden
a 20 pulg/min y 30 pulg/min, es decir “𝑎” es igual a dos. En el caso de “𝑏” corresponde a las
subclasificaciones de la profundidad de corte, las que corresponden a 0.025 pulgadas y 0.040 pulgadas,
por lo que “𝑏” también toma el valor de dos.
El valor de “𝑛” corresponde a dos, debido a que la cantidad de datos por bloque es igual a dos. Al
reemplazar dichos valores se obtiene la siguiente tabla:
Luego se determina el factor observado, el cual se encuentra dado por la división entre la media de
cuadrados de cada factor y la media de cuadrados del error, lo que se muestra a continuación:
Fuente de Suma de Grados de Media de
F0
Variación Cuadrados Libertad Cuadrados
A 𝑆𝑆𝑎 = 52.5625 1 52.5625 52.5625/7.0625 = 7.44
B 𝑆𝑆𝑏 = 60.0625 1 60.0625 60.0625/7.0625 = 8.50
C 𝑆𝑆𝑐 = 0.5625 1 0.5625 0.5625/7.0625 = 0.80
AB 𝑆𝑆𝑎𝑏 = 3.0625 1 3.0625 3.0625/7.0625 = 0.43
AC 𝑆𝑆𝑎𝑐 = 5.0625 1 5.0625 5.0625/7.0625 = 0.72
BC 𝑆𝑆𝑏𝑐 = 27.5625 1 27.5625 27.5625/7.0625 = 3.90
ABC 𝑆𝑆𝑎𝑏𝑐 = 1.5625 1 1.5625 1.5625/7.0625 = 0.22
ERROR 𝑆𝑆𝐸 = 56.5000 8 7.0625
TOTAL 𝑆𝑆𝑇 = 206.9375 15
Finalmente se define el coeficiente de correlación (𝑟), el que se encuentra dado por la raíz cuadrada del
coeficiente de determinación (𝑟 2 ). Lo anterior se muestra en la siguiente tabla resumen:
Coeficiente de Porcentaje de Coeficiente de
Suma de Cuadrados
Determinación (𝑟2 ) Determinación (% 𝑟 2 ) Correlación (𝑟2 )
𝑆𝑆𝑎 = 52.5625 0.2540 25.40 0.5040
𝑆𝑆𝑏 = 60.0625 0.2902 29.02 0.5387
𝑆𝑆𝑐 = 0.5625 0.0027 0.27 0.0521
𝑆𝑆𝑎𝑏 = 3.0625 0.0148 1.48 0.1217
𝑆𝑆𝑎𝑐 = 5.0625 0.0245 2.45 0.1564
𝑆𝑆𝑏𝑐 = 27.5625 0.1332 13.32 0.3650
𝑆𝑆𝑎𝑏𝑐 = 1.5625 0.0076 0.76 0.0869
𝑆𝑆𝐸 = 56.5000 0.2730 27.30 0.5225
𝑆𝑆𝑇 = 206.9375
𝑆𝑆𝑆𝑢𝑚𝑎 = 150.4375 0.7270 72.70 0.8526
Reemplazando, se tiene:
3.625 3.875 0.375 0.875 1.125 2.625 0.625
𝑦 = 17.9375 + ∙ 𝑥𝑎 + ∙ 𝑥𝑏 + ∙ 𝑥𝑐 − ∙ 𝑥𝑎𝑏 + ∙ 𝑥𝑎𝑐 − ∙ 𝑥𝑏𝑐 + ∙ 𝑥𝑎𝑏𝑐
2 2 2 2 2 2 2
𝑦 = 17.9375 + 1.8125𝑥𝑎 + 1.9375𝑥𝑏 + 0.1875𝑥𝑐 − 0.4375𝑥𝑎𝑏 + 0.5625𝑥𝑎𝑐 − 1.3125𝑥𝑏𝑐 + 0.3125𝑥𝑎𝑏𝑐
16.- Preguntas
e) Si tiene dos grupos de muestras como determina la que tiene mayor dispersión
La forma de determinar que muestra presenta mayor dispersión, es necesario calcular el coeficiente
de variación, el cual está dado por la división de la desviación estándar por el promedio.
𝑆𝑥
𝐶𝑉(𝑥) =
𝑥̅
17.- En un yacimiento de Cobre que tiene 1 millón de TM, cuya distribución de leyes tiene los
siguientes datos:
Limite Superior Frecuencia %
1.22 2.5
1.27 5.0
1.32 5.0
1.37 10.0
1.42 15.0
1.47 20.0
1.52 12.5
1.57 10.0
1.62 5.0
1.67 7.5
1.72 2.5
1.80 5.0
a) Determine la calidad del mineral realizando la Curva Tonelaje Ley , para leyes de corte 1.3;
1.4; 1.5; 1.6, 1.7 determinando tonelajes y ley media
b) Determinar entre - y la ley corte anterior el tonelaje y la ley media, que representa el
lastre.
a) Solución: Lo primero que determinar las frecuencias acumuladas a partir de la información entregada
por el ejercicio, lo cual se muestra a continuación:
Luego se grafica en el papel probabilístico se deben disponer los valores de los límite superior en el eje
y, con la respectiva frecuencia acumulada en el eje x (puntos de color azul), para luego crear una recta
que siga la tendencia de los datos (línea de color rojo)
En seguida como se considera que la variable se distribuye de forma normal, se tiene que la mediana y
la media presentan el mismo valor. Sabiendo que la mediana es el valor que deja el 50% de los datos
debajo, al entrar al papel probabilístico con el 50% se obtiene que el promedio es de 1.47 (punto de color
rojo)
1.7
1.9
1.3
1.6
1.5
1.4
1.2
𝑉𝑎𝑙𝑜𝑟 𝑑𝑒𝑙 95% − 𝑉𝑎𝑙𝑜𝑟 𝐼𝑛𝑖𝑐𝑖𝑎𝑙 𝑉𝑎𝑙𝑜𝑟 𝑑𝑒𝑙 50% − 𝑉𝑎𝑙𝑜𝑟 𝐼𝑛𝑖𝑐𝑖𝑎𝑙
𝜎= ≈
4 2
Al reemplazar se obtiene:
1.47 − 1.22
𝜎≈ ≈ 0.13
2
En conclusión las leyes presentan una distribución normal con una media de 1.47% y desviación
estándar de 0.13%
Luego se plantearán las ecuaciones que serán utilizadas para determinar la ley media en la curva
tonelaje v/s ley:
𝑥−𝜇
𝑌𝑐 =
𝜎
Donde 𝑥 corresponde a la ley de corte, 𝜇 la media de las leyes obtenidas a partir del papel probabilístico,
y 𝜎 la desviación estándar obtenida de la misma forma.
%(+𝑋𝑐 ) = 1 − 𝑃(𝑌𝑐 )
Donde %(+𝑋𝑐 ) es la proporción sobre la ley de corte, y 𝑃(𝑌𝑐 ) es la probabilidad obtenida a partir del 𝑌𝑐
determinado anteriormente.
Finalmente para confeccionar la curva tonelaje v/s ley se utilizan los valores de las leyes de corte (𝑋𝑐 ),
Tonelaje sobre la ley de corte (𝑇𝑜𝑛(+𝑋𝑐 )) y la ley media sobre la ley de corte (𝜇 (+𝑋𝑐 )), lo cual se muestra
a continuación:
b) Solución: Para determinar el tonelaje (𝑇𝑜𝑛(−𝑋𝑐 )) y la ley media (𝜇(−𝑋𝑐 )) que representan el lastre se
deben realizar los mismos cálculos, pero considerando los valores bajo las leyes de corte, lo que se
muestra en la siguiente tabla:
18.- En un yacimiento de Cobre que tiene 1 millón de TM, cuya distribución de leyes tiene los
siguientes datos:
Limite Superior Frecuencia Acum%
31.5 0.1
35.5 1.5
39.5 7.1
43.5 24.3
47.5 48.8
51.5 75.1
55.5 90.7
59.5 97.4
63.5 99.7
67.5 100.0
Preguntas
a) Determine la calidad del mineral realizando la Curva Tonelaje Ley, para leyes de corte 35, 40,
45, 50, 55, 60 determinando tonelajes y ley media.
b) Determinar entre - y la ley corte anterior el tonelaje y la ley media, que representa el lastre.
a) Solución: Se grafica en el papel probabilístico los valores de los límite superior en el eje y, con la
respectiva frecuencia acumulada en el eje x (puntos de color azul), para luego crear una recta que siga
la tendencia de los datos (línea de color rojo)
En seguida como se considera que la variable se distribuye de forma normal, se tiene que la mediana y
la media presentan el mismo valor. Sabiendo que la mediana es el valor que deja el 50% de los datos
debajo, al entrar al papel probabilístico con el 50% se obtiene que el promedio es de 48 (punto de color
rojo)
𝑉𝑎𝑙𝑜𝑟 𝑑𝑒𝑙 95% − 𝑉𝑎𝑙𝑜𝑟 𝐼𝑛𝑖𝑐𝑖𝑎𝑙 𝑉𝑎𝑙𝑜𝑟 𝑑𝑒𝑙 50% − 𝑉𝑎𝑙𝑜𝑟 𝐼𝑛𝑖𝑐𝑖𝑎𝑙
𝜎= ≈
4 2
Al reemplazar se obtiene:
48 − 31.5
𝜎≈ ≈ 8.25
2
En conclusión las leyes presentan una distribución normal con una media de 48% y desviación estándar
de 8.25%
90
80
40
70
60
50
30
Finalmente para confeccionar la curva tonelaje v/s ley se utilizan los valores de las leyes de corte (𝑋𝑐 ),
Tonelaje sobre la ley de corte (𝑇𝑜𝑛(+𝑋𝑐 )) y la ley media sobre la ley de corte (𝜇 (+𝑋𝑐 )), lo cual se muestra
a continuación:
b) Solución: Para determinar el tonelaje (𝑇𝑜𝑛(−𝑋𝑐 )) y la ley media (𝜇(−𝑋𝑐 )) que representan el lastre se
deben realizar los mismos cálculos, pero considerando los valores bajo las leyes de corte, lo que se
muestra en la siguiente tabla: