Kdtyjb

Retract an already submitted recommendation letter (written for an undergrad student)

What's the difference between using dependency injection with a container and using a service locator?

Are these square matrices always diagonalisable?

How do I check if a string is entirely made of the same substring?

Is Bran literally the world's memory?

Split coins into combinations of different denominations

Is Electric Central Heating worth it if using Solar Panels?

"My boss was furious with me and I have been fired" vs. "My boss was furious with me and I was fired"

Putting Ant-Man on house arrest

Are all CP/M-80 implementations binary compatible?

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

Map material from china not allowed to leave the country

Why didn't the Space Shuttle bounce back into space as many times as possible so as to lose a lot of kinetic energy up there?

My admission is revoked after accepting the admission offer

How to not starve gigantic beasts

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

Justification for leaving new position after a short time

Co-worker works way more than he should

Trumpet valves, lengths, and pitch

Align column where each cell has two decimals with siunitx

Arriving in Atlanta after US Preclearance in Dublin. Will I go through TSA security in Atlanta to transfer to a connecting flight?

Identify story/novel: Tribe on colonized planet, not aware of this. "Taboo," altitude sickness, robot guardian (60s? Young Adult?)

Why does the Cisco show run command not show the full version, while the show version command does?

How to keep bees out of canned beverages?

Does Mathematica have an implementation of the Poisson Binomial Distribution?

7

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?

probability-or-statistics

share|improve this question

edited 3 hours ago

Chris K

7,32722143

asked 3 hours ago

user120911user120911

81338

$endgroup$

add a comment | 

7

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?

probability-or-statistics

share|improve this question

edited 3 hours ago

Chris K

7,32722143

asked 3 hours ago

user120911user120911

81338

$endgroup$

add a comment | 

$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?

probability-or-statistics

share|improve this question

edited 3 hours ago

Chris K

7,32722143

asked 3 hours ago

user120911user120911

81338

$endgroup$

share|improve this question

edited 3 hours ago

Chris K

7,32722143

asked 3 hours ago

user120911user120911

81338

share|improve this question

edited 3 hours ago

Chris K

7,32722143

asked 3 hours ago

user120911user120911

81338

edited 3 hours ago

Chris K

7,32722143

edited 3 hours ago

Chris K

7,32722143

asked 3 hours ago

user120911user120911

81338

asked 3 hours ago

user120911user120911

81338

1 Answer
1

active

oldest

votes

9

$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

edited 2 hours ago

answered 2 hours ago

gwrgwr

8,73322861

$endgroup$

add a comment | 

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

Sign up or log in

StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});

Sign up using Google

Sign up using Facebook

Sign up using Email and Password

Post as a guest

Name

Email

Required, but never shown

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

Name

Email

Required, but never shown

9

$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

edited 2 hours ago

answered 2 hours ago

gwrgwr

8,73322861

$endgroup$

add a comment | 

9

$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

edited 2 hours ago

answered 2 hours ago

gwrgwr

8,73322861

$endgroup$

add a comment | 

$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

edited 2 hours ago

answered 2 hours ago

gwrgwr

8,73322861

$endgroup$

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}

    ]

 ]

dist = PoissonBinomialDistribution[ {0.04, 0.07, 0.07} ];

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

share|improve this answer

edited 2 hours ago

answered 2 hours ago

gwrgwr

8,73322861

answered 2 hours ago

gwrgwr

8,73322861

answered 2 hours ago

gwrgwr

8,73322861