# Expectation of CauchyDistribution

I've noticed this strange behavior and I'm wondering if it's a bug.

I define a Cauchy distribution:

c = CauchyDistribution[0, 1];


If I evaluate Mean[c], I get Indeterminate, as expected.

If I evaluate Expectation[x, x \[Distributed] c], I get Expectation[x, x \[Distributed] CauchyDistribution[0, 1]]. I would have expected Indeterminate, too, but that's ok.

However, if I evaluate Expectation[x + y, {x \[Distributed] c, y \[Distributed] c}], I get 0.

The sum of two independent Cauchy distributed random variables should be another Cauchy distributed random variable, right? Why is the expectation 0?

-

You can check that :

TransformedDistribution[ x + y, {x \[Distributed] CauchyDistribution[0, 1],
y \[Distributed] CauchyDistribution[0, 1]}]

(* CauchyDistribution[0, 2] *)

Expectation[z, z \[Distributed] CauchyDistribution[0, 2]]
(* Expectation[z, z \[Distributed] CauchyDistribution[0, 2]] *)

Mean[TransformedDistribution[ x + y, {x \[Distributed] c, y \[Distributed] c}]]
(* Indeterminate *)

-
So, the conclusion is it's a bug? – Sjoerd C. de Vries Jan 11 '13 at 12:19
@SjoerdC.deVries Yes it looks like it. – b.gatessucks Jan 11 '13 at 12:42
@AndyRoss can you confirm this diagnosis? – rcollyer Jan 11 '13 at 15:48
@rcollyer I was hoping for news from him too. – b.gatessucks Jan 11 '13 at 16:31