Sesión 1. Medidas de tendencia central

data=c(266.63,163.41,219.41,162.64,187.16,289.17,306.55,335.48,343.50,226.80,208.99,230.46)

#Media
media=mean(data)
media
## [1] 245.0167
#Mediana
mediana= median (data)
mediana
## [1] 228.63
#Moda
#En R no hay una funcion directa para la moda

#Rango 
rango = max(data)-min(data)
rango
## [1] 180.86
#Varianza
data1=data-media
data1
##  [1]  21.61333 -81.60667 -25.60667 -82.37667 -57.85667  44.15333  61.53333
##  [8]  90.46333  98.48333 -18.21667 -36.02667 -14.55667
data2= data1*data1

data3=sum(data2)
data3
## [1] 43375.91
varianza=data3/12
desviacion=sqrt(varianza)
desviacion
## [1] 60.12203

Sesión 2. Distribución normal

#Ejemplo 1

#a 
a=pnorm(600,1300,600)*100
a
## [1] 12.16725
#b
b=(pnorm(1500,1300,600)-pnorm(1000,1300,600))*100
b
## [1] 32.20211
#c
c= (1-pnorm(2200,1300,600))*100
c
## [1] 6.68072
#Problema 2

p=pnorm(21,18.7,5)*100
p
## [1] 67.72419
p2=100-p
p2
## [1] 32.27581
#Problema 3

a3=(1-pnorm(90,80,4))
a3
## [1] 0.006209665
b3=(pnorm(85,80,4) - pnorm(70,80,4))
b3
## [1] 0.8881406
c3=(1-pnorm(100,80,4))  * 1000 #numero de pilas
d3=(1-pnorm(90,80,4)) *1000 #numero de pilas

c3
## [1] 0.0002866516
d3
## [1] 6.209665

Sesión 3. Pruebas de hipótesis

Paso 1: Plantear la hipótesis
Paso 2: Nivel de significancia
Paso 3: Zona de aceptación/rechazo
Paso 4: Función pivotal
Paso 5: Conclusión

Sesión 4. Ejercicios del mundo real

Capítulo 3: Medidas de Tendencia central y dispersión

3-84

¿Qué reacción tendría usted si un aficionado al fútbol americano le dijera lo siguiente? “Los Raiders de Rockland tienen un promedio de 3.6 yardas de recorrido por tierra. Como sólo necesitan 10 yardas para anotar y tienen cuatro oportunidades para lograrlo, la anotación es segura, siempre y cuando mantengan su forma de jugar por tierra.”

Respuesta: La afirmación es incorrecta porque ignora completamente la variabilidad en yardas ganadas por carrera.

3-86

A continuación, se tienen 3 partes del presupuesto de defensa de un año, a cada una de éstas se le asignó, por parte del Congreso mexicano, la misma cantidad de financiamiento:
* Salario de oficiales (total).
* Mantenimiento de la flota aérea.
* Adquisiciones de alimentos (total).
Tomando en cuenta la distribución de posibles resultados para los gastos reales en cada una de éstas áreas, haga corresponder cada sección a una de las curvas de la siguiente figura, fundamente su respuesta.

Respuesta: Salarios de funcionarios: A; mantenimiento de flota: C; adquisiciones de alimentos: B.

3-92

El 30 de junio de 1992, la capitalización de nueve mercados de valores del Pacífico y Asia fue:
País Capitalización (En miles millones de dólares)
Filipinas 17
Indonesia 21
Tailandia 44
Singapur 50
Malasia 79
Corea del Sur 86
Taiwan 140
Hong Kong 178
Australia 203
a) Encuentre la media aritmética de los datos.
capitalizacion = c(17,21,44,50,79,86,140,178,203)
media = mean(capitalizacion)
media
## [1] 90.88889
b) Encuentre la mediana de los datos.
mediana = median(capitalizacion)
mediana
## [1] 79
c) Encuentre la moda de los datos.
#No hay función de moda y los datos al ser variables muy distintas, no tienen moda.
d) ¿Cuál es la mejor medida de la tendencia central del conjunto de datos?
histograma= hist(capitalizacion)

#La mejor medida es la mediana ya que hay un sesgo positivo como se aprecia en el histograma de los datos

