As we saw with the cumulative distribution FX for a single continuous random variable, there is a function FX,Y called the joint probability distribution function.
FX,Y(x,y)=P(X≤x,Y≤y)
Just as there is the probability density function to the cumulative distribution, we also have a joint probability density function.
fX,Y(x,y)=δxδyδ2FX,Y(x,y)
FX,Y(x,y)=∫−∞x∫−∞yfX,Y(s,t),dt,ds
This can also be thought of as a sum of discrete joint probabilities:
FX,Y(xi,yj)=P(X≤xi,Y≤yj)=k,l∑P(X=xk,Y=yl)
For k,l:xk≤xi,yl≤yj.
Double integrals
While double integrals might appear more complicated than single ones we have looked at before, they are actually very similar just with more steps. For example, looking at the following example:
FX,Y(x,y)=∫−∞x∫−∞yfX,Y(s,t),dt,ds
We would take the first integral, the one closest to the function, with respect to t.
FX,Y(x,y)=∫−∞x[∫−∞yfX,Y(s,t),dt],ds
Then take the outer one of that answer with respect to s.
When each integral is calculated, all not respecting variables are treated as constants. The two integrals can be calculated in any order as long as the variables they respect are correct. You can find a reminder of some useful integration techniques in the calculus notes.
Graphical area
Problems in this topic often take the form of the probability of two random variables, X & Y, each between two bounds.
P(a≤X≤b,c≤Y≤d)=∫cd∫abfX,Y(s,t),ds,dt
If you are able to work out the distribution and get your integration all worked out from the question then you are sorted. However, there is another way to think about it that can help get your head around a confusing context. Suppose we draw a 2D graph with the values of Y on the y-axis and the values of X on the x-axis, then:
Each bound of the integration draws a side of the coloured rectangle that holds the area that represents the probability we are trying to find. This method also extends to when we have x as a bound and we simply draw a line at y=x, it can even do infinite bounds as we will see in the following example: