5
$\begingroup$

Let's consider two balanced dice. I am modelling each of them using

d1 = ProbabilityDistribution[1/6, {x, 1, 6, 1}]
d2 = ProbabilityDistribution[1/6, {x, 1, 6, 1}]

Call $x$ the result of the first dice and $y$ the result of the second dice after they are thrown, let $z=2x+y$. How can I make a probability distribution corresponding to $z$ using $d_1$ and $d_2$ ? I tried adding them

d3=2*d1+d2

but this didn't work.

Improve this question
$\endgroup$
1
  • 2
    $\begingroup$ I think you want TransformedDistribution. $\endgroup$
    – Eric Brown
    Nov 26, 2015 at 18:57

1 Answer 1

Reset to default
11
$\begingroup$

As Eric Brown suggested you should use TransformedDistribution[]. Also I'd suggest using built in distributions if possible, here I mean DiscreteUniformDistribution[]. 

tf = TransformedDistribution[
  A + B, {A \[Distributed] DiscreteUniformDistribution[{1, 6}], 
   B \[Distributed] DiscreteUniformDistribution[{1, 6}]}]
PDF[tf, x]
DiscretePlot[PDF[tf, x], {x, 2, 12}]

enter image description here

Improve this answer
$\endgroup$
1
  • 3
    $\begingroup$ Just as a tiny comment: it would be better to first evaluate: ee = PDF[tf, x], and then DiscretePlot[ee {x, 2, 12}] . By contrast, DiscretePlot[PDF[tf, x], {x, 2, 12}] might take about 100 times longer to plot, because it tells Mma to calculate PDF[t,x] at each value of x (rather than just solving it once). $\endgroup$
    – wolfies
    Nov 28, 2015 at 17:37

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge that you have read and understand our privacy policy and code of conduct.

Not the answer you're looking for? Browse other questions tagged or ask your own question.