A U-statistic for the tail index of a multivariate stable random vector is given as an extension of the univariate case introduced by Fan (2006). Asymptotic normality and consistency of the proposed U-statistic for the tail index are proved theoretically. The proposed estimator is used to estimate the spectral measure. The performance of both introduced tail index and spectral measure estimators is compared with the known estimators by comprehensive simulations and real datasets.
1. Introduction
In recent years, stable distributions have received extensive use in a vast number of fields including physics, economics, finance, insurance, and telecommunications. Different sorts of data found in applications arise from heavy tailed or asymmetric distribution, where normal models are clearly inappropriate. In fact, stable distributions have theoretical underpinnings to accurately model a wide variety of processes. Stable distribution has originated with the work of Lévy [1]. There are a variety of ways to introduce a stable random vector. In the following, two definitions are proposed for a stable random vector; see Samorodnitsky and Taqqu [2].
Definition 1.
A random vector X=(X1,…,Xd)T is said to be stable in Rd if for any positive numbers A and B there are a positive number C and a vector D∈Rd such that(1)AX1+BX2=dCX+D,where X1 and X2 are independent and identical copies of X and C=Aα+Bα1/α.
Definition 2.
Let 0<α<2. Then X is a non-Gaussian α-stable random vector in Rd if there exist a finite measure Γ on the unit sphere Sd=x=x1,…,xdT∈Rd∣x,x=1 and a vector μ=μ1,…,μdT∈Rd such that(2)φXt=logEexpit,X=-∫Sdt,sα1-isgnt,stanπα2Γds+it,μ,α≠1,-∫Sdt,s1+isgnt,s2πlogt,sΓds+it,μ,α=1,where t,s=∑i=1dtisi for t=(t1,…,td)T, s=(s1,…,sd)T, i2=-1, and sgn(·) denotes the sign function. The pair Γ,μ is unique.
The parameter α, in Definitions 1 and 2, is called tail index. A random vector X is said to be a strictly α-stable random vector in Rd if μ=0 for α≠1; see Samorodnitsky and Taqqu [2]. We note that X is strictly α-stable, in the sense of Definition 1, if D=0. Throughout we assume that X is strictly α-stable and α≠1. The probability density function of a stable distribution has no closed-form expression and moments with orders greater than or equal to α are not finite for the members of this class. The two aforementioned difficulties make statistical inference about the parameters of a stable distribution hard. However, a series of contributions has permitted inference about the parameters of univariate and multivariate stable distributions. For example, in the univariate case, maximum likelihood (ML) estimation was studied first by DuMouchel (1971) and then by Nolan [3]. Although the ML approach leads to an efficient estimate for samples of large size, it involves numerical complexities. A program, called STABLE, uses a cubic spline interpolation of stable densities for this purpose; see Nolan [4]. STABLE estimates all four parameters of a stable distribution for α≥0.4. Sample quantile (SQ) technique is another approach proposed by McCulloch [5]. The results are simple and consistent estimators of all four parameters based on five sample quantiles. The empirical characteristic function (CF) is suggested by Kogon and Williams [6]. The CF and SQ methods work well but are not as efficient as the ML method. As the last approach considered here, U-statistics for the tail index and scale parameters of a univariate strictly stable distribution are introduced by Fan [7]. In multivariate case, the focus of interest is the spectral measure estimation. Among them, we refer to Nolan et al. [8], Pivato and Seco [9], Ogata [10], and Mohammadi et al. [11].
The structure of the paper is as follows. In Section 2, new estimators for the tail index and spectral measure of a strictly stable distribution are presented which is an extension of the U-statistic proposed by Fan [7] for the univariate case. A comprehensive simulation study is performed in Section 3 to compare the performance of the introduced estimators and the known estimators. Two real data sets are analyzed in this section to illustrate the performance of the proposed method.
2. New Estimators
This section consists of two subsections. Firstly, we propose an estimator for the tail index. Secondly, an estimator for the spectral measure is given.
2.1. Estimation of Tail Index
The main result of this section is given in Theorem 4, which gives U-statistic for the inverse of tail index of a strictly stable distribution. We present the main result in the light of Lemma 3 given as follows. The proofs are given in the Appendix.
Lemma 3.
Let X=X1,…,XdT be a d-dimensional strictly stable random vector. Then, VarlogX is finite, where · denotes the Euclidean norm.
Theorem 4.
Let x1,…,xn be a sequence of n observations from a d-dimensional strictly stable random vector. Then (3)Un=n2-1∑1≤i<j≤nHxi,xj,where(4)Hxi,xj=logxi+xjlog2-logxi+logxj2log2is the U-statistic for 1/α.
As it is seen, from Theorem 4, the introduced U-statistic is an unbiased estimator for 1/α. Hereafter, we write α^MU=1/Un as introduced estimator for α. Here, subscript MU indicates that α^MU is constructed based on multivariate U-statistic defined in Theorem 4. It should be noted that when the true value of α is near two, the kernel given in (4) could be less than 0.5. So, α^MU is greater than two. In this case, we set α^MU=2.
2.2. Spectral Measure Estimation
We use α^MU to estimate an m-point discrete approximation to the exact spectral measure of the form(5)Γ·=∑j=1mγjIsj·,where γj is a mass at point sj in the unit sphere Sd and Isj(·) is an indicator function at point sj; for j=1,…,m, see Byczkowski et al. [12]. To estimate Γ(·), we replace Definition 2 for a strictly d-dimensional stable random vector with(6)φXt=-∑j=1mt,sjα1-isgnt,sjGα,t,sjγj,=-∑j=1mψt,sj,αγj,where G(α,t,sj)=tan(πα/2) for α≠1 and -2/πlog|t,sj| for α=1. Define(7)Λ=ψt1,s1,αψt1,s2,α⋯ψt1,sm,αψt2,s1,αψt2,s2,α⋯ψt2,sm,α⋮⋮⋱⋮ψtm,s1,αψtm,s2,α⋯ψtm,sm,α,(8)V=-logφXt1,…,-logφXtmT,where tj=(tj1,…,tjd)T∈Sd, for j=1,…,m. Using (7) and (8), both sides of (6) are connected together through the following linear system:(9)V=Λγ,where γ=(γ1,…,γm)T. Assuming that Λ in (9) is nonsingular, then γ=Λ-1V. Hence, we estimate the vector of the masses as(10)γ^=Λ^-1V^,where(11)Λ^=ψt1,s1,α^MUψt1,s2,α^MU⋯ψt1,sm,α^MUψt2,s1,α^MUψt2,s2,α^MU⋯ψt2,sm,α^MU⋮⋮⋱⋮ψtm,s1,α^MUψtm,s2,α^MU⋯ψtm,sm,α^MU,V^=-log1n∑i=1nexpit1,xi,…,-log1n∑i=1nexpitm,xiT,in which xi is i-th vector observation in random sample of size n.
Due to the standard error of V^, we have two problems with direct use of (10). Firstly, γ^ may be complex, and secondly, its real part may be quite negative. Since Λ^ and V^ are complex while gamma is constrained to be real (and nonnegative), the Euclidean norm used by McCulloch [13] and Nolan et al. [8] must be replaced with the complex modulus to solve both problems in a novel way. For this, we use the nnls(·) library in the R package. In the next section, the estimated spectral measure γ^, based on α^MU, is shown by γ^MU. We note that another estimator of γ can be constructed by separating both of the real and imaginary parts in the structure of V^. But simulation results show that constructed estimator gives the same performance.
3. Simulation Study
This section is in three parts. Firstly, we study the performance of the proposed estimator with the known ones for estimating the tail index. Secondly, we compare the performance of the spectral measure estimator developed through the introduced tail index estimator with the known approaches. In the last subsection, we give a real data example to illustrate the efficiency of the proposed estimators.
3.1. Performance Analysis of the Tail Index Estimators
Here, we perform a simulation study to compare the performance of α^MU and four other estimators for α, including (1) α^ML, (2) α^SQ, (3) α^CF, and (4) α^MM. The first three competitors are ML, SQ, and CF estimations for the tail index, respectively. Each of three competitors is obtained as α^PROJ=1/m∑j=1mα^(uj) after projecting the d-dimensional stable random vector using uj,X. Here, m is the number of masses, uj=uj1,…,ujdT is an arbitrary unit vector, and X is the d-dimensional stable random vector. It is worth noting that the first three competitors are computed by the help of STABLE software after projecting. The fourth estimator, that is, α^MM, is the second estimator for tail index proposed by Mohammadi et al. [11]. We compare both the bias and root mean-squared error (RMSE) of estimators for 500 replications of samples of size n=500 and 5000 of a bivariate stable random vector generated by the method given in Modarres and Nolan [14]. We use two settings for discrete spectral measure with m=8 masses, including γ1=(0.1,0.2,0.3,0.4,0.1,0.2,0.3,0.4)T and γ2=(0,0.1,0.7,0.3,0.7,0.3,0.7,0.1)T. In both cases, masses are concentrated on points sj = cos(2πj-1/m),sin2πj-1/mT for j=1,…,8. In the first case that data are coming from a stable distribution with γ1, we generate tj from a uniform distribution on the unit sphere Sd. For the second case, we set tj=sj. Biases and RMSEs for α = (0.1 : 0.1 : 0.9, 0.95, 1.05, 1.1 : 0.1 : 1.9, 1.95, 2) are shown in Figures 1 and 2. As Figure 1 shows, when γ1=0.1,0.2,0.3,0.4,0.1,0.2,0.3,0.4T, we observe that α^MU is more efficient than α^SQ for n=5000. Also, it works better than α^MM in terms of RMSE (for α≤1.8). Based on Figure 2, when γ2=(0,0.1,0.7,0.3,0.7,0.3,0.7,0.1)T, we observe that α^MU is more efficient than other methods when α<1.4 and n=5000 in the sense of RMSE. Also, when α<1.7 and n=500, α^MU is more efficient than α^SQ, α^CF, and α^MM with respect to RMSE.
Biases and RMSEs of estimators when data are generated from a strictly stable distribution with discrete spectral measure γ1=0.1,0.2,0.3,0.4,0.1,0.2,0.3,0.4T. (a) Bias when n=500, (b) RMSE when n=500, (c) bias when n=5000, and (d) RMSE when n=5000.
Biases and RMSEs of estimators when data are generated from a strictly stable distribution with discrete spectral measure γ2=(0,0.1,0.7,0.3,0.7,0.3,0.7,0.1)T. (a) Bias when n=500, (b) RMSE when n=500, (c) bias when n=5000, and (d) RMSE when n=5000.
3.2. Performance Analysis of the Spectral Measure Estimators
Here, we compare the performance of the estimator for masses of spectral measure γ=(γ1,…,γm)T constructed based on U-statistic, γ^MU with the other four known estimators for the spectral measure. The competitors are three types of estimators for γ based on empirical characteristic function method: (1) γ^MLE-cf; (2) γ^SQ-cf; (3) γ^CF-cf; and (4) Mohammadi et al. [11] estimator for γ, γ^MM. For computing γ^MLE-cf, γ^SQ-cf, and γ^CF-cf, we use command mvstable.fit(x,nspectral,method1d,method2d,param) in the STABLE program, where x is data vector, nspectral is number of spectral measure masses, method1d is the method to use for estimating parameters of univariate stable distribution, that is, MLE, SQ, and CF (corresponding codes in STABLE are 1, 2, and 3, respectively), method2d is the method to use for estimating parameters of bivariate stable distribution (we set method2d=2 which corresponds to empirical characteristic function approach, cf), and param refers to kind of parameterization. Here, we set param=1 since we are using the characteristic function in (2). More information about the first three competitors is given in Robust Analysis Inc. [15]. The estimators γ^MLE-cf, γ^SQ-cf, γ^CF-cf, and γ^MM are obtained by substituting α^ML, α^SQ, α^CF, and α^MM into (7) and then solving linear system (10), respectively. Comparisons are based on the RMSE of γ^j, for j=1,…,m, which is defined as 1/N∑i=1Nγ^ij-γij2, where N is the number of iterations and γ^ij is the estimation of jth component of γ^ at ith iteration. We consider five scenarios for the structure of discrete spectral measure as follows.
We note that the first and the third scenarios above are similar to Examples 2 and 1 of Nolan et al. [8], respectively. The fourth scenario is called Triangle since corresponding density contour plot is similar to a triangle. For each of the above five scenarios, we arrange the settings of simulation as m=8, α = 1.25; 1.75; n = 2000; 5000 (n is sample size), and N=500. It should be noted that masses are located at sj=cos(2π(j-1)/m),sin(2π(j-1)/m)T, for j=1,…,m, and components of tj=(tj1,…,tjd)T are generated from a uniform distribution on the unit sphere Sd. The results of simulations are given in Figures 3–6. As it is seen, γ^MU shows better performance than γ^MM.
RMSEs of γ^ under different scenarios when α=1.25. We use the following symbol scheme: ⋄ for γ^UM, ○ for γ^MM, + for γ^ML-cf, × for γ^SQ-cf, and ∆ for γ^CF-cf.
RMSEs of γ^ under different scenarios when α=1.75. We use the following symbol scheme: ⋄ for γ^UM, ○ for γ^MM, + for γ^ML-cf, × for γ^SQ-cf, and ∆ for γ^CF-cf.
RMSEs of γ^ under different scenarios when α=1.25. We use the following symbol scheme: ⋄ for γ^UM, ○ for γ^MM, + for γ^ML-cf, × for γ^SQ-cf, and ∆ for γ^CF-cf.
RMSEs of γ^ under different scenarios when α=1.75. We use the following symbol scheme: ⋄ for γ^UM, ○ for γ^MM, + for γ^ML-cf, × for γ^SQ-cf, and ∆ for γ^CF-cf.
3.3. Real Data Analysis
Here, we give two examples. In the first example, adjusted daily log-return (in percent) for the 30 stocks at the Dow Jones index is collected between January 3, 2000, and December 31, 2004. The log-return percent of 1247 closing prices has been computed for AXP (American Express Company) and MRK (Merck & Co. Inc.) stocks after multiplying the daily log-return by 100; see Nolan [16]. The scatter plot of AXP and MRK stocks log-return percent values, X=(AXP,MRK)T, is shown in Figure 7. We use a bivariate α-stable distribution with m=12 points of masses for spectral measure addressed by sj=cos(2π(j-1)/m),sin(2π(j-1)/m)T, for j=1,…,12. We estimate the location parameter as μ^ML-cf=(-3.438E-07,-2.402E-07)T. So, a strictly α-stable distribution is fitted to the Y=(X-μ^ML-cf)T. For this, we set tj=sj, for j=1,…,12. Table 1 shows the results for modelling data through five methods. We note that estimated tail indices are α^MU=1.581, α^MM=1.734, α^ML-cf=1.618, α^SQ-cf=1.493, and α^CF-cf=1.723. As it is seen, estimated tail indices through estimators α^MU and α^ML-cf are closer together than the other estimators. In the second example, we focus on the cubic-root of the monthly average of river discharge. We choose discharge of the Odra and Wisla rivers in Poland during 1901 to 1986 (raw data are in m3/s. They are available at https://nelson.wisc.edu/sage/data-and-models/riverdata/). The scatter plot for cubic-root of Odra river discharge versus cubic-root of Wisla river discharge is shown in Figure 8. Setting m=8, sj=cos(2π(j-1)/m),sin(2π(j-1)/m)T, and tj=sj, for j=1,…,8, we obtain μ^ML-cf=(7.9837,9.9947)T. After fitting a strictly α-stable distribution to the shifted data, results for estimating spectral measure are given in Table 2. Estimated tail indices are α^MU=1.860, α^MM=1.312, α^ML-cf=1.813, α^SQ-cf=1.936, and α^CF-cf=1.962. Based on results given in Table 2, estimated masses through estimators γ^MU, γ^ML-cf, and γ^SQ-cf are closer together than the other estimators. We compare here γ^MU with γ^ML-cf and γ^SQ-cf since the latter estimators are among the best estimators for the masses as shown in the previous subsection.
Estimation results after fitting a strictly bivariate α-stable distribution to AXP and MRK stocks data.
Estimator
γ^1
γ^2
γ^3
γ^4
γ^5
γ^6
γ^7
γ^8
γ^9
γ^10
γ^11
γ^12
γ^MU
0.396
0.070
0.077
0.433
0.087
0
0.231
0.378
0
0.439
0.078
0
γ^MM
0.373
0
0
0.463
0.021
0
0.333
0.479
0
0.563
0.100
0
γ^ML-cf
0.338
0.156
0.162
0.412
0
0
0.490
0.206
0.055
0.425
0.122
0
γ^SQ-cf
0.232
0.216
0.089
0.404
0.012
0
0.480
0.115
0.177
0.254
0.156
0
γ^CF-cf
0.421
0.129
0.237
0.434
0
0
0.537
0.209
0
0.527
0.089
0
Estimation results after fitting a strictly bivariate α-stable distribution to Odra and Wisla discharge data.
Estimator
γ^1
γ^2
γ^3
γ^4
γ^5
γ^6
γ^7
γ^8
γ^MU
0.006
1.459
0.498
0
0
0
0
0
γ^MM
0
0.800
0.257
0
0.026
0.483
0.165
0
γ^ML-cf
0
1.380
0.508
0
0
0
0
0
γ^SQ-cf
0
1.490
0.532
0
0
0
0
0
γ^CF-cf
0
0.196
0.199
0
0.009
1.283
0.339
0
Scatter plot for AXP versus MRK daily log-return percent, X=(AXP,MRK)T.
Scatter plot for cubic-root of Odra and Wisla rivers discharge.
4. Conclusion
We compare the performance of the introduced U-statistic for the tail index with the well-known methods, including maximum likelihood, empirical characteristic function, sample quantile, and that introduced in Mohammadi et al. [11] through a simulation study. In the sense of root mean-squared error, it is proved that proposed tail index estimator always outperforms Mohammadi et al. [11] and SQ methods when α≤1.4. This is while ML and CF methods show better performance than the proposed estimator for large α, say α>1.4 in terms of root mean-squared error. Simulation studies for estimating the discrete spectral measure γ under five scenarios prove that estimator of γ based on introduced U-statistic shows, in terms of root mean-squared error, better performance than Mohammadi et al. [11] estimator. Analysis of two sets of real data reveals that estimator of the tail index and γ based on U-statistic shows expedient performance. As some possible future works, firstly, we aim to introduce a U-statistic for the case of a nonzero location parameter. Secondly, we look for methodology possibly based on a U-statistic, to estimate tail, masses, and location parameters simultaneously. Finally, recalling that the approach employed in this work is based on characteristic function, the discrete spectral measure using α^MU can be estimated through projection approach.
AppendixProof of Lemma 3.
We show that Elog2X<∞. Suppose d=2 and p-=PX≤1, p+=1-p-, p--=PX1<0,X2<0, p-+=PX1<0,X2>0, p+-=PX1>0,X2<0, and p++=PX1>0,X2>0. So(A.1)Elog2X=Elog2X∣X<1p-+Elog2X∣X≥1p+≤Elog2X∣X<1p-+Elog2X∣X≥1p+≤Elog2X+Elog2X∣X≥1,where X, in the above, denotes one of the components of vector X. It should be noted that inequality Elog2X∣X<1≤Elog2X∣X<1 holds irrespective of d. On the other hand, (A.2)Elog2X∣X≥1=Elog2X∣X≥1,X1<0,X2<0p--+Elog2X∣X≥1,X1<0,X2>0p-++Elog2X∣X≥1,X1>0,X2<0p+-+Elog2X∣X≥1,X1>0,X2>0p++≤Elog2X1+X2∣X≥1,X1<0,X2<0p--+Elog2-X1+X2∣X≥1,X1<0,X2>0p-++Elog2X1-X2∣X≥1,X1>0,X2<0p+-+Elog2X1+X2∣X≥1,X1>0,X2>0p++≤2Elog2X1+X2+2Elog2X1+X2.Thus, (A.3)Elog2X≤Elog2X+2Elog2X1+X2+2Elog2X1-X2.Generally, for d≥2, one can write (A.4)Elog2X≤Elog2X+∑i=0ddiElog2-∑j=1iXj+Xi+1+⋯+Xd,where we adopt this convention that ∑j=10Xj=0. Let S(α,β,σ,μ=0) stands for a univariate strictly stable random variable with tail index α, scale parameter σ, and skewness parameter β. It is well known that if X=(X1,…,Xd)T is an α-stable random vector, then any linear combination of its components such as 〈bi,X〉=∑j=1dbijXj, for bi=-1,…,-1︷i,1,…,1︷d-iT, follows a stable distribution with tail index α,(A.5)σi=∫Sdbi,sαΓds1/α,βi=1σiα∫Sdbi,sαΓds1/α,where Γ(·) is spectral measure and |x|[r]=sgn(x)|x|r; see Samorodnitsky and Taqqu [2]. It follows, from Kuruoglu [17], that if X~S(α,β,σ,0), then (A.6)eα,β,σ=Elog2X=7π2-6θ2+6logσ/cosθ+γ1-α2-π2α26α2,where θ=arctan(βtan(πα/2)). Also, (A.7)eα,βi,σi=Elog2Z=7π2-6θi2+6logσi/cosθi+γ1-α2-π2α26α2,where Z=-∑j=1iXj+Xi+1+⋯+Xd, for i=0,…,d and θi=arctan(βitan(πα/2)). Parameters σi and βi are defined in (A.5). Finally,(A.8)Elog2X≤eα,β,σ+∑i=0ddieα,βi,σi.The proof is complete since all terms on the right-hand side of (A.8) are finite.
Proof of Theorem 4.
We rewrite Definition 1 as(A.9)AX1+BX2=dCX1+D.Setting A=1, B=1, and D=0 in (A.9), it yields(A.10)X1+X2=d21/αX1.By applying log-transformation, after taking the Euclidean norm, to both sides of (A.10), we have(A.11)1α=logX1+X2-logX1log2.The right-hand side of (A.11) can be used to define a symmetric kernel of the form(A.12)HX1,X2=logX1+X2log2-logX1+logX22log2.To guarantee the asymptotic normality of the introduced U-statistics for 1/α with kernel (A.12), we need to check that EH(X1,X2)2<∞. It suffices to show that VarH(X1,X2)<∞. For this, the result of Lemma 3 shows that VarlogX1 is finite. On the other hand, from (A.10) it turns out that(A.13)VarlogX1+X2=VarlogX1=VarlogX2.We use property (A.13) to calculate variance of the right-hand side of (A.12) as (A.14)VarHX1,X2=VarlogX1log22+VarlogX12log22-Klog22,where(A.15)K=CovlogX1+X2,logX1+logX2≤2VarlogX1+X2VarlogX1.Applying property (A.13) again on the right-hand side of (A.15), we have(A.16)VarHX1,X2≤7VarlogX12log22,where we used the result of Lemma 3 to get the right-hand side of (A.16). This means that(A.17)λ=VarEHX1,X2∣X1≤VarHX1,X2<∞.Therefore, EH(X1,X2)2<∞. Now, we define U-statistic for 1/α with kernel given in (A.12) as(A.18)Un=n2-1∑1≤i<j≤nHxi,xj.By definition, given U-statistic in (A.18) is unbiased estimator for 1/α.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
LévyP.Théorie des erreurs. La loi de Gauss et les lois exceptionnelles1924524985MR1504839SamorodnitskyG.TaqquM. S.1994New York, NY, USAChapman & HallStochastic ModelingMR1280932NolanJ. P.Barndorff-NielsenO. E.MikoschT.ResnickI.Maximum likelihood estimation of stable parameters2001Boston, Mass, USABirkhöuser379400NolanJ. P.Numerical calculation of stable densities and distribution functions199713475977410.1080/15326349708807450MR1482292ZBL0899.60012McCullochJ. H.Simple consistent estimators of stable distribution parameters19861541109113610.1080/03610918608812563MR8767832-s2.0-38049084243KogonS. M.WilliamsD. B.AdlerR.FeldmanR.TaqquM.Characteristic function based estimation of stable parameters1998Boston, Mass, USABirkhäuser311338FanZ.Parameter estimation of stable distributions2006351-324525510.1080/03610920500439992MR2274047ZBL1084.620192-s2.0-32444450412NolanJ. P.PanorskaA. K.McCullochJ. H.Estimation of stable spectral measures2001349–111113112210.1016/s0895-7177(01)00119-4MR18588392-s2.0-0035943852PivatoM.SecoL.Estimating the spectral measure of a multivariate stable distribution via spherical harmonic analysis200387221924010.1016/S0047-259X(03)00052-6MR2016936ZBL1041.600192-s2.0-0345492357OgataH.Estimation for multivariate stable distributions with generalized empirical likelihood2013172224825410.1016/j.jeconom.2012.08.017MR30106152-s2.0-84872043401MohammadiM.MohammadpourA.OgataH.On estimating the tail index and the spectral measure of multivariate αα-stable distributions201578554956110.1007/s00184-014-0515-7MR33554612-s2.0-84930823558ByczkowskiT.NolanJ. P.RajputB.Approximation of multidimensional stable densities1993461133110.1006/jmva.1993.1044MR1231241ZBL0790.600202-s2.0-0001858020McCullochJ. H.Estimation of the bivariate stable spectral representation by the projection method2000161-2476210.1023/A:1008797318867ZBL0964.621082-s2.0-0003188008ModarresR.NolanJ. P.A method for simulating stable random vectors1994911119MR1278554ZBL0944.65004Robust Analysis Inc2010Software and User Manualhttp://www.robustanalysis.comNolanJ. P.Multivariate elliptically contoured stable distributions: theory and estimation20132852067208910.1007/s00180-013-0396-7MR3107292ZBL1306.651182-s2.0-84884701033KuruogluE. E.Density parameter estimation of skewed /spl alpha/-stable distributions200149102192220110.1109/78.950775MR18758192-s2.0-0035485846