G OODNESS - OF - FIT T EST - 相依結構對多資產選擇權定價影響之實證分析

As we obtain the marginal distributions of the underlying asset returns and

combine them into the joint distribution with copula, next we suppose to test whether

the specific distribution for a random variable accurately fits the corresponding

observations. We propose to use Kolmogorov-Smirnov (K-S) test and Anderson-Darling

(A-D) test as our measurement, which belong to the large class of goodness-of-fit tests.

They test the null hypothesis that the populations drawn from two or more

independent samples of data are identical. This test therefore can be used to decide

which dependence structure best fit the empirical data, namely the null hypothesis H

₀

of

same population distributions cannot be rejected. Then we can compare the P-value

with the desired significance level to facilitate a decision about this null hypothesis.

Let x be a set of sample observations. According to the Glivenko–Cantelli theorem,

the empirical cumulative distribution F_E will converge to the hypothetical cumulative

distribution F_H almost surely. The test statistics of K-S test and A-D test are described

as follows:

D_KS = max

t F_E x_t − F_H(x_t) (3.13) D_AD = _F^[F^E^{x −F}^H^{x ]}²

H x [1−F_H x ]

x dF_H x (3.14) The empirical cumulative distribution F_E is defined as

F_E(x) =¹_T ^T_t=1 ^N_j=1I(x_j,t≤ x_j^u^j^∙T) (3.15)

where I(∙) is the indicator function, u is a vector of marginal probabilities, and x_j^u^j^∙T is the u_j∙ T th order statistic.

The K-S test is distribution free in the sense that the critical values do not depend

on the specific distribution being tested, and it is more sensitive to deviations in the

center of the distribution. The A-D test introduced by Anderson and Darling (1952,

1954) is a modification of the K-S test that gives more weight to deviations in the tails

than K-S test. By using these two goodness-of-fit tests simultaneously, we could have

more evidence to prove our model selection.

15 3.5 MONTE CARLO SIMULATION

Here we propose to apply a conditional sampling method to simulate draws from a

specific bivariate copula. Given the parameters of each copula function, our task is to

generate pairs of observations (𝑢, 𝑣) of the random vector (𝑈, 𝑉), where 𝑈, 𝑉 are

both independent uniform distribution and their joint distribution function is C. The

main idea of conditional sampling is based on the conditional distribution 𝑐_𝑢 𝑣 = Pr(𝐹₂ ≤ 𝑣 𝐹₁ = 𝑢) = lim

Δ𝑢→∞

𝐶 𝑢+Δ𝑢,𝑣 −𝐶(𝑢,𝑣)

Δ𝑢 =^𝜕𝐶_𝜕𝑢 = 𝐶_𝑢 𝑣 (3.16) where 𝐶_𝑢(𝑣) is the partial derivative of the copula. The conditional distribution 𝑐_𝑢 𝑣

is a non-decreasing function and exists for all 𝑣 ∈ [0,1].

The general procedure to simulate the desired return process is stated as follows:

(i) Generate two independent random variables (𝑢, 𝑢′) ∈ Uniform(0,1), where 𝑢

is the first desired draw.

(ii) Compute the quasi-inverse function of 𝑐_𝑢 𝑣 , which depend on the parameters

of the selected copula and on 𝑢. Set 𝑣 = 𝑐_𝑢⁻¹(𝑢′) to obtain the second desired

draw. The following Table 2 reports the quasi-inverse functions of copulas.

(iii) Let 𝐹_𝑖 ∙ be the cumulative density function of the standardized marginal

innovations z_i,t~skewed − t z_i η_i, ϕ_i . We take the inverse functions 𝑆𝑇_1,𝑡, 𝑆𝑇_2,𝑡 = 𝐹⁻¹ 𝑢₁ , 𝐹⁻¹ 𝑢₂ as the observation of innovations at time 𝑡.

Then calculate the daily asset returns 𝑟_1,𝑡, 𝑟_2,𝑡 by the following formula:

𝑟_𝑖,𝑡 = 𝑟_𝑓 −¹₂ℎ_{𝑖,𝑡−1}+ 𝑆𝑇_𝑖,𝑡, for 𝑖 ∈ 1,2 , 𝑡 ∈ 1, … , 𝑇 (3.17)

Table 2: Quasi-inverse function of copula function 𝒄_𝑢(𝒗)

Gaussian Φ 1 − 𝜌²Φ⁻¹ 𝑣 + 𝜌Φ⁻¹ 𝑢

Stud𝑒𝑛𝑡’s t

t_𝜈 𝜌 t_𝜈⁻¹ 𝑢 + 𝜈 + t_𝜈⁻¹ 𝑢 ² × 1 − 𝜌²

𝜈 − 1 t_{𝜈 +1}⁻¹ (𝑣)

Gumbel No closed-form solution

Clayton 𝑣 = u^−∝ u^′−^∝+1^∝ − 1 + 1

−1

∝

Frank 𝑣 = −1

∝ln 1 + u^′(1 − e^−∝) u^′ e^−∝u− 1 − e^−∝u

Note: The inverse function of 𝑐_𝑢(𝑣) for Gumbel copula is calculated by solving an equation by numerical way. The details are discussed in the book Copula Methods in Finance, written by Cherubini, Luciano, and Vecchiato, 2004.

(iv) After we have all daily asset returns (𝑟_1,𝑡, 𝑟_2,𝑡) for 𝑡 = {1, … , 𝑇}, we obtain the

percentage growth (𝑅_1,𝑡, 𝑅_2,𝑡) of underlying assets defined as:

𝑅_i,𝑡 = exp ^𝑡_{𝑗 =1}𝑟_𝑖,𝑗 , for 𝑖 ∈ 1,2 , 𝑡 ∈ 1, … , 𝑇 (3.18)

Finally, we can calculate the bivariate option prices under different payoff functions

with underlying asset return (𝑅_1,𝑇, 𝑅_2,𝑇) and count it in one simulation sample.

4. D

ATASETS AND

O

PTION

P

AYOFF

S

TRUCTURES

4.1DATA AND DIAGNOSTIC ANALYSIS

We use several types of indices as proxy for the equity, commodity and foreign

exchange returns with different regions or industries: MSCI Pacific Index and MSCI Far

East Index (cross regions), MSCI World Energy Index and MSCI World Bank Index

(cross industries), Reuters/Jefferies CRB Index and MSCI Latin America Index (cross

assets), S&P 500 Index and DXY currency Index (cross assets). The linear correlations

of each pair of underlying financial assets are shown in Table 3. We would like to

compare the impact of different types of dependence structure to the bivariate options

pricing.

Table 3: The linear correlation of each pair assets

Linear Correlation Pair Composition

Pair 1 0.9828 MSCI Pacific Index and MSCI Far East Index

Pair 2 0.6330 MSCI World Energy Index and MSCI World Bank Index Pair 3 0.4160 Reuters/Jefferies CRB Index and MSCI Latin America Index Pair 4 0.0317 S&P 500 Index and DXY currency Index

All data are collected from Bloomberg over the period from July 1, 2002 to

December 31, 2008. The analysis is based on the logarithm returns of daily closing

prices which exclude non-trading days to ensure that enough observations of the tail of

18 the distributions are available.

The following Table 4 reports the summary statistics on the returns in the samples.

All of them have apparently skewness and leptokurtotic, which implies the

unconditional distributions of asset returns are asymmetric and fat-tailed. Moreover, the

Jarque-Bera test also shows obviously non-Gaussian characteristic at the 5%

significance level for all underlying assets.

Table 4: Summary Statistics Asset

Statistics MSCI Pacific

MSCI Far East

MSCI

Energy MSCI Bank CRB MSCI Latin

