Suppose we have a joint distribution $P(D,X,L) = P(D|X)P(X|L)P(L)$. Here, D and L are discrete but X is a continuous random variable. I want to compute $P(D=d)$. How do I do this numerically? The fact that $X$ is continuous confuses me. 



Marginalizing over discrete and continuous random variables
I'm assuming that your "continuous" is actually "absolutely continuous", i.e. $X$ has a density.


$$P(D=d) = E[P(D=d|X)] = \int dx \ P(D=d|X=x) f\_X(x) = \int dx \sum\_\ell \ P(D=d|X=x) f\_{X|L}(x|\ell) P(L=\ell)$$


where the sum is over the possible values of $L$, and $f\_{X|L}(x|\ell)$ is the conditional density of $X$ given $L=\ell$.