e) Encuentre la desviación estándar de los datos. (La población completa está incluida en ellos.)
suma=capitalizacion-media
suma2=suma*suma
suma3=sum(suma2)
suma4=suma3/9
suma5=sqrt(suma4)
suma5
## [1] 64.07365

3-100

Matthews, Young y Asociados, una agencia de consultorías de Chapell Hill, tiene los siguientes registros que indican el número de días que cada uno de sus 10 consultores de planta cobró el último año:
212 220 230 210 228 229 231 219 221 222
a) Sin calcular el valor de las medidas, ¿cuál de ellas cree usted que le daría una mayor información acerca de esta distribución: el rango (alcance) o la desviación estándar?

Respuesta: Desviación estándar.

b) Calcule: Rango, Varianza y Desviación Estándar.

Respuesta

dias= c(212, 220, 230, 210, 228, 229, 231, 219, 221, 222)
rango_dias= max(dias)-min(dias)
media_dias= mean(dias) 
#Desviación estandar
dias2= dias= dias - media_dias
dias3= dias2*dias2
dias4=sum(dias3)
dias5= dias4/10
dias6_desv=sqrt(dias5)

#Respuestas
rango_dias
## [1] 21
media_dias
## [1] 222.2
#Desviación estandar
dias6_desv
## [1] 6.982836
c) Tomando en cuenta la dificultad y el tiempo para calcular cada medida que revisó en el inciso a), ¿cuál sugeriría como la mejor?

Respuesta: La desviación estandar

d) ¿Qué haría que usted cambiara su opinión al respecto?

Respuesta: Nada

3-106

Allison Barett realiza análisis estadísticos para un equipo de carreras automovilísticas. A continuación, se presentan las cifras en kilómetros por litro del gasto de combustible de sus automóviles en las carreras recientes:
4.77 6.11 6.11 5.05 5.99 4.91 5.27 6.01 5.75 4.89 6.05 5.22 6.02 5.24 6.11 5.02
a) Calcule la mediana del consumo de combustible.

Respuesta:

km = c(4.77, 6.11,  6.11,   5.05,   5.99,   4.91,   5.27,   6.01, 5.75, 4.89,   6.05,   5.22,   6.02,   5.24,   6.11,   5.02)
mediana_km = median(km)
mediana_km
## [1] 5.51
b) Calcule la media del mismo consumo.

Respuesta:

media_km = mean(km)
media_km
## [1] 5.5325
c) Agrupe los datos en 5 clases de igual tamaño. ¿Cuál es el intervalo del valor de consumo de combustible para la clase modal?

Respuesta:

clases <- cut(km, breaks=5)
clases2= table(clases)
clases2
## clases
## (4.77,5.04] (5.04,5.31] (5.31,5.57] (5.57,5.84] (5.84,6.11] 
##           4           4           0           1           7
d) ¿Cuál de las 3 medidas de tendencia central es la que mejor puede servirle a Allison cuando haga un pedido de combustible? Explique su respuesta.

Respuesta:

histograma_km= hist(km)

histograma_km
## $breaks
## [1] 4.6 4.8 5.0 5.2 5.4 5.6 5.8 6.0 6.2
## 
## $counts
## [1] 1 2 2 3 0 1 1 6
## 
## $density
## [1] 0.3125 0.6250 0.6250 0.9375 0.0000 0.3125 0.3125 1.8750
## 
## $mids
## [1] 4.7 4.9 5.1 5.3 5.5 5.7 5.9 6.1
## 
## $xname
## [1] "km"
## 
## $equidist
## [1] TRUE
## 
## attr(,"class")
## [1] "histogram"
#Depende 
e) ¿Cuál es el rango?

Respuesta:

rango= max(km) - min(km)
rango
## [1] 1.34

8-64

Inspectores del gobierno, al investigar los cargos levantados contra una embotelladora de bebidas no alcohólicas, de Texas, que no llenaba adecuadamente sus productos, han muestreado 200 botellas y encontraron que el promedio de llenado es de 31.7 onzas líquidas (aproximadamente 930 ml.) Se anuncia que las botellas contienen 32 onzas líquidas (946.33 ml.). Se sabe que la desviación estándar de la población es de 1.5 onzas líquidas (44.36 ml.). ¿Deberían concluir los inspectores, al nivel de significancia de 2%, que las botellas están siendo llenadas con menos contenido, y emitir una sanción?

Respuesta:

