Ejercicios Diseño Multifactorial
Ejercicios Diseño Multifactorial
Ejercicios Diseño Multifactorial
Semestre 2023-1
Índice
Ejemplo 5.2…………………………………………………………………………………………………………………………………….2
Ejercicio 19……………………………………………………………………………………………………………………………………13
Ejercicio 20……………………………………………………………………………………………………………………………………19
Ejercicio 21……………………………………………………………………………………………………………………………………24
Ejercicio 23……………………………………………………………………………………………………………………………………30
Ejercicio 25……………………………………………………………………………………………………………………………………33
Ejercicio 26……………………………………………………………………………………………………………………………………38
Ejercicio 27……………………………………………………………………………………………………………………………………43
1
Ejemplo 5.2
El acabado (Y) está en unidades de gramos e interesa minimizar su valor. La representación gráfica
del diseño 4 x 3 se muestra en la figura 5.3.
2
Para la resolución del problema primero se realizó la suma de cada valor de respuesta de las
pruebas Ai - Bj junto con cada una de sus repeticiones tal como se muestra en la tabla de arriba.
Se determina que los factores A y B respectivamente corresponden a la Profundidad y la
velocidad de cada una de las pruebas, siendo el valor del acabado la respuesta obtenida.
3
Posteriormente, se realizan otras dos tablas dónde se registran los cuadrados de las sumas de
los datos de acabado obtenidos, así como los valores individuales del acabado para corrida y
repetición.
4
Se calculan las Sumas de cuadrados, obteniendo los siguientes resultados:
− SCa: 2,125.11
− SCb: 3,160.5
5
Se continua el cálculo de las Sumas de cuadrados, obteniendo los siguientes resultados:
− SCab: 557.06
− SCt: 6,532
− SCe: 689.33
6
Se calculan los cuadrados medios, obteniendo los siguientes resultados:
− CMa: 708.37
− CMb: 1,580.25
− CMab: 92.84
− CMe: 28.72
7
Posteriormente, se construye la tabla ANOVA y se evalúan las hipótesis comparando los valores
calculados de F con los que se leen en las tablas, resultando en que se rechaza H0 para todos
los casos evaluados.
8
Se realiza la comparación de medias por el método LSD obteniéndose que para el factor A
existen diferencias significativas entre cada una de las medias a excepción del caso – .
Mientras que para el caso del factor B, se obtiene que existen diferencias significativas entre
cada par de medias a excepción del caso –
9
De acuerdo con las gráficas de comparación de medias para ambos factores junto con las gráficas de
interacción, se determina que los niveles que arrojan un menor valor para la respuesta de acabado
son la Profundidad A1= 0.15 y la Velocidad B1= 0.2, ya que se observa que es en estas dónde se
obtienen los valores más bajos para el acabado. Así mismo, la combinación de factores que presenta
el valor de acabado más bajo es Profundidad A1 – Velocidad B1, con lo que se concluye que estos
son las mejores condiciones para obtener el mejor valor de acabado.
10
De acuerdo a lo observado en la gráfica de Probabilidad Normal, se determina que si se cumple el
supuesto de Normalidad con únicamente un punto aberrante presente.
11
De acuerdo con lo observado en la gráfica de residuos vs orden, se determina que si se cumple el
supuesto de independencia ya que no se observa un patrón sobre los puntos graficados.
12
Ejercicio 19.- Se corre un diseño factorial 3 x 2 con 10 réplicas para investigar el hinchamiento del
catalizador después de la extrusión en la fabricación de botellas de polietileno de alta densidad. EI
catalizador se utiliza en la obtención de dicho polietileno. Los factores investigados son: molde (con
dos niveles) y B: catalizador (con tres niveles). Los datos obtenidos se muestran en la siguiente tabla:
c) Dibuje las gráficas de medias para los dos efectos principales con los métodos LSD y de Tukey.
Compare los resultados de ambos métodos.
g) Utilice la gráfica de residuos contra factores para detectar posibles efectos sobre la dispersión del
hinchamiento. ¿En cuál molde parece que es menor la dispersión?
Se realiza la prueba sobre un modelo de diseño factorial 3x2 y se formulan las siguientes hipótesis:
13
𝐻0 : 𝐵1 = 𝐵2 = 𝐵3 = 0
𝐻1 : 𝐵𝑗 ≠ 0 𝑝𝑎𝑟𝑎 𝑎𝑙𝑔ú𝑛 𝑣𝑎𝑙𝑜𝑟 𝑑𝑒 𝑗
c) Gráfico de medias
94
93
Hinchamiento
92
91
90
89
88
1 2 3
Catalizador
14
Medias y 95.0% de Tukey HSD
94
93
Hinchamiento
92
91
90
89
88
1 2 3
Catalizador
93
92
Hinchamiento
91
90
89
88
1 2
Molde
93
92
Hinchamiento
91
90
89
88
1 2
Molde
15
d) Gráfica de interacción con intervalos de confianza sobrepuestos
96 Catalizador
1
2
94 3
Hinchamiento
92
90
88
86
1 2
Molde
16
f) Verificación de los supuestos de normalidad y varianza constante
99.9
99
95
porcentaje
80
50
20
0.1
-2.2 -0.2 1.8 3.8 5.8
RESIDUOS
17
g) Detección de posibles efectos sobre la dispersión del hinchamiento
2
residuos
-2
-4
-6
1 2
Molde
2
residuos
-2
-4
-6
1 2 3
Catalizador
Dados los datos observados en las gráficas de Residuos vs. Factores, se determina que la menor
dispersión se obtiene con el tipo de molde 1
18
Ejercicio 20.- Para mejorar la resistencia a la torsión de las adhesiones de componentes electrónicos
sobre placas, se estudiaron dos tipos de pegamentos (A1 y A2) y tres temperaturas de curado (60,
80 y 100 °C). En cada combinación se analizaron dos componentes y los
e) Verifique residuos.
19
Seguidamente se le pide al software que realice un Análisis de Varianza multifactorial con un nivel
de confianza del 95% y se obtiene la siguiente tabla ANOVA:
c) Gráficas de efectos
Se realiza la comparación de medias por método LSD con un nivel de confianza del 95% y se obtienen
las siguientes gráficas:
20
− Para la temperatura de curado (Factor B):
21
− Para la interacción de factores:
De acuerdo con lo obtenido en las gráficas se determina que el mejor efecto se obtiene al utilizar el
tratamiento con el tipo de pegamento A2 aplicando la temperatura T3 = 100°C, ya que con esta
combinación de factores se puede observar un incremento significativo en las medias de tensión
obtenidas.
e) Verificación de residuos
a)
22
b)
c)
d)
23
Ejercicio 21.- Se desea investigar de qué manera afecta el tiempo de curado y el tipo del acelerante
a la resistencia de caucho vulcanizado. Se realiza un experimento y se obtienen los siguientes datos:
c) Realice el análisis estadístico apropiado para probar las hipótesis que formuló.
d) ¿Hay algún tiempo de cura que es mejor para aumentar la resistencia? Argumente su respuesta.
g) Explique de manera gráfica cómo se obtuvo en la computadora el valor-p para tiempo de cura.
h) Verifique que se cumplan los supuestos. En caso de que no se cumpliera el supuesto de varianza
constante para el tiempo de cura, ¿qué significaría eso y cómo pudiera corregirse?
Para este problema se utilizó un Diseño Factorial del tipo 23 y se analizó por medio de un ANOVA
multifactorial.
b) Formulación de hipótesis
24
• Para la interacción de tipo de acelerante y tiempo de cura:
𝐻0 : (𝐴𝑇)𝑖𝑗 = 0 𝑝𝑎𝑟𝑎 𝑡𝑜𝑑𝑜 𝑣𝑎𝑙𝑜𝑟 𝑑𝑒 𝑖𝑗
𝐻1 : (𝐴𝑇)𝑖𝑗 ≠ 0 𝑝𝑎𝑟𝑎 𝑎𝑙𝑔ú𝑛 𝑣𝑎𝑙𝑜𝑟 𝑑𝑒 𝑖𝑗
c) Análisis estadístico
Para el análisis estadístico se realizó un análisis de varianza (ANOVA) multifactorial, para lo cual se
realizó un diseño factorial en el software estadístico Statgraphics y se ingresaron los datos
correspondientes a las respuestas de tensión obtenidas tal como se muestra a continuación:
Seguidamente se le pidió al software que realice un Análisis de Varianza multifactorial con un nivel
de confianza del 95% y se generó la siguiente tabla ANOVA:
Se realizó la comparación de medias por método LSD con un nivel de confianza del 95% y se
obtuvieron las siguientes gráficas:
25
− Para el tipo de acelerante (Factor A):
26
− Para la interacción de factores:
De acuerdo con lo obtenido en las gráficas se determina que el mejor tiempo de cura para obtener
una mayor resistencia del caucho es el tiempo a 60 minutos, ya que durante las evaluaciones con
este tiempo se obtuvieron las medias de resistencia más elevadas en las pruebas con los tres tipos
de acelerantes, sin embargo, no existen diferencias estadísticamente significativas entre alguna de
las tres medias obtenidas.
De la misma manera, se determina que el mejor acelerante para obtener una mayor resistencia del
caucho es el acelerante A2, ya que con este se obtuvieron las medias de resistencia más elevadas
27
durante las evaluaciones, sin embargo, tampoco existen diferencias estadísticamente significativas
entre las medias obtenidas.
Por último, se determina que el mejor efecto en cuanto a la resistencia de caucho se obtiene al
utilizar el tratamiento con el tipo de acelerante A2 aplicando el tiempo T2 = 60 min, ya que con esta
combinación de factores se puede observar un incremento significativo en las medias de resistencia
del caucho obtenidas.
Se realiza el cálculo del valor p de Levene por medio del software Excel utilizando la fórmula
=DISTR.F(x, grados_de_libertad_1, grados_de_libertad_2)
- Para el tipo de acelerante se toma el valor calculado de F=0.12, y los grados de libertad del
Factor A = 2 y de los residuos = 9, obteniendo un valor p = 0.89
- Para el tiempo de curado se toma el valor calculado de F=0.63, y los grados de libertad del
Factor B = 2 y de los residuos = 9, obteniendo un valor p = 0.55
28
- Para la interacción de factores se toma el valor calculado de F=0.23, y los grados de libertad
de la interacción AxB = 4 y de los residuos = 9, obteniendo un valor p = 0.91
29
Ejercicio 23.- En un laboratorio de microbiología se realiza un experimento para investigar si influye
el tipo de verdura (lechuga-L, cilantro-C, zanahoria-Z) y la temperatura (8 y 200C) de
almacenamiento en la sobrevivencia del vidrio colerae. Se hicieron varias réplicas. El porcentaje de
sobrevivencia obtenido después de 24 horas de inoculado el alimento se muestra a continuación:
a) Señale el nombre del diseño empleado y formule las hipótesis que pueden ser probadas.
c) Verifique el supuesto de igual varianza entre los tratamientos (varianza constante). ¿se cumple
satisfactoriamente?
d) En caso de que no se cumpla el supuesto anterior, ¿cómo afecta esto a sus conclusiones?
Se realiza la prueba sobre un modelo de diseño factorial 3x2 y se formulan las siguientes hipótesis:
• Para la temperatura:
𝐻0 : 𝑇1 = 𝑇2 = 0
𝐻1 : 𝑇1 ≠ 𝑇2
30
b) Construcción de tabla ANOVA
Se realiza una prueba de rangos múltiples para conocer cuales son los valores que presentan
diferencias significativas entre sí.
Para los tipos de vegetales, se observa que el alimento 3 es diferente del 1 y 2 siendo este el que
presenta una menor supervivencia.
Para las temperaturas, se observa que son diferentes entre sí y que a los 20 grados se presenta una
menor supervivencia de los microorganismos.
31
c) Verificación del supuesto de igual varianza entre tratamientos (varianza constante).
d) En caso de que no se cumpla el supuesto anterior, ¿cómo afecta esto a sus conclusiones?
La varianza constate al no cumplirse nos indica que los valores no son precisos o confiables.
32
Ejercicio 25.- Con el objetivo de estudiar la producción de huitlacoche (hongo comestible del maíz)
se decide correr un experimento con tres variedades de maíz en dos localidades maiceras. Las
variables de interés fueron: el porcentaje de cobertura de la mazorca por el hongo, el peso total de
la mazorca y el peso del huitlacoche. Se hicieron cuatro réplicas.
33
a) Escriba el nombre y modelo estadístico del diseño que se está empleando.
b) ¿Hay un efecto significativo de los factores variedad y localidad en las tres variables de respuesta?
c) ¿Existe claramente una localidad y variedad de maíz donde se produce más huitlacoche? Apóyese
en gráficas y pruebas estadísticas.
d) ¿El que haya mayor cobertura del hongo garantiza mayor producción de huitlacoche? Vea de
manera simultánea las gráficas de interacción para ambas variables de respuesta y/o haga un análisis
de correlación entre estas dos variables.
f) Haga los análisis de residuos para verificar los supuestos del modelo.
b) Efecto significativo de los factores variedad y localidad en las tres variables de respuesta
34
c) Determinación de localidad y variedad de maíz donde se produce más huitlacoche
35
d) Análisis de la interacción entre la cobertura del hongo y la producción de huitlacoche.
36
De acuerdo con la prueba de rangos múltiples, es posible obtener una producción promedio del
34.63% de huitlacoche en la localidad A
De acuerdo con la gráfica de residuos vs. predichos se confirma que no se cumple el supuesto de
varianza constante.
37
Ejercicio 26.- Los siguientes datos corresponden a diseño 3 x 3 con tres réplicas. Interesa investigar
el efecto de ambos factores sobre Y, para encontrar las condiciones adecuadas para maximizar.
d) Verifique supuestos.
Se utiliza un diseño factorial 3x3 para el cual se manejan las siguientes hipótesis
• Para el factor A
𝐻0 : 𝐴1 = 𝐴2 = 𝐴3 = 0
𝐻1 : 𝐴𝑖 ≠ 0 𝑝𝑎𝑟𝑎 𝑎𝑙𝑔ú𝑛 𝑣𝑎𝑙𝑜𝑟 𝑑𝑒 𝑖
• Para el Factor B
𝐻0 : 𝐵1 = 𝐵2 = 𝐵3 = 0
𝐻1 : 𝐵𝑗 ≠ 0 𝑝𝑎𝑟𝑎 𝑎𝑙𝑔ú𝑛 𝑣𝑎𝑙𝑜𝑟 𝑑𝑒 𝑗
38
Dados los resultados obtenidos en la tabla ANOVA, se concluye que únicamente existe un efecto
significativo dentro del Factor A, mientras que el efecto del factor B y el de interacción no influyen
en la respuesta del experimento.
39
d) Verificación de supuestos
De acuerdo con lo observado en la gráfica de Probabilidad Normal para el efecto, se determina que
si se cumple el supuesto de normalidad con únicamente un punto aberrante presente.
De acuerdo con lo observado en la gráfica de Predichos vs. Residuos, se determina que no se cumple
el supuesto de varianza constante.
40
De acuerdo con lo observado en la gráfica de Orden vs. Residuos, se determina que si se cumple el
supuesto de independencia al no observarse claramente un patrón no aleatorio entre los puntos.
41
De acuerdo con las pruebas de rangos múltiples, se puede determinar que para el caso del factor A
se tiene que todas las medias son significativamente diferentes entre sí, por lo que el nivel que
muestra un mejor efecto es el A2 debido a que posee una media más alta.
Para el caso del factor B se tiene que las medias de B1=B2 y B2=B3 pero B1≠B3, por lo que se podrían
utilizar los niveles B1 o B2 para obtener un mejor efecto ya que son los que poseen medias más altas
y no presentan diferencias significativas entre sí.
42
Ejercicio 27.- En una empresa alimenticia se desean evaluar cuatro antioxidantes, a través de su
efecto en un aceite vegetal. El propósito es seleccionar el producto que retrase más la oxidación. Las
pruebas se hacen en condiciones de estrés, midiendo como variable de respuesta al índice de
peróxidos. Se evalúan diferentes unidades experimentales a diferentes tiempos. Los datos obtenidos
se muestran a continuación (en el control no se agrega ningún antioxidante). Dado que uno de los
factores es el tiempo, y éste no se puede aleatorizar, entonces se le puede ver como un factor de
bloques.
b) Formule el modelo estadístico más apropiado al problema y las hipótesis estadísticas que se
pueden probar.
e) Considerando que a menor índice de peróxidos mejor es el producto, ¿hay algún producto que
sea mejor estadísticamente?
• Factores controlables:
− Producto (5 niveles)
− Tiempo (3 niveles)
• Variable respuesta:
− Índice de peróxidos
43
𝑌𝑖 𝑗 𝑘 = µ + 𝛼𝑖 + 𝛽𝑗 + (𝛼𝛽)𝑖𝑗 + 𝜀𝑖 𝑗 𝑘
𝑆𝑖𝑒𝑛𝑑𝑜: 𝑖 = 1,2,3,4,5; 𝑗 = 1,2,3; 𝑘 = 1,2
• Para el Producto
𝐻0 : 𝑃𝑐 = 𝑃1 = 𝑃2 = 𝑃3 = 𝑃4 = 0
𝐻1 : 𝑃𝑖 ≠ 0 𝑝𝑎𝑟𝑎 𝑎𝑙𝑔ú𝑛 𝑣𝑎𝑙𝑜𝑟 𝑑𝑒 𝑖
• Para el Tiempo
𝐻0 : 𝑇1 = 𝑇2 = 𝑇3 = 0
𝐻1 : 𝑇𝑗 ≠ 0 𝑝𝑎𝑟𝑎 𝑎𝑙𝑔ú𝑛 𝑣𝑎𝑙𝑜𝑟 𝑑𝑒 𝑗
De acuerdo con la tabla ANOVA, se puede concluir con un nivel de confianza del 95%, que tanto el
tiempo como la interacción producto X tiempo presentan un efecto significativo sobre el índice de
peróxidos en los antioxidantes.
Por otro lado, el tipo de producto no muestra un efecto significativo sobre el índice de peróxidos
en los antioxidantes.
Dado a que la interacción AB muestra efecto significativo, se puede concluir que el tiempo afecta la
presencia de peróxidos de acuerdo al tipo de producto utilizado, es decir que entre más tiempo
transcurra mayor será el retraso en la oxidación de los productos.
d) Verificación de supuestos.
44
De acuerdo con la gráfica de Probabilidad Normal, se determina que si se cumple el supuesto de
Normalidad con solamente un punto aberrante presente.
45
De acuerdo con la gráfica de Orden vs. Residuos, se determina que no se cumple el supuesto de
independencia, ya que se observa un patrón de comportamiento sobre los puntos graficados.
46