Does Mathematica have an implementation of the Poisson Binomial Distribution? Announcing the arrival of Valued Associate #679: Cesar Manara Unicorn Meta Zoo #1: Why another podcast?RandomVariate from 2-dimensional probability distributionBayesian Inference with Continuous prior distributionGrain(Particle) Size Distribution (PSD) Analysis with MathematicaHow to get probabilities for multinomial & hypergeometric distribution ranges more quickly?Given an exact formula, how can Mathematica find a probability distribution whose PDF matches it?Probability of an eventLikelihood for BetaBinomialDistribution with variable number of trialsInverse Fourier Transform of Poisson Characteristic Function?How to use EventData correctly to model a trial sequenceComputing the Expectation

Do I need to protect SFP ports and optics from dust/contaminants? If so, how?

Second order approximation of the loss function (Deep learning book, 7.33)

Has a Nobel Peace laureate ever been accused of war crimes?

Split coins into combinations of different denominations

A strange hotel

Does Feeblemind produce an ongoing magical effect that can be dispelled?

Are all CP/M-80 implementations binary compatible?

Where did Arya get these scars?

Could Neutrino technically as side-effect, incentivize centralization of the bitcoin network?

Why did C use the -> operator instead of reusing the . operator?

Align column where each cell has two decimals with siunitx

What is a 'Key' in computer science?

Married in secret, can marital status in passport be changed at a later date?

How to get even lighting when using flash for group photos near wall?

What do you call the part of a novel that is not dialog?

Book with legacy programming code on a space ship that the main character hacks to escape

Did the Roman Empire have penal colonies?

Why did Israel vote against lifting the American embargo on Cuba?

Multiple options vs single option UI

What is it called when you ride around on your front wheel?

I preordered a game on my Xbox while on the home screen of my friend's account. Which of us owns the game?

What is the term for a person whose job is to place products on shelves in stores?

Israeli soda type drink

With indentation set to `0em`, when using a line break, there is still an indentation of a size of a space



Does Mathematica have an implementation of the Poisson Binomial Distribution?



Announcing the arrival of Valued Associate #679: Cesar Manara
Unicorn Meta Zoo #1: Why another podcast?RandomVariate from 2-dimensional probability distributionBayesian Inference with Continuous prior distributionGrain(Particle) Size Distribution (PSD) Analysis with MathematicaHow to get probabilities for multinomial & hypergeometric distribution ranges more quickly?Given an exact formula, how can Mathematica find a probability distribution whose PDF matches it?Probability of an eventLikelihood for BetaBinomialDistribution with variable number of trialsInverse Fourier Transform of Poisson Characteristic Function?How to use EventData correctly to model a trial sequenceComputing the Expectation










8












$begingroup$


I need to work out the probability of having $k$ successful trials out of a total of $n$ when success probabilities are heterogeneous. This calculation relates to the Poisson Binomial Distribution. Does Mathematica, or perhaps the Mathstatica add-on, have an implementation for that?










share|improve this question











