5 votos

Mejores prácticas para la estimación de densidad de discretos y continuos de variables aleatorias

Actualmente estoy tratando de estimar la densidad de una distribución conjunta de K unidimensional, vehículos recreativos. Tengo a mi disposición un conjunto de N puntos de muestreo, cada uno de los cuales representa un resultado de la K RVs.

Algunos detalles acerca de mi problema:

  • la RVs son independientes

  • la RVs necesidad no pertenecen a la misma familia de distribuciones

  • la Rv puede todos ser discretos o de todo ser continua, pero no tanto

  • Sé que los límites superior e inferior para cada RV; en el caso discreto, conozco los valores que la RV puede tomar.

Ahora, yo soy la estimación de la densidad de cada uno de RV por separado, utilice la ksdensity función en MATLAB. La independencia de la asunción, a continuación, me permite producir un conjunto de densidad mediante el producto de las densidades individuales. Estoy esperando para mejorar la precisión de mi estimación, esto ya sea por el uso de otro método (que puedo código en MATLAB) o por jugar con las opciones en ksdensity (tales como el núcleo del tipo, el apoyo, la anchura de la densidad de la ventana).

Estoy específicamente con la esperanza de que la gente puede arrojar luz sobre:

  • ¿Cuál es el método a utilizar para el caso discreto vs el caso continuo. En el caso discreto, vale la pena especificar los límites y los valores? En el caso continuo,

  • Si alguna vez hace sentido para olvidar la independencia de la asunción y estimar la distribución conjunta como una distribución conjunta

  • Si alguien sabe de algún simple material de lectura sobre el tema.

11voto

brian buck Puntos 1103

Sólo voy a hablar sobre el caso continuo.

Primero de todo yo no uso de MATLAB (en lugar de R o código de mí lo que necesito), pero he mirado en su enlace :

  • el ksdensity función de los usos normales del kernel que está bien (otra suave kernel que podían probar el kernel de Epanechnikov, otras dos buenas opciones pueden ser biweight kernel o Silverman kernel pero tendrá que poner en práctica por sí mismo).

  • pero esta función utiliza una forma muy extraña elección de anchura, a menos que intente estimar normal RV, si no la manera más fácil de hacer el trabajo es el uso de la validación cruzada para elegir su ancho

  • si su RV son de forma compacta compatible el kernel smoothing va a generar artefactos cerca de los puntos finales del apoyo que usted desea, para corregir (el estimado de RV puede ir sobre el apoyo de la "real")

Si está seguro de que su RV son independientes, se puede estimar cada uno por separado. y hacer el producto (y no, usted NO debe usar cúpulas puesto que su propósito es el modelo de dependencia)

Voy a terminar con dos de las cosas que usted desea, puede leer :

  • una muy buena (pero bastante teórico) de referencia es Tsybakov del libro "Introducción a los no paramétricos de estimación"

  • muy fácil de leer, si se puede acceder a ella de referencia es el "Kernel smoothing" por Chapman & Hall.

2voto

Akira Puntos 1061

Es posible que desee probar el uso de una Cópula. Hay algunas versiones gratuitas de este en la web, el más prometedor de lo que parece ser Andrew Patton Cópula del cuadro de herramientas (tengo no se utiliza este, en los que la mente, pero se ve sobre la derecha).

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