America S&P 500 DXY Mean× 10³ 0.0273 0.0154 0.1149 -0.1114 0.0245 0.2477 -0.0270 -0.0750 Std. Dev. 0.0059 0.0062 0.0071 0.0059 0.0048 0.0086 0.0060 0.0023 Skewness -0.4306 -0.2449 -0.6977 -0.0771 -0.3822 -0.5251 -0.1874 -0.0164 Kurtosis 9.2944 8.3391 16.2382 16.1062 7.7685 14.3414 14.5425 4.7099 J-B 2854.33* 2032.39* 12294.54* 11919.17* 1577.92* 9027.30* 8980.17* 196.21*

Q(20) 33.40* 40.31* 96.83* 136.82* 60.52* 76.20* 120.07* 28.11 Q²(20) 2809.89* 1942.96* 3225.62* 2411.35* 1989.51* 3194.73* 2954.04* 587.34*

ARCH(5) 615.74* 452.75* 589.36* 303.39* 318.75* 658.89* 465.85* 104.25*

Note: The sample period for the daily logarithm returns runs from July 1, 2002 to December 31, 2008 excluding the holidays. J-B represents the Jarque-Bera test for normality; Q(20) is the Ljung-Box lack-of-fit hypothesis test for up to the 20^th order serial correlation in the returns; Q²(20) is the Ljung-Box test for the serial correlation in the squared returns; and ARCH(5) is the LM test for up to the fifth-order ARCH effects. *indicates significance at the 5% level.

We also implement the Ljung-Box lack-of-fit hypothesis test, computing the

Q-statistic for autocorrelation lags 20 at the 5% significance level. There are significant

serial correlations in the returns and squared returns for most of samples, but the DXY

Index have no serial correlation in the returns. Moreover, the Engle's ARCH test for up

to the fifth-order shows significant evidence in support of the GARCH effects, or

namely the heteroscedasticity. These statistics support the preliminary hypothesis of

skewed-t GARCH model settings that we discussed before.

4.2OPTION PAYOFF STRUCTURES

We select three types of exotic bivariate options commonly seen in practice, the

digital options, the spread options, and the rainbow options, as our proxy to examine the impact of each dependence structure. Let X_i be the logarithm return at maturity of asset

i ∈ {1,2}, K is the pre-determined strike price, and I ∙ is the indicator function, which

equals 1 if the statement in parentheses holds and zero otherwise. The payoff structures

are shown in Table 5.

Table 5: Bivariate options payoff structures

Option Payoff Structure

Digital Options I(X₁> 𝐾 , X₂> 𝐾)

Spread Options max⁡(X₁− X₂− K , 0)

Rainbow Options max max⁡(X₁ , X₂ − K, 0]

5. E

MPIRICAL

R

ESULTS

We implement the skewed-t GARCH model to the above dataset and derive the

empirical copula of the standardized GARCH innovations displayed in Figure 1. The

empirical copula would converge to the true copula function under the regularity

condition; for instance, Gaenssler and Stute (1987) or van der Vaart and Wallner (1996).

Figure 1: The empirical copula of the standardized GARCH innovations Pair 1: MSCI Pacific vs MSCI Far East Pair 2: MSCI Energy vs MSCI Bank

Pair 3: CRB vs MSCI Latin America Pair 4: S&P500 vs DXY

Note: The skewed-t GARCH model is described as follows:

r_i,t= μ_i+ ε_i,t,

21 hi,t= ωi+ βihi,t−1+ αiε_i,t−1² ,

εi,t|ℱt −1 = hi,tzi,t , zi,t~skewed − t zi η_i, ϕi

Then the standardized GARCH innovations for each pair of underlying assets are εi,t hi,t, for 𝑖 ∈ 1,2 and 𝑡 ∈ 1, … , 𝑇 .

We observe that there exists strong positive dependence between the innovations of

first two pairs of assets, particularly in the tails. On the other hand, the innovations of

last two pairs of underlying assets are quite irrelevant and scattered all over the figure.

The dependence structures are dissimilar among these four pairs of assets.

The following Table 6 reports the maximum likelihood estimation results of the

copula-based GARCH model. The estimates of parameters for marginal skewed-t

GARCH process are presented along with the P-values in parentheses in the panel A,

