capítulo 5 resultados -...

69

Capítulo 5

RESULTADOS En este capitulo, se presentan las aplicaciones numéricas correspondientes con el fin de

obtener resultados concretos en la comparación de los dos tipos de modelos, para tal se

corren ocho ejemplos diferentes en las dos modalidades de bloques al azar. Así como la

aplicación de la metodología de la prueba DE programada en Fortran.

5.1 CORRIDA EN EL LENGUAJE FORTRAN (EJEMPLO 1).

En esta sección se pretende instruir a lector sobre los pasos necesarios que habrá de

considerar para correr los ocho ejemplos que permiten hacer la comparación de los dos

modelos bootstrap paramétricos en bloques al azar.

5.1.1 EJEMPLO DE UNA CORRIDA PARA EL MODELO N(0,1)

Se tomo el ejemplo 5-1 de la página 149 de Montgomery (1999).

Una química desea probar el efecto que tienen cuatro agentes químicos sobre la

resistencia de un tipo particular de tela. Como puede haber variabilidad entre un rollo de

tela y otro, decide utilizar un diseño aleatorizado por bloques, considerando los rollos de

tela como bloques. Ella selecciona 5 rollos y les aplica los cuatro agentes químicos en

orden aleatorio. A continuación en la tabla 5.1 se proporcionan los resultados de la

resistencia a la tensión.

CAPÍTULO 5. RESULTADOS

70

Rollo de tela

Agente químico 1 2 3 4 5 1 73 68 74 71 67 2 73 67 75 72 70 3 75 68 78 73 68 4 73 71 75 75 69

Tabla 5.1

A continuación se enumeran los pasos necesarios para la resolución del problema,

tomando en cuenta que para esta primera corrida se decide generar números aleatorios que

se distribuyan N(0,1).

1. Crear un archivo externo cuya extensión sea *.DAT, este archivo contendrá

los datos del diseño a analizar, la configuración de los datos será como sigue: primeramente

se introduce el número de tratamientos, enseguida el de bloques o repeticiones por

tratamiento, y finalmente las observaciones comenzando por le primer tratamiento, ver

para este ejemplo la tabla 5.2.

Tabla 5.2

4 No. de tratamientos 5 No. de bloques 73 68 74 71 67 73 67 75 72 70 Observaciones 75 68 78 73 68 73 71 75 75 69

CAPÍTULO 5. RESULTADOS

71

2. Abrir el archivo con la extensión *.FOR, que contiene el programa fuente.

Para correr el programa es necesario ejecutar los comandos: Compile file, Build y Execute,

del menú Project en la pantalla de dialogo de Fortran. Esto con el fin de verificar posibles

errores.

3. Con los pasos 1 y 2 se despliega una pantalla de ejecución, donde el

programa ha comenzado a desplegarse, la primera instrucción es: “DESEAS ANALIZAR

UN EXPERIMENTO? (SI=1)”. Si la repuesta fue un uno, el programa pregunta el

nombre del archivo de entrada que contiene los datos: “INTRODUCE EL NOMBRE DEL

ARCHIVO DE ENTRADA (INPUT)”. Después se pide el nombre del archivo de salida

donde se desplegaran los resultados con la extensión *.LIS: “INTRODUCE EL NOMBRE

DEL ARCHIVO DE SALIDA (OUTPUT)”. Después pregunta la semilla del simulador:

“INTRODUCE LA SEMILLA DEBE SER UN ENTERO MAYOR QUE 0 Y MENOR

QUE 2147483647”. Posteriormente se pide el número de simulaciones y el número de

iteraciones por simulación: “INTRODUCE EL NUMERO DE SIMULACIONES (IB) Y

EL DE ITERACIONES POR SIMULACION (IREPB)”. Finalmente se pregunta por el

nivel de confianza: “INTRODUCE EL NIVEL DE CONFIANZA DEBE SER UN

NUMERO ENTRE 0 Y 100”.

4. Ahora el programa ha obtenido los datos necesarios para comenzar la

simulación. Sin embargo, necesita saber que tipo de modelo es el que se correrá si es el

modelo N(0,1) introduce un uno , o si es le modelo N(0, CME), se introduce cualquier otro

número: “SI DESEAS CORRER EL MODELO N(0,1) INTRODUCE UNO”.

5. Un vez determinado el tipo de diseño experimental, el programa pregunta si

se desea un probar un contraste: “QUIERES PROBAR UN CONTRASTE?(SI =1)”, si la

respuesta es afirmativa se introduce un uno, para que después el programa pregunte:

“METE LOS COEFICIENTES DEL CONTRASTE”.

CAPÍTULO 5. RESULTADOS

72

6. Una ves introducidos los coeficientes 2 1 -1 -2, se da inicio a la simulación,

desplegándose en pantalla:

SIMULACIÓN 1

SIMULACIÓN 2

SIMULACIÓN 3

…

SIMULACIÓN 30

El numero de simulaciones para cada uno de los ejemplos que se expondrán este

capitulo, será de 30, para facilitar la comparación entre modelos con los mismos

parámetros.

7. Desplegándose igualmente y por último en pantalla, el intervalo de

confianza para el valor P calculado del contraste basado en el modelo N(0,1) y se

pregunta si se desea probar otro contraste, sea:

“INTERVALO DE CONFIANZA DEL PVALUE DEL CONTRASTE:

.0991< .0995< .1000

DESEAS ANALIZAR UN EXPERIMENTO? (SI=1)”

Para responder a esta pregunta, un número diferente de uno termina con la corrida

del programa.

8. Para conocer los resultados es necesario abrir el archivo de salida que se

genero con la corrida del programa este es:

PROGRAMA PARA INVESTIGAR LA COMPARACION DE DOS MODELOS DE BOOTSTRAP. PROGRAMADO POR: ESTEBAN BURGUETE HERNANDEZ,

CAPÍTULO 5. RESULTADOS

73

JOSE FCO. TAMBORERO ARNAL. MANUEL ALEJANDRO PEREZ MIRELES. ARCHIVO DE ENTRADA: mon51.dat ARCHIVO DE SALIDA: mon56.lis SEMILLA PARA SIMULACION: 123412414 NUMERO TOTAL DE SIMULACIONES: 30 ITERACIONES POR SIMULACION: 1000000 NIVEL DE CONFIANZA: 95.00 Z DE TABLAS: 1.96 DISEÑO: BLOQUES AL AZAR NUMERO DE TRATAMIENTOS: 4 NUMERO DE BLOQUES: 5 DATOS: TRAT( 1) BLQ( 1) 73.0000 TRAT( 1) BLQ( 2) 68.0000 TRAT( 1) BLQ( 3) 74.0000 TRAT( 1) BLQ( 4) 71.0000 TRAT( 1) BLQ( 5) 67.0000 TRAT( 2) BLQ( 1) 73.0000 TRAT( 2) BLQ( 2) 67.0000 TRAT( 2) BLQ( 3) 75.0000 TRAT( 2) BLQ( 4) 72.0000 TRAT( 2) BLQ( 5) 70.0000 TRAT( 3) BLQ( 1) 75.0000 TRAT( 3) BLQ( 2) 68.0000 TRAT( 3) BLQ( 3) 78.0000 TRAT( 3) BLQ( 4) 73.0000 TRAT( 3) BLQ( 5) 68.0000 TRAT( 4) BLQ( 1) 73.0000 TRAT( 4) BLQ( 2) 71.0000 TRAT( 4) BLQ( 3) 75.0000 TRAT( 4) BLQ( 4) 75.0000 TRAT( 4) BLQ( 5) 69.0000 MODELO= 1 SI MODELO=1, EL MODELO ES EL N(0,1) TOTALES DE LOS TRATAMIENTOS. TOTAL DEL TRATAMIENTO: 1 353.0000 TOTAL DEL TRATAMIENTO: 2 357.0000 TOTAL DEL TRATAMIENTO: 3 362.0000 TOTAL DEL TRATAMIENTO: 4 363.0000