$endgroup$
















    8












    $begingroup$


    I need to work out the probability of having $k$ successful trials out of a total of $n$ when success probabilities are heterogeneous. This calculation relates to the Poisson Binomial Distribution. Does Mathematica, or perhaps the Mathstatica add-on, have an implementation for that?










    share|improve this question











    $endgroup$














      8












      8








      8


      2



      $begingroup$


      I need to work out the probability of having $k$ successful trials out of a total of $n$ when success probabilities are heterogeneous. This calculation relates to the Poisson Binomial Distribution. Does Mathematica, or perhaps the Mathstatica add-on, have an implementation for that?










      share|improve this question











      $endgroup$




      I need to work out the probability of having $k$ successful trials out of a total of $n$ when success probabilities are heterogeneous. This calculation relates to the Poisson Binomial Distribution. Does Mathematica, or perhaps the Mathstatica add-on, have an implementation for that?







      probability-or-statistics






      share|improve this question















      share|improve this question













      share|improve this question




      share|improve this question








      edited 6 hours ago









      Chris K

      7,32722143




      7,32722143










      asked 6 hours ago









      user120911user120911

      81838




      81838




















          1 Answer
          1






          active

          oldest

          votes


















          10












          $begingroup$

          Mathematica does not know about the PoissonBinomialDistribution, but you can use the formula given for the PDF on Wikipedia:



          PoissonBinomialDistribution[ plist : __Real ] := With[

          n = Length @ plist,
          c = Exp[(2 I [Pi])/(Length@plist + 1)]

          ,
          ProbabilityDistribution[
          1/(n + 1) Sum[c^(-l k) Product[1 + (c^l - 1) plist[[m]], m, 1, n ], l, 0, n]
          ,
          k, 0, n, 1
          ]
          ]


          Now we may model a quality control where fault type 1 has a prob of 4% and fault types 2 and 3 have a prob of 7%:



          dist = PoissonBinomialDistribution[ 0.04, 0.07, 0.07 ];


          With this we find the probability for 3 faults:



          Probability[ k == 3, k [Distributed] dist ]// PercentForm



          0.0196 %







          share|improve this answer











          $endgroup$













            Your Answer








            StackExchange.ready(function()
            var channelOptions =
            tags: "".split(" "),
            id: "387"
            ;
            initTagRenderer("".split(" "), "".split(" "), channelOptions);

            StackExchange.using("externalEditor", function()
            // Have to fire editor after snippets, if snippets enabled
            if (StackExchange.settings.snippets.snippetsEnabled)
            StackExchange.using("snippets", function()
            createEditor();
            );

            else
            createEditor();

            );

            function createEditor()
            StackExchange.prepareEditor(
            heartbeatType: 'answer',
            autoActivateHeartbeat: false,
            convertImagesToLinks: false,
            noModals: true,
            showLowRepImageUploadWarning: true,
            reputationToPostImages: null,
            bindNavPrevention: true,
            postfix: "",
            imageUploader:
            brandingHtml: "Powered by u003ca class="icon-imgur-white" href="https://imgur.com/"u003eu003c/au003e",
            contentPolicyHtml: "User contributions licensed under u003ca href="https://creativecommons.org/licenses/by-sa/3.0/"u003ecc by-sa 3.0 with attribution requiredu003c/au003e u003ca href="https://stackoverflow.com/legal/content-policy"u003e(content policy)u003c/au003e",
            allowUrls: true
            ,
            onDemand: true,
            discardSelector: ".discard-answer"
            ,immediatelyShowMarkdownHelp:true
            );



            );













            draft saved

            draft discarded


















            StackExchange.ready(
            function ()
            StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fmathematica.stackexchange.com%2fquestions%2f196962%2fdoes-mathematica-have-an-implementation-of-the-poisson-binomial-distribution%23new-answer', 'question_page');

            );

            Post as a guest















            Required, but never shown

























            1 Answer
            1






            active

            oldest

            votes








            1 Answer
            1






            active

            oldest

            votes









            active

            oldest

            votes






            active

            oldest

            votes









            10












            $begingroup$

            Mathematica does not know about the PoissonBinomialDistribution, but you can use the formula given for the PDF on Wikipedia:



            PoissonBinomialDistribution[ plist : __Real ] := With[

            n = Length @ plist,
            c = Exp[(2 I [Pi])/(Length@plist + 1)]

            ,
            ProbabilityDistribution[
            1/(n + 1) Sum[c^(-l k) Product[1 + (c^l - 1) plist[[m]], m, 1, n ], l, 0, n]
            ,
            k, 0, n, 1
            ]
            ]


            Now we may model a quality control where fault type 1 has a prob of 4% and fault types 2 and 3 have a prob of 7%:



            dist = PoissonBinomialDistribution[ 0.04, 0.07, 0.07 ];


            With this we find the probability for 3 faults:



            Probability[ k == 3, k [Distributed] dist ]// PercentForm



            0.0196 %







            share|improve this answer











            $endgroup$

















              10












              $begingroup$

              Mathematica does not know about the PoissonBinomialDistribution, but you can use the formula given for the PDF on Wikipedia:



              PoissonBinomialDistribution[ plist : __Real ] := With[

              n = Length @ plist,
              c = Exp[(2 I [Pi])/(Length@plist + 1)]

              ,
              ProbabilityDistribution[
              1/(n + 1) Sum[c^(-l k) Product[1 + (c^l - 1) plist[[m]], m, 1, n ], l, 0, n]
              ,
              k, 0, n, 1
              ]
              ]


              Now we may model a quality control where fault type 1 has a prob of 4% and fault types 2 and 3 have a prob of 7%:



              dist = PoissonBinomialDistribution[ 0.04, 0.07, 0.07 ];


              With this we find the probability for 3 faults:



              Probability[ k == 3, k [Distributed] dist ]// PercentForm



              0.0196 %







              share|improve this answer











              $endgroup$















                10












                10








                10





                $begingroup$

                Mathematica does not know about the PoissonBinomialDistribution, but you can use the formula given for the PDF on Wikipedia:



                PoissonBinomialDistribution[ plist : __Real ] := With[

                n = Length @ plist,
                c = Exp[(2 I [Pi])/(Length@plist + 1)]

                ,
                ProbabilityDistribution[
                1/(n + 1) Sum[c^(-l k) Product[1 + (c^l - 1) plist[[m]], m, 1, n ], l, 0, n]
                ,
                k, 0, n, 1
                ]
                ]


                Now we may model a quality control where fault type 1 has a prob of 4% and fault types 2 and 3 have a prob of 7%:



                dist = PoissonBinomialDistribution[ 0.04, 0.07, 0.07 ];


                With this we find the probability for 3 faults:



                Probability[ k == 3, k [Distributed] dist ]// PercentForm



                0.0196 %







                share|improve this answer











                $endgroup$



                Mathematica does not know about the PoissonBinomialDistribution, but you can use the formula given for the PDF on Wikipedia:



                PoissonBinomialDistribution[ plist : __Real ] := With[

                n = Length @ plist,
                c = Exp[(2 I [Pi])/(Length@plist + 1)]

                ,
                ProbabilityDistribution[
                1/(n + 1) Sum[c^(-l k) Product[1 + (c^l - 1) plist[[m]], m, 1, n ], l, 0, n]
                ,
                k, 0, n, 1
                ]
                ]


                Now we may model a quality control where fault type 1 has a prob of 4% and fault types 2 and 3 have a prob of 7%:



                dist = PoissonBinomialDistribution[ 0.04, 0.07, 0.07 ];


                With this we find the probability for 3 faults:



                Probability[ k == 3, k [Distributed] dist ]// PercentForm



                0.0196 %








                share|improve this answer














                share|improve this answer



                share|improve this answer








                edited 5 hours ago

























                answered 5 hours ago









                gwrgwr

                8,74322862




                8,74322862



























                    draft saved

                    draft discarded
















































                    Thanks for contributing an answer to Mathematica Stack Exchange!


                    • Please be sure to answer the question. Provide details and share your research!

                    But avoid


                    • Asking for help, clarification, or responding to other answers.

                    • Making statements based on opinion; back them up with references or personal experience.

                    Use MathJax to format equations. MathJax reference.


                    To learn more, see our tips on writing great answers.




                    draft saved


                    draft discarded














                    StackExchange.ready(
                    function ()
                    StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fmathematica.stackexchange.com%2fquestions%2f196962%2fdoes-mathematica-have-an-implementation-of-the-poisson-binomial-distribution%23new-answer', 'question_page');

                    );

                    Post as a guest















                    Required, but never shown





















































                    Required, but never shown














                    Required, but never shown












                    Required, but never shown







                    Required, but never shown

































                    Required, but never shown














                    Required, but never shown












                    Required, but never shown







                    Required, but never shown







                    Popular posts from this blog

                    How to create a command for the “strange m” symbol in latex? Announcing the arrival of Valued Associate #679: Cesar Manara Planned maintenance scheduled April 23, 2019 at 23:30 UTC (7:30pm US/Eastern)How do you make your own symbol when Detexify fails?Writing bold small caps with mathpazo packageplus-minus symbol with parenthesis around the minus signGreek character in Beamer document titleHow to create dashed right arrow over symbol?Currency symbol: Turkish LiraDouble prec as a single symbol?Plus Sign Too Big; How to Call adfbullet?Is there a TeX macro for three-legged pi?How do I get my integral-like symbol to align like the integral?How to selectively substitute a letter with another symbol representing the same letterHow do I generate a less than symbol and vertical bar that are the same height?

                    Българска екзархия Съдържание История | Български екзарси | Вижте също | Външни препратки | Литература | Бележки | НавигацияУстав за управлението на българската екзархия. Цариград, 1870Слово на Ловешкия митрополит Иларион при откриването на Българския народен събор в Цариград на 23. II. 1870 г.Българската правда и гръцката кривда. От С. М. (= Софийски Мелетий). Цариград, 1872Предстоятели на Българската екзархияПодмененият ВеликденИнформационна агенция „Фокус“Димитър Ризов. Българите в техните исторически, етнографически и политически граници (Атлас съдържащ 40 карти). Berlin, Königliche Hoflithographie, Hof-Buch- und -Steindruckerei Wilhelm Greve, 1917Report of the International Commission to Inquire into the Causes and Conduct of the Balkan Wars

                    Category:Tremithousa Media in category "Tremithousa"Navigation menuUpload media34° 49′ 02.7″ N, 32° 26′ 37.32″ EOpenStreetMapGoogle EarthProximityramaReasonatorScholiaStatisticsWikiShootMe