where 0.0000 indicates that the P-value is less than 0.00005.

The value of ω, α and β are all nonnegative that ensure the conditional variance h_i,t is positive, and the sum of α and β is close to one for all pairs of assets, which

implies the shocks between the underlying assets have high persistence in volatility and

stationary process. For the estimation of skewed-t distribution, the kurtosis parameter η

and the asymmetry parameter φ show that the conditional process of all assets is

fat-tailed and asymmetric, and this is consistence with previous studies.

The panel B of Table 6 reports the estimates of parameters for different types of

copula functions. All of these estimates are significant under 5% level, except for the

degrees of freedom ν of the last three pairs of Student’s t copulas.

Table 6: Copula-based GARCH Model Estimations Asset

Panel A: Estimates of marginal process

ω × 10⁵ 0.0179 0.0187 0.0115 0.0029 0.0100 0.0231 0.0030 0.0025

Panel B: Estimates of dependence process

Gaussian ρ 0.9842 0.5146 0.2801 0.0777

Note: The values in parentheses are P-values, where 0.0000 indicates that the value is less than 0.00005.

We compare the model goodness-of-fit through the Kolmogorov-Smirnov (K-S)

test and Anderson-Darling (A-D) test along with the P-values in parentheses in Panel A

and B of Table 7. The log-likelihood estimator and Akaike Information Criterion (AIC)

of the above copula estimates are also reported in the Panel C and D of Table 7. The

lower values of the K-S test statistics, A-D test statistics and AIC information criterion

represent better goodness-of-fit. On the contrary, the larger of log-likelihood estimates

indicate better fit.

According to both of the K-S test and A-D test, the first pair of assets, which

components are MSCI Pacific and MSCI Far East Index, has the largest P-value that

means the best model fit under the Gumbel copula function. The Frank copula is the

most ideal model for the second pair of assets, which components are MSCI Energy and

MSCI Bank Index. Although the P-values of these five dependence structures are not

significant enough for the third and fourth pair of assets, the Frank copula is still the

model that relatively describes the dependence structure best based on the empirical

distribution under the K-S and A-D test statistics.

Nevertheless, the inferred conclusion from log-likelihood estimator and AIC

information criterion indicates that the Student’s t is the best dependence structure for

all of these four pairs of asset returns, which obviously contradicted to our above

induction by the K-S and A-D tests. We still choose to apply the K-S and A-D test as our

basis to the following bivariate option pricing comparison.

Table 7: Goodness-of-fit test to Copula-based GARCH Models Asset Panel A: Kolmogorov-Smirnov (K-S) test

Gaussian 0.0316 0.0597 0.0886 0.1297

Gaussian 1.4626 9.0281 19.2397 43.9457

(0.1865) (0.0000) (0.0000) (0.0000)

Student’s t 1.4293 8.3982 18.8446 43.9345

(0.1951) (0.0000) (0.0000) (0.0000)

Gumbel 1.2187 10.7238 21.6726 44.4410

(0.2599) (0.0000) (0.0000) (0.0000)

Clayton 4.9203 15.9543 26.2062 46.3411

(0.0029) (0.0000) (0.0000) (0.0000)

Frank 2.3892 8.0226 17.7837 42.2574

(0.0566) (0.0000) (0.0000) (0.0000)

Panel C: Maximum Likelihood Estimation

Gaussian 2958.28 258.44 67.11 4.89

Student’s t 2999.23 268.76 68.08 13.49

Gumbel 2940.42 238.37 52.00 8.17

Clayton 2322.66 187.30 52.07 5.52

Frank 2854.21 244.85 62.38 5.30

Panel D: Akaike Information Criterion (AIC)

Gaussian -5914.57 -514.89 -132.23 -7.77

Student’s t -5996.46 -535.52 -134.16 -24.98

Gumbel -5878.84 -474.74 -102.00 -14.33

Clayton -4643.33 -372.59 -102.14 -9.03

Frank -5706.42 -487.69 -122.76 -8.60