CAPÍTULO 5. RESULTADOS

74

CUADRADO MEDIO DEL ERROR: 1.8167 RESULTADOS FINALES. COEFICIENTES DEL CONTRASTE: COEF. DEL TRAT( 1) 2.00 COEF. DEL TRAT( 2) 1.00 COEF. DEL TRAT( 3) -1.00 COEF. DEL TRAT( 4) -2.00 F CALCULADA DEL CONTRASTE: 6.8807 INTERVALO DE CONFIANZA DEL PVALUE DEL CONTRASTE: .0991< .0995< .1000 GRACIAS POR USAR ESTE PROGRAMA

5.1.2 CORRIDA PARA EL MODELO N(0,CME)

Se considera el mismo ejemplo de la sección 5.1.1, considere igualmente el mismo

contraste. Dado que el propósito del presente documento es comparar el comportamiento

de cada uno de los modelos, se toma exactamente los mismos parámetros; datos, número de

simulaciones, iteraciones por simulación y contraste. Al correr el programa para este

segundo ejemplo se tiene:

PROGRAMA PARA INVESTIGAR LA COMPARACION DE DOS MODELOS DE BOOTSTRAP. PROGRAMADO POR: ESTEBAN BURGUETE HERNANDEZ, JOSE FCO. TAMBORERO ARNAL. MANUEL ALEJANDRO PEREZ MIRELES. ARCHIVO DE ENTRADA: mon51.dat ARCHIVO DE SALIDA: mon58.lis SEMILLA PARA SIMULACION: 123124142 NUMERO TOTAL DE SIMULACIONES: 30 ITERACIONES POR SIMULACION: 1000000 NIVEL DE CONFIANZA: 95.00

CAPÍTULO 5. RESULTADOS

75

Z DE TABLAS: 1.96 DISEÑO: BLOQUES AL AZAR NUMERO DE TRATAMIENTOS: 4 NUMERO DE BLOQUES: 5 DATOS: TRAT( 1) BLQ( 1) 73.0000 TRAT( 1) BLQ( 2) 68.0000 TRAT( 1) BLQ( 3) 74.0000 TRAT( 1) BLQ( 4) 71.0000 TRAT( 1) BLQ( 5) 67.0000 TRAT( 2) BLQ( 1) 73.0000 TRAT( 2) BLQ( 2) 67.0000 TRAT( 2) BLQ( 3) 75.0000 TRAT( 2) BLQ( 4) 72.0000 TRAT( 2) BLQ( 5) 70.0000 TRAT( 3) BLQ( 1) 75.0000 TRAT( 3) BLQ( 2) 68.0000 TRAT( 3) BLQ( 3) 78.0000 TRAT( 3) BLQ( 4) 73.0000 TRAT( 3) BLQ( 5) 68.0000 TRAT( 4) BLQ( 1) 73.0000 TRAT( 4) BLQ( 2) 71.0000 TRAT( 4) BLQ( 3) 75.0000 TRAT( 4) BLQ( 4) 75.0000 TRAT( 4) BLQ( 5) 69.0000 MODELO= 2 SI MODELO=1, EL MODELO ES EL N(0,1) TOTALES DE LOS TRATAMIENTOS. TOTAL DEL TRATAMIENTO: 1 353.0000 TOTAL DEL TRATAMIENTO: 2 357.0000 TOTAL DEL TRATAMIENTO: 3 362.0000 TOTAL DEL TRATAMIENTO: 4 363.0000 CUADRADO MEDIO DEL ERROR: 1.8167 RESULTADOS FINALES. COEFICIENTES DEL CONTRASTE: COEF. DEL TRAT( 1) 2.00 COEF. DEL TRAT( 2) 1.00 COEF. DEL TRAT( 3) -1.00 COEF. DEL TRAT( 4) -2.00 F CALCULADA DEL CONTRASTE: 6.8807

CAPÍTULO 5. RESULTADOS

76

INTERVALO DE CONFIANZA DEL PVALUE DEL CONTRASTE: .0988< .0994< .1001 GRACIAS POR USAR ESTE PROGRAMA 5.2 CORRIDA EN EL LENGUAJE FORTRAN (EJEMPLO 2).

Ahora que el lector está familiarizado con la manera de correr el programa, introducir los

datos, etc. Se expondrá un segundo ejemplo tomado de Montgomery (1999) ejemplo 5-13

de la pagina 153.

Un ingeniero industrial realiza un experimento para estudiar el tiempo en que tarda

un ojo en enfocar. Esta interesado en la relación que existe entre la distancia del objeto al

ojo y el tiempo que el objeto tarda en enfocar. Cuatro diferentes distancias resultan de

interés. Hay cinco sujetos disponibles para el experimento. Como puede haber diferencias

entre los sujetos, él decide efectuar un diseño aleatorizado por bloques. Los datos

recopilados se muestran en la tabla 5.3.

Sujeto

Distancia (ft) 1 2 3 4 5 4 10 6 6 6 6 6 7 6 6 1 6 8 5 3 3 2 5

10 6 4 4 2 3 Tabla 5.3

El contraste que se desea probar es:

CAPÍTULO 5. RESULTADOS

77

4321 µµµµ +=+

5.2.1 CORRIDA PARA EL MODELO N (0,1).

Siguiendo a pie de la letra las etapas descritas en la sección 5.1, se obtiene como archivo de

salida:

PROGRAMA PARA INVESTIGAR LA COMPARACION DE DOS MODELOS DE BOOTSTRAP. PROGRAMADO POR: ESTEBAN BURGUETE HERNANDEZ, JOSE FCO. TAMBORERO ARNAL. MANUEL ALEJANDRO PEREZ MIRELES. ARCHIVO DE ENTRADA: MON513.DAT ARCHIVO DE SALIDA: MON513.LIS SEMILLA PARA SIMULACION: 1231232142 NUMERO TOTAL DE SIMULACIONES: 30 ITERACIONES POR SIMULACION: 1000000 NIVEL DE CONFIANZA: 95.00 Z DE TABLAS: 1.96 DISEÑO: BLOQUES AL AZAR NUMERO DE TRATAMIENTOS: 4 NUMERO DE BLOQUES: 5 DATOS: TRAT( 1) BLQ( 1) 10.0000 TRAT( 1) BLQ( 2) 6.0000 TRAT( 1) BLQ( 3) 6.0000 TRAT( 1) BLQ( 4) 6.0000 TRAT( 1) BLQ( 5) 6.0000 TRAT( 2) BLQ( 1) 7.0000 TRAT( 2) BLQ( 2) 6.0000 TRAT( 2) BLQ( 3) 6.0000 TRAT( 2) BLQ( 4) 1.0000 TRAT( 2) BLQ( 5) 6.0000

CAPÍTULO 5. RESULTADOS

