The Math Forum

Search All of the Math Forum:

Views expressed in these public forums are not endorsed by NCTM or The Math Forum.

Math Forum » Discussions » Inactive » comp.soft-sys.math.mathematica

Notice: We are no longer accepting new posts, but the forums will continue to be readable.

Topic: MarcumQ and Speeding-up Computations
Replies: 0  

Advanced Search

Back to Topic List Back to Topic List

Posts: 5
Registered: 5/4/13
MarcumQ and Speeding-up Computations
Posted: Jul 16, 2013 5:55 AM
  Click to see the message monospaced in plain text Plain Text   Click to reply to this topic Reply

Re 9.0.1 under Windows 7.

Here I find the expected value of a MarcumQ function over a 3D-domain whose joint probability density function is f.

How can I more quickly compute the answer 0.930445?

I think the domain causes the long computation time (if I had only one computation, this would be okay--but I have many). I had hoped simulation would give a good estimate but it didn't. I can't figure out how to use Compile with Expectation.

My hunch says that any increase in speed will come from using some bizarre fact that applies the MarcumQ function to chi-square tests. I've looked for such results but have come up dry.

Any ideas to speed up Mathematica computations? Thanks.


f = ProbabilityDistribution[Piecewise@{{6, x1 + x2 + x3 <= 1}},
{x1, 0, 1}, {x2, 0, 1}, {x3, 0, 1}];

q = {0.1, 0.2, 0.3, 0.4};
t = Array[s, 3];
t = Append[t, 1 - Total@t];

MarcumQ[3/2, Sqrt[ 50 Total@(((q - t)^2)/q)], Sqrt[7.8]],
Most@t \[Distributed] f]

{555.613162, 0.930445}

Point your RSS reader here for a feed of the latest messages in this topic.

[Privacy Policy] [Terms of Use]

© The Math Forum at NCTM 1994-2018. All Rights Reserved.