Note: The test statistics of K-S test and A-D test are described as follows:

DKS= max

t F_E x_t − F_H(xt) D_AD = _F^[F^E^{x −F}^H^{x ]}²

H x [1−F_H x ]

x dF_H x .

The copula estimates in Table 5.1 are implemented in the hypothetical cumulative distribution F_H and the empirical cumulative distribution F_E is defined as

F_E(x) =¹

T ^T_t=1 ^N_j=1I(x_j,t≤ x_j^u^j^∙T)

where I(∙) is the indicator function, u is a vector of marginal probabilities, and x_j^u^j^∙T is the u_j∙ T th order statistic. Besides, the values in parentheses are P-values, where 0.0000 indicates that the value is less than 0.00005.

After we derive all desired parameter estimates, we simulate the asset return

process that given the parameters of skewed-t GARCH model and copula functions.

Especially we apply the parameter estimates of each marginal skew-t GARCH process

with initial conditional volatility h_i,0 equals the long-term variance level ω_i 1 − α_i− β_i . The estimates of Gaussian, Student’s t, Gumbel, Clayton and Frank

copula functions are also implemented to simulate different types of dependence

structure. To ensure the convergence of simulated option prices, we simulate 10,000

times and assume that the initial underlying asset return equals one for each return

process.

The following Figure 2 reports the simulation results of digital option pricing of

MSCI Pacific Index and MSCI Far East Index with different strike prices K from 0.80 to

1.20 and different time-to-maturity T from one to six months. We also fix the strike

price K to divide them into three groups of in-the-money (ITM), at-the-money (ATM)

and out-of-the-money (OTM) respectively in Figure 3.

The differences of digital option price between each type of dependence structure

tend to become wider when the maturity increases and the initial bivariate option prices

are around at-the-money. However, the tendency of price change is not very obvious,

and the same situation happens both for the spread options and rainbow options

displayed in Figure 4 to 6.

In order to confirm whether the differences among these bivariate option prices are

significant or not, we perform a paired t-test of the null hypothesis that the 10,000 times

simulated asset returns in the difference between each pairs of dependence structures are

a random sample from a normal distribution with mean zero and unknown variance,

against the alternative that the mean is not zero.

Figure 1:Digital Option Pricing Result Under Different Copulas vs. Different Strike Prices - Pair 1- MSCI Pacific Index and MSCI Far East Index

A: 1 Month Maturity B: 2 Month Maturity

C: 3 Month Maturity D: 4 Month Maturity

E: 5 Month Maturity F: 6 Month Maturity

Note: The payoff function of digital option is I(X₁> 𝐾 , X₂> 𝐾), where I(∙) is the indicator function that equals to one if the statement in parentheses holds and zero otherwise. We assume that the number of trading days for a month is twenty days.

0.00

0.80 0.84 0.88 0.92 0.96 1.00 1.04 1.08 1.12 1.16 1.20

0.00

0.80 0.84 0.88 0.92 0.96 1.00 1.04 1.08 1.12 1.16 1.20

0.00

0.80 0.84 0.88 0.92 0.96 1.00 1.04 1.08 1.12 1.16 1.20

0.00

0.80 0.84 0.88 0.92 0.96 1.00 1.04 1.08 1.12 1.16 1.20

0.00

0.80 0.84 0.88 0.92 0.96 1.00 1.04 1.08 1.12 1.16 1.20

0.00

0.80 0.84 0.88 0.92 0.96 1.00 1.04 1.08 1.12 1.16 1.20

Figure 2: Digital Option Price Result Under Different Copulas vs. Different Time to Maturities - Pair 1- MSCI Pacific Index and MSCI Far East Index

ITM (𝐾 = 0.9) ATM (𝐾 = 1) OTM (𝐾 = 1.1)

Note: ITM, ATM and OTM represent that the digital option is in-the-money, at-the-money and out-of-the-money respectively. Since both the underlying asset returns are set to start at the same value 1 in the simulation process, ATM indicates that 𝐾 equals to 1, ITM indicates that K is less than 1 and OTM indicates that K is larger than 1.

