I'm trying to model the distribution of effects of mutations (let's call it s) in evolution but I'm stuck in generating the probability distribution function (pdf) for my model. So, my model is a mixture of 3 components. As far as I understood, if I assume independence between all components, my final pdf is a weighted sum of the pdf of each component.
My problem starts in defining the pdf of each component. I have 2 major problems. 1: In two of my components I have a function F() that relates some biochemical value (let's call it E) with the effect I want to model. I assume that E is distributed as a Gaussian (or a mixture of gaussians). My solution was to find the function G(s) = E (receive a value of s and convert in the value of E) and substitute that functions in the Gaussian pdf. Is this correct?
2: For this I have no solution yet. So I assume that s depend on 2 different variables (E and K). In my previous point I got the pdf for the cases in which only one is changing (assuming the changes are Gaussian). Know I want to build a component where I allow the 2 variables to change. So in the end I know the distribution of E and K (I assume they are independent) and the function H(E, K) = s. I can also get the function G(s, K), = E and J(s, E) = K. How can I get the pdf for this component? I found somewhere that I can use the method of the cdf (integrate the joint probability distribution for each variable), but the Gaussian pdf is not integrable!!