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

                    The Calvary Singular or Plural The 2019 Stack Overflow Developer Survey Results Are InAre collective nouns always plural, or are certain ones singular?Is “audience” singular or plural?“Wasn't” vs. “weren't” in a vernacular sentence“My last couple of years” — singular or plural?Is 'rest' singular or plural?Is “all but one” singular or plural?Whether to use the singular or plural form of basis?Singular and Plural for numbersIs there a plural form of teeth?performance: plural vs singular?singular or plural nouns?Singular and Plural

                    How does one intimidate enemies without having the capacity for violence?Ideas for how aliens would approach this fight?How to convey the scale of my humanoid without science or units?How does the “space drive” conserve momentum?Planet Vanishes - How does this affect the orbiting starships?How does a community of a Universe Simulator have the same language as its creator?How would US Presidential elections be affected if voters could choose the state their vote for President was counted in?How Does One Ensures the Immortality of Their ConsciousnessHow do I retain national independence while also having a one world government?How can Ganymede have an Earth-like gravity without us having realized it?How would one make a lion mount for a fantasy world?

                    Output visual diagram of pictureASCII-art logic gate diagramBooks on a ShelfDetermine the Dimensions of a Rotated RectangleDraw a Houndstooth PatternDraw and label an ASCII hexagonal gridGolf me an ASCII AlphabetASCII Jigsaw PuzzleOutput a pretty boxASCII-Art Venn DiagramASCII Exact Cover with Rectangles