Mryntu

How to not starve gigantic beasts

What is ls Largest Number Formed by only moving two sticks in 508?

What *exactly* is electrical current, voltage, and resistance?

Is Bran literally the world's memory?

"Whatever a Russian does, they end up making the Kalashnikov gun"? Are there any similar proverbs in English?

What is the ongoing value of the Kanban board to the developers as opposed to management

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

SQL Query not selecting all points that it should?

Will I lose my paid in full property

Protagonist's race is hidden - should I reveal it?

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

c++ diamond problem - How to call base method only once

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

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

Can you stand up from being prone using Skirmisher outside of your turn?

How to keep bees out of canned beverages?

My admission is revoked after accepting the admission offer

A strange hotel

Would reducing the reference voltage of an ADC have any effect on accuracy?

Expansion//Explosion and Siren Stormtamer

What if Force was not Mass times Acceleration?

What was Apollo 13's "Little Jolt" after MECO?

How to find the right literary agent in the USA?

"Rubric" as meaning "signature" or "personal mark" -- is this accepted usage?

Does Mathematica have an implementation of the Poisson Binomial Distribution?

6

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 2 hours ago

Chris K

7,32722143

asked 2 hours ago

user120911user120911

80838

$endgroup$

add a comment | 

6

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 2 hours ago

Chris K

7,32722143

asked 2 hours ago

user120911user120911

80838

$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 2 hours ago

Chris K

7,32722143

asked 2 hours ago

user120911user120911

80838

$endgroup$

share|improve this question

edited 2 hours ago

Chris K

7,32722143

asked 2 hours ago

user120911user120911

80838

share|improve this question

edited 2 hours ago

Chris K

7,32722143

asked 2 hours ago

user120911user120911

80838

edited 2 hours ago

Chris K

7,32722143

edited 2 hours ago

Chris K

7,32722143

asked 2 hours ago

user120911user120911

80838

asked 2 hours ago

user120911user120911

80838

1 Answer
1

active

oldest

votes

5

$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 1 hour ago

answered 1 hour ago

gwrgwr

8,69322861

$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

5

$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 1 hour ago

answered 1 hour ago

gwrgwr

8,69322861

$endgroup$

add a comment | 

5

$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 1 hour ago

answered 1 hour ago

gwrgwr

8,69322861

$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 1 hour ago

answered 1 hour ago

gwrgwr

8,69322861

$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 1 hour ago

answered 1 hour ago

gwrgwr

8,69322861

answered 1 hour ago

gwrgwr

8,69322861

answered 1 hour ago

gwrgwr

8,69322861