Ejercicio de Diseño Experimental

Descargar como docx, pdf o txt
Descargar como docx, pdf o txt
Está en la página 1de 8

Enunciado

Considérese una investigación llevada a cabo para estudiar el efecto que tiene la
concentración de un reactivo y la presencia de un catalizador sobre la conversión del
rendimiento un proceso químico.  Sea la concentración del reactivo el factor A con dos
niveles de interés, 15% y 20%.  El catalizador constituye el factor B; el nivel alto o
superior denota el uso de dos sacos de catalizador y el nivel bajo o inferior denota el uso
de un solo saco.  El experimento se realiza tres veces con los siguientes datos:

Resultados

Tabla 1. Datos y promedios de la variable de respuesta


Tota
factor A factor B Replicas l Yates
Concentració Catalizado Combinación de
n del reactivo r tratamientos x1 x2 R1 R2 R3
15% 1 A: Bajo ¿ B: Bajo ¿ -1 -1 28 25 27 80 1
20% 1 A: Alto¿ B: Bajo ¿ 1 -1 36 32 32 100 a
15% 2 A: Bajo ¿ B: Alto¿ -1 1 18 19 23 60 b
20% 2 A: Alto¿ B: Alto¿ 1 1 31 30 29 90 ab
330

Tabla 2. Factores y niveles

Factores Código Niveles


Concentración del A a 1: 15% 1: Alto
a 2: 25% -1: Bajo
reactivo
Catalizador B b1: 2 1: Alto
b 2: 1 -1: Bajo
Tabla 3. Tratamiento

Tratamient
Combinaciones
o
1 a1b1 Concentración del reactivo 15% - Catalizador 1
2 a2b1 Concentración del reactivo 20% - Catalizador 1
3 a1b2 Concentración del reactivo 15% - Catalizador 2
4 a2b2 Concentración del reactivo 20% - Catalizador 2

Tabla 4.Características del Experimento


Factor de experimental 2
Variable de respuesta 1
Replicas 3
Unidades experimentales 12
Aleatorizados No
Factor en estudio Tiempo de reacción de un proceso
químico
Factor A Concentración del reactivo
Factor B Catalizador

Se elaboró un diseño factorial 22aleatorio con 12 datos

Hipótesis

H 0 :μ 1=μ2=⋯ … . μ n

H 0 :Todos los tratamientos son iguales.

H 0 : Efecto A=0

H 0 : Efecto B=0

H 0 : Efecto AB=0

H 1: μ1≠ μ2 ⋯ … . μn

H 1 : Al menos un tratamiento es diferente

H 1 : Efecto A ≠ 0

H 1 : Efecto B ≠0

H 1 : Efecto AB≠ 0

Modelo estadístico del diseño factorial 22

Y ijk =μ+∝i+ β j+ ¿

Dónde:

μ → Efecto media global

∝i →Efecto incremental sobre la media causada por el nivel i del factor A.

β j →Efecto incremental sobre la media causada por el nivel j del factor B.


¿ Efecto incremental sobre la media causada por la interacción del niveli del factor A y
del nivel jdel factor B.

ε ijk → Termino del error.

Tabla 5. Signos de un diseño factorial 22.

A B AB
- - +
+ - -
- + +
+ + -

Tabla 6. Datos del corridos en el programa statgraphics.

factor A factor B
Concentración del Catalizador
Bloque Tiempo de reacción
reactivo
1 15 1 28
1 20 2 18
1 15 1 36
1 20 2 31
2 15 1 25
2 20 2 19
2 15 1 32
2 20 2 30
3 15 1 27
3 20 2 23
3 15 1 32
3 20 2 29

Tabla 7. Análisis de varianza ANOVA

Fuente de Suma de cuadrados G.L Cuadrado F0 Valor-p


