6 votos

Generar la cola de la distribución de una muestra determinada en R

Tengo una muestra de mediciones a partir de la vida real del dispositivo que pierde todas las medidas que están a menos de un cierto umbral (dado dispositivo no es lo suficientemente precisa).

De la teoría y también mediciones más precisas de los dispositivos sé que distribuyen los datos muy de cerca a la distribución de Gumbel.

He encontrado los parámetros de una distribución de Gumbel para una muestra determinada, con el MLE y ahora quiero generar "perfecto" falta de valores y agregar a la muestra.

No estoy seguro de cuál es la manera más fácil de hacer que en R. Tengo mi muestra y tamaño de la muestra, obviamente. Quiero crear una copia de la misma y, a continuación, agregue algunos de los valores de uso teórico PDF.

4voto

AdamSane Puntos 1825

Estoy suponiendo que los valores son truncados por debajo del umbral, $t$ más que censurado por debajo de $t$ (es decir, no se sabe cuántos hay por debajo del umbral).

Deje que el número de puntos que se observa por encima del punto de truncamiento de ser $n_o$.

Un enfoque simple podría ir de la siguiente manera:

  1. Estimación de $p$, la proporción de la distribución por debajo del umbral de la equipada con los parámetros y el punto de truncamiento (esto ignora la incertidumbre en las estimaciones de los parámetros; sin embargo, si la muestra es grande que no puede ser un terrible aproximación)

  2. Simular $N_u$, el número de valores perdidos de la binomial negativa $NB(n_o,p)$, obteniendo el valor específico de la $n_u$.

  3. Simular $n_u$ valores de la truncada Gumbel en (0,$t$); por ejemplo, usted podría considerar la posibilidad de alguna forma de aceptar-rechazar.

i-Ciencias.com

I-Ciencias es una comunidad de estudiantes y amantes de la ciencia en la que puedes resolver tus problemas y dudas.
Puedes consultar las preguntas de otros usuarios, hacer tus propias preguntas o resolver las de los demás.

Powered by:

X