78

TRAT( 3) BLQ( 1) 5.0000 TRAT( 3) BLQ( 2) 3.0000 TRAT( 3) BLQ( 3) 3.0000 TRAT( 3) BLQ( 4) 2.0000 TRAT( 3) BLQ( 5) 5.0000 TRAT( 4) BLQ( 1) 6.0000 TRAT( 4) BLQ( 2) 4.0000 TRAT( 4) BLQ( 3) 4.0000 TRAT( 4) BLQ( 4) 2.0000 TRAT( 4) BLQ( 5) 3.0000 MODELO= 1 SI MODELO=1, EL MODELO ES EL N(0,1) TOTALES DE LOS TRATAMIENTOS. TOTAL DEL TRATAMIENTO: 1 34.0000 TOTAL DEL TRATAMIENTO: 2 26.0000 TOTAL DEL TRATAMIENTO: 3 18.0000 TOTAL DEL TRATAMIENTO: 4 19.0000 CUADRADO MEDIO DEL ERROR: 1.2750 RESULTADOS FINALES. COEFICIENTES DEL CONTRASTE: COEF. DEL TRAT( 1) 1.00 COEF. DEL TRAT( 2) 1.00 COEF. DEL TRAT( 3) -1.00 COEF. DEL TRAT( 4) -1.00 F CALCULADA DEL CONTRASTE: 20.7451 INTERVALO DE CONFIANZA DEL PVALUE DEL CONTRASTE: .0019< .0019< .0020 GRACIAS POR USAR ESTE PROGRAMA 5.2.2 CORRIDA PARA EL MODELO N (0, CME).

De igual manera se sigue a pie de la letra las etapas descritas en la sección 5.1, se obtiene

como archivo de salida para este modelo:

PROGRAMA PARA INVESTIGAR LA COMPARACION DE DOS MODELOS DE BOOTSTRAP.

CAPÍTULO 5. RESULTADOS

79

PROGRAMADO POR: ESTEBAN BURGUETE HERNANDEZ, JOSE FCO. TAMBORERO ARNAL. MANUEL ALEJANDRO PEREZ MIRELES. ARCHIVO DE ENTRADA: MON513.DAT ARCHIVO DE SALIDA: MON515.LIS SEMILLA PARA SIMULACION: 2132443243 NUMERO TOTAL DE SIMULACIONES: 30 ITERACIONES POR SIMULACION: 1000000 NIVEL DE CONFIANZA: 95.00 Z DE TABLAS: 1.96 DISEÑO: BLOQUES AL AZAR NUMERO DE TRATAMIENTOS: 4 NUMERO DE BLOQUES: 5 DATOS: TRAT( 1) BLQ( 1) 10.0000 TRAT( 1) BLQ( 2) 6.0000 TRAT( 1) BLQ( 3) 6.0000 TRAT( 1) BLQ( 4) 6.0000 TRAT( 1) BLQ( 5) 6.0000 TRAT( 2) BLQ( 1) 7.0000 TRAT( 2) BLQ( 2) 6.0000 TRAT( 2) BLQ( 3) 6.0000 TRAT( 2) BLQ( 4) 1.0000 TRAT( 2) BLQ( 5) 6.0000 TRAT( 3) BLQ( 1) 5.0000 TRAT( 3) BLQ( 2) 3.0000 TRAT( 3) BLQ( 3) 3.0000 TRAT( 3) BLQ( 4) 2.0000 TRAT( 3) BLQ( 5) 5.0000 TRAT( 4) BLQ( 1) 6.0000 TRAT( 4) BLQ( 2) 4.0000 TRAT( 4) BLQ( 3) 4.0000 TRAT( 4) BLQ( 4) 2.0000 TRAT( 4) BLQ( 5) 3.0000 MODELO= 2 SI MODELO=1, EL MODELO ES EL N(0,1) TOTALES DE LOS TRATAMIENTOS. TOTAL DEL TRATAMIENTO: 1 34.0000 TOTAL DEL TRATAMIENTO: 2 26.0000 TOTAL DEL TRATAMIENTO: 3 18.0000 TOTAL DEL TRATAMIENTO: 4 19.0000

CAPÍTULO 5. RESULTADOS

80

CUADRADO MEDIO DEL ERROR: 1.2750 RESULTADOS FINALES. COEFICIENTES DEL CONTRASTE: COEF. DEL TRAT( 1) 1.00 COEF. DEL TRAT( 2) 1.00 COEF. DEL TRAT( 3) -1.00 COEF. DEL TRAT( 4) -1.00 F CALCULADA DEL CONTRASTE: 20.7451 INTERVALO DE CONFIANZA DEL PVALUE DEL CONTRASTE: .0019< .0019< .0020 GRACIAS POR USAR ESTE PROGRAMA

5.3 CORRIDA EN EL LENGUAJE FORTRAN (EJEMPLO 3).

Sea el ejemplo 9 de la página 142 de Gutiérrez y de la Vara (2003), se tiene únicamente la

tabla 5.4, sea:

Tratamiento

Bloque A B C 1 3 7 4 2 4 9 6 3 2 3 3 4 6 10 7

Tabla 5.4

El contraste a probar es:

321 2µµµ =+

5.3.1 CORRIDA PARA EL MODELO N(0,1).

El archivo de salida que arroja Fortran para este ejemplo se presenta a continuación,

recalcando que se toma el mismo número de repeticiones por simulación, un millón de

iteraciones y treinta simulaciones para estimar el intervalo de confianza del valor P, sea:

CAPÍTULO 5. RESULTADOS

81

PROGRAMA PARA INVESTIGAR LA COMPARACION DE DOS MODELOS DE BOOTSTRAP. PROGRAMADO POR: ESTEBAN BURGUETE HERNANDEZ, JOSE FCO. TAMBORERO ARNAL. MANUEL ALEJANDRO PEREZ MIRELES. ARCHIVO DE ENTRADA: gut9.dat ARCHIVO DE SALIDA: gut91.lis SEMILLA PARA SIMULACION: 1232323 NUMERO TOTAL DE SIMULACIONES: 30 ITERACIONES POR SIMULACION: 1000000 NIVEL DE CONFIANZA: 95.00 Z DE TABLAS: 1.96 DISEÑO: BLOQUES AL AZAR NUMERO DE TRATAMIENTOS: 3 NUMERO DE BLOQUES: 4 DATOS: TRAT( 1) BLQ( 1) 3.0000 TRAT( 1) BLQ( 2) 4.0000 TRAT( 1) BLQ( 3) 2.0000 TRAT( 1) BLQ( 4) 6.0000 TRAT( 2) BLQ( 1) 7.0000 TRAT( 2) BLQ( 2) 9.0000 TRAT( 2) BLQ( 3) 3.0000 TRAT( 2) BLQ( 4) 10.0000 TRAT( 3) BLQ( 1) 4.0000 TRAT( 3) BLQ( 2) 6.0000 TRAT( 3) BLQ( 3) 3.0000 TRAT( 3) BLQ( 4) 7.0000 MODELO= 1 SI MODELO=1, EL MODELO ES EL N(0,1) TOTALES DE LOS TRATAMIENTOS. TOTAL DEL TRATAMIENTO: 1 15.0000 TOTAL DEL TRATAMIENTO: 2 29.0000

CAPÍTULO 5. RESULTADOS

82