variabilidad medio
A [ a+ab−b−(1)]
2
A× A SC A CM A P (F¿ F 0 ¿
SC A = 2 Gl A CM E
n2

B [ b+ab−a−(1)]
2
B× B SC B CM B P (F¿ F 0 ¿
SC B= Gl B CM E
n 22
AB [ ab+ ( 1 )−a−b ]
2
A×B SC AB CM AB P (F¿ F 0 ¿
SC AB = Gl AB CM E
n 22

Error SC E=SC A −SC B−SC AB −SC T 4 (n−1) SC E


Gl E

Total SCT =SC A + SC B + SC AB n 22−1

Factor de Correlación
2
∑ (Y i . )
FC =
r∗t

(330 )2 108900
F C= = =9075
4 ×3 12
Contrastes
A=[ a+ab−b−(1)] =[ 100+90−60−80 ] =50

B=[ b+ ab−a−(1) ] =[ 60+90−100−80 ] =30

AB=[ ab+ (1 )−a−b ] = [ 90+80−100−60 ]=10

Efectos
1 1 1
A= [ a+ab−b−(1)] = [ 100+ 90−60−80 ] = ( 50 )=8,33
2n 2(3) 6
1 1 1
B= [ b+ ab−a−(1) ]= [ 60+ 90−100−80 ] = (30 )=5
2n 2(3) 6
1 1 1
AB= [ ab+ ( 1 ) −a−b ] = [ 90+80−100−60 ] = ( 10 ) =1,67
2n 2(3) 6

Suma de cuadrados
Suma de cuadrados de A
2
[ a+ab−b−( 1 ) ]
SC A = 2
n ( 2)

[ 100+ 90−60−80 ]2 (50)2


SC A = 2
= =208,33
n ( 2) 3(4 )

Suma de cuadrados de B
2
[ b+ab−a−( 1 ) ]
SC B= 2
n ( 2)

[ 60+ 90−100−80 ]2 (30)2


SC A = 2
= =75
n ( 2) 3(4 )

Suma de cuadrados de AB
2
[ ab+ ( 1 )−a−b ]
SC B= 2
n ( 2)

[ 90+80−100−60 ]2 (10)2
SC A = 2
= =8,33
n ( 2) 3(4 )

Suma de cuadrados Totales


SCT =SC A + SC B + SC AB

SCT =( 208,33 ) + ( 75 ) + ( 8,33 ) =291,66

SCT =[ ( 28 )2− (36 )2+ (18 )2 + ( 31 )2 + ( 25 )2 + (32 )2 + ( 19 )2 + ( 30 )2 + ( 27 )2 + ( 32 )2+ ( 23 )2 + ( 29 )2 ] −9075

SCT =[ 784+1296+324 +961+625+1024 +361+900+729+1024+ 529+841 ] −9075


SCT =9398−9075=323

Suma del Error


SC E=SC T −SC A −SC B−SC AB

SC E=323−208,33−75−8,33=31,34

Grados de libertad
Gl A = A × A=1 ×1=1

Gl B=B × B=1 ×1=1

Gl AB =A × B=1 ×1=1

Gl Error =4 ( n−1 )=4 ( 3−1 ) =4 ( 2 ) =8

GlTotal =n2 2−1=3 ¿

Cuadrados Medios

SC A 208,33
CM A = = =208,33
Gl A 1
SC B 75
CM B= = =75
Gl B 1

SC AB 8,33
CM AB = = =8,33
Gl AB 1

SC Error 31,34
CM Error = = =3,92
Gl Error 8

Razón de F

CM A 208,33
F A= = =53,15
CM Error 3,92

CM B 75
F B= = =19,13
CM Error 3,92

CM AB 8,33
F AB= = =2,13
CM Error 3,92

Tabla 8. Resultado de análisis de varianza ANOVA

Fuente Suma de Gl Cuadrado Razón-F Valor-P


Cuadrados Medio
A:
Concentració 208,33 1 208,33 53,15 0,0001
n del reactivo
B: 75 1 75 19,13 0,0024
Catalizador
AB 8,33 1 8,33 2,13 0,1826
Error total 31,34 8 3,92
Total 323 11
(corregidos)
La tabla ANOVA nos indica que los valores de p en dos factores que son el factor A
(concentraccion del reactivo) y el factor B (Catalizador) son menores que 0,05, el cual
es estadísticamente significativo sobre el Tiempo de reacción de un proceso químico del
producto con un 95,0% de nivel de confianza.

Diagrama de Pareto
En el diagrama de Pareto nos indica que el efecto del factor A que es la concentración
de reactivo se acerca más a la variable de respuesta que es el Tiempo de un proceso
qímico.
Grafica de interacciones

Tiempo de reaccion

Concentracion del reactivo

También podría gustarte