Simulación del Monte CarloHola comunidad, creé una tabla de meses (x) por frecuencia de falla (y). Después analicé la distribución y creé otra columna con aleatorio para la distribución que mejor se ajusta (AIC más pequeños). Pero necesito calcular la probabilidad final de unión de los cinco eventos (EVENTO SUPERIOR) y la simulación (10000 iteraciones) me dan en frecuencia de falla. ¿Cómo obtuve la probabilidad de fallar X meses en el EVENTO SUPERIOR con 10000 iteraciones?
Re: Simulación de MontecarloHola comunidad,
Hice la simulación de Monte Carlo con 10.000 iteraciones, ahora necesito saber cuáles son los parámetros estadísticos de evaluación de la simulación. Cuál es el mejor parámetro estadístico para evaluar la simulación.
Atentamente
mayre
Re: Simulación de MontecarloHola Kevin,
Hice la simulación considerando la probabilidad como la frecuencia de las fallas y el número total de muestras generadas, luego obtuve la probabilidad de fallas. El algoritmo utilizado fue la fórmula de probabilidad para eventos independientes. Consideré mejor la distribución exponencial porque los valores de falla no pueden ser negativos.
¡¡Muchas gracias!!
Re: Simulación de MontecarloHola Mayre! Algunas personas se obsesionan con la frase "Monte Carlo" porque significa algo muy específico para ellas, pero hablemos de simulación por un momento. Para que las simulaciones sean exitosas, se requieren algunas cosas: 1) el generador de números produce valores que pasan las pruebas de aleatoriedad, 2) se generan suficientes muestras para producir resultados precisos, 3) se utilizan técnicas de muestreo adecuadas, 4) el algoritmo utilizado es válido para lo que se está modelando, y 5) simula el fenómeno en cuestión. JMP puede manejar bien los números 1 y 2, pero 3, 4 y 5 dependen del modelador (usted) para especificarlos correctamente. Según mi experiencia, 4 y 5 son los más difíciles y propensos a errores. Le aconsejo que dedique suficiente tiempo a estudiar esos requisitos en relación con su situación. Por ejemplo, podría ser mejor utilizar distribuciones elegidas por teoría que por mejor ajuste. Por cierto, una definición de probabilidad es una frecuencia relativa de ocurrencia. Dado que conoce la frecuencia de las fallas y el número total de muestras generadas, ¿por qué no sería una estimación de la probabilidad de fallas?Re: Simulación de MontecarloPrimero, le agradezco su atención.
Supongo que todos los subsistemas son completamente independientes. No existe ningún subsistema F2, sólo F1, F3, F4, F5 y F6, todos independientes.
Estoy siguiendo los pasos del enlace anterior: "Uso de JMP para evaluar el riesgo en predicciones financieras mediante simulaciones de Monte Carlo
Pero en mi trabajo estoy usando el FTA (Análisis de árbol de fallas) para calcular la probabilidad del EVENTO SUPERIOR. Estoy usando los conjuntos de cortes y la fórmula final fue
mi = F1 +F3+F4+F5+F6
Como los eventos son independientes, se debe calcular la probabilidad. Pero mi resultado en la frecuencia de fallas no es probabilidad.
Espero ayudarte!!
Re: Simulación de MontecarloHola Mayre!
Supongo que si falla un subsistema, falla todo el sistema. ¿Es eso correcto?
¿Son los subsistemas completamente independientes?
¿Cometiste un error en tu fórmula al omitir p(F2)? ¿O esa era tu intención?
Más información hace que sea más fácil ayudar.
Re: Simulación de MontecarloHola Kevin,
Necesito simular la probabilidad de falla de un evento específico (Evento principal de FTA - ANÁLISIS DE ÁRBOL DE FALLAS).
Realicé los siguientes pasos:
Creé la tabla con los datos:
1) X (filas) - meses e Y (columna) la frecuencia que son fechas históricas de 23 meses de observaciones.
ejemplo: 1 1
2 1
3 0 ......
2) Ajustar distribuciones al subconjunto histórico (5 subconjuntos en total)
3) Se creó un nuevo conjunto de datos (con una nueva columna) con extracciones aleatorias de las distribuciones para los cinco subconjuntos.
4) Creó la nueva columna (EVENTO SUPERIOR), pero la columna debe ser la probabilidad de falla del EVENTO SUPERIOR (EVENTO TOPO = PROBABILIDAD TOTAL)
Pregunta: No puedo sumar las frecuencias de fallas. Necesito calcular la probabilidad. La fórmula final es:
P(T) = P(F1)+P(F3)+P(F4)+P(F5)+P(F6) - (P(F1).P(F2) - P(F3).P(F4) - P(F4).P(F5) - P(F5).P(F6)) - P(F1).P(F3).P(F4).P(F5).P(F6)
¿Cómo transformo frecuencia en probabilidad en este caso?
Re: Simulación de MontecarloHola Mayre!
No entiendo tu petición.
JMP y JSL tienen grandes capacidades de simulación.
Si pudiera reformular su problema, sus datos y su enfoque previsto con mucho más detalle, ¿tal vez podríamos ayudarlo?