TOTAL DEL TRATAMIENTO: 3 20.0000 CUADRADO MEDIO DEL ERROR: .9167 RESULTADOS FINALES. COEFICIENTES DEL CONTRASTE: COEF. DEL TRAT( 1) 1.00 COEF. DEL TRAT( 2) 1.00 COEF. DEL TRAT( 3) -2.00 F CALCULADA DEL CONTRASTE: .7273 INTERVALO DE CONFIANZA DEL PVALUE DEL CONTRASTE: .0875< .0880< .0886 GRACIAS POR USAR ESTE PROGRAMA 5.3.2 CORRIDA PARA EL MODELO N (0, CME).

Se presenta igualmente los resultados contenidos en el archivo de salida para el modelo

bloques al zar con observaciones simuladas que se distribuyen N(0, CME), sea entonces:

PROGRAMA PARA INVESTIGAR LA COMPARACION DE DOS MODELOS DE BOOTSTRAP. PROGRAMADO POR: ESTEBAN BURGUETE HERNANDEZ, JOSE FCO. TAMBORERO ARNAL. MANUEL ALEJANDRO PEREZ MIRELES. ARCHIVO DE ENTRADA: gut9.dat ARCHIVO DE SALIDA: gut93.lis SEMILLA PARA SIMULACION: 1223234 NUMERO TOTAL DE SIMULACIONES: 30 ITERACIONES POR SIMULACION: 1000000 NIVEL DE CONFIANZA: 95.00 Z DE TABLAS: 1.96 DISEÑO: BLOQUES AL AZAR NUMERO DE TRATAMIENTOS: 3 NUMERO DE BLOQUES: 4

CAPÍTULO 5. RESULTADOS

83

DATOS: TRAT( 1) BLQ( 1) 3.0000 TRAT( 1) BLQ( 2) 4.0000 TRAT( 1) BLQ( 3) 2.0000 TRAT( 1) BLQ( 4) 6.0000 TRAT( 2) BLQ( 1) 7.0000 TRAT( 2) BLQ( 2) 9.0000 TRAT( 2) BLQ( 3) 3.0000 TRAT( 2) BLQ( 4) 10.0000 TRAT( 3) BLQ( 1) 4.0000 TRAT( 3) BLQ( 2) 6.0000 TRAT( 3) BLQ( 3) 3.0000 TRAT( 3) BLQ( 4) 7.0000 MODELO= 2 SI MODELO=1, EL MODELO ES EL N(0,1) TOTALES DE LOS TRATAMIENTOS. TOTAL DEL TRATAMIENTO: 1 15.0000 TOTAL DEL TRATAMIENTO: 2 29.0000 TOTAL DEL TRATAMIENTO: 3 20.0000 CUADRADO MEDIO DEL ERROR: .9167 RESULTADOS FINALES. COEFICIENTES DEL CONTRASTE: COEF. DEL TRAT( 1) 1.00 COEF. DEL TRAT( 2) 1.00 COEF. DEL TRAT( 3) -2.00 F CALCULADA DEL CONTRASTE: .7273 INTERVALO DE CONFIANZA DEL PVALUE DEL CONTRASTE: .0876< .0881< .0885 GRACIAS POR USAR ESTE PROGRAMA 5.4 CORRIDA EN EL LENGUAJE FORTRAN (EJEMPLO 4).

Sea el ejercicio 10 de la página 143 de Gutiérrez y de la Vara (2003).

Se hace un estudio sobre la efectividad de tres marcas de atomizador para matar

moscas, Para ello, cada atomizador se aplica a un grupo de 100 moscas, y se cuenta el

número de moscas muertas, expresado en porcentajes. Se hicieron seis réplicas, pero éstas

CAPÍTULO 5. RESULTADOS

84

se hicieron en días diferentes, por ello se sospecha que puede haber algún efecto

importante debido a esta fuente de vacación. Los datos obtenidos se muestran en la tabla

5.5.

Numero de réplica

Marca de

atomizador 1 2 3 4 5 6

1 72 65 67 75 62 73 2 55 59 68 70 53 50 3 64 74 61 58 51 69

Tabla 5.5

Se desea probar el contraste: 3212 µµµ +=

5.4.1 CORRIDA PARA EL MODELO N(0,1).

Los resultados son los siguientes:

PROGRAMA PARA INVESTIGAR LA COMPARACION DE DOS MODELOS DE BOOTSTRAP. PROGRAMADO POR: ESTEBAN BURGUETE HERNANDEZ, JOSE FCO. TAMBORERO ARNAL. MANUEL ALEJANDRO PEREZ MIRELES. ARCHIVO DE ENTRADA: gut10.dat ARCHIVO DE SALIDA: gut101.lis SEMILLA PARA SIMULACION: 1232344 NUMERO TOTAL DE SIMULACIONES: 30 ITERACIONES POR SIMULACION: 1000000 NIVEL DE CONFIANZA: 95.00 Z DE TABLAS: 1.96

CAPÍTULO 5. RESULTADOS

85

DISEÑO: BLOQUES AL AZAR NUMERO DE TRATAMIENTOS: 3 NUMERO DE BLOQUES: 6 DATOS: TRAT( 1) BLQ( 1) 72.0000 TRAT( 1) BLQ( 2) 65.0000 TRAT( 1) BLQ( 3) 67.0000 TRAT( 1) BLQ( 4) 75.0000 TRAT( 1) BLQ( 5) 62.0000 TRAT( 1) BLQ( 6) 73.0000 TRAT( 2) BLQ( 1) 55.0000 TRAT( 2) BLQ( 2) 59.0000 TRAT( 2) BLQ( 3) 68.0000 TRAT( 2) BLQ( 4) 70.0000 TRAT( 2) BLQ( 5) 53.0000 TRAT( 2) BLQ( 6) 50.0000 TRAT( 3) BLQ( 1) 64.0000 TRAT( 3) BLQ( 2) 74.0000 TRAT( 3) BLQ( 3) 61.0000 TRAT( 3) BLQ( 4) 58.0000 TRAT( 3) BLQ( 5) 51.0000 TRAT( 3) BLQ( 6) 69.0000 MODELO= 1 SI MODELO=1, EL MODELO ES EL N(0,1) TOTALES DE LOS TRATAMIENTOS. TOTAL DEL TRATAMIENTO: 1 414.0000 TOTAL DEL TRATAMIENTO: 2 355.0000 TOTAL DEL TRATAMIENTO: 3 377.0000 CUADRADO MEDIO DEL ERROR: 51.4333 RESULTADOS FINALES. COEFICIENTES DEL CONTRASTE: COEF. DEL TRAT( 1) -2.00 COEF. DEL TRAT( 2) 1.00 COEF. DEL TRAT( 3) 1.00 F CALCULADA DEL CONTRASTE: 4.9773 INTERVALO DE CONFIANZA DEL PVALUE DEL CONTRASTE: .0739< .0744< .0749 GRACIAS POR USAR ESTE PROGRAMA

CAPÍTULO 5. RESULTADOS

86

5.4.2 CORRIDA PARA EL MODELO N (0, CME).

Se presenta el archivo de salida para el modelo N(0, CME), tomamos el mismo ejemplo y

los mismos parámetros que el modelo del apartado anterior, sea entonces:

