Conférencier | Transcription |
| Notre exposé d'aujourd'hui montre comment |
| utiliser JMP pour faire du bayésien |
| estimation. Voici un aperçu |
| de mon discours. Je vais commencer |
| avec une brève introduction à |
| Méthodes statistiques bayésiennes. |
| Alors je vais passer par |
| quatre exemples différents qui |
| viennent de la fiabilité, |
| mais les méthodes que nous présentons |
| sont vraiment beaucoup plus généraux et |
| peut être appliqué dans d'autres domaines de |
| application. Alors je vais |
| cède la parole à Peng et il est |
| va te montrer à quel point c'est facile |
| faire réellement ces choses dans |
| JMP. Techniquement, fiabilité |
| est une probabilité. Le |
| probabilité d'un système, |
| véhicule, machine ou quoi que ce soit |
| est-ce que c'est intéressant, sera |
| remplir sa fonction prévue |
| sous fonctionnement rencontré |
| conditions pour une période déterminée |
| de temps. Je souligne rencontré |
| ici pour souligner que |
| la fiabilité dépend essentiellement |
| sur l'environnement dans lequel un |
| le produit est utilisé. Condra |
| définition de la fiabilité comme qualité |
| heures supplémentaires. Et de nombreux ingénieurs |
| penser que la fiabilité est |
| éviter les échecs, c'est-à-dire |
| concevoir et fabriquer un produit |
| cela n'échouera pas. Fiabilité |
| est un très quantitatif |
| discipline d'ingénierie, mais |
| nécessite souvent sophistiqué |
| statistique et probabiliste |
| des idées. Au cours des 30 dernières années, |
| il y a eu un virtuel |
| révolution où les méthodes Bayes |
| sont maintenant couramment utilisés et dans |
| de nombreux domaines différents de |
| application. Cette révolution |
| commencé par la redécouverte de |
| Chaînes de Markov Méthodes de Monte Carlo |
| et a été accéléré par |
| améliorations spectaculaires dans |
| la puissance de calcul dont nous disposons |
| aujourd'hui, ainsi que le |
| développement de relativement facile |
| utiliser un logiciel pour implémenter |
| Méthodes Bayes. Dans les années 1990, nous |
| avait des BUGS. Aujourd'hui Stan et d'autres |
| des paquets similaires ont largement |
| remplacé BUGS, mais l'autre |
| ce qui se passe, c'est que nous |
| commence à voir plus bayésien |
| méthodes mises en œuvre dans |
| logiciel commercial. Donc pour |
| exemple, SAS a PROC MCMC. |
| Et maintenant JMP a quelques très |
| des outils puissants qui étaient |
| développé pour la fiabilité, mais |
| comme je l'ai dit, ils peuvent être appliqués |
| dans d'autres domaines également, et |
| il y a une forte motivation pour |
| l'utilisation des méthodes bayésiennes. |
| D'une part, il fournit un |
| des moyens pour combiner avant |
| informations avec des données limitées |
| pouvoir rendre utile |
| inférences. En outre, il existe de nombreux |
| situations, en particulier avec |
| complications à effets aléatoires |
| comme les données de censure, où le maximum |
| la probabilité est difficile à |
| mettre en œuvre, mais où Bayes |
| méthodes sont relativement faciles à |
| mettre en place. Il y a un petit |
| inconvénient de l'utilisation de Bayes |
| méthodes. Il faut réfléchir un peu |
| plus difficile sur certaines choses, |
| en particulier sur |
| paramétrage et comment |
| spécifier les distributions antérieures. |
| Mon premier exemple concerne un |
| cage de roulement de moteur d'avion. |
| Ce sont des données de défaillance sur le terrain |
| où il y avait 1703 avions |
| moteurs qui contenaient ce |
| cage de roulement. Les plus anciens |
| avait 2 220 heures de fonctionnement. Le |
| spécification de la durée de vie de |
| cette cage de roulement était que non |
| plus de 10% des unités |
| échouent après 8 000 heures de fonctionnement. |
| Cependant, 6 unités avaient échoué et |
| cela a soulevé la question | est-ce que nous |
| avez un sérieux problème ici? Fais |
| nous devons repenser ce roulement |
| cage pour répondre à cette fiabilité |
| état? Ceci est un graphique d'événement |
| des données. L'intrigue de l'événement |
| illustre la structure |
| des données, et en particulier, |
| on peut voir les six échecs |
| ici. En plus de cela, nous |
| avoir droit censuré |
| observations, indiquées ici par |
| les flèches. Donc ce sont des unités |
| qui sont toujours en service et |
| ils n'ont pas encore échoué, et |
| la flèche droite indique que |
| tout ce que nous savons c'est que si nous attendons longtemps |
| assez à droite, les unités |
| finira par échouer. Voici un |
| analyse du maximum de vraisemblance |
| ces données, donc la probabilité |
| l'intrigue suggère ici que le |
| La distribution Weibull fournit un |
| bonne description de ces données. |
| Cependant, lorsque nous utilisons le |
| profileur de distribution pour |
| estimer la fraction échouant à |
| 8000 heures, on voit que le |
| l'intervalle de confiance est énorme, |
| compris entre environ 3% de tous les |
| jusqu'à 100%. Ce n'est pas très |
| utile. Donc méthodes de vraisemblance |
| fonctionne comme ça. Nous spécifions le |
| modèle et les données. |
| Cela définit la probabilité et |
| alors nous utilisons la probabilité de |
| faire des inférences. Bayésien |
| les méthodes sont similaires, sauf que nous |
| avoir également des informations préalables |
| spécifié. Théorème de Bayes |
| combine la vraisemblance avec le |
| informations préalables, fournissant un |
| distribution postérieure, puis |
| nous utilisons le postérieur |
| distribution pour faire des inférences. |
| Voici l'analyse Bayes du |
| cage de roulement. Les priors sont |
| spécifié ici pour le B10 |
| ou moment auquel 10% échoueraient. |
| Nous avons un intervalle très large |
| ici. La gamme, effectivement 1000 |
| heures jusqu'à 50 000 heures. |
| Tout le monde conviendrait que B10 |
| est quelque part dans cette gamme. Pour |
| le paramètre de forme de Weibull, |
| cependant, nous allons utiliser un |
| distribution préalable informative |
| basé sur les ingénieurs |
| connaissance de l'échec |
| mécanisme et leur vaste |
| expérience précédente avec ça |
| mécanisme. Ils peuvent dire avec |
| peu de doute que le Weibull |
| le paramètre de forme doit être |
| entre 1,5 et 3, et |
| voici où nous spécifions que |
| informations. Donc au lieu de |
| spécifiant des informations pour |
| paramètres traditionnels de Weibull, |
| nous avons reparameterized, où maintenant |
| le B10 est l'un des |
| paramètres, et voici les |
| plage spécifiée. Et puis nous |
| avoir l'information préalable pour |
| le paramètre de forme de Weibull |
| spécifié ici. Et puis JMP |
| générera des échantillons à partir du |
| articulation postérieure, conduisant à |
| ces estimations de paramètres et |
| intervalles de confiance indiqués ici. |
| Voici une représentation graphique de |
| l'analyse de Bayes. Le noir |
| points voici un échantillon de |
| la distribution précédente, donc encore |
| très large pour le quantile .1 |
| et quelque peu contraint pour le |
| Paramètre de forme Weibull beta. Sur |
| le droit ici, nous avons le |
| articulation postérieure, qui en effet |
| est l'endroit où la vraisemblance contourne |
| et l'échantillon précédent |
| chevaucher. Et puis ces tirages |
| de l'articulation postérieure sont |
| utilisé pour calculer des estimations et |
| Intervalles crédibles bayésiens. Alors |
| voici le même profileur que nous |
| vu précédemment où le |
| l'intervalle de confiance n'était pas |
| utile. Après avoir apporté le |
| informations sur le Weibull |
| paramètre de forme, maintenant nous pouvons voir |
| que l'intervalle de confiance |
| varie entre 12% et 83%, |
| illustrant clairement que nous |
| ont manqué l’objectif de 10%. |
| Alors qu'avons-nous appris ici? |
| Avec un petit nombre d'échecs, |
| il n'y a pas beaucoup d'informations |
| sur la fiabilité. Mais les ingénieurs |
| ont souvent des informations qui peuvent |
| être utilisé, et en utilisant ce |
| information, nous pouvons nous améliorer |
| précision et plus utile |
| inférences. Et méthodes bayésiennes |
| fournir une méthode formelle pour |
| combinant ces informations antérieures |
| avec nos données limitées. Voici |
| un autre exemple. Le moteur de la fusée est |
| l'un des cinq composants critiques |
| dans un missile. Dans ce particulier |
| application, il y avait |
| environ 20 000 missiles |
| dans l'inventaire. Au fil du temps, 1940 |
| de ces missiles avaient été tirés |
| et ils ont tous travaillé, sauf en |
| trois cas, où il y avait |
| Échec catastrophique. Et ceux-ci |
| étaient des missiles plus anciens, et donc |
| il y avait une certaine inquiétude que |
| il peut y avoir un échec d'usure |
| mécanisme qui mettrait en |
| mettre en péril les quelque 20 000 |
| missiles restant en inventaire. |
| On pensait que les échecs étaient |
| en raison du cyclage thermique, mais |
| il n'y avait aucune information sur |
| le nombre de cycles thermiques. Nous |
| ont seulement l'âge du |
| missile quand il a été tiré. |
| C'est un substitut utile, mais |
| l'effet de l'utilisation d'un substitut |
| comme ça c'est que tu en as plus |
| la variabilité de vos données. À présent |
| dans ce cas, il n'y avait pas |
| temps de défaillance directement observés. |
| Quand une fusée est appelée à |
| fonctionner et il fonctionne |
| avec succès, tout ce que nous savons est |
| qu'ils n'avaient pas échoué à la |
| l'âge de ces unités lorsqu'elles |
| ont été invités à opérer. Et pour |
| les unités qui ont échoué |
| de façon catastrophique, encore une fois, nous |
| je ne sais pas à quelle heure ces |
| les unités ont échoué. À un moment donné |
| avant qu'ils ne soient appelés à |
| fonctionner. Ils avaient échoué, donc tous |
| nous savons que l'échec était |
| avant l'âge auquel il était |
| mis à la porte. Donc, comme je l'ai dit, il y avait |
| crainte qu'il y ait une usure |
| mécanisme de défaillance qui se déclenche |
| ici qui mettrait en |
| mettre en péril le montant restant |
| vie pour les unités dans le |
| stocker. Alors voici le tableau |
| des données. Ici nous avons le |
| unités qui ont fonctionné |
| avec succès, et donc ce sont |
| observations censurées à droite, mais |
| ces observations ici |
| sont ceux qui ont échoué et |
| J'ai dit, à un niveau relativement plus élevé |
| âge. Ceci est le tracé des événements de |
| les données et encore une fois, nous pouvons voir |
| les bonnes observations censurées |
| ici avec la flèche pointant vers |
| la droite, et nous pouvons voir le |
| a laissé des observations censurées avec |
| la flèche pointant vers la gauche |
| indiquant la région de |
| incertitude. Mais même avec ceux |
| données que nous pouvons toujours adapter à un Weibull |
| Distribution. Et voici le |
| diagramme de probabilité montrant le |
| estimation du maximum de vraisemblance et |
| bandes de confiance. En savoir plus |
| informations du maximum |
| analyse de vraisemblance. Et ici nous |
| avoir l'estimation de la fraction |
| échec à 20 ans, ce qui était |
| la durée de vie de ces fusées |
| moteurs, et encore l'intervalle |
| est énorme, compris entre 3% et |
| 100%. Encore une fois, pas très utile. |
| Mais les ingénieurs, |
| savoir ce qu'est l'échec |
| mécanisme était, encore une fois |
| informations sur le Weibull |
| paramètre de forme. Le maximum |
| l'estimation de la probabilité était |
| extrêmement grand et le |
| les ingénieurs étaient à peu près sûrs que |
| c'était faux, surtout avec |
| la variabilité supplémentaire du |
| des données qui auraient tendance à conduire |
| le paramètre de forme Weibull à un |
| valeur inférieure. Comme je vous l'ai montré sur |
| la diapositive précédente, confiance |
| intervalle pour la fraction échouant à |
| 20 ans, c'était énorme. Donc encore une fois, |
| nous allons spécifier un a priori |
| distribution, puis utilisez cela |
| dans une analyse Bayes. Encore une fois, le |
| avant pour B10, l'heure à laquelle |
| 10% échouera, est choisi pour être |
| extrêmement large. Nous n'avons pas vraiment |
| veux assumer quoi que ce soit là-bas, |
| et tout le monde conviendrait que |
| cette quantité est quelque part |
| entre cinq ans et 400 |
| années. Mais pour la forme Weibull |
| paramètre, nous allons supposer |
| que c'est entre un et cinq. |
| Encore une fois, nous savons que c'est plus grand que |
| un parce que c'est une usure |
| mécanisme de défaillance, et le |
| les ingénieurs étaient sûrs qu'il |
| ne ressemblait en rien au nombre |
| 8 que nous avions vu dans le |
| estimation maximale probable. |
| Et en effet, cinq est aussi un très |
| grand paramètre de forme Weibull. |
| Une fois de plus, JMP est sollicité |
| pour générer des tirages à partir du |
| distribution postérieure. Et voici |
| intrigue similaire à celles que nous |
| scie dans l'exemple de la cage de roulement. |
| Les points noirs ici, encore une fois, sont |
| un échantillon du prieur |
| Distribution. Encore une fois très large dans |
| termes du B10, mais un peu |
| contraint pour la version bêta, de sorte que le |
| bêta est un préalable informatif |
| Distribution. Et encore une fois le |
| les courbes de niveau représentent les |
| informations dans les données limitées. |
| Dans notre postérieur, une fois de plus, est |
| où nous obtenons un chevauchement entre les |
| avant et la probabilité et nous |
| peut le voir ici. Donc encore une fois |
| nous avons une comparaison entre le |
| intervalle de vraisemblance maximale, |
| qui est extrêmement large, et le |
| intervalle que nous obtenons pour le |
| même quantité en utilisant les Bayes |
| inférence, qui a incorporé le |
| informations préalables sur le Weibull |
| paramètre de forme. Et maintenant le |
| l'intervalle est compris entre .002 et |
| 0,98 ou environ 0,1, environ 10% |
| échouer, donc cela pourrait être |
| acceptable. Certaines des choses |
| que nous avons appris ici, même |
| bien qu'il n'y ait pas eu de réel |
| temps d'échec, nous pouvons encore obtenir |
| informations de fiabilité provenant du |
| données, mais avec très peu de pannes |
| il n'y a pas beaucoup d'informations |
| là. Mais nous pouvons utiliser le |
| les connaissances de l'ingénieur sur le |
| Paramètre de forme Weibull à |
| compléter les données pour obtenir |
| inférences utiles et JMP fait |
| c'est vraiment facile à faire. Mon dernier |
| deux exemples concernent |
| tests accélérés. Accéléré |
| le test est un |
| technique pour obtenir des informations |
| sur la fiabilité des composants |
| rapidement lors de la conception d'un |
| produit. L'idée de base est de |
| unités de test à des niveaux élevés de |
| des variables comme la température ou |
| tension pour faire échouer les choses |
| rapidement et ensuite utiliser un modèle |
| pour extrapoler vers le bas |
| conditions d'utilisation. L'extrapolation est |
| toujours dangereux et nous devons |
| garde cela à l'esprit. C'est le |
| raison pour laquelle nous aimerions avoir notre |
| modèle être motivé physiquement. |
| Voici donc un exemple de |
| test de durée de vie accéléré |
| sur un laser. Les unités ont été testées à |
| 40, 60 et 80 degrés C, mais le |
| la condition d'utilisation était de 10 degrés C. |
| C'est la température nominale |
| au fond de l'Atlantique |
| Océan, là où se trouvaient ces lasers |
| va être utilisé dans un nouveau |
| système de télécommunications. Le |
| le test a duré 5000 heures, un peu |
| peu plus de six mois. Le |
| les ingénieurs voulaient estimer |
| fraction échouant à environ 30000 |
| les heures. Cela fait environ 3,5 ans |
| et encore, à 10 degrés C. |
| Voici les résultats de la |
| analysi. Pour |
| tester et construire de manière appropriée |
| modèle, JMP utilise ces trois |
| différentes analyses. La première |
| on ajuste des distributions séparées |
| à chaque niveau de température. |
| Le modèle suivant fait de même |
| chose, sauf qu'elle contraint |
| le paramètre de forme Sigma doit être |
| la même chose à tous les niveaux de |
| Température. C'est analogue |
| à l'hypothèse de Sigma constante |
| que nous fabriquons généralement |
| analyse de régression. Puis |
| enfin, nous ajustons la régression |
| modèle, qui en fait, est un |
| régression linéaire simple |
| connexion à vie à |
| Température. Et pour compléter |
| cette visualisation de ces |
| trois modèles, JMP fait |
| tests de rapport de vraisemblance à tester |
| s'il y a des preuves que |
| le Sigma pourrait dépendre de |
| température puis pour tester |
| s'il y a des preuves de manque |
| d'ajustement dans le modèle de régression. |
| Et à partir des grandes valeurs P |
| ici, on voit qu'il n'y a pas |
| preuves contre ce modèle. |
| Une autre façon de tracer les résultats |
| de monter ce modèle |
| consiste à tracer la durée de vie par rapport à |
| température sur des échelles spéciales. UNE |
| échelle rythmique logarithmique pour des heures |
| opération dans ce que l'on appelle une |
| Échelle d'Arrhenius pour la température. |
| Correspondant à l'Arrhenius |
| modèle, qui décrit comment |
| la température affecte la réaction |
| taux, et donc à vie. Et |
| ce sont les résultats de la |
| Estimation de vraisemblance maximale |
| pour notre modèle. Le JMP |
| le profileur de distribution nous donne |
| une estimation de la fraction |
| échec à 30 000 heures. |
| Et nous pouvons voir que cela varie entre |
| .002 et environ .12, soit 12% |
| échouer. Les ingénieurs de |
| des applications comme celle-ci, cependant, |
| ont souvent des informations sur |
| ce qu'on appelle l'efficace |
| énergie d'activation de la panne |
| mécanisme, et cela correspond |
| à la pente de la régression |
| ligne dans le modèle Arrhenius. Alors |
| nous avons fait une analyse Bayes et en |
| cette analyse, nous avons fait un |
| hypothèse sur l'efficacité |
| énergie d'activation. Et c'est |
| va fournir plus de précision |
| pour nous. Donc ce que nous avons ici est |
| un nuage de points matriciel du |
| distribution postérieure commune |
| après avoir spécifié au préalable |
| distributions pour |
| paramètres, faiblement informatifs |
| pour le quantile .1 à 40 |
| degrés C. Encore une fois, tout le monde |
| serait d'accord que ce nombre est |
| quelque part entre 100 et |
| 32 000. Aussi faiblement informatif |
| pour la forme log-normale |
| paramètre. Encore une fois, tout le monde |
| serait d'accord que ce nombre est |
| quelque part entre 0,05 |
| et 20. Mais pour la pente de la |
| ligne de régression, nous avons un |
| informatif avant de télémétrie |
| entre 0,6 et 0,8, basé sur |
| expérience antérieure avec le |
| mécanisme de défaillance. Et cela mène |
| à cette comparaison, où maintenant |
| sur le côté droit ici, le |
| intervalle pour la fraction échouant à |
| 30000 heures, c'est beaucoup plus restreint |
| qu'il ne l'était avec le maximum |
| estimation de la vraisemblance. Dans |
| en particulier, la limite supérieure maintenant |
| est seulement d'environ 4% par rapport à |
| 12% pour le maximum de vraisemblance |
| estimations. Donc les leçons apprises. |
| Les tests accélérés fournissent |
| informations de fiabilité rapidement, |
| et les ingénieurs ont souvent |
| des informations sur l'effet de |
| énergie d'activation. Et cela peut |
| être utilisé pour améliorer |
| précision ou pour réduire les coûts en |
| pas besoin d'en tester autant |
| unités. Et encore une fois, bayésien |
| les méthodes fournissent un |
| méthode pour combiner les ingénieurs |
| connaissances avec les données limitées. |
| Mon dernier exemple concerne un |
| test de durée de vie accéléré de |
| un dispositif à circuit intégré. Unités |
| ont été testés à haute température |
| et les données résultantes étaient |
| intervalle censuré. C'est |
| parce que des échecs ont été découverts |
| seulement lors des inspections |
| ont été menées périodiquement. Dans |
| ce test, cependant, il y avait |
| seulement des échecs aux deux hauts |
| niveaux de température. Le but |
| du test consistait à estimer le |
| .01 quantile à 100 degrés C. |
| Ceci est un tableau des données |
| où nous pouvons voir les échecs à |
| 250 et 300 degrés C. |
| Et pas d'échecs, d'accord |
| censuré aux trois inférieurs |
| niveaux de température. Maintenant, quand |
| nous avons fait le maximum de vraisemblance |
| estimation, dans ce cas, nous avons vu |
| preuve solide que le Weibull |
| le paramètre de forme dépendait de |
| Température. Donc, la valeur P est |
| environ 0,03. Cela s'avère |
| être une preuve contre le |
| Modèle d'Arrhenius, et c'est |
| car le modèle d'Arrhenius devrait |
| seulement le temps d'échelle. Mais si tu |
| changer le paramètre de forme en |
| augmentation de la température, vous êtes |
| faire plus que réduire le temps. |
| Et c'est donc un problème, et ça |
| a suggéré qu'à 300 degrés C, |
| il y a eu un échec différent |
| mécanisme. Et en effet, quand le |
| les ingénieurs ont suivi et |
| déterminé la cause de l'échec |
| des unités à 250 et 300, |
| ils ont vu qu'il y avait un |
| mécanisme différent à 300. Quoi |
| cela signifiait que nous devions |
| jetez ces données. Et alors |
| faisons-nous alors? Maintenant nous avons seulement |
| eu des échecs à 250 degrés C |
| et JMP ne fait pas très bien |
| avec ça. C'est étonnamment, |
| court réellement et donne |
| réponses, mais la confiance |
| les intervalles sont extrêmement larges |
| ici, comme on pouvait s'y attendre. Mais |
| les ingénieurs savaient ce que |
| mécanisme de défaillance était et ils |
| avait eu une expérience antérieure et |
| pour qu'ils puissent apporter ça |
| informations sur la pente vers |
| l'analyse à l'aide de Bayes |
| méthodes. Alors encore une fois, voici le |
| articulation postérieure et la largeur de |
| la distribution dans le |
| postérieur pour bêta 1 est |
| effectivement ce que nous avons supposé |
| quand on met en avant |
| distribution pour ce paramètre. |
| Alors encore une fois, voici le |
| spécification du prieur |
| distributions, où nous avons utilisé |
| faiblement informatif pour le |
| quantile et pour Sigma, mais |
| distribution préalable informative |
| pour la pente beta 1. Et je peux |
| obtenir une estimation du temps à |
| dont 1% échouent. Donc l'extrémité inférieure |
| point de l'intervalle de confiance |
| pour le moment auquel 1% échouera |
| est plus de 140 000 heures. |
| Donc ça fait environ 20 ans, beaucoup |
| plus long que le technologique |
| la vie de ces produits dans laquelle |
| ce circuit intégré sera |
| utilisé. Alors qu'avons-nous appris ici? |
| Eh bien, dans certaines applications, nous |
| avoir une censure par intervalles car |
| les échecs sont découverts uniquement |
| quand il y a une inspection. Nous |
| besoin de statistiques appropriées |
| les méthodes de traitement de ces données, |
| et JMP a ces méthodes. Si |
| vous utilisez des niveaux excessifs d'un |
| accélération variable comme |
| température, vous pouvez générer |
| de nouveaux modes de défaillance qui rendent le |
| informations trompeuses. Alors on |
| a dû jeter ces unités. |
| Mais même avec des échecs à seulement |
| un niveau de température, si nous |
| avoir des informations préalables |
| à propos de l'activation effective |
| l'énergie, on peut combiner ça |
| informations avec le limité |
| données pour faire des inférences utiles. |
| Enfin, certains concluent |
| remarques, améliorations de |
| matériel informatique et logiciels |
| ont grandement amélioré notre capacité |
| pour analyser la fiabilité et autres |
| Les données. Maintenant, nous pouvons également utiliser Bayes |
| méthodes, fournissant un autre ensemble de |
| outils pour combiner des informations |
| avec des données limitées et JMP a |
| des outils puissants pour ce faire. |
| Donc, bien que ces bayésiens |
| des capacités ont été développées pour |
| la partie fiabilité de JMP, |
| ils peuvent certainement être utilisés dans |
| autres domaines d'application. Et |
| voici quelques références, |
| dont la 2ème édition de |
| Fiabilité des méthodes statistiques, |
| qui devrait être probablement dans |
| Juin 2021. OK, alors maintenant je suis |
| va céder la parole à Peng |
| et il va vous montrer comment |
| il est facile de faire ces analyses. |
| Merci professeur. |
| Avant de commencer mon |
| démonstration, je voudrais |
| pour montrer cette diapositive sur |
| Flux de travail d'analyse bayésienne |
| dans la distribution de la vie et |
| Fit Life par X. |
| Tout d'abord, vous devez adapter un |
| modèle paramétrique utilisant le maximum |
| probabilité. Je t'assume déjà |
| savoir comment faire cela dans ces deux |
| plates-formes. Ensuite, vous devez |
| examiner ou trouver un modèle |
| utilisateur graphique de spécification |
| interface pour bayésien |
| estimation dans le rapport |
| de l'étape précédente. Pour |
| exemple, ceci est une capture d'écran de |
| Modèle Weibull dans la vie |
| Distribution. Vous devez aller au |
| menu triangle rouge |
| et choisissez des estimations bayésiennes |
| pour révéler le graphique |
| interface utilisateur pour le |
| Analyse bayésienne. |
| Dans Fit Life by X, veuillez consulter |
| la capture d'écran d'un lognormal |
| résultat. Et l'utilisateur graphique |
| interface pour le bayésien |
| l'inférence est sur la dernière étape. |
| Après avoir trouvé l'utilisateur graphique |
| interface pour l'analyse bayésienne, |
| vous devrez fournir le |
| informations sur les prieurs. |
| Vous devez décider des priors |
| dispersion individuelle |
| paramètres. Vous devez fournir |
| les informations pour le |
| hyperparamètres et autres |
| des informations telles que |
| probabilité pour le quantile. |
| En plus de cela, nous |
| besoin de fournir le numéro |
| d'échantillons postérieurs. |
| Besoin également de fournir un aléatoire |
| graine au cas où vous voudriez |
| reproduisez votre résultat dans le |
| avenir. Ensuite, vous pouvez cliquer sur Fit |
| Modèle. Cela générera un |
| rapport du modèle. |
| Vous pouvez adapter plusieurs modèles |
| au cas où vous voudriez étudier le |
| sensibilité de différents |
| Modèles bayésiens donnés |
| distribution antérieure différente. |
| Le résultat d'un modèle bayésien, |
| y compris les choses suivantes, |
| la première est une méthode d'échantillonnage. |
| Et puis est une copie des priors. |
| Puis a eu une estimation postérieure |
| des paramètres. |
| Et puis il y a des nuages de points, |
| soit pour le prieur, soit |
| les postérieurs. |
| En fin de compte, nous avons deux |
| profileurs. Un pour la distribution |
| et celui du quantile. |
| En utilisant ces résultats, vous |
| peut faire d'autres inférences |
| comme la prédiction des pannes. |
| Regardez maintenant la démonstration. |
| Nous démontrerons avec le |
| dernier exemple que le professeur |
| mentionné que ??? présentation. |
| C'est le périphérique IC là-bas. |
| Nous avons deux colonnes pour l'heure |
| à l'événement | HeuresL et HeuresU à |
| représentent la situation de censure. |
| Nous avons un compte |
| pour l'observation individuelle et |
| température individuelle |
| observation. Nous excluons le dernier |
| quatre observations parce qu'elles |
| sont associés à un autre |
| mode de défaillance. Nous voulons exclure |
| ces observations de la |
| Analyse. Maintenant, commencez à spécifier. |
| ???. |
| Nous mettons des heures dans Y |
| Nous mettons Count en fréquence. |
| Nous mettons les degrés C dans X. |
| Nous utilisons l'Arrhenius Celsius |
| pour notre relation. |
| Nous utilisons lognormal pour notre |
| Distribution. |
| Cliquez ensuite sur OK. |
| Le résultat est le maximum |
| influence de la probabilité |
| pour le lognormal. Nous allons à la |
| Estimations bayésiennes, et |
| commencez à spécifier nos priors comme |
| Le professeur a fait dans son |
| présentation. Nous choisissons un |
| lognormal pour le quantile. |
| Il fait 250 degrés C. |
| Et sa vie B10. Alors |
| la probabilité est de 0,1. |
| Les deux extrémités de la |
| la distribution log-normale est de 100 |
| et 10 000. |
| Spécifiez maintenant la pente. |
| La distribution est log-normale. |
| Et deux extrémités de la |
| la distribution est de 0,65 |
| et .85 parce que c'est informatif |
| pour exiger la gamme est étroite. |
| Maintenant, nous spécifions le prieur |
| distribution pour Sigma, |
| qui est un lognormal |
| et il avait une large gamme; |
| c'est 0,05 et 5. |
| Nous avons décidé de dessiner 5000 |
| échantillons postérieurs. Et attribuez un |
| graine aléatoire arbitraire. Puis |
| nous cliquons ... cliquez sur Fit Model. |
| Et ce que le rapport génère pour cela |
| spécification. La méthode est |
| rejet simple. Et voici un |
| copie de notre spécification proir. |
| Les estimations postérieures |
| résumer nos échantillons postérieurs. |
| Vous pouvez exporter les échantillons postérieurs |
| soit en cliquant sur cette exportation |
| Échantillons Monte Carlo |
| ou choisissez-le dans le menu, |
| qui est ici | Exporter Monte |
| Échantillons de Carlo. |
| Des échantillons postérieurs sont illustrés |
| dans ces nuages de points. |
| Nous avons ici deux diagrammes de dispersion. |
| Le premier à utiliser |
| la même paramatérisation que la précédente |
| spécification, qui est ... qui |
| utiliser le quantile, la pente et le signal. |
| L'autre nuage de points ... le |
| le deuxième nuage de points utilise un |
| paramétrage traditionnel, |
| qui comprend l'interception de |
| régression, une pente de la |
| régression et Sigma. |
| En fin de compte faire ???, nous |
| pouvez-vous nous pouvons regarder le |
| profileur. Regardons ici |
| au deuxième profileur, |
| profil quantile, afin que nous puissions |
| trouver le même résultat |
| comme ce que le professeur avait montré dans |
| l'une des diapositives précédentes. |
| Entrez 0,1 ... |
| 0,01 pour la probabilité. Alors ceci est |
| 1%. Nous entrons 100 degrés C |
| pour DegreesC. |
| Et nous ajustons |
| les axes. |
| Alors maintenant, nous voyons un semblable |
| profileur. |
| Il a déjà été dans |
| la diapositive précédente. |
| Et nous pouvons lire sur l'axe Y |
| pour obtenir le résultat souhaité, qui |
| est le temps pendant lequel 1% du |
| l'appareil échouera à 100 degrés |
| C. Donc, cela conclut mon |
| manifestation. Et laissez |
| je passe à autre chose. |
| Cette diapositive explique |
| about ... explique l'implémentation JMP |
| d'algorithme d'échantillonnage. |
| Nous avons vu que le simple |
| le rejet est apparu dans le |
| exemple précédent et c'est le |
| première étape de notre |
| la mise en oeuvre. Le simple |
| l'algorithme de rejet est essayé et |
| vraie méthode à vos échantillons, |
| mais cela peut être peu pratique |
| si le taux de rejet est élevé. |
| Donc si le rejet |
| le taux est élevé, nous ... |
| les implémentations passent à |
| deuxième étape, qui est un aléatoire |
| Marcher Metropolis-Hastings |
| algorithme. Le deuxième algorithme |
| est efficace, mais au cas où ... dans |
| situations ça peut échouer |
| indétectable si la probabilité |
| est irrégulière. Par exemple, le |
| la probabilité est plutôt plate. Nous |
| conçu cette implémentation |
| parce que nous avons une situation, |
| il y a très peu de pannes ou |
| même pas d'échecs. En cela |
| situation, la probabilité est |
| relativement élevé, mais ??? |
| situation, nous utilisons simple |
| l'algorithme de rejet et le |
| le taux de rejet n'est pas si mal |
| et cette méthode suffira. |
| Quand on a de plus en plus |
| échecs, la probabilité |
| devient plus régulier. Donc il a |
| un pic au milieu. En cela |
| situation, le simple |
| taux de rejet ... le simple |
| la méthode de rejet devient |
| peu pratique à cause de la haute |
| taux de réjection. Mais le hasard |
| L'algorithme de marche devient de plus en plus |
| plus prometteur pour réussir |
| sans échec. Alors c'est notre |
| mise en œuvre et explication de |
| pourquoi nous faisons cela. |
| Cette diapositive explique comment nous |
| spécifier un précédent normal tronqué |
| dans ces deux plates-formes. Parce que |
| la normale tronquée n'est pas un |
| construction de la distribution préalable dans |
| ces deux plates-formes. |
| Tout d'abord, regardez ce que |
| normal tronqué. Ici nous |
| donner un exemple de normale tronquée |
| avec deux extrémités à 5 et 400. Le |
| deux extrémités sont illustrées par |
| ce L et ce R. |
| Mais la normale tronquée n'est rien |
| mais une normale en rejetant |
| toutes les valeurs qui sont |
| négatif, qui est représenté |
| par cette courbe traversante comme le |
| distribution normale équivalente |
| pour ce tronqué particulier |
| distribution normale. |
| Si nous voulons le spécifier |
| tronqué normal ou nous devons |
| définir est équivalent ... équivalent |
| distribution normale à deux extrémités |
| qui avait eu le même nouveau |
| et les paramètres Sigma comme ceux de |
| ces normaux tronqués |
| distributions. Nous fournissons donc un |
| script pour ce faire. |
| Dans ce script, le |
| le calcul est le |
| Suivant. Nous trouvons le nouveau |
| et Sigma du tronqué |
| Ordinaire. |
| Donc, nous pouvons obtenir le |
| équivalent normal |
| Distribution. |
| Et nous ça ??? Sigma de cela |
| distribution normale que nous pouvons trouver |
| les deux extrémités de cette normale |
| Distribution. |
| Donc pour spécifier le tronqué |
| normal avec ces deux extrémités |
| valeur pour spécifier l'équivalent |
| distribution normale avec deux |
| fin ... ces deux valeurs finales. |
| C'est ainsi que nous spécifions |
| normal tronqué dans ces |
| deux plates-formes utilisant |
| équivalent normal |
| Distribution. Et voici le |
| contenu du script. |
| Tout ce que vous avez à faire est d'appeler |
| une fonction, qui est ici le |
| inverser le paramètre tnorm à |
| valeur normale. De quoi avez-vous besoin |
| à fournir sont les deux extrémités de |
| la normale tronquée |
| distribution et ça donnera |
| les deux extrémités de l'équivalent |
| distribution normale et vous pouvez |
| utilisez ces deux nombres pour spécifier |
| la distribution antérieure. |
| Donc, cela conclut mon |
| manifestation. Dans mon |
| démonstration j'ai montré comment |
| démarrer l'analyse bayésienne dans Life |
| distribution et Fix Life par |
| X, comment entrer avant |
| informations, |
| et quel est le contenu de |
| Résultat bayésien. J'ai aussi expliqué |
| ce que notre mise en œuvre de la |
| échantillonnage et pourquoi faites-vous |
| cette. Et à la fin j'explique |
| comment spécifier un tronqué |
| normal avant d'utiliser un |
| équivalent normal antérieur dans |
| ces deux ... dans ces deux |
| plates-formes. Merci. |
