Re: Area Under Curve (Min Length Interval)
- To: mathgroup at smc.vnet.net
- Subject: [mg53901] Re: Area Under Curve (Min Length Interval)
- From: "Carl K. Woll" <carlw at u.washington.edu>
- Date: Wed, 2 Feb 2005 18:10:43 -0500 (EST)
- Organization: University of Washington
- References: <ctqdm2$sd3$1@smc.vnet.net>
- Sender: owner-wri-mathgroup at wolfram.com
Bruce,
One approach is to use a little calculus and FindRoot. Loading the
appropriate package, and defining F,
In[1]:=
<<Statistics`ContinuousDistributions`
In[2]:=
F[x_] := CDF[ChiSquareDistribution[5], x]
we find that
In[5]:=
FindRoot[{F'[b]-F'[a]==0,F[b]-F[a]==93/100},{a,1},{b,5}]
Out[5]=
{a -> 0.37253, b -> 10.3441}
Carl Woll
"Bruce Colletti" <vze269bv at verizon.net> wrote in message
news:ctqdm2$sd3$1 at smc.vnet.net...
> Re Mathematica 5.1.
>
> How would I compute the minimum length interval over which the area under
> f(x) is given?
>
> For instance, as shown below, f(x) is the PDF of a chi-square distributed
> random variable whose CDF is F[x]. Seeking the minimum length
> 93%-interval [a,b], the code returns "Obtained solution does not satisfy
> the following constraints within Tolerance -> 0.001..." Fiddling with
> options has been futile.
>
> Any ideas? Thankx.
>
> Bruce
>
> F[x_] := CDF[ChiSquareDistribution[5], x]
>
> Minimize[{b - a, F[b] - F[a] == 0.93, b > a > 0}, {a, b}]
>
> NMinimize[{b - a, F[b] - F[a] == 0.93, b > a > 0}, {a, b}]
>
>