PROGRAMA PARA INVESTIGAR LA COMPARACION DE DOS MODELOS DE BOOTSTRAP. PROGRAMADO POR: ESTEBAN BURGUETE HERNANDEZ, JOSE FCO. TAMBORERO ARNAL. MANUEL ALEJANDRO PEREZ MIRELES. ARCHIVO DE ENTRADA: gut10.dat ARCHIVO DE SALIDA: gut102.lis SEMILLA PARA SIMULACION: 12323123 NUMERO TOTAL DE SIMULACIONES: 30 ITERACIONES POR SIMULACION: 1000000 NIVEL DE CONFIANZA: 95.00 Z DE TABLAS: 1.96 DISEÑO: BLOQUES AL AZAR NUMERO DE TRATAMIENTOS: 3 NUMERO DE BLOQUES: 6 DATOS: TRAT( 1) BLQ( 1) 72.0000 TRAT( 1) BLQ( 2) 65.0000 TRAT( 1) BLQ( 3) 67.0000 TRAT( 1) BLQ( 4) 75.0000 TRAT( 1) BLQ( 5) 62.0000 TRAT( 1) BLQ( 6) 73.0000 TRAT( 2) BLQ( 1) 55.0000 TRAT( 2) BLQ( 2) 59.0000 TRAT( 2) BLQ( 3) 68.0000 TRAT( 2) BLQ( 4) 70.0000 TRAT( 2) BLQ( 5) 53.0000 TRAT( 2) BLQ( 6) 50.0000 TRAT( 3) BLQ( 1) 64.0000 TRAT( 3) BLQ( 2) 74.0000 TRAT( 3) BLQ( 3) 61.0000 TRAT( 3) BLQ( 4) 58.0000 TRAT( 3) BLQ( 5) 51.0000 TRAT( 3) BLQ( 6) 69.0000 MODELO= 2

CAPÍTULO 5. RESULTADOS

87

SI MODELO=1, EL MODELO ES EL N(0,1) TOTALES DE LOS TRATAMIENTOS. TOTAL DEL TRATAMIENTO: 1 414.0000 TOTAL DEL TRATAMIENTO: 2 355.0000 TOTAL DEL TRATAMIENTO: 3 377.0000 CUADRADO MEDIO DEL ERROR: 51.4333 RESULTADOS FINALES. COEFICIENTES DEL CONTRASTE: COEF. DEL TRAT( 1) -2.00 COEF. DEL TRAT( 2) 1.00 COEF. DEL TRAT( 3) 1.00 F CALCULADA DEL CONTRASTE: 4.9773 INTERVALO DE CONFIANZA DEL PVALUE DEL CONTRASTE: .0739< .0744< .0749 GRACIAS POR USAR ESTE PROGRAMA

5.5 CORRIDA EN EL LENGUAJE FORTRAN (EJEMPLO 5).

Sea el ejercicio 11 de la página 143 de Gutiérrez y de la Vara (2003).

En una empresa le tienen varios silos para almacenar lache (cisternas de 60000L).

Un aspecto crítico para conservar la leche es la temperatura de almacenamiento. Se

sospecha que en algunos sitios hay problemas, por ello durante cinco días se decide

registrar la temperatura a cierta hora crìtica. Obviamente la temperatura se un día a otro es

una variabilidad que podría impactar la variabilidad total. Ver tabla 5.6

Dìa

Silo Lun Mar Mier Ju Vi A 4 4 5 .5 3 B 5 6 2 4 4 C 4.5 4 3.5 2 3 D 2.5 4 6.5 4.5 4 E 4 4 3.5 2 4

Tabla 5.6

CAPÍTULO 5. RESULTADOS

88

Se desea probar el contraste: 543212 µµµµµ ++=+

5.5.1 CORRIDA PARA EL MODELO N (0, 1).

PROGRAMA PARA INVESTIGAR LA COMPARACION DE DOS MODELOS DE BOOTSTRAP. PROGRAMADO POR: ESTEBAN BURGUETE HERNANDEZ, JOSE FCO. TAMBORERO ARNAL. MANUEL ALEJANDRO PEREZ MIRELES. ARCHIVO DE ENTRADA: gut11.dat ARCHIVO DE SALIDA: gut12.lis SEMILLA PARA SIMULACION: 1223233 NUMERO TOTAL DE SIMULACIONES: 30 ITERACIONES POR SIMULACION: 1000000 NIVEL DE CONFIANZA: 95.00 Z DE TABLAS: 1.96 DISEÑO: BLOQUES AL AZAR NUMERO DE TRATAMIENTOS: 5 NUMERO DE BLOQUES: 5 DATOS: TRAT( 1) BLQ( 1) 4.0000 TRAT( 1) BLQ( 2) 4.0000 TRAT( 1) BLQ( 3) 5.0000 TRAT( 1) BLQ( 4) .5000 TRAT( 1) BLQ( 5) 3.0000 TRAT( 2) BLQ( 1) 5.0000 TRAT( 2) BLQ( 2) 6.0000 TRAT( 2) BLQ( 3) 2.0000 TRAT( 2) BLQ( 4) 4.0000 TRAT( 2) BLQ( 5) 4.0000 TRAT( 3) BLQ( 1) 4.5000 TRAT( 3) BLQ( 2) 4.0000 TRAT( 3) BLQ( 3) 3.5000 TRAT( 3) BLQ( 4) 2.0000 TRAT( 3) BLQ( 5) 3.0000 TRAT( 4) BLQ( 1) 2.5000 TRAT( 4) BLQ( 2) 4.0000 TRAT( 4) BLQ( 3) 6.5000 TRAT( 4) BLQ( 4) 4.5000 TRAT( 4) BLQ( 5) 4.0000 TRAT( 5) BLQ( 1) 4.0000 TRAT( 5) BLQ( 2) 4.0000

CAPÍTULO 5. RESULTADOS

89

TRAT( 5) BLQ( 3) 3.5000 TRAT( 5) BLQ( 4) 2.0000 TRAT( 5) BLQ( 5) 4.0000 MODELO= 1 SI MODELO=1, EL MODELO ES EL N(0,1) TOTALES DE LOS TRATAMIENTOS. TOTAL DEL TRATAMIENTO: 1 16.5000 TOTAL DEL TRATAMIENTO: 2 21.0000 TOTAL DEL TRATAMIENTO: 3 17.0000 TOTAL DEL TRATAMIENTO: 4 21.5000 TOTAL DEL TRATAMIENTO: 5 17.5000 CUADRADO MEDIO DEL ERROR: 1.6150 RESULTADOS FINALES. COEFICIENTES DEL CONTRASTE: COEF. DEL TRAT( 1) -2.00 COEF. DEL TRAT( 2) -1.00 COEF. DEL TRAT( 3) 1.00 COEF. DEL TRAT( 4) 1.00 COEF. DEL TRAT( 5) 1.00 F CALCULADA DEL CONTRASTE: .0619 INTERVALO DE CONFIANZA DEL PVALUE DEL CONTRASTE: .9637< .9645< .9652 GRACIAS POR USAR ESTE PROGRAMA 5.5.2 CORRIDA PARA EL MODELO N (0, CME).

PROGRAMA PARA INVESTIGAR LA COMPARACION DE DOS MODELOS DE BOOTSTRAP. PROGRAMADO POR: ESTEBAN BURGUETE HERNANDEZ, JOSE FCO. TAMBORERO ARNAL. MANUEL ALEJANDRO PEREZ MIRELES. ARCHIVO DE ENTRADA: gut11.dat ARCHIVO DE SALIDA: gut13.lis