Figure 3: Spread Option Pricing Result Under Different Copulas vs. Different Strike Prices - Pair 1- MSCI Pacific Index and MSCI Far East Index

A: 1 Month Maturity B: 2 Month Maturity

C: 3 Month Maturity D: 4 Month Maturity

0.8

0.80 0.84 0.88 0.92 0.96 1.00 1.04 1.08 1.12 1.16 1.20

0.00

0.80 0.84 0.88 0.92 0.96 1.00 1.04 1.08 1.12 1.16 1.20

0.00

0.80 0.84 0.88 0.92 0.96 1.00 1.04 1.08 1.12 1.16 1.20

0.00

0.80 0.84 0.88 0.92 0.96 1.00 1.04 1.08 1.12 1.16 1.20

E: 5 Month Maturity F: 6 Month Maturity

Note: The payoff function of spread option is denoted by max⁡(X₁− X₂− K , 0). We assume that the number of trading days for a month is twenty days.

Figure 4: Spread Option Price Result Under Different Copulas vs. Different Time to Maturities - Pair 1- MSCI Pacific Index and MSCI Far East Index

ITM (𝐾 = −0.1) ATM (𝐾 = 0) OTM (𝐾 = 0.1)

Figure 5: Rainbow Option Pricing Result Under Different Copulas vs. Different Strike Prices - Pair 1- MSCI Pacific Index and MSCI Far East Index

A: 1 Month Maturity B: 2 Month Maturity

0.00

0.80 0.84 0.88 0.92 0.96 1.00 1.04 1.08 1.12 1.16 1.20

0.00

0.80 0.84 0.88 0.92 0.96 1.00 1.04 1.08 1.12 1.16 1.20

0.09

0.80 0.84 0.88 0.92 0.96 1.00 1.04 1.08 1.12 1.16 1.20

0.00

0.80 0.84 0.88 0.92 0.96 1.00 1.04 1.08 1.12 1.16 1.20

C: 3 Month Maturity D: 4 Month Maturity

E: 5 Month Maturity F: 6 Month Maturity

Note: The payoff function of rainbow option is max max⁡(X₁ , X₂ − K, 0]. We assume that the number of trading days for a month is twenty days.

Figure 6: Rainbow Option Price Result Under Different Copulas vs. Different Time to Maturities - Pair 1- MSCI Pacific Index and MSCI Far East Index

ITM (𝐾 = 0.9) ATM (𝐾 = 1) OTM (𝐾 = 1.1)

0.80 0.84 0.88 0.92 0.96 1.00 1.04 1.08 1.12 1.16 1.20

0.00

0.80 0.84 0.88 0.92 0.96 1.00 1.04 1.08 1.12 1.16 1.20

0.00

0.80 0.84 0.88 0.92 0.96 1.00 1.04 1.08 1.12 1.16 1.20

0.00

0.80 0.84 0.88 0.92 0.96 1.00 1.04 1.08 1.12 1.16 1.20

0.10

dependence structure between the returns of MSCI Pacific and MSCI Far East Index as

our basis model for the comparison.

In the scenario of one month time-to-maturity in Table 8, the difference between

each pairs of dependence structures are significant under 5% level, except that the

digital and rainbow option prices between Gumbel copula and Gaussian copula are not.

However, the price differences among each pairs of dependence structure become more

significant except for the rainbow option price between Gumbel copula and Frank

copula in the scenario of six months time-to-maturity. The results of paired t-test

support our previous observation that the differences between the bivariate options

prices are more significant as the time-to-maturity increases.

Table 8: T test for the distance between the ATM (K=1) option prices of different dependence structures with 1 month time-to-maturity - Pair 1- MSCI Pacific Index and MSCI Far East Index

Panel 1: Digital Option

Student’s t Clayton Gumbel Frank

Gaussian -0.0009 0.0267 0.0014 0.0138

(0.0833) (0.0000) (0.1221) (0.0000)

