Mdtryhht

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?

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?

probability-or-statistics

share|improve this question

edited 6 hours ago

Chris K

7,32722143

asked 6 hours ago

user120911user120911

81838

$endgroup$

add a comment | 

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?

probability-or-statistics

share|improve this question

edited 6 hours ago

Chris K

7,32722143

asked 6 hours ago

user120911user120911

81838

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

Chris K

7,32722143

asked 6 hours ago

user120911user120911

81838

$endgroup$

share|improve this question

edited 6 hours ago

Chris K

7,32722143

asked 6 hours ago

user120911user120911

81838

share|improve this question

edited 6 hours ago

Chris K

7,32722143

asked 6 hours ago

user120911user120911

81838

edited 6 hours ago

Chris K

7,32722143

edited 6 hours ago

Chris K

7,32722143

asked 6 hours ago

user120911user120911

81838

asked 6 hours ago

user120911user120911

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

edited 5 hours ago

answered 5 hours ago

gwrgwr

8,74322862

$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

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

edited 5 hours ago

answered 5 hours ago

gwrgwr

8,74322862

$endgroup$

add a comment | 

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

edited 5 hours ago

answered 5 hours ago

gwrgwr

8,74322862

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

answered 5 hours ago

gwrgwr

8,74322862

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

answered 5 hours ago

gwrgwr

8,74322862

answered 5 hours ago

gwrgwr

8,74322862

answered 5 hours ago

gwrgwr

8,74322862