CAPÍTULO 5. RESULTADOS

90

SEMILLA PARA SIMULACION: 1232214 NUMERO TOTAL DE SIMULACIONES: 30 ITERACIONES POR SIMULACION: 1000000 NIVEL DE CONFIANZA: 95.00 Z DE TABLAS: 1.96 DISEÑO: BLOQUES AL AZAR NUMERO DE TRATAMIENTOS: 5 NUMERO DE BLOQUES: 5 DATOS: TRAT( 1) BLQ( 1) 4.0000 TRAT( 1) BLQ( 2) 4.0000 TRAT( 1) BLQ( 3) 5.0000 TRAT( 1) BLQ( 4) .5000 TRAT( 1) BLQ( 5) 3.0000 TRAT( 2) BLQ( 1) 5.0000 TRAT( 2) BLQ( 2) 6.0000 TRAT( 2) BLQ( 3) 2.0000 TRAT( 2) BLQ( 4) 4.0000 TRAT( 2) BLQ( 5) 4.0000 TRAT( 3) BLQ( 1) 4.5000 RAT( 3) BLQ( 2) 4.0000 TRAT( 3) BLQ( 3) 3.5000 TRAT( 3) BLQ( 4) 2.0000 TRAT( 3) BLQ( 5) 3.0000 TRAT( 4) BLQ( 1) 2.5000 TRAT( 4) BLQ( 2) 4.0000 TRAT( 4) BLQ( 3) 6.5000 TRAT( 4) BLQ( 4) 4.5000 TRAT( 4) BLQ( 5) 4.0000 TRAT( 5) BLQ( 1) 4.0000 TRAT( 5) BLQ( 2) 4.0000 TRAT( 5) BLQ( 3) 3.5000 TRAT( 5) BLQ( 4) 2.0000 TRAT( 5) BLQ( 5) 4.0000 MODELO= 2 SI MODELO=1, EL MODELO ES EL N(0,1) TOTALES DE LOS TRATAMIENTOS. TOTAL DEL TRATAMIENTO: 1 16.5000 TOTAL DEL TRATAMIENTO: 2 21.0000 TOTAL DEL TRATAMIENTO: 3 17.0000 TOTAL DEL TRATAMIENTO: 4 21.5000 TOTAL DEL TRATAMIENTO: 5 17.5000

CAPÍTULO 5. RESULTADOS

91

CUADRADO MEDIO DEL ERROR: 1.6150 RESULTADOS FINALES. COEFICIENTES DEL CONTRASTE: COEF. DEL TRAT( 1) -2.00 COEF. DEL TRAT( 2) -1.00 COEF. DEL TRAT( 3) 1.00 COEF. DEL TRAT( 4) 1.00 COEF. DEL TRAT( 5) 1.00 F CALCULADA DEL CONTRASTE: .0619 INTERVALO DE CONFIANZA DEL PVALUE DEL CONTRASTE: .9642< .9644< .9647 GRACIAS POR USAR ESTE PROGRAMA 5.6 CORRIDA EN EL LENGUAJE FORTRAN (EJEMPLO 6).

Sea el ejercicio 12 de la página 144 de Gutiérrez y de la Vara (2003).

Se diseño un experimento para estudiar el rendimiento de cuatro detergentes. Las

siguientes lecturas de “blancura” se obtuvieron con un equipo especial diseñado para 12

cargas de lavado distribuidas en tres modelos de lavadoras. Ver tabla 5.7

Tratamiento

Detergente Lavadora 1 Lavadora 2 Lavadora 3 A 45 43 51 B 47 44 52 C 50 49 57 D 42 37 49

Tabla 5.7

Se desea probar el contraste: 5321 3µµµµ =++

5.6.1 CORRIDA PARA EL MODELO N (0, 1).

CAPÍTULO 5. RESULTADOS

92

PROGRAMA PARA INVESTIGAR LA COMPARACION DE DOS MODELOS DE BOOTSTRAP. PROGRAMADO POR: ESTEBAN BURGUETE HERNANDEZ, JOSE FCO. TAMBORERO ARNAL. MANUEL ALEJANDRO PEREZ MIRELES. ARCHIVO DE ENTRADA: gute12.dat ARCHIVO DE SALIDA: gute12.lis SEMILLA PARA SIMULACION: 12323213 NUMERO TOTAL DE SIMULACIONES: 30 ITERACIONES POR SIMULACION: 1000000 NIVEL DE CONFIANZA: 95.00 Z DE TABLAS: 1.96 DISEÑO: BLOQUES AL AZAR NUMERO DE TRATAMIENTOS: 4 NUMERO DE BLOQUES: 3 DATOS: TRAT( 1) BLQ( 1) 45.0000 TRAT( 1) BLQ( 2) 43.0000 TRAT( 1) BLQ( 3) 51.0000 TRAT( 2) BLQ( 1) 47.0000 TRAT( 2) BLQ( 2) 44.0000 TRAT( 2) BLQ( 3) 52.0000 TRAT( 3) BLQ( 1) 50.0000 TRAT( 3) BLQ( 2) 49.0000 TRAT( 3) BLQ( 3) 57.0000 TRAT( 4) BLQ( 1) 42.0000 TRAT( 4) BLQ( 2) 37.0000 TRAT( 4) BLQ( 3) 49.0000 MODELO= 1 SI MODELO=1, EL MODELO ES EL N(0,1) TOTALES DE LOS TRATAMIENTOS. TOTAL DEL TRATAMIENTO: 1 139.0000 TOTAL DEL TRATAMIENTO: 2 143.0000 TOTAL DEL TRATAMIENTO: 3 156.0000 TOTAL DEL TRATAMIENTO: 4 128.0000 CUADRADO MEDIO DEL ERROR: 1.3056

CAPÍTULO 5. RESULTADOS

93

RESULTADOS FINALES. COEFICIENTES DEL CONTRASTE: COEF. DEL TRAT( 1) 1.00 COEF. DEL TRAT( 2) 1.00 COEF. DEL TRAT( 3) 1.00 COEF. DEL TRAT( 4) -3.00 F CALCULADA DEL CONTRASTE: 62.0426 INTERVALO DE CONFIANZA DEL PVALUE DEL CONTRASTE: .0004< .0004< .0005 GRACIAS POR USAR ESTE PROGRAMA

5.6.2 CORRIDA PARA EL MODELO N (0, CME).

PROGRAMA PARA INVESTIGAR LA COMPARACION DE DOS MODELOS DE BOOTSTRAP. PROGRAMADO POR: ESTEBAN BURGUETE HERNANDEZ, JOSE FCO. TAMBORERO ARNAL. MANUEL ALEJANDRO PEREZ MIRELES. ARCHIVO DE ENTRADA: gute12.dat ARCHIVO DE SALIDA: gute13.lis SEMILLA PARA SIMULACION: 1212312 NUMERO TOTAL DE SIMULACIONES: 30 ITERACIONES POR SIMULACION: 1000000 NIVEL DE CONFIANZA: 95.00 Z DE TABLAS: 1.96 DISEÑO: BLOQUES AL AZAR NUMERO DE TRATAMIENTOS: 4 NUMERO DE BLOQUES: 3 DATOS: TRAT( 1) BLQ( 1) 45.0000 TRAT( 1) BLQ( 2) 43.0000 TRAT( 1) BLQ( 3) 51.0000 TRAT( 2) BLQ( 1) 47.0000 TRAT( 2) BLQ( 2) 44.0000 TRAT( 2) BLQ( 3) 52.0000 TRAT( 3) BLQ( 1) 50.0000 TRAT( 3) BLQ( 2) 49.0000 TRAT( 3) BLQ( 3) 57.0000

