Cinturones de confianza por Monte Carlo

Parte A
1. Generar n números aleatorios xi con distribución exponencial con un
valor dado de τ y calcular un estadístico t igual al promedio de los xi.
2. Repetir el punto anterior, un numero NExperimentos de veces y guardar
el valor del estadístico t obtenido cada vez en un histograma.
3. Para ese histograma, encontrar un valor de tmin y tmax tal que entre
ellos se encuentre una fracción CL del total de los eventos.
4. Repetir los tres puntos anteriores para 100 valores de tau entre 0 y 10
y graficar  el cinturón de confianza en el plano τ vs t.
Parte B
Utilizar la herramienta generada en la Parte A para otros dos estadísticos:
1. la mediana
2. t=Σi  (xi+xi**2+xi**3+xi**4)

Parte C
Discutir que pasa si como “estadístico” t utilizamos el q de Wilks visto en clase.

Parte D
Utilizar Wilks para calcular el intervalo de confianza para el parámetro τ de la  exponencial con un dado CL.
Parte E
(Opcional) Calcular la cobertura de los intervalos calculados en los ítems anteriores en función del parámetro τ
Print Friendly, PDF & Email