##### Paso 1: Plantear la hipótesis
#h0: x = 32
#h1: x ≠ 31.7

##### Paso 2: Nivel de significancia
# 99% (de cada lado) = 2.718 (z) 
#Alfa = 0.02 /2

##### Paso 3: Zona de aceptación/rechazo

# En el rango de -2.33 a 2.33

##### Paso 4: Función pivotal

z= (31.7-32) / (1.5/sqrt(200))
z
## [1] -2.828427
##### Paso 5: Conclusión
# Se rechaza H0 porque está afuera del rango, por lo tanto se concluye con que las botellas se llenan con menos contenido. 
LS0tCnRpdGxlOiAiV29ya3Nob3AxIgphdXRob3I6ICJFbGkgR2FicmllbCBIZXJuw6FuZGV6IE1lZGluYSAtIEEwMTI3NTc2MyIKZGF0ZTogIjIwMjMtMDUtMTgiCm91dHB1dDogCiAgaHRtbF9kb2N1bWVudDoKICAgIHRvYzogVFJVRQogICAgdG9jX2Zsb2F0OiBUUlVFCiAgICBjb2RlX2Rvd25sb2FkOiBUUlVFCi0tLQoKIVtdKC9Vc2Vycy9nYWJyaWVsbWVkaW5hL0Rvd25sb2Fkcy9yZWNpYm8tZGUtbHV6LTg3MHg2MDAuanBnKQoKIyMjIFNlc2nDs24gMS4gTWVkaWRhcyBkZSB0ZW5kZW5jaWEgY2VudHJhbAoKYGBge3J9CmRhdGE9YygyNjYuNjMsMTYzLjQxLDIxOS40MSwxNjIuNjQsMTg3LjE2LDI4OS4xNywzMDYuNTUsMzM1LjQ4LDM0My41MCwyMjYuODAsMjA4Ljk5LDIzMC40NikKCiNNZWRpYQptZWRpYT1tZWFuKGRhdGEpCm1lZGlhCgojTWVkaWFuYQptZWRpYW5hPSBtZWRpYW4gKGRhdGEpCm1lZGlhbmEKCiNNb2RhCiNFbiBSIG5vIGhheSB1bmEgZnVuY2lvbiBkaXJlY3RhIHBhcmEgbGEgbW9kYQoKI1JhbmdvIApyYW5nbyA9IG1heChkYXRhKS1taW4oZGF0YSkKcmFuZ28KCiNWYXJpYW56YQpkYXRhMT1kYXRhLW1lZGlhCmRhdGExCgpkYXRhMj0gZGF0YTEqZGF0YTEKCmRhdGEzPXN1bShkYXRhMikKZGF0YTMKCnZhcmlhbnphPWRhdGEzLzEyCmRlc3ZpYWNpb249c3FydCh2YXJpYW56YSkKZGVzdmlhY2lvbgoKYGBgCgohW10oL1VzZXJzL2dhYnJpZWxtZWRpbmEvRG93bmxvYWRzLzYzMTc1Yzc1ZDE2YjkyMjI3Zjc3MDYzZi5qcGcpCgojIyMgU2VzacOzbiAyLiBEaXN0cmlidWNpw7NuIG5vcm1hbAoKYGBge3J9CiNFamVtcGxvIDEKCiNhIAphPXBub3JtKDYwMCwxMzAwLDYwMCkqMTAwCmEKCiNiCmI9KHBub3JtKDE1MDAsMTMwMCw2MDApLXBub3JtKDEwMDAsMTMwMCw2MDApKSoxMDAKYgoKI2MKYz0gKDEtcG5vcm0oMjIwMCwxMzAwLDYwMCkpKjEwMApjCgojUHJvYmxlbWEgMgoKcD1wbm9ybSgyMSwxOC43LDUpKjEwMApwCgpwMj0xMDAtcApwMgoKI1Byb2JsZW1hIDMKCmEzPSgxLXBub3JtKDkwLDgwLDQpKQphMwoKYjM9KHBub3JtKDg1LDgwLDQpIC0gcG5vcm0oNzAsODAsNCkpCmIzCgpjMz0oMS1wbm9ybSgxMDAsODAsNCkpICAqIDEwMDAgI251bWVybyBkZSBwaWxhcwpkMz0oMS1wbm9ybSg5MCw4MCw0KSkgKjEwMDAgI251bWVybyBkZSBwaWxhcwoKYzMKZDMKCmBgYAoKIVtdKC9Vc2Vycy9nYWJyaWVsbWVkaW5hL0Rvd25sb2Fkcy9FU1RBLVVOSTMtTFBILTAwMS5wbmcpCgojIyMgU2VzacOzbiAzLiBQcnVlYmFzIGRlIGhpcMOzdGVzaXMKCiMjIyMjIFBhc28gMTogUGxhbnRlYXIgbGEgaGlww7N0ZXNpcwojIyMjIyBQYXNvIDI6IE5pdmVsIGRlIHNpZ25pZmljYW5jaWEKIyMjIyMgUGFzbyAzOiBab25hIGRlIGFjZXB0YWNpw7NuL3JlY2hhem8KIyMjIyMgUGFzbyA0OiBGdW5jacOzbiBwaXZvdGFsCiMjIyMjIFBhc28gNTogQ29uY2x1c2nDs24KCiFbXSgvVXNlcnMvZ2FicmllbG1lZGluYS9Eb3dubG9hZHMvZXN0YWRpc3RpY2EuanBnKQoKIyMjIFNlc2nDs24gNC4gRWplcmNpY2lvcyBkZWwgbXVuZG8gcmVhbAoKIyMjIyMgQ2Fww610dWxvIDM6IE1lZGlkYXMgZGUgVGVuZGVuY2lhIGNlbnRyYWwgeSBkaXNwZXJzacOzbgoKIyMjIyAqKjMtODQqKgojIyMjIyDCv1F1ZcyBIHJlYWNjaW/MgW4gdGVuZHJpzIFhIHVzdGVkIHNpIHVuIGFmaWNpb25hZG8gYWwgZnXMgXRib2wgYW1lcmljYW5vIGxlIGRpamVyYSBsbyBzaWd1aWVudGU/IOKAnExvcyBSYWlkZXJzIGRlIFJvY2tsYW5kIHRpZW5lbiB1biBwcm9tZWRpbyBkZSAzLjYgeWFyZGFzIGRlIHJlY29ycmlkbyBwb3IgdGllcnJhLiBDb21vIHNvzIFsbyBuZWNlc2l0YW4gMTAgeWFyZGFzIHBhcmEgYW5vdGFyIHkgdGllbmVuIGN1YXRybyBvcG9ydHVuaWRhZGVzIHBhcmEgbG9ncmFybG8sIGxhIGFub3RhY2lvzIFuIGVzIHNlZ3VyYSwgc2llbXByZSB5IGN1YW5kbyBtYW50ZW5nYW4gc3UgZm9ybWEgZGUganVnYXIgcG9yIHRpZXJyYS7igJ0KIyMjIyBSZXNwdWVzdGE6IExhIGFmaXJtYWNpb8yBbiBlcyBpbmNvcnJlY3RhIHBvcnF1ZSBpZ25vcmEgY29tcGxldGFtZW50ZSBsYSB2YXJpYWJpbGlkYWQgZW4geWFyZGFzIGdhbmFkYXMgcG9yIGNhcnJlcmEuCgojIyMjICoqMy04NioqCgojIyMjIyBBIGNvbnRpbnVhY2nDs24sIHNlIHRpZW5lbiAzIHBhcnRlcyBkZWwgcHJlc3VwdWVzdG8gZGUgZGVmZW5zYSBkZSB1biBhw7FvLCBhIGNhZGEgdW5hIGRlIMOpc3RhcyBzZSBsZSBhc2lnbsOzLCBwb3IgcGFydGUgZGVsIENvbmdyZXNvIG1leGljYW5vLCBsYSBtaXNtYSBjYW50aWRhZCBkZSBmaW5hbmNpYW1pZW50bzoKIyMjIyMgKglTYWxhcmlvIGRlIG9maWNpYWxlcyAodG90YWwpLgojIyMjIyAqCU1hbnRlbmltaWVudG8gZGUgbGEgZmxvdGEgYcOpcmVhLgojIyMjIyAqCUFkcXVpc2ljaW9uZXMgZGUgYWxpbWVudG9zICh0b3RhbCkuCgojIyMjIyBUb21hbmRvIGVuIGN1ZW50YSBsYSBkaXN0cmlidWNpw7NuIGRlIHBvc2libGVzIHJlc3VsdGFkb3MgcGFyYSBsb3MgZ2FzdG9zIHJlYWxlcyBlbiBjYWRhIHVuYSBkZSDDqXN0YXMgw6FyZWFzLCBoYWdhIGNvcnJlc3BvbmRlciBjYWRhIHNlY2Npw7NuIGEgdW5hIGRlIGxhcyBjdXJ2YXMgZGUgbGEgc2lndWllbnRlIGZpZ3VyYSwgZnVuZGFtZW50ZSBzdSByZXNwdWVzdGEuIAoKIVtdKC9Vc2Vycy9nYWJyaWVsbWVkaW5hL0Rvd25sb2Fkcy9DYXB0dXJhIGRlIFBhbnRhbGxhIDIwMjMtMDUtMTggYSBsYShzKSA5LjU4LjE0LnBuZykKCiMjIyMgUmVzcHVlc3RhOiBTYWxhcmlvcyBkZSBmdW5jaW9uYXJpb3M6IEE7IG1hbnRlbmltaWVudG8gZGUgZmxvdGE6IEM7IGFkcXVpc2ljaW9uZXMgZGUgYWxpbWVudG9zOiBCLgoKIyMjIyAqKjMtOTIqKgojIyMjIyBFbCAzMCBkZSBqdW5pbyBkZSAxOTkyLCBsYSBjYXBpdGFsaXphY2nDs24gZGUgbnVldmUgbWVyY2Fkb3MgZGUgdmFsb3JlcyBkZWwgUGFjw61maWNvIHkgQXNpYSBmdWU6CgpQYcOtcyAgICAgICAgICB8IENhcGl0YWxpemFjacOzbiAoRW4gbWlsZXMgbWlsbG9uZXMgZGUgZMOzbGFyZXMpCi0tLS0tLS0tLS0tLS0gfCAtLS0tLS0tLS0tLS0tCkZpbGlwaW5hcyAgIAl8MTcKSW5kb25lc2lhICAgCXwyMQpUYWlsYW5kaWEJICAgIHw0NApTaW5nYXB1cgkgICAgfDUwCk1hbGFzaWEJICAgICAgfCA3OQpDb3JlYSBkZWwgU3VyIHwgODYKVGFpd2FuCSAgICAgIHwgMTQwCkhvbmcgS29uZwkgICAgfCAxNzgKQXVzdHJhbGlhCSAgICB8IDIwMwoKIyMjIyMgYSkJRW5jdWVudHJlIGxhIG1lZGlhIGFyaXRtw6l0aWNhIGRlIGxvcyBkYXRvcy4KCmBgYHtyfQpjYXBpdGFsaXphY2lvbiA9IGMoMTcsMjEsNDQsNTAsNzksODYsMTQwLDE3OCwyMDMpCm1lZGlhID0gbWVhbihjYXBpdGFsaXphY2lvbikKbWVkaWEKYGBgCgojIyMjIyBiKQlFbmN1ZW50cmUgbGEgbWVkaWFuYSBkZSBsb3MgZGF0b3MuCmBgYHtyfQptZWRpYW5hID0gbWVkaWFuKGNhcGl0YWxpemFjaW9uKQptZWRpYW5hCmBgYAoKIyMjIyMgYykJIEVuY3VlbnRyZSBsYSBtb2RhIGRlIGxvcyBkYXRvcy4gCmBgYHtyfQojTm8gaGF5IGZ1bmNpw7NuIGRlIG1vZGEgeSBsb3MgZGF0b3MgYWwgc2VyIHZhcmlhYmxlcyBtdXkgZGlzdGludGFzLCBubyB0aWVuZW4gbW9kYS4KYGBgCgojIyMjIyBkKQkgwr9DdcOhbCBlcyBsYSBtZWpvciBtZWRpZGEgZGUgbGEgdGVuZGVuY2lhIGNlbnRyYWwgZGVsIGNvbmp1bnRvIGRlIGRhdG9zPwoKYGBge3J9Cmhpc3RvZ3JhbWE9IGhpc3QoY2FwaXRhbGl6YWNpb24pCgpgYGAKIVtdKC9Vc2Vycy9nYWJyaWVsbWVkaW5hL0Rvd25sb2Fkcy9HcmFmTWVkQ29ycjIucG5nKQoKI0xhIG1lam9yIG1lZGlkYSBlcyBsYSBtZWRpYW5hIHlhIHF1ZSBoYXkgdW4gc2VzZ28gcG9zaXRpdm8gY29tbyBzZSBhcHJlY2lhIGVuIGVsIGhpc3RvZ3JhbWEgZGUgbG9zIGRhdG9zCgojIyMjIyBlKQkgRW5jdWVudHJlIGxhIGRlc3ZpYWNpw7NuIGVzdMOhbmRhciBkZSBsb3MgZGF0b3MuIChMYSBwb2JsYWNpw7NuIGNvbXBsZXRhIGVzdMOhIGluY2x1aWRhIGVuIGVsbG9zLikKCmBgYHtyfQpzdW1hPWNhcGl0YWxpemFjaW9uLW1lZGlhCnN1bWEyPXN1bWEqc3VtYQpzdW1hMz1zdW0oc3VtYTIpCnN1bWE0PXN1bWEzLzkKc3VtYTU9c3FydChzdW1hNCkKc3VtYTUKCmBgYAoKIyMjIyAqKjMtMTAwKioKCiMjIyMjIE1hdHRoZXdzLCBZb3VuZyB5IEFzb2NpYWRvcywgdW5hIGFnZW5jaWEgZGUgY29uc3VsdG9yw61hcyBkZSBDaGFwZWxsIEhpbGwsIHRpZW5lIGxvcyBzaWd1aWVudGVzIHJlZ2lzdHJvcyBxdWUgaW5kaWNhbiBlbCBuw7ptZXJvIGRlIGTDrWFzIHF1ZSBjYWRhIHVubyBkZSBzdXMgMTAgY29uc3VsdG9yZXMgZGUgcGxhbnRhIGNvYnLDsyBlbCDDumx0aW1vIGHDsW86IAojIyMjIyAyMTIgICAgICAyMjAgICAgICAyMzAgICAgICAyMTAgICAgICAyMjggICAgICAyMjkgICAgICAyMzEgICAgICAyMTkgICAgICAyMjEgICAgICAyMjIKIyMjIyMgYSkJU2luIGNhbGN1bGFyIGVsIHZhbG9yIGRlIGxhcyBtZWRpZGFzLCDCv2N1w6FsIGRlIGVsbGFzIGNyZWUgdXN0ZWQgcXVlIGxlIGRhcsOtYSB1bmEgbWF5b3IgaW5mb3JtYWNpw7NuIGFjZXJjYSBkZSBlc3RhIGRpc3RyaWJ1Y2nDs246IGVsIHJhbmdvIChhbGNhbmNlKSBvIGxhIGRlc3ZpYWNpw7NuIGVzdMOhbmRhcj8KIyMjIyAqKlJlc3B1ZXN0YSoqOiBEZXN2aWFjaW/MgW4gZXN0YcyBbmRhci4KIyMjIyMgYikJQ2FsY3VsZTogUmFuZ28sIFZhcmlhbnphIHkgRGVzdmlhY2nDs24gRXN0w6FuZGFyLgojIyMjICoqUmVzcHVlc3RhKioKYGBge3J9CmRpYXM9IGMoMjEyLCAyMjAsIDIzMCwgMjEwLCAyMjgsIDIyOSwgMjMxLCAyMTksIDIyMSwgMjIyKQpyYW5nb19kaWFzPSBtYXgoZGlhcyktbWluKGRpYXMpCm1lZGlhX2RpYXM9IG1lYW4oZGlhcykgCiNEZXN2aWFjacOzbiBlc3RhbmRhcgpkaWFzMj0gZGlhcz0gZGlhcyAtIG1lZGlhX2RpYXMKZGlhczM9IGRpYXMyKmRpYXMyCmRpYXM0PXN1bShkaWFzMykKZGlhczU9IGRpYXM0LzEwCmRpYXM2X2Rlc3Y9c3FydChkaWFzNSkKCiNSZXNwdWVzdGFzCnJhbmdvX2RpYXMKbWVkaWFfZGlhcwojRGVzdmlhY2nDs24gZXN0YW5kYXIKZGlhczZfZGVzdgpgYGAKCiMjIyMjIGMpCSBUb21hbmRvIGVuIGN1ZW50YSBsYSBkaWZpY3VsdGFkIHkgZWwgdGllbXBvIHBhcmEgY2FsY3VsYXIgY2FkYSBtZWRpZGEgcXVlIHJldmlzw7MgZW4gZWwgaW5jaXNvIGEpLCDCv2N1w6FsIHN1Z2VyaXLDrWEgY29tbyBsYSBtZWpvcj8KIyMjIyAqKlJlc3B1ZXN0YToqKiBMYSBkZXN2aWFjacOzbiBlc3RhbmRhciAKIyMjIyMgZCkJwr9RdcOpIGhhcsOtYSBxdWUgdXN0ZWQgY2FtYmlhcmEgc3Ugb3BpbmnDs24gYWwgcmVzcGVjdG8/CiMjIyMgKipSZXNwdWVzdGE6KiogTmFkYQoKIyMjIyAqKjMtMTA2KioKCiMjIyMjIEFsbGlzb24gQmFyZXR0IHJlYWxpemEgYW7DoWxpc2lzIGVzdGFkw61zdGljb3MgcGFyYSB1biBlcXVpcG8gZGUgY2FycmVyYXMgYXV0b21vdmlsw61zdGljYXMuIEEgY29udGludWFjacOzbiwgc2UgcHJlc2VudGFuIGxhcyBjaWZyYXMgZW4ga2lsw7NtZXRyb3MgcG9yIGxpdHJvIGRlbCBnYXN0byBkZSBjb21idXN0aWJsZSBkZSBzdXMgYXV0b23Ds3ZpbGVzIGVuIGxhcyBjYXJyZXJhcyByZWNpZW50ZXM6CiMjIyMjIDQuNzcJNi4xMQk2LjExCTUuMDUJNS45OQk0LjkxCTUuMjcJNi4wMSA1Ljc1CTQuODkJNi4wNQk1LjIyCTYuMDIJNS4yNAk2LjExCTUuMDIKIyMjIyMgYSkJQ2FsY3VsZSBsYSBtZWRpYW5hIGRlbCBjb25zdW1vIGRlIGNvbWJ1c3RpYmxlLgoKIyMjIyAqKlJlc3B1ZXN0YToqKgpgYGB7cn0Ka20gPSBjKDQuNzcsIDYuMTEsCTYuMTEsCTUuMDUsCTUuOTksCTQuOTEsCTUuMjcsCTYuMDEsIDUuNzUsCTQuODksCTYuMDUsCTUuMjIsCTYuMDIsCTUuMjQsCTYuMTEsCTUuMDIpCm1lZGlhbmFfa20gPSBtZWRpYW4oa20pCm1lZGlhbmFfa20KYGBgCgoKIyMjIyMgYikJQ2FsY3VsZSBsYSBtZWRpYSBkZWwgbWlzbW8gY29uc3Vtby4KCiMjIyMgKipSZXNwdWVzdGE6KioKYGBge3J9Cm1lZGlhX2ttID0gbWVhbihrbSkKbWVkaWFfa20KYGBgCgojIyMjIyBjKQlBZ3J1cGUgbG9zIGRhdG9zIGVuIDUgY2xhc2VzIGRlIGlndWFsIHRhbWHDsW8uIMK/Q3XDoWwgZXMgZWwgaW50ZXJ2YWxvIGRlbCB2YWxvciBkZSBjb25zdW1vIGRlIGNvbWJ1c3RpYmxlIHBhcmEgbGEgY2xhc2UgbW9kYWw/CgojIyMjICoqUmVzcHVlc3RhOioqCmBgYHtyfQpjbGFzZXMgPC0gY3V0KGttLCBicmVha3M9NSkKY2xhc2VzMj0gdGFibGUoY2xhc2VzKQpjbGFzZXMyCmBgYAoKCiMjIyMjIGQpCcK/Q3XDoWwgZGUgbGFzIDMgbWVkaWRhcyBkZSB0ZW5kZW5jaWEgY2VudHJhbCBlcyBsYSBxdWUgbWVqb3IgcHVlZGUgc2VydmlybGUgYSBBbGxpc29uIGN1YW5kbyBoYWdhIHVuIHBlZGlkbyBkZSBjb21idXN0aWJsZT8gRXhwbGlxdWUgc3UgcmVzcHVlc3RhLgoKIyMjIyAqKlJlc3B1ZXN0YToqKgoKYGBge3J9Cmhpc3RvZ3JhbWFfa209IGhpc3Qoa20pCmhpc3RvZ3JhbWFfa20KI0RlcGVuZGUgCmBgYAoKCiMjIyMjIGUpCcK/Q3XDoWwgZXMgZWwgcmFuZ28/CiMjIyMgKipSZXNwdWVzdGE6KioKYGBge3J9CnJhbmdvPSBtYXgoa20pIC0gbWluKGttKQpyYW5nbwpgYGAKCiMjIyMgKio4LTY0KioKCgojIyMjIyBJbnNwZWN0b3JlcyBkZWwgZ29iaWVybm8sIGFsIGludmVzdGlnYXIgbG9zIGNhcmdvcyBsZXZhbnRhZG9zIGNvbnRyYSB1bmEgZW1ib3RlbGxhZG9yYSBkZSBiZWJpZGFzIG5vIGFsY29ow7NsaWNhcywgZGUgVGV4YXMsIHF1ZSBubyBsbGVuYWJhIGFkZWN1YWRhbWVudGUgc3VzIHByb2R1Y3RvcywgaGFuIG11ZXN0cmVhZG8gMjAwIGJvdGVsbGFzIHkgZW5jb250cmFyb24gcXVlIGVsIHByb21lZGlvIGRlIGxsZW5hZG8gZXMgZGUgMzEuNyBvbnphcyBsw61xdWlkYXMgKGFwcm94aW1hZGFtZW50ZSA5MzAgbWwuKSAgU2UgYW51bmNpYSBxdWUgbGFzIGJvdGVsbGFzIGNvbnRpZW5lbiAzMiBvbnphcyBsw61xdWlkYXMgKDk0Ni4zMyBtbC4pLiBTZSBzYWJlIHF1ZSBsYSBkZXN2aWFjacOzbiBlc3TDoW5kYXIgZGUgbGEgcG9ibGFjacOzbiBlcyBkZSAxLjUgb256YXMgbMOtcXVpZGFzICg0NC4zNiBtbC4pLiDCv0RlYmVyw61hbiBjb25jbHVpciBsb3MgaW5zcGVjdG9yZXMsIGFsIG5pdmVsIGRlIHNpZ25pZmljYW5jaWEgZGUgMiUsIHF1ZSBsYXMgYm90ZWxsYXMgZXN0w6FuIHNpZW5kbyBsbGVuYWRhcyBjb24gbWVub3MgY29udGVuaWRvLCB5IGVtaXRpciB1bmEgc2FuY2nDs24/CgojIyMjICoqUmVzcHVlc3RhOioqCgohW10oL1VzZXJzL2dhYnJpZWxtZWRpbmEvRG93bmxvYWRzL0NhcHR1cmEgZGUgUGFudGFsbGEgMjAyMy0wNS0xOSBhIGxhKHMpIDEwLjA3LjM5LnBuZykKCmBgYHtyfQojIyMjIyBQYXNvIDE6IFBsYW50ZWFyIGxhIGhpcMOzdGVzaXMKI2gwOiB4ID0gMzIKI2gxOiB4IOKJoCAzMS43CgojIyMjIyBQYXNvIDI6IE5pdmVsIGRlIHNpZ25pZmljYW5jaWEKIyA5OSUgKGRlIGNhZGEgbGFkbykgPSAyLjcxOCAoeikgCiNBbGZhID0gMC4wMiAvMgoKIyMjIyMgUGFzbyAzOiBab25hIGRlIGFjZXB0YWNpw7NuL3JlY2hhem8KCiMgRW4gZWwgcmFuZ28gZGUgLTIuMzMgYSAyLjMzCgojIyMjIyBQYXNvIDQ6IEZ1bmNpw7NuIHBpdm90YWwKCno9ICgzMS43LTMyKSAvICgxLjUvc3FydCgyMDApKQp6CiMjIyMjIFBhc28gNTogQ29uY2x1c2nDs24KIyBTZSByZWNoYXphIEgwIHBvcnF1ZSBlc3TDoSBhZnVlcmEgZGVsIHJhbmdvLCBwb3IgbG8gdGFudG8gc2UgY29uY2x1eWUgY29uIHF1ZSBsYXMgYm90ZWxsYXMgc2UgbGxlbmFuIGNvbiBtZW5vcyBjb250ZW5pZG8uIAoKYGBgCgoKCgoK