CAPÍTULO 5. RESULTADOS

94

TRAT( 4) BLQ( 1) 42.0000 TRAT( 4) BLQ( 2) 37.0000 TRAT( 4) BLQ( 3) 49.0000 MODELO= 2 SI MODELO=1, EL MODELO ES EL N(0,1) TOTALES DE LOS TRATAMIENTOS. TOTAL DEL TRATAMIENTO: 1 139.0000 TOTAL DEL TRATAMIENTO: 2 143.0000 TOTAL DEL TRATAMIENTO: 3 156.0000 TOTAL DEL TRATAMIENTO: 4 128.0000 CUADRADO MEDIO DEL ERROR: 1.3056 RESULTADOS FINALES. COEFICIENTES DEL CONTRASTE: COEF. DEL TRAT( 1) 1.00 COEF. DEL TRAT( 2) 1.00 COEF. DEL TRAT( 3) 1.00 COEF. DEL TRAT( 4) -3.00 F CALCULADA DEL CONTRASTE: 62.0426 INTERVALO DE CONFIANZA DEL PVALUE DEL CONTRASTE: .0004< .0004< .0005 GRACIAS POR USAR ESTE PROGRAMA

5.7 CORRIDA EN EL LENGUAJE FORTRAN (EJEMPLO 7)

Sea el ejemplo 14.10 del Mendenhall y Sincich (1955) página 832:

Prior to submitting a bid for a construction job, cost engineers prepare a detailed

analysis of the estimated labor and materials costs require to complete the job. This

estimate will depend on the engineer who performs the analysis. An overly large estimate

Hill reduce the chance of acceptance of a company’s bid price, whereas an estimate that is

too low Hill reduce the profit or even cause the company to lose money on the job. A

company that employs three job cost engineers wanted to compare the mean level of the

CAPÍTULO 5. RESULTADOS

95

engineer’s estimates. This was done by having each engineer estimate the cost of the same

four jobs. The data (in hundreds of thousands of dollars) are shown in Table 5.8

Job

Engineer 1 2 3 4

1 4.6 6.2 5 6.6 2 4.9 6.3 5.4 6.8 3 4.4 5.9 5.4 6.3

Tabla 5.8

Se desea probar el contraste: 3212 µµµ += 5.7.1 CORRIDA PARA EL MODELO N (0, 1).

PROGRAMA PARA INVESTIGAR LA COMPARACION DE DOS MODELOS DE BOOTSTRAP. PROGRAMADO POR: ESTEBAN BURGUETE HERNANDEZ, JOSE FCO. TAMBORERO ARNAL. MANUEL ALEJANDRO PEREZ MIRELES. ARCHIVO DE ENTRADA: mende14.dat ARCHIVO DE SALIDA: mende14.lis SEMILLA PARA SIMULACION: 123213214 NUMERO TOTAL DE SIMULACIONES: 30 ITERACIONES POR SIMULACION: 1000000 NIVEL DE CONFIANZA: 95.00 Z DE TABLAS: 1.96 DISEÑO: BLOQUES AL AZAR NUMERO DE TRATAMIENTOS: 3 NUMERO DE BLOQUES: 4 DATOS: TRAT( 1) BLQ( 1) 4.6000 TRAT( 1) BLQ( 2) 6.2000 TRAT( 1) BLQ( 3) 5.0000 TRAT( 1) BLQ( 4) 6.6000 TRAT( 2) BLQ( 1) 4.9000 TRAT( 2) BLQ( 2) 6.3000 TRAT( 2) BLQ( 3) 5.4000 TRAT( 2) BLQ( 4) 6.8000

CAPÍTULO 5. RESULTADOS

96

TRAT( 3) BLQ( 1) 4.4000 TRAT( 3) BLQ( 2) 5.9000 TRAT( 3) BLQ( 3) 5.4000 TRAT( 3) BLQ( 4) 6.3000 MODELO= 1 SI MODELO=1, EL MODELO ES EL N(0,1) TOTALES DE LOS TRATAMIENTOS. TOTAL DEL TRATAMIENTO: 1 22.4000 TOTAL DEL TRATAMIENTO: 2 23.4000 TOTAL DEL TRATAMIENTO: 3 22.0000 CUADRADO MEDIO DEL ERROR: .0311 RESULTADOS FINALES. COEFICIENTES DEL CONTRASTE: COEF. DEL TRAT( 1) -2.00 COEF. DEL TRAT( 2) 1.00 COEF. DEL TRAT( 3) 1.00 F CALCULADA DEL CONTRASTE: .4820 INTERVALO DE CONFIANZA DEL PVALUE DEL CONTRASTE: .1428< .1434< .1440 GRACIAS POR USAR ESTE PROGRAMA 5.7.2 CORRIDA PARA EL MODELO N (0, CME).

PROGRAMA PARA INVESTIGAR LA COMPARACION DE DOS MODELOS DE BOOTSTRAP. PROGRAMADO POR: ESTEBAN BURGUETE HERNANDEZ, JOSE FCO. TAMBORERO ARNAL. MANUEL ALEJANDRO PEREZ MIRELES. ARCHIVO DE ENTRADA: mende14.dat ARCHIVO DE SALIDA: mende141.lis SEMILLA PARA SIMULACION: 123123123 NUMERO TOTAL DE SIMULACIONES: 30 ITERACIONES POR SIMULACION: 1000000

CAPÍTULO 5. RESULTADOS

97

NIVEL DE CONFIANZA: 95.00 Z DE TABLAS: 1.96 DISEÑO: BLOQUES AL AZAR NUMERO DE TRATAMIENTOS: 3 NUMERO DE BLOQUES: 4 DATOS: TRAT( 1) BLQ( 1) 4.6000 TRAT( 1) BLQ( 2) 6.2000 TRAT( 1) BLQ( 3) 5.0000 TRAT( 1) BLQ( 4) 6.6000 TRAT( 2) BLQ( 1) 4.9000 TRAT( 2) BLQ( 2) 6.3000 TRAT( 2) BLQ( 3) 5.4000 TRAT( 2) BLQ( 4) 6.8000 TRAT( 3) BLQ( 1) 4.4000 TRAT( 3) BLQ( 2) 5.9000 TRAT( 3) BLQ( 3) 5.4000 TRAT( 3) BLQ( 4) 6.3000 MODELO= 2 SI MODELO=1, EL MODELO ES EL N(0,1) TOTALES DE LOS TRATAMIENTOS. TOTAL DEL TRATAMIENTO: 1 22.4000 TOTAL DEL TRATAMIENTO: 2 23.4000 TOTAL DEL TRATAMIENTO: 3 22.0000 CUADRADO MEDIO DEL ERROR: .0311 RESULTADOS FINALES. COEFICIENTES DEL CONTRASTE: COEF. DEL TRAT( 1) -2.00 COEF. DEL TRAT( 2) 1.00 COEF. DEL TRAT( 3) 1.00 F CALCULADA DEL CONTRASTE: .4820 INTERVALO DE CONFIANZA DEL PVALUE DEL CONTRASTE: .1426< .1433< .1440 GRACIAS POR USAR ESTE PROGRAMA

