3
\title{Mode of the Stable Distribution Function}
5
Computes the mode of the stable distribution function.
7
Finds the maximum of \code{\link{dstable}} numerically, using
11
stableMode(alpha, beta,
13
tol = .Machine$double.eps^0.25)
16
\item{alpha, beta}{numeric parameters:
17
value of the index parameter \code{alpha} in the range \eqn{(0,2]},
18
and the skewness parameter \code{beta}, in the range \eqn{[-1, 1]}.}
19
\item{beta.max}{for numerical purposes, values of beta too close to 1,
20
are set to \code{beta.max}. Do not modify unless you know what
22
\item{tol}{numerical tolerance for \code{\link{optimize}()}.}
25
returns a numeric value, the location of the stable mode.
29
Diethelm Wuertz for the Rmetrics \R-port;
30
minor cleanup by Martin Maechler.
34
Chambers J.M., Mallows, C.L. and Stuck, B.W. (1976);
35
\emph{A Method for Simulating Stable Random Variables},
36
J. Amer. Statist. Assoc. 71, 340--344.
39
\emph{Stable Distributions},
40
Preprint, University Washington DC, 30 pages.
43
\emph{Numerical Calculation of Stable Densities and Distribution
45
Preprint, University Washington DC, 16 pages.
47
Samoridnitsky G., Taqqu M.S. (1994);
48
\emph{Stable Non-Gaussian Random Processes, Stochastic Models
49
with Infinite Variance},
50
Chapman and Hall, New York, 632 pages.
52
Weron, A., Weron R. (1999);
53
\emph{Computer Simulation of Levy alpha-Stable Variables and
55
Preprint Technical Univeristy of Wroclaw, 13 pages.
60
stableMode(alpha=1, beta=0)
61
stableMode(alpha=1.2, beta=0.1)
63
\keyword{distribution}