MathGroup Archive 2011

[Date Index] [Thread Index] [Author Index]

Search the Archive

Re: TransformedDistribution -- odd problem

Thanks -- that fixes it!

Bonus question: if we don't specify v2>2, why doesn't Mathematica return a two part solution:
 k v2 / (v2-2) v2>2

 Indeterminate True

That's what it does if we request the mean of F -- why doesn't it do the same if we request the mean of k*F.

From: Barrie Stokes <Barrie.Stokes at>
To: mathgroup at; paulvonhippel at yahoo <paulvonhippel at>
Sent: Monday, July 25, 2011 8:09 PM
Subject: [mg120497] Re: TransformedDistribution -- odd problem

Hi Paul

There are conditions on the v1 and v2, the degrees of freedom of the F distribution:

Assuming[v2 > 2, 
Mean[TransformedDistribution[F , 
  F \[Distributed] FRatioDistribution[v1, v2]]]]

{Assuming[v2 > 2, 
 Mean[TransformedDistribution[k*F , 
  F \[Distributed] FRatioDistribution[v1, v2]]]], k*v2/(-2 + v2)}

{Assuming[v2 > 2, 
  Mean[TransformedDistribution[k + F , 
   F \[Distributed] FRatioDistribution[v1, v2]]]], 
 k + v2/(-2 + v2)} // FullSimplify

which shows precisely what you expect for k*F and k+F.



>>> On 25/07/2011 at 9:29 pm, in message <201107251129.HAA25540 at>,
paulvonhippel at yahoo <paulvonhippel at> wrote:
> A little more experimenting shows that the TransformedDistribution
> function will also not provide the mean of k+F where k is a constant
> and F has an F distribution -- i.e.,
>  Mean[TransformedDistribution[k*F ,  F \[Distributed]
> FRatioDistribution[v, v]]]
> If I changce the distribution of F to NormalDistribution or
> ChiSquareDistribution, I can get a mean for k*F or k+F. So the problem
> only occurs when I define a simple function of an F variable using the
> TransformedDistribution function.
> This all strikes me as very strange, and I'd be curious to know if
> others can reproduce my results. If you can't reproduce my results,
> I'd be interested in theories about why my results differ from yours.
> E.g., is there a setting I should change in the software?
> I am using version and looking to upgrade to 8.0.1, if that
> makes a difference.
> Many thanks for any pointers.
> On Jul 24, 2:22 am, paulvonhippel at yahoo <paulvonhip... at>
> wrote:
>> I'm having a very strange problem with TransformedDistribution, where
>> I can calculate the mean of an F distribution but I cannot calculate
>> the mean of a constant multiplied by an F distribution. That is, if I
>> type
>>  Mean[TransformedDistribution[F, F \[Distributed]
>> FRatioDistribution[v, v]]]
>> Mathematica gives me an answer. But if I type
>>  Mean[TransformedDistribution[k*F ,  F \[Distributed]
>> FRatioDistribution[v, v]]]
>> Mathematica just echoes the input. I swear I got an answer for the
>> second expression earlier today. What am I doing wrong?

  • Prev by Date: Roots of a Jacobi polynomial
  • Next by Date: Re: And now for something completely different
  • Previous by thread: Re: TransformedDistribution -- odd problem
  • Next by thread: Re: TransformedDistribution -- odd problem