Gsyfuy

Implementing 3DES algorithm in Java: is my code secure?

How to translate "red flag" into Spanish?

Co-worker works way more than he should

Is Diceware more secure than a long passphrase?

Is Electric Central Heating worth it if using Solar Panels?

Multiple fireplaces in an apartment building?

Error: Syntax error. Missing ')' for CASE Statement

Is a 5 watt UHF/VHF handheld considered QRP?

What is /etc/mtab in Linux?

All ASCII characters with a given bit count

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

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

Function to calculate red-edgeNDVI in Google Earth Engine

Justification for leaving new position after a short time

Will I lose my paid in full property

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

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

Need of separate security plugins for both root and subfolder sites Wordpress?

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

My bank got bought out, am I now going to have to start filing tax returns in a different state?

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

Is Bran literally the world's memory?

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

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

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