CAPÍTULO 5. RESULTADOS

98

5.8 CORRIDA EN EL LENGUAJE FORTRAN (EJEMPLO 8).

Sea el ejercisio 14.32 del Mendenhall y Sincich (1955) página 839: Refer to Microcomputers in Civil Engineering study of lateral drift in a building.

The lateral displacements (in inches) estimated by three different computer programs are

recorded in the table for each of five building levels: 1, 5. 10, 15, and 21. Use the ANOVA

formulas to analyze the data. Recall that the goal is to compare the mean drift ratios

estimated by the three computer programs, STAAD-III (run1), STAAD (run 2), and

DRIFT. Use α=0.05. See table 5.9

Tratamiento

Level STAAD-III STAAD-III (2) Drift 1 .17 .16 .16 2 1.35 1.26 1.27 3 3.04 2.76 2.77 4 4.54 3.98 3.99 5 5.94 4.99 5

Tabla 5.9

Se desea probar el contraste: 321 2µµµ =+

5.8.1 CORRIDA PARA EL MODELO N (0, 1).

PROGRAMA PARA INVESTIGAR LA COMPARACION DE DOS MODELOS DE BOOTSTRAP. PROGRAMADO POR: ESTEBAN BURGUETE HERNANDEZ, JOSE FCO. TAMBORERO ARNAL. MANUEL ALEJANDRO PEREZ MIRELES. ARCHIVO DE ENTRADA: mend32.dat ARCHIVO DE SALIDA: mend33.lis

CAPÍTULO 5. RESULTADOS

99

SEMILLA PARA SIMULACION: 12321312 NUMERO TOTAL DE SIMULACIONES: 30 ITERACIONES POR SIMULACION: 1000000 NIVEL DE CONFIANZA: 95.00 Z DE TABLAS: 1.96 DISEÑO: BLOQUES AL AZAR NUMERO DE TRATAMIENTOS: 3 NUMERO DE BLOQUES: 5 DATOS: TRAT( 1) BLQ( 1) .1700 TRAT( 1) BLQ( 2) 1.3500 TRAT( 1) BLQ( 3) 3.0400 TRAT( 1) BLQ( 4) 4.5400 TRAT( 1) BLQ( 5) 5.9400 TRAT( 2) BLQ( 1) .1600 TRAT( 2) BLQ( 2) 1.2600 TRAT( 2) BLQ( 3) 2.7600 TRAT( 2) BLQ( 4) 3.9800 TRAT( 2) BLQ( 5) 4.9900 TRAT( 3) BLQ( 1) .1600 TRAT( 3) BLQ( 2) 1.2700 TRAT( 3) BLQ( 3) 2.7700 TRAT( 3) BLQ( 4) 3.9900 TRAT( 3) BLQ( 5) 5.0000 MODELO= 1 SI MODELO=1, EL MODELO ES EL N(0,1) TOTALES DE LOS TRATAMIENTOS. TOTAL DEL TRATAMIENTO: 1 15.0400 TOTAL DEL TRATAMIENTO: 2 13.1500 TOTAL DEL TRATAMIENTO: 3 13.1900 CUADRADO MEDIO DEL ERROR: .0487 RESULTADOS FINALES. COEFICIENTES DEL CONTRASTE: COEF. DEL TRAT( 1) 1.00 COEF. DEL TRAT( 2) 1.00 COEF. DEL TRAT( 3) -2.00 F CALCULADA DEL CONTRASTE: 2.2414 INTERVALO DE CONFIANZA DEL PVALUE DEL CONTRASTE: .0088< .0090< .0092

CAPÍTULO 5. RESULTADOS

100

5.8.2 CORRIDA PARA EL MODELO N (0, CME)

ROGRAMA PARA INVESTIGAR LA COMPARACION DE OS MODELOS DE BOOTSTRAP. ROGRAMADO POR: ESTEBAN BURGUETE HERNANDEZ, JOSE FCO. TAMBORERO ARNAL. MANUEL ALEJANDRO PEREZ MIRELES. ARCHIVO DE ENTRADA: mend32.dat ARCHIVO DE SALIDA: mend36.lis SEMILLA PARA SIMULACION: 12213123 NUMERO TOTAL DE SIMULACIONES: 30 ITERACIONES POR SIMULACION: 1000000 NIVEL DE CONFIANZA: 95.00 Z DE TABLAS: 1.96 DISEÑO: BLOQUES AL AZAR NUMERO DE TRATAMIENTOS: 3 NUMERO DE BLOQUES: 5 DATOS: TRAT( 1) BLQ( 1) .1700 TRAT( 1) BLQ( 2) 1.3500 TRAT( 1) BLQ( 3) 3.0400 TRAT( 1) BLQ( 4) 4.5400 TRAT( 1) BLQ( 5) 5.9400 TRAT( 2) BLQ( 1) .1600 TRAT( 2) BLQ( 2) 1.2600 TRAT( 2) BLQ( 3) 2.7600 TRAT( 2) BLQ( 4) 3.9800 TRAT( 2) BLQ( 5) 4.9900 TRAT( 3) BLQ( 1) .1600 TRAT( 3) BLQ( 2) 1.2700 TRAT( 3) BLQ( 3) 2.7700 TRAT( 3) BLQ( 4) 3.9900 TRAT( 3) BLQ( 5) 5.0000 MODELO= 2 SI MODELO=1, EL MODELO ES EL N(0,1) TOTALES DE LOS TRATAMIENTOS. TOTAL DEL TRATAMIENTO: 1 15.0400

CAPÍTULO 5. RESULTADOS

101

TOTAL DEL TRATAMIENTO: 2 13.1500 TOTAL DEL TRATAMIENTO: 3 13.1900 CUADRADO MEDIO DEL ERROR: .0487 RESULTADOS FINALES. COEFICIENTES DEL CONTRASTE: COEF. DEL TRAT( 1) 1.00 COEF. DEL TRAT( 2) 1.00 COEF. DEL TRAT( 3) -2.00 F CALCULADA DEL CONTRASTE: 2.2414 INTERVALO DE CONFIANZA DEL PVALUE DEL CONTRASTE: .0089< .0090< .0091 5.9 RESUMEN

A continuación se muestra la tabla 5.10 resumen donde se obtiene el valor P, el α real, con

su respectivo intervalo de confianza de cada ejemplo, con los dos tipos de modelos

simulados, sea:

Valor P Intervalo de confianza Ejemplo

N(0,1) N(0,CME) N(0,1) N(0,CME)

1 0.0995 0.0994 (0.0991; 0.1000) (0.0988; 0 .1001)

2 0.0019 0.0019 (0.0019; 0.0020) (0.0019; 0.0020)

3 0.0880 0.0881 (0.0875; 0.0886) (0.0876; 0.0885)

4 0.0744 0.0744 (0.0739; 0.0749) (0.0739; 0.0749)

5 0.9645 0.9644 (0.9637; 0.9652) (0.9642; 0.9647)

6 0.0004 0.0004 (0.0004; 0.0005) (0.0004; 0.0005)

7 0.1434 0.1433 (0.1428; 0.1440) (0.1426; 0.1440)

8 0.0090 0.0090 (0.0088; 0.0092) (0.0089; 0.0091)

Tabla 5.10