Wolfgang Hitzl wrote:
>
> Hallo Mathematica Freunde,
>
> ich hab Mathematica 3.0.1, im Package
> Statistics`ContinuousDistributions`
> ist die cumulative density function (CDF) der nichtzentralen
> F-Verteilung zur Freude jedes Statistikers gegeben.
>
> Allerdings:
> Wenn ich den Wert der CDF der nichtzentralen FVerteilung mit 1 und 2000
> Freiheitsgraden berechnen will an der Stelle 3.88, mit
>
> CDF[NoncentralFRatioDistribution[1,2000,1],3.88]
>
> dann erhalte ich die Fehlermeldung:
>
> NIntegrate::"ploss":
> "Numerical integration stopping due to loss of precision. Achieved \
>
> neither the requested PrecisionGoal nor AccuracyGoal; suspect highly \
> oscillatory integrand, or the true value of the integral is 0. If your \
>
> integrand is oscillatory try using the option Method->Oscillatory in \
> NIntegrate."
>
> Wer hat einen konstruktiven Vorschlag, dass ich den Wert dennoch
> berechnen kann?
>
> Diese Problem tritt für einen Statistiker bei praktischen
> Fragestellungen leicht auf.
>
> Viele Grüße aus Österreich
>
> Dr. Wolfgang Hitzl
Eine Mö?glichkeit ist:
Mathematica 3.0 for Linux
Copyright 1988-97 Wolfram Research, Inc.
-- Motif graphics initialized --
In[1]:= <<Statistics`ContinuousDistributions`
In[2]:= SetOptions[NIntegrate, Compiled->False];
In[3]:= Timing[CDF[NoncentralFRatioDistribution[1,2000.,1],3.88] ]
Out[3]= {3.96 Second, 0.8322429605}
--
Rolf Mertig
Mertig Research & Consulting
http://www.mertig.com
http://www.feyncalc.com
|