Student’s t 0.0276 0.0023 0.0147

(0.0000) (0.0116) (0.0000)

Clayton -0.0253 -0.0129

(0.0000) (0.0000)

Gumbel 0.0124

(0.0000)

32 Panel 2: Spread Option

Student t Clayton Gumbel Frank

Gaussian 0.0001 -0.0018 -0.0003 -0.0009

(0.0000) (0.0000) (0.0000) (0.0000)

Student’s t -0.0019 -0.0003 -0.0009

(0.0000) (0.0000) (0.0000)

Clayton 0.0016 0.0010

(0.0000) (0.0000)

Gumbel -0.0006

(0.0000) Panel 3: Rainbow Option

Student t Clayton Gumbel Frank

Gaussian 0.0000 -0.0014 0.0000 -0.0002

(0.0000) (0.0000) (0.4577) (0.0000)

Student’s t -0.0014 0.0000 -0.0002

(0.0000) (0.0001) (0.0000)

Clayton 0.0014 0.0012

(0.0000) (0.0000)

Gumbel -0.0002

(0.0000)

Note: The table reports the distance between the ATM (K=1) bivariate option prices of different dependence structures with 1 months time-to-maturity for the first pair of assets returns, which components are MSCI Pacific and MSCI Far East Index. The values in parentheses are p-values under the paired t tests.

Table 9: T test for the distance between the ATM (K=1) option prices of different dependence structures with 6 months time-to-maturity - Pair 1- MSCI Pacific Index and MSCI Far East Index

Panel 1: Digital Option

Student’s t Clayton Gumbel Frank

Gaussian -0.0014 0.0352 0.0048 0.0333

(0.0163) (0.0000) (0.0000) (0.0000)

Student’s t 0.0366 0.0062 0.0347

(0.0000) (0.0000) (0.0000)

Clayton -0.0304 -0.0019

(0.0000) (0.4016)

Gumbel 0.0285

(0.0000) Panel 2: Spread Option

Student t Clayton Gumbel Frank

Gaussian 0.0002 -0.0051 -0.0007 -0.0019

(0.0000) (0.0000) (0.0000) (0.0000)

Student’s t -0.0053 -0.0009 -0.0021

(0.0000) (0.0000) (0.0000)

Clayton 0.0044 0.0031

(0.0000) (0.0000)

Gumbel -0.0012

(0.0000) Panel 3: Rainbow Option

Student t Clayton Gumbel Frank

Gaussian 0.0001 -0.0035 -0.0001 -0.0002

(0.0000) (0.0000) (0.0000) (0.0009)

Student’s t -0.0036 -0.0002 -0.0003

(0.0000) (0.0000) (0.0000)

Clayton 0.0033 0.0033

(0.0000) (0.0000)

Gumbel -0.0001

(0.2549) Note: The table reports the distance between the ATM (K=1) bivariate option prices of different dependence structures with 6 months time-to-maturity for the first pair of assets returns, which components are MSCI Pacific and MSCI Far East Index. The values in parentheses are p-values under the paired t tests.

We also compare the options price distance with different strike prices at maturity

in Table 10 and 11. The significance of option price distance decreases when the strike

price is set to be in-the-money or out-of-the-money, expect for the rainbow options. The

differences of rainbow option prices among these copulas are almost all significant

under 5% level, which represents that the option payoff functions also determine the

importance of dependence structure of underlying asset returns.

Hence, we can prove the hypothesis that the differences of bivariate option prices

under distinct dependence structures are more significant when the option is

at-the-money. Besides, the option payoff function is also a key factor to affect the

importance of dependence structure selection.

Table 10: T test for the distance between the ITM (K=0.9) option prices of different dependence structures with 6 months time-to-maturity - Pair 1- MSCI Pacific Index and MSCI Far East Index

Panel 1: Digital Option

Student’s t Clayton Gumbel Frank

Gaussian 0.0001 0.0039 0.0014 0.0092

(0.7815) (0.0002) (0.0433) (0.0000)

