5
$\begingroup$

How to analytically express cov(X,Y), when:

X=C*A/(A+B) 

and

Y=C*B/(A+B)

Here C, A and B are independent variables with normal distributions.

More specifically I would like to express cov(X,Y) using the expected means of each of the input variables C, A and B, as well as their standard deviations or variances.

A is a normal distribution with an expected mean E(A) and a variance var(A)
B is a normal distribution with an expected mean E(B) and a variance var(B)
C is a normal distribution with an expected mean E(C) and a variance var(C)

Thank you so much in advance!

$\endgroup$
  • $\begingroup$ if this is HW please add the self-study tag $\endgroup$ – Antoine Jun 9 '16 at 9:42
  • $\begingroup$ Thank you for asking but it is for my research at work. I have really searched a lot the last few days for the answer, but I can only find examples where X and Y are linear combinations of random variables like X=aA+bB+cC (with small a, b and c as constants) but I cannot find solutions where X is a function that is the product of a combination of A, B and C. $\endgroup$ – Astrid Marie Jun 9 '16 at 10:18
  • $\begingroup$ I think that would help if you could show some of those examples you refer to and then explain how it is different in your case. $\endgroup$ – Antoine Jun 9 '16 at 10:27
  • $\begingroup$ Thanks for the tip. The question ressembles a bit this stats.stackexchange.com/questions/137571/covariance-of-product , but the example in the link is in matrix notation and it is not clear for me if Xi and Xj in the example are correlated, so I cannot used this directly. $\endgroup$ – Astrid Marie Jun 9 '16 at 10:36
  • 1
    $\begingroup$ Given that this involves the ratio of Normals, convergence may be an issue $\endgroup$ – wolfies Jun 9 '16 at 11:24
3
$\begingroup$

One can find a general solution, without assuming Normality. In particular, if $A$, $B$ and $C$ are independent, and noting that the Covariance operator is the {1,1} central moment, then $\text{Cov}( \frac{C A}{A+B},\frac{C B}{A+B})$ is:

enter image description here

where I am using a developmental version of the CentralMomentToCentral function in mathStatica (alas, not in any public release yet).

The reason the solution does not converge (assuming Normal parents) is because $A+B$ is Normal, and thus $E[\frac{1}{A+B}]$ is the expectation of the inverse of a Normal variable, which is known not to converge (see, for instance, Mean and variance of inverse of a normal RV)

$\endgroup$
  • $\begingroup$ Thank you so much Wolfies! It is extremely helpful to see your answer, and I really appreciate that you took your time to work with this problem. Now I try to implement the result in my calculations and compare with the simulations :-) The mathStatica seems like a really powerful tool. $\endgroup$ – Astrid Marie Jun 10 '16 at 12:58
  • $\begingroup$ Just one little question, I see that the C in your code has a different colour than A and B. Is there any reason for this? $\endgroup$ – Astrid Marie Jun 10 '16 at 13:21

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

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