Ingeniería Agrícola y Uso Integral del Agua

Proyecto:

Diseño Experimental Completamente Aleatorio:

Tiempo de combustión de cerillos

Semestre: 2

Profesor:

Dr. Mario Alberto Vázquez Peña

Presenta:

Angel Garduño García

CHAPINGO, MÉXICO. MAYO DE 2011

Introducción:

El diseño experimental completamente aleatorio consiste en la asignación de los

tratamientos en forma completamente aleatoria a las unidades experimentales (individuos,
grupos, parcelas, jaulas, animales, insectos, etc.). Debido a su aleatorización irrestricta, es
conveniente que se utilicen unidades experimentales de lo más homogéneas posibles: animales de
la misma edad, del mismo peso, similar estado fisiológico; parcelas de igual tamaño, etc., de
manera de disminuir la magnitud del error experimental, ocasionado por la variación intrínseca de
las unidades experimentales. Este diseño es apropiado para experimentos de laboratorio,
invernadero, etc., es decir, situaciones experimentales en las que las variables rodean el
experimento puedan ser controlables o lo más homogéneas posibles.

Objetivos:

 Implementar un diseño experimental completamente aleatorio, para analizar las

diferencias existentes en el tiempo de combustión de fósforos con distintos tratamientos.
 Realizar un análisis de varianza (A de V) a los datos obtenidos con los distintos
tratamientos.
 Realizar la prueba de Bartlett con la finalidad de probar la igualdad de varianzas en los
tratamientos y dar sustento al A de V.
 Realizar alguna prueba para determinar el supuesto de normalidad en las poblaciones.
 Realizar alguna prueba de separación de medias, para identificar el mejor tratamiento.

Materiales:

 Plastilina
 Cronómetro
 Cerillos de 3 marcas diferentes (En este caso, “Monarca, Manola, Clásicos”)
 Superficie plana

Metodología:

Para la realización del experimento se utilizaron 3 marcas diferentes de cerillos, para los cuales
se deseaba comparar el tiempo de combustión, sin embargo, para cada marca se realizó el
experimento para dos posiciones distintas del cerillo, es decir; en posición horizontal y en posición
vertical, esto debido a que el tiempo de combustión para las distintas posiciones era diferente.
Con lo anterior, quedaron especificados 6 tratamientos, de la siguiente forma:

 Monarca- Vertical
 Monarca- Horizontal
 Manola-Vertical
 Manola-Horizontal
 Clásicos-Vertical
 Clásicos-Horizontal
 Para cada uno de los tratamientos anteriores se realizaron 10 repeticiones, obteniendo de
esta forma una tabla con 60 datos agrupados para su análisis.

En la figura 1, puede observarse la disposición de los cerillos tanto en forma vertical como
en forma horizontal:

Figura 1. “Disposición de los cerillo en forma horizontal y vertical”

Como puede observarse, los cerillos fueron sujetados con plastilina, con la finalidad de
evitar que existiera un error grande en la medición, ya sea porque el fósforo se apagara debido al
movimiento de mano o tuviera un tiempo diferente de combustión debido a la inclinación del
mismo.

De la misma forma en la figura 1, se muestran dos tratamientos: el de Clásicos-Vertical y

Clásicos horizontal, es decir que este paso se repitió para los cerillos de marca Manola y Monarca.

Figura la 2. “medida del cerillo”

En la figura 2 se muestra la medida que tiene el cerillo, lo cual sirve para referenciar el
hecho, de que al momento de clavarse en la plastilina, se cuidó que 3 cm quedarán expuestos al
ambiente, mientras que 0.2 cm, eran la parte que quedaba clavada.

En la figura 3, se muestra la combustión de un cerillo en posición vertical y horizontal, así

como la marca que dejan en la plastilina al apagarse. Con esto es posible explicar, que la medición
del tiempo se realizó justo en el momento de encendido del cerillo, hasta que este mismo se
apagó sólo (en el momento de que la flama tocaba la plastilina). Cabe resaltar que el encendido de
lo cerillos se realizó con una flama externa para evitar el movimiento del mismo.
 Figura 3. “Combustión de los cerillos en posición vertical y horizontal”

De la misma forma estos pasos fueron repetidos para las marcas Manola y Monarca, para
finalmente registrar los datos en una tabla de resultados que contiene los distintos tratamientos
con sus respectivas repeticiones y tiempos de combustión de los cerillos.

Cabe resaltar que para llevar a cabo el experimento se tuvo que tomar medidas, que
hicieran posible que el error experimental fuera el mínimo, entre estas se pueden mencionar las
siguientes:

 El experimento fue realizado en un cuarto cerrado para evitar corrientes de aire que