Student’s t 0.0038 0.0013 0.0091

(0.0005) (0.0741) (0.0000)

Clayton -0.0025 0.0053

(0.0433) (0.0003)

Gumbel 0.0078

(0.0000) Panel 2: Spread Option

Student t Clayton Gumbel Frank

Gaussian 0.0000 0.0001 -0.0001 0.0032

(0.3348) (0.7802) (0.1778) (0.0000)

Student’s t 0.0001 -0.0001 0.0032

(0.7284) (0.3229) (0.0000)

Clayton -0.0001 0.0032

(0.5916) (0.0000)

Gumbel 0.0033

(0.0000) Panel 3: Rainbow Option

Student t Clayton Gumbel Frank

Gaussian 0.0002 -0.0050 -0.0006 -0.0015

(0.0000) (0.0000) (0.0000) (0.0000)

Student’s t -0.0052 -0.0008 -0.0017

(0.0000) (0.0000) (0.0000)

Clayton 0.0044 0.0035

(0.0000) (0.0000)

Gumbel -0.0010

(0.0000) Note: The table reports the distance between the ITM (K=0.9) bivariate option prices of different dependence structures with 6 months time-to-maturity for the first pair of assets returns, which components are MSCI Pacific and MSCI Far East Index. The values in parentheses are p-values under the paired t tests.

Table 11: T test for the distance between the OTM (K=1.1) option prices of different dependence structures with 6 months time-to-maturity - Pair 1- MSCI Pacific Index and MSCI Far East Index

Panel 1: Digital Option

Student’s t Clayton Gumbel Frank

Gaussian -0.0005 0.0151 -0.0009 0.0116

(0.0588) (0.0000) (0.0495) (0.0000)

Student’s t 0.0156 -0.0004 0.0121

(0.0000) (0.3173) (0.0000)

Clayton -0.0160 -0.0035

(0.0000) (0.0001)

Gumbel 0.0125

(0.0000) Panel 2: Spread Option

Student t Clayton Gumbel Frank

Gaussian 0.0000 -0.0001 0.0000 0.0000

(1.0000) (0.0010) (1.0000) (0.0971)

Student’s t -0.0001 0.0000 0.0000

(0.0010) (1.0000) (0.0971)

Clayton 0.0001 0.0001

(0.0010) (0.0032)

Gumbel 0.0000

(0.0971) Panel 3: Rainbow Option

Student t Clayton Gumbel Frank

Gaussian 0.0000 -0.0005 0.0000 0.0001

(0.0008) (0.0000) (0.3102) (0.0000)

Student’s t -0.0005 0.0000 0.0001

(0.0000) (0.6609) (0.0001)

Clayton 0.0005 0.0007

(0.0000) (0.0000)

Gumbel 0.0001

(0.0001) Note: The table reports the distance between the OTM (K=1.1) bivariate option prices of different dependence structures with 6 months time-to-maturity for the first pair of assets returns, which components are MSCI Pacific and MSCI Far East Index. The values in parentheses are p-values under the paired t tests.

The simulation results for the last three pairs of underlying assets are shown in the

Appendix A to C. The Frank copula is the best fit model based on our previous K-D test

and A-D test, which has neither lower nor upper tail dependency. After we compare the

distance of simulated bivariate option prices between Frank copula and other four types

of copula functions by paired t-test, the results show that the option price distinguishes

between copulas are not very significant in Table 12. Especially for the fourth pair of

asset returns that composed by S&P 500 Index and DXY currency Index, the P-values

from the paired t tests in parentheses are the highest in these four pairs of assets.

From the linear correlation estimates and the scatter plotted in Figure 1, we can

observe that the dependence is weakest for the fourth pair. Accordingly, the strength of

dependency to the underlying assets also determines the impact to bivariate option

prices. As the dependency is weaker, the distance between option prices are less

significant.

Table 12: T test for the distance between the ATM (K=1) option prices of different dependence structures

在文檔中相依結構對多資產選擇權定價影響之實證分析 (頁 20-0)