pudieran apagar el fósforo.

 Debido a que el cerillo necesitaba oxígeno para la combustión, cada que se realizaba una
repetición de experimento, se esperaba 1 minuto para realizar una nueva, esto con la
finalidad de que la cantidad de oxígeno en la atmosfera se normalizara.

 Los cerillos eran fijados en plastilina para evitar el error humano al estar agarrándolo
durante la combustión y que el pulso la afectara.

 La misma persona que encendía el cerillo era quién hacia iniciar el conteo en el
cronómetro, esto con el fin de evitar un desfasamiento entre el tiempo de encendido y el
inicio del conteo.
Resultados:

Se obtuvo el siguiente diseño experimental completamente aleatorio donde se ensayaron

diferentes marcas y posición para la combustión de cerillos expresada en tiempo (segundos).

Tratamientos
Monarca Monarca Manola Manola Clásicos Clásicos
Vertical Horizontal Vertical Horizontal Vertical Horizontal
Número de
muestra Tiempo de combustión del cerillo (segundos)
44.4 20.2 40.2 24.7 41.2 23.8
1
43.1 21.9 42.1 23.2 42.3 21.5
2
44.3 22.1 43.2 25.0 42.4 20.1
3
44.5 19.8 41.3 24.3 44.7 22.4
4
44.0 19.6 40.3 21.5 43.0 23.8
5
43.7 18.5 40.3 22.9 41.3 21.9
6
43.1 17.6 39.5 22.3 43.2 23.4
7
43.6 22.3 40.1 21.3 42.2 21.4
8
45.2 23.1 38.3 22.1 40.8 21.9
9
45.1 19.8 40.2 21.6 39.7 22.3
10
0.5467 3.2010 1.8583 1.8788 1.9796 1.3672

9 9 9 9 9 9

Análisis de Resultados:

ANALISIS DE VARIANZA:

En Microsoft Excel se realizó un análisis de varianza de un factor, de lo cual surgieron las siguientes
tablas:

RESUMEN
Grupos Cuenta Suma Promedio Varianza

Monarca Vertical 10 441 44.1 0.546666667

Monarca Horizontal 10 204.9 20.49 3.201
Manola Vertical 10 405.5 40.55 1.858333333
 Manola Horizontal 10 228.9 22.89 1.878777778
Clásicos Vertical 10 420.8 42.08 1.979555556
Clásicos Horizontal 10 222.5 22.25 1.367222222

ANÁLISIS DE VARIANZA
Origen de las Suma de Grados de Promedio de los Valor crítico
variaciones cuadrados libertad cuadrados Fo Probabilidad para Ft
Tratamientos 6316.32 5 1263.264 699.7687 8.959E-48 3.3769
Error 97.484 54 1.805259259
Total 6413.804 59

Se desea probar la hipótesis:

Ho: vs. Ha: Al menos un tratamiento es diferente a los demás.

Como Fo>Ft para α=0.01, se procede a rechazar Ho y se concluye que existen diferencias
significativas en el tiempo de combustión de los cerillos al utilizar diferentes tratamientos respecto
a la marca y posición de combustión del cerillo.

PRUEBA DE BARTLETT PARA PROBAR IGUALDAD DE VARIANZAS

Se desea probar la hipótesis:

Ho:
en oposición a
Ha: al menos una es diferente de las demás.

De la tabla de datos se obtuvo lo siguiente:

1.9796 1.3672

Con los datos anteriores se forma una tabla de cálculos para obtener

0.5467 9 4.9203 -0.6039 -5.4351 0.1111

3.2010 9 28.809 1.1635 10.4715 0.1111
1.8583 9 16.7247 0.6197 5.5773 0.1111
1.8788 9 16.9092 0.6306 5.9454 0.1111
1.9796 9 17.8164 0.6829 6.1461 0.1111
1.3672 9 12.3048 0.3128 2.8152 0.1111
Sumatoria l=54 97.4844 25.5207 0.6667
De la tabla anterior es fácil obtener:

∑
∑

(∑ ) ( )

Estadística utilizada para la prueba:

( ∑ )

Criterio de rechazo:

Rechazar Ho si , con

Con

Por lo que no se rechaza Ho, lo cual quiere decir que se

cumple la homogeneidad en las varianzas .

VERIFICACIÓN DE SUPUESTO DE NORMALIDAD:

Debido a que se está trabajando con n=60, y las tablas de los coeficientes para la prueba de
Shapiro-Wilk tiene datos hasta n=50, entonces se verificará el supuesto de normalidad con la
prueba Gráfica. Para esto que se debe obtener la gráfica a partir de las siguientes parejas de
valores ( )( ) ( ) donde ( ) son los residuales ordenados
de menor a mayor y Zp son los percentiles de la distribución Normal estándar de nivel ,
i=1, 2, …, n respectivamente:

Tratamiento Dato Media Residuo ri residuo ordenado i pi zpi

Monarca Vertical 44.4 0.300 -2.89 1 0.0083 -2.394
Monarca Vertical 43.1 -1.000 -2.38 2 0.0250 -1.960
Monarca Vertical 44.3 0.200 -2.25 3 0.0417 -1.732
Monarca Vertical 44.5 0.400 -2.15 4 0.0583 -1.569
Monarca Vertical 44.0 44.1 -0.100 -1.99 5 0.0750 -1.440
Monarca Vertical 43.7 -0.400 -1.59 6 0.0917 -1.331
Monarca Vertical 43.1 -1.000 -1.39 7 0.1083 -1.235
Monarca Vertical 43.6 -0.500 -1.29 8 0.1250 -1.150
Monarca Vertical 45.2 1.100 -1.28 9 0.1417 -1.073
 Monarca Vertical 45.1 1.000 -1.05 10 0.1583 -1.001
Monarca Horizontal 20.2 -0.290 -1 11 0.1750 -0.935
Monarca Horizontal 21.9 1.410 -1 12 0.1917 -0.872
Monarca Horizontal 22.1 1.610 -0.89 13 0.2083 -0.812
Monarca Horizontal 19.8 -0.690 -0.88 14 0.2250 -0.755
Monarca Horizontal 19.6 20.49 -0.890 -0.85 15 0.2417 -0.701
Monarca Horizontal 18.5 -1.990 -0.79 16 0.2583 -0.648
Monarca Horizontal 17.6 -2.890 -0.78 17 0.2750 -0.598
Monarca Horizontal 22.3 1.810 -0.75 18 0.2917 -0.549
Monarca Horizontal 23.1 2.610 -0.69 19 0.3083 -0.501
Monarca Horizontal 19.8 -0.690 -0.69 20 0.3250 -0.454
Manola Vertical 40.2 -0.350 -0.59 21 0.3417 -0.408
Manola Vertical 42.1 1.550 -0.5 22 0.3583 -0.363
Manola Vertical 43.2 2.650 -0.45 23 0.3750 -0.319
Manola Vertical 41.3 0.750 -0.4 24 0.3917 -0.275
Manola Vertical 40.3 40.55 -0.250 -0.35 25 0.4083 -0.232
Manola Vertical 40.3 -0.250 -0.35 26 0.4250 -0.189
Manola Vertical 39.5 -1.050 -0.35 27 0.4417 -0.147
Manola Vertical 40.1 -0.450 -0.35 28 0.4583 -0.105
Manola Vertical 38.3 -2.250 -0.29 29 0.4750 -0.063
Manola Vertical 40.2 -0.350 -0.25 30 0.4917 -0.021
Manola Horizontal 24.7 1.810 -0.25 31 0.5083 0.021
Manola Horizontal 23.2 0.310 -0.1 32 0.5250 0.063
Manola Horizontal 25.0 2.110 0.01 33 0.5417 0.105
Manola Horizontal 24.3 1.410 0.05 34 0.5583 0.147
Manola Horizontal 21.5 -1.390 0.12 35 0.5750 0.189
Manola Horizontal 22.9 22.89 0.010 0.15 36 0.5917 0.232
Manola Horizontal 22.3 -0.590 0.2 37 0.6083 0.275
Manola Horizontal 21.3 -1.590 0.22 38 0.6250 0.319
Manola Horizontal 22.1 -0.790 0.3 39 0.6417 0.363
Manola Horizontal 21.6 -1.290 0.31 40 0.6583 0.408
Clásicos Vertical 41.2 -0.880 0.32 41 0.6750 0.454
Clásicos Vertical 42.3 0.220 0.4 42 0.6917 0.501
Clásicos Vertical 42.4 0.320 0.75 43 0.7083 0.549
Clásicos Vertical 44.7 2.620 0.92 44 0.7250 0.598
Clásicos Vertical 43.0 42.08 0.920 1 45 0.7417 0.648
Clásicos Vertical 41.3 -0.780 1.1 46 0.7583 0.701
Clásicos Vertical 43.2 1.120 1.12 47 0.7750 0.755
Clásicos Vertical 42.2 0.120 1.15 48 0.7917 0.812
Clásicos Vertical 40.8 -1.280 1.41 49 0.8083 0.872
Clásicos Vertical 39.7 -2.380 1.41 50 0.8250 0.935
 Clásicos Horizontal 23.8 1.550 1.55 51 0.8417 1.001
Clásicos Horizontal 21.5 -0.750 1.55 52 0.8583 1.073
Clásicos Horizontal 20.1 -2.150 1.55 53 0.8750 1.150
Clásicos Horizontal 22.4 0.150 1.61 54 0.8917 1.235
Clásicos Horizontal 23.8 1.550 1.81 55 0.9083 1.331
Clásicos Horizontal 21.9 22.25 -0.350 1.81 56 0.9250 1.440
Clásicos Horizontal 23.4 1.150 2.11 57 0.9417 1.569
Clásicos Horizontal 21.4 -0.850 2.61 58 0.9583 1.732
Clásicos Horizontal 21.9 -0.350 2.62 59 0.9750 1.960
Clásicos Horizontal 22.3 0.050 2.65 60 0.9917 2.394

Gráfica de Residuales vs. Zpi

3

1
Residuales

0
0.021
0.105
0.189
0.275
0.363
0.454
0.549
0.648
0.755
0.872
1.001
1.150
1.331
1.569
1.960
-0.701
-2.394
-1.732
-1.440
-1.235
-1.073
-0.935
-0.812

-0.598
-0.501
-0.408
-0.319
-0.232
-0.147
-0.063

-1

-2
Comportamiento Lineal
-3

-4
Zpi

La gráfica mostrada fue realizada en Excel y se puede apreciar el comportamiento lineal de los
datos en la parte central que es la más importante, lo cual muestra una aproximación para
verificar el supuesto de Normalidad.
SEPARACIÓN DE MEDIAS:

Una vez que se rechaza la hipótesis nula para los resultados obtenidos en el diseño experimental
anteriormente y se ha comprobado la igualdad de varianzas entre tratamientos, así como el
supuesto de normalidad, es posible proceder a la separación de medias, para lo cual se realizará
mediante la prueba de Tukey:

Comparaciones posibles de medias:

Posteriormente se procede a la obtención del valor crítico de la prueba aplicando la siguiente

fórmula:
̅

Donde los valores de se obtiene de la tabla VIII del libro Montgomery:

Fijando el valor , con t=6 y glError=54 se tiene:

En la tabla de valores críticos para la prueba de Tukey, se obtuvieron los siguientes valores por lo
que se realiza una interpolación:

Grados de Libertad t=6

glError=40 4.232
glError=54 x=4.184
glError=60 4.163

Por otra parte se tiene que ̅ √ √ , por lo tanto:

De tal manera que para aquella diferencia entre medias que superen este valor se consideran
estadísticamente diferentes a nivel de significancia de
̅ ̅ -20.49 ̅ -22.25 ̅ -22.89 ̅ =-40.55 ̅ -42.08
Monarca - Vertical 44.10 23.61* 21.85* 21.21* 3.55* 2.02*
Clásicos - Vertical 42.08 21.59* 19.83* 19.19* 1.53
Manola - Vertical 40.55 20.06* 18.30* 17.66*
Manola - Horizontal 22.89 2.40* 0.64
Clásicos - Horizontal 22.25 1.76
Monarca -Horizontal 20.49

Puede verse que las diferencias que superan el valor de Tα son las indicadas en asterisco, lo que
permite concluir que el tiempo de combustión promedio en las marcas Manola, Clásicos y
Monarca en posición horizontal se diferencian estadísticamente del tiempo de combustión
promedio obtenido en las marcas Manola, Clásicos y Monarca en posición vertical a un nivel de
significancia de 0.05.

Conclusiones:

De este experimento se puede deducir fácilmente que el tiempo de combustión en un cerillo

cuando está en posición vertical, es mayor que cuando se encuentra en posición horizontal, y del
análisis de datos es posible concluir que entre los 6 tratamientos, el que un mayor tiempo de
combustión presenta es el Monarca-Vertical, sin embargo, no existe una gran diferencia, con el
Clásico-Vertical y el Manola-Vertical, por lo que a la hora de recomendar utilizar cualquiera de las
tres marcas sería conveniente verificar otros factores tales como el precio, la disponibilidad,
cantidad de cerillos por caja, etc.

Bibliografía:

 Said, Infante. Gil; Guillermo, P. Sárate Lara. “Métodos Estadísticos”. Segunda Edición.
Editorial Trillas. Sexta Reimpresión, 2000.

 Montgomery. “Diseño y análisis de Experimentos”. Segunda Edición. Editorial Limusa.

2004.

 Análisis de Diseños Experimentales Básicos:

Fuente: http://www.galeon.com/colposfesz/est501/dca/dca.htm

 Combustión y cmbustibles:
Fuente: http://www.energia.inf.cu/iee-mep/SyT/CDG/Taller1BURE/COMBUSTION.PDF