Analysis of US Sector of Services with a New Fama-French 5-Factor Model

Show more

1. Introduction

After the Capital Asset Pricing Model (CAPM) was created by Sharpe (1964) [2] and Lintner (1965) [3] , it makes a fundamental contribution to understand the relationship between expected returns and market risks. Fama and French (1993) added size and book-to-market factors into the CAPM, many empirical results show it’s capable to explain the stock returns better than the CAPM. After that, many new factor models are developed. Panel A of Table 1 documents the development of the factor model in stock market. For example, Carhart (1997) [4] introduced a Carhart 4-factor (C) model by augumenting the Fama-French 3-factor (FF3) model with momentum factor which can explain the short-term persistence in expected returns. Chan and Faff (2005) [5] advocated a liquidity-agumented FF3 model by using Australian data and find the liquidity factor is very robust to sensitivity checks. Connor, Hagmann and Linton (2012)

Table 1. Researches about the Fama-French 5-factor Model.

Notes: “-” means that no information is available in this paper; CAPM = Capital Asset Pricing Model; FF3 = Fama and French (1993) 3-factor model; FF4 = Fama and French 4-factor model (2013); FF5 = Fama and French (2015) 5-factor model; C = Carhart (1997) 4-factor; q-factor = Hou, Xue, and Zhang (2012) q-factor model; 14-factor = Harvay and Liu (2015) 14-factor model; Mkt = Market; SMB = Size; HML = Book-to-market; WML = Momentum; IMV = liquidity; Vol = Own-volatility; RMW = Profitability; CMA = Investment; PMU = Profitable Minus Unprofitable; HML = High Minus Low; HAC-adjusted OLS = Newey-West heteroskedasticity; and autocorrelation-adjusted OLS. WLS = Weighted least squares.

[6] considered a five-factor extension of the C model which suggests an own-volatility factor.

In 2015, Fama and French proposed 5 factor model(FF5), it adds profitability and investment factors into their 3-factor model proposed in 1993. Since then, many studies about Fama-French 5-factor (FF5) model have been done. Panel B of Table 1 presents the researches for the FF5 Model. And these researches mainly apply the FF5 model to empirical stock markets and compare the FF5 model with others.

For example, Hou, Xue and Zhang (2015) [7] found that the 4-factor q-model created by Hou, Xue and Zhang (2014) [8] performs better than the FF5 model in US market. Harshita et al. (2015) [9] pointed out that the FF5 model works better in India than CAPM and FF3 model. Fama and French (2015) [9] also showed that the FF5 model can explain quite well for North America and other 3 regions. Mardy et al. (2016) [10] empirically investigated the FF5 Model in Australia, finding after adding the profitability and investment factors, FF5 model is really able to explain more asset pricing anomalies than other competing asset pricing models (like Fama-French 3-factor model and Carhart 4-factor model).

Although FF5 model has better performance in many case, it’s not adapted to every situation. Fama and French (2017) [11] analyzed the international market and found that the investment factor CMA is redundant for Europe, Japan and Asia Pacific. Meanwhile, Fama and French also found the new factors’ performance are different for small and big stock market. And for different regions, factors’ performance also exist difference. Besides, Guo et al. (2017) [12] found that the profitability factor significantly improves the description of average return, and investment pattern in average returns is weak in China stock market.

In 2017, Li et al. [13] added non-normal errors of SSAEPD proposed by Zhu and Zinde-Walsh (2009) [14] and the EGARCH-type volatilities suggested in Nelson (1991) [15] to extend the 5 factor model in Fama and French (2015). They called this new model as FF5-SSAEPD-EGARCH. Both EGARCH equation and SSAEPD can be used to capture the fat-tailedness. SSAEPD can be used to capture the asymmetric kurtosis of data. Thus, in this paper we use the data of US services industry to empirically test the new model and compare it with Fama-French 5 factors (FF5). In this paper, following two hypotheses will be tested:

1) With EGARCH-type volatilities and SSAEPD errors, are Fama-French 5 factors still alive?

2) Can this new model explain services industry better than the 5 factor model in Fama and French (2015)?

To answer these questions, we run simulation to test the validity of MatLab program used in this paper. Then, the industry of services in US are analyzed. Data are downloaded from the French’s Data Library, and the sample period is from Jul. 1990 to Feb. 2017. Method of Maximum Likelihood Estimation (MLE) is used to estimate the parameters. Likelihood Ratio test (LR) and Kolmogorov-Smirnov test (KS) are used for model diagnostics. Akaike Information Criterion(AIC) is used for model comparsion.

We find out the Fama-French 5 factors are still alive. The new model has better in-sample fit than the 5-factor model in Fama and French (2015). The industry of services can earn extra Alpha returns since the constant term in the new model is statistically significant. The Beta ( ${\beta}_{1}$ ) coefficient (for US, North American) is very close to 1. We also find out models with GARCH-typed volatility fit data better than those with EGARCH-typed volatility. To capture fat-tailedness, GARCH equation is better than non-normal error terms of SSAEPD.

The organization of this paper is as follows: The model and methodology are discussed in Section 2; Empirical results and the model comparisons will be presented in Section 3; Section 4 is the conclusions and future extensions.

2. Model and Methodology

2.1. Models

2.1.1. Fama-French 5-Factor Model (FF5-Normal)

Fama and French(2015) propose a 5-factor model (denoted as FF5) to explain market, size, value, profitability, and investment patterns in expected stock returns, and show this model empirically outperforms their 3 factor model. The 5-factor model is:

${R}_{t}-{R}_{ft}={\beta}_{0}+{\beta}_{1}\ast \left({R}_{mt}-{R}_{ft}\right)+{\beta}_{2}\ast SM{B}_{t}+{\beta}_{3}\ast HML{O}_{t}$ (1)

$+{\beta}_{4}\ast RM{W}_{t}+{\beta}_{5}\ast CM{A}_{t}+{u}_{t},\text{\hspace{0.17em}}{u}_{t}~\text{Normal}\left(\mu ,{\sigma}^{2}\right).$ (2)

where
$\theta =\left({\beta}_{0},{\beta}_{1},{\beta}_{2},{\beta}_{3},{\beta}_{4},{\beta}_{5},\mu ,\sigma \right)$ are parameters to be estimated in this model.
${R}_{t}$ is the rate of return for stock portfolio.
${R}_{ft}$ is the rate of return for the risk-free asset.
${R}_{mt}$ is the rate of return for the market.
$SM{B}_{t}$ stands for small market capitalization minus big market capitalization.
$HML{O}_{t}$ is the high book-to-market ratio minus low book-to-market ratio orthogonalized^{1}. RMW_{t} stands for robust operating profitability portfolios minus weak operating profitability portfolios.
$CM{A}_{t}$ stands for conservative investment portfolios minus aggressive investment portfolios. The error term
${u}_{t}$ is distributed as the Normal.
$\text{\hspace{0.05em}}t=1,2,\cdots ,T$ .

2.1.2. FF5-SSAEPD-EGARCH Model

Li et al. (2016) extend Fama-French(2015) five-factor model by introducing a Standardized Standard Asymmetric Exponential Power Distribution (SSAEPD) errors and the EGARCH -type volatilites. The new model we proposed is (denoted as the FF5-SSAEPD-EGARCH model):

${R}_{t}-{R}_{ft}={\beta}_{0}+{\beta}_{1}\left({R}_{mt}-{R}_{ft}\right)+{\beta}_{2}SM{B}_{t}+{\beta}_{3}HML{O}_{t}$ (5)

$+{\beta}_{4}RM{W}_{t}+{\beta}_{5}CM{A}_{t}+{u}_{t}\mathrm{,}$

${u}_{t}={\sigma}_{t}{z}_{t},\text{\hspace{0.17em}}{z}_{t}~SSAEPD\left(\alpha ,{p}_{1},{p}_{2}\right),$ (6)

$\mathrm{ln}\left({\sigma}_{t}^{2}\right)=a+{\displaystyle \underset{i=1}{\overset{s}{\sum}}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}g\left({z}_{t-i}\right)+{\displaystyle \underset{j=1}{\overset{m}{\sum}}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}{b}_{j}\mathrm{ln}\left({\sigma}_{t-j}^{2}\right),$ (7)

$\begin{array}{c}g\left({z}_{t-i}\right)={c}_{i}{z}_{t-i}+{d}_{i}\left[\left|{z}_{t-i}\right|-E\left(\left|{z}_{t-i}\right|\right)\right]\\ =(\begin{array}{ll}\left({c}_{i}+{d}_{i}\right){z}_{t-i}-{d}_{i}E\left(\left|{z}_{t-i}\right|\right)\mathrm{,}\hfill & \text{if}\text{\hspace{0.17em}}{z}_{t-i}\ge \mathrm{0,}\hfill \\ \left({c}_{i}-{d}_{i}\right){z}_{t-i}-{d}_{i}E\left(\left|{z}_{t-i}\right|\right)\mathrm{,}\hfill & \text{else}\mathrm{.}\hfill \end{array}\end{array}$ (8)

where $\theta =\left({\beta}_{0},{\beta}_{1},{\beta}_{2},{\beta}_{3},{\beta}_{4},{\beta}_{5},\alpha ,{p}_{1},{p}_{2},a,{\left\{{b}_{j}\right\}}_{j=1}^{m},{\left\{{c}_{i}\right\}}_{i=1}^{s},{\left\{{d}_{i}\right\}}_{i=1}^{s}\right)$ are the parameters to be estimated. Definitions of variables are the same as before. ${\sigma}_{t}$ is the conditional standard deviation, i.e., volatility. The error term ${z}_{t}$ is distributed as the Standardized Standard Asymmetric Exponential Power Distribution (SSAEPD) proposed in Zhu and Zinde-Walsh (2009).

•Standardized Standard AEPD (SSAEPD)

According to Zhu and Zinde-Walsh (2009), the AEPD density has following form^{2}:

${f}_{\text{AEPD}}\left(x\right)=\{\begin{array}{l}\left(\frac{\alpha}{{\alpha}^{*}}\right)\frac{1}{\sigma}K\left({p}_{1}\right)\mathrm{exp}\left(-\frac{1}{{p}_{1}}{\left|\frac{x-\mu}{2{\alpha}^{*}\sigma}\right|}^{{p}_{1}}\right),\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{if}\text{\hspace{0.17em}}x\le \mu ,\hfill \\ \left(\frac{1-\alpha}{1-{\alpha}^{*}}\right)\frac{1}{\sigma}K\left({p}_{2}\right)\mathrm{exp}\left(-\frac{1}{{p}_{2}}{\left|\frac{x-\mu}{2\left(1-{\alpha}^{*}\right)\sigma}\right|}^{{p}_{2}}\right),\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{if}\text{\hspace{0.17em}}x>\mu .\hfill \end{array}$ (10)

${f}_{\text{AEPD}}\left(x\right)=\{\begin{array}{l}\frac{1}{\sigma}\mathrm{exp}\left(-\frac{1}{{p}_{1}}{\left|\frac{x-\mu}{2\alpha \sigma K\left({p}_{1}\right)}\right|}^{{p}_{1}}\right),\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{if}\text{\hspace{0.17em}}x\le \mu ,\hfill \\ \frac{1}{\sigma}\mathrm{exp}\left(-\frac{1}{{p}_{2}}{\left|\frac{x-\mu}{2\left(1-\alpha \right)\sigma K\left({p}_{2}\right)}\right|}^{{p}_{2}}\right),\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{if}\text{\hspace{0.17em}}x>\mu .\hfill \end{array}$ (9)

$\theta =\left(\alpha ,{p}_{1},{p}_{2},\mu ,\sigma \right)$

where
$\theta =\left(\alpha \mathrm{,}{p}_{1}\mathrm{,}{p}_{2}\mathrm{,}\mu \mathrm{,}\sigma \right)$ is the parameter vector.
$\mu \in R$ and
$\sigma \mathrm{>0}$ represent location and scale, respectively^{3}.
$\alpha \in \left(\mathrm{0,1}\right)$ is the skewness parameter.
${p}_{1}>0$ and
${p}_{2}>0$ are the left and the right tail parameters, respectively.
$K\left(p\right)$ and
${\alpha}^{\mathrm{*}}$ are defined as

$K\left(p\right)=\frac{1}{2{p}^{1/p}\Gamma \left(1+1/p\right)}\mathrm{,}$ (11)

${\alpha}^{*}=\frac{\alpha K\left({p}_{1}\right)}{\alpha K\left({p}_{1}\right)+\left(1-\alpha \right)K\left({p}_{2}\right)}.$ (12)

If we set the location parameter
$\mu =0$ and the scale parameter
$\sigma =1$ , then we say X is a random variable distributed as Standard AEPD, denote it as
$X~\text{SAEPD}\left(\alpha \mathrm{,}{p}_{1}\mathrm{,}{p}_{2}\mathrm{,0,1}\right)$ . Its PDF^{4}, mean and variance are

${f}_{\text{SAEPD}}\left(x\right)=\{\begin{array}{l}\left(\frac{\alpha}{{\alpha}^{*}}\right)K\left({p}_{1}\right)\mathrm{exp}\left(-\frac{1}{{p}_{1}}{\left|\frac{x}{2{\alpha}^{*}}\right|}^{{p}_{1}}\right),\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{if}\text{\hspace{0.17em}}x\le 0,\hfill \\ \left(\frac{1-\alpha}{1-{\alpha}^{*}}\right)K\left({p}_{2}\right)\mathrm{exp}\left(-\frac{1}{{p}_{2}}{\left|\frac{x}{2\left(1-{\alpha}^{*}\right)}\right|}^{{p}_{2}}\right),\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{if}\text{\hspace{0.17em}}x>0,\hfill \end{array}$ (14)

$E\left(X\right)=\frac{1}{B}\left[{\left(1-\alpha \right)}^{2}\frac{{p}_{2}\Gamma \left(2/{p}_{2}\right)}{{\Gamma}^{2}\left(1/{p}_{2}\right)}-{\alpha}^{2}\frac{{p}_{1}\Gamma \left(2/{p}_{1}\right)}{{\Gamma}^{2}\left(1/{p}_{1}\right)}\right]\mathrm{,}$ (15)

$\begin{array}{c}Var\left(X\right)=\frac{1}{{B}^{2}}\{{\left(1-\alpha \right)}^{3}\frac{{p}_{2}^{2}\Gamma \left(3/{p}_{2}\right)}{{\Gamma}^{3}\left(1/{p}_{2}\right)}+{\alpha}^{2}\frac{{p}_{1}^{2}\Gamma \left(3/{p}_{1}\right)}{{\Gamma}^{3}\left(1/{p}_{1}\right)}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}-{\left[\left(1-\alpha \right)\frac{{p}_{2}\Gamma \left(2/{p}_{2}\right)}{{\Gamma}^{2}\left(1/{p}_{2}\right)}-{\alpha}^{2}\frac{{p}_{1}\Gamma \left(2/{p}_{1}\right)}{{\Gamma}^{2}\left(1/{p}_{1}\right)}\right]}^{2}\}.\end{array}$ (16)

Then, if we standardize X with its mean and standard deviation, we can get

$\theta =\left(\alpha ,{p}_{1},{p}_{2},0,1\right)$

$Z=\frac{X-E\left(X\right)}{\sqrt{Var\left(X\right)}}$ , which we call Standardized Standard AEPD (SSAEPD). The

PDF of Z can be got by transformation.

${f}_{\text{SSAEPD}}\left(Z\right)=\left|J\right|{f}_{\text{SAEPD}}\left(E\left(X\right)+Z\sqrt{Var\left(X\right)}\right)$ (17)

$=\delta {f}_{\text{SAEPD}}\left(\omega +Z\delta \right)$ (18)

where $\omega =E\left(X\right)$ , $\left|J\right|=\delta $ and $\delta =\sqrt{Var\left(X\right)}$ , we can get the probability density function (PDF) of the SSAEPD

${f}_{\text{SSAEPD}}\left(z\right)=\{\begin{array}{l}\delta \left(\frac{\alpha}{{\alpha}^{*}}\right)K\left({p}_{1}\right)\mathrm{exp}\left(-\frac{1}{{p}_{1}}{\left|\frac{\omega +z\delta}{2{\alpha}^{*}}\right|}^{{p}_{1}}\right),\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{if}\text{\hspace{0.17em}}z\le -\frac{\omega}{\delta},\hfill \\ \delta \left(\frac{1-\alpha}{1-{\alpha}^{*}}\right)K\left({p}_{2}\right)\mathrm{exp}\left(-\frac{1}{{p}_{2}}{\left|\frac{\omega +z\delta}{2\left(1-{\alpha}^{*}\right)}\right|}^{{p}_{2}}\right),\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{if}\text{\hspace{0.17em}}z>-\frac{\omega}{\delta}.\hfill \end{array}$ (19)

$E\left(z\right)=0,Var\left(z\right)=1.$ With $\alpha =0.5$ , ${p}_{1}={p}_{2}=2$ , SSAEPD reduces to Normal (0,1).

2.2. Method of Maximum Likelihood Estimation (MLE)

We estimate the FF5-SSAEPD-EGARCH model with the method of Maximum Likelihood Estimation (MLE). The likelihood function is

$L\left({\left\{{R}_{t}-{R}_{ft},{R}_{mt}-{R}_{ft},SM{B}_{t},HML{O}_{t},RM{W}_{t},CM{A}_{t}\right\}}_{t=1}^{T};\theta \right)$ (20)

$={\displaystyle \underset{t=1}{\overset{T}{\prod}}}f\left({R}_{t}-{R}_{ft}\right)$

$={\displaystyle \underset{i=1}{\overset{n}{\prod}}}\{\begin{array}{ll}\frac{\delta}{\eta}\left(\frac{\alpha}{{\alpha}^{\ast}}\right)K\left({p}_{1}\right)\mathrm{exp}\left(-\frac{1}{{p}_{1}}{\left|\frac{\omega +\delta {z}_{t}}{2{\alpha}^{\ast}}\right|}^{{p}_{1}}\right),\hfill & {z}_{t}\le -\frac{\omega}{\delta},\hfill \\ \frac{\delta}{\eta}\left(\frac{1-\alpha}{1-{\alpha}^{\ast}}\right)K\left({p}_{2}\right)\mathrm{exp}\left(-\frac{1}{{p}_{2}}{\left|\frac{\omega +\delta {z}_{t}}{2\left(1-{\alpha}^{\ast}\right)}\right|}^{{p}_{2}}\right),\hfill & {z}_{t}>-\frac{\omega}{\delta}.\hfill \end{array}$ (21)

where

${z}_{t}=\frac{{R}_{t}-{R}_{ft}-{\beta}_{0}-{\beta}_{1}\left({R}_{mt}-{R}_{ft}\right)-{\beta}_{2}SM{B}_{t}-{\beta}_{3}HML{O}_{t}-{\beta}_{4}RM{W}_{t}-{\beta}_{5}CM{A}_{t}}{{\sigma}_{t}},$ (22)

$\mathrm{ln}\left({\sigma}_{t}^{2}\right)=a+{\displaystyle \underset{i=1}{\overset{s}{\sum}}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}g\left({z}_{t-i}\right)+{\displaystyle \underset{j=1}{\overset{m}{\sum}}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}{b}_{j}\mathrm{ln}\left({\sigma}_{t-j}^{2}\right),$ (23)

$\begin{array}{c}g\left({z}_{t-i}\right)=c{z}_{t-i}+{d}_{i}\left[\left|{z}_{t-i}\right|-E\left(\left|{z}_{t-i}\right|\right)\right],\\ =(\begin{array}{ll}\left({c}_{i}+{d}_{i}\right){z}_{t-i}-{d}_{i}E\left(\left|{z}_{t-i}\right|\right),\hfill & \text{if}\text{\hspace{0.17em}}{z}_{t-i}\ge 0,\hfill \\ \left({c}_{i}-{d}_{i}\right){z}_{t-i}-{d}_{i}E\left(\left|{z}_{t-i}\right|\right),\hfill & \text{else}\text{.}\hfill \end{array}\end{array}$ (24)

3. Empirical Analysis

3.1. Data

In this paper, the sector of services in US is analyzed. Monthly return and 3 types of 5 factors (US 5 factors, North American 5 factors and Global 5 factors)^{5} are downloaded from French’s Data Library^{6}. Sample period is from July 1990 to Feb. 2017. 3 types of 5 factors (US, north American, global) are compared.

North American 5 factors are construted by Canada and United States.

US 5 factors are constructed by United States.

Table 2 lists the descriptive statistics calculated by Matlab^{7}. The values of skewness are not equal to 0 and those of Kurtosis are not 3. Especially, kurtosis values are all greater than 3. P-values of JB tests are 0, which are smaller than

Table 2. Descriptive statistics (1990:7 to 2017:2, monthly).

Notes: Med. = Median; Max. = Maximum; Min. = Minimum; St.De. = Standard Devistion; Ske. = Skewness; Kur. = Kurtosis; P = englishP-value of Jarque-Bera Test; ME = Market Excess Return; SMB = Small minus Big; HML = High minus Low; RMW = Robust minus Weak; CMA = Conservation minus Aggressive; The null hypothesis of JB test is H_{0}: Data are distributed as Normal(0,1).

0.05. That means, under 5% significance level, we can reject the null hypothesis and conclude that data do not follow Normal distribution. Hence, non-Normal error of SSAEPD may be proper. And from Table 2 we can see US 5 factors are very similar to North American 5 factors.

3.2. Estimation Results

The estimates are listed in Table 3. For FF5-SSAEPD-EGARCH, the Alpha returns for Global five factors is 0.78, bigger than the ones calculated from both US five factors and North American five factors ( 0.35 and 0.27, respectively). And the values of Beta ( ${\beta}_{1}$ coefficient) for US five factors is close to that from the North American five factors, which is very close to 1. Meanwhile, the value of Beta ( ${\beta}_{1}$ ) for Global five factors is 0.81, which is the smallest. It is interesting to find the coefficient of ${\beta}_{2}$ is negative for global five factors, which means the small-size effect documented in US market can not be found in the global market. Similar conclusions can be found from model of FF5-Normal.

Table 3. Estimates.

Notes: FF5-Normal is the model used in Fama-French (2015); FF5-SSAEPD-EGARCH is the new 5-factor model suggested by Zhu and Li (2016) supposing the error term meet the EGARCH-type volatilities and SSAEPD errors.

3.2.1. Fama-French 5 Factors Still Alive

• Parameter Restriction Tests

Likelihood Ratio test (LR)^{8} is used to test the significance of regressors in these models. The P-values for Likelihood Ratio tests are listed in Table 4. We find out with non-Normal errors such as SSAEPD and EGARCH-type volatilities, the Fama-French 5 factors are still alive for the sector of services.

Panel A of Table 4 lists the test results for the FF5-SSAEPD-EGARCH model. For example, the P-values of the joint significance test (see column T1)^{9} for all 3 types of 5 factors (US, North Ameirican and Global) are approximately equal to 0, which means the coefficient of
${\beta}_{1}\mathrm{,}{\beta}_{2}\mathrm{,}{\beta}_{3},{\beta}_{4}$ and
${\beta}_{5}$ are statistically joint significant under 5% significance level.

For 3 kinds of five factors, the individual significance tests show ${\beta}_{1}$ is statistically significant (see column T3). That is, market returns have significant effect on this sector returns. Same is true for ${\beta}_{3}\mathrm{,}{\beta}_{4}$ and ${\beta}_{5}$ . For 2/3 types of 5 factors, ${\beta}_{0}$ and ${\beta}_{2}$ are statistically significant (see column T2 and T4, respectively).

Panel B of Table 4 lists the test results for the FF5-Normal model. For this model, this sector doesn’t have a statistically significant coefficient
${\beta}_{0}$ under 5% significance level (see column T2) which means they can not earn statistically signicant Alpha returns. But with FF5-SSAEPD-EGARCH model, this sector in both US and Global market have a statistically significant coefficient
${\beta}_{0}$ under 5% significance level, expecially we can earn more Alpha return from Global market because
${\beta}_{0}$ in this market is 0.78 (see Table 3, column 2), the highest among 3 types of 5 factors. Furthermore, the size factor seems statistically significant for 2/3 kinds of 5 factors and is not significant in Global Market^{10}.

Table 4. P-values of likelihood ratio test (LR).

For the FF5-SSAEPD-EGARCH model, among both US 5 factors and Global 5 factors, all individual coefficients in the mean equation are statistically significant. Hence, we conclude that Fama-French 5 factors are alive even if EGARCH and SSAEPD considered.

In this part, some restrictions on the parameters in the EGARCH equation are also tested with Likelihood Ratio test (LR). And the results are also listed in Table 4. Results show the EGARCH-type volatility should be included in Fama-French 5 factor model. For instance, we do the joint significance test for hypothesis ${H}_{0}:b=c=d=0$ . The P-value of the LR are all smaller than the significance level 5%, which means our EGARCH-type volatilities is necessary, english ARCH and GARCH terms should be added into Fama-French 5-factor model since they are all statistically signicant (see column T13).

Next, we test the parameters in the SSAEPD with same method and the results show englishparameter α is statistically significant equal to 0.5, so skewnessenglish is not documentedenglish. Non-Normality is confirmed (see column T8, T10, T11, T12). The left and the right tail parametersenglish ( ${p}_{1}$ and ${p}_{2}$ ) are jointly statistically different from 2 (see column T10). And leftenglish tail is statistically different from 2 in all markets (see column T11) but right tail is only statistically

different from 2 in North American market (see column T12). i.e., strong left fat-tailedness is documented. Therefore. this new 5-factor model can capture the fat-tailedness better than FF5-Normal model.

• Kolmogorov-Smirnov Test for Residuals

We check the residuals for models with Kolmogorov-Smirnov test (KS). The P-values of KS test are listed in Table 5, the P-value of the Global five factors is 0.07, greater than 5%, which means under 5% significance level, the null hypothesis is not rejected and the residuals from FF5-SSAEPD-EGARCH do follow the SSAEPD. 2/3 markets support this result. For the FF5-Normal model^{11}, the P-values of the KS test are also listed in Table 5. All of them have smaller P-values than 0.05, which means reject the nulls. Hence, the residuals of the FF5-Normal model don’t follow Normal distribution. And the FF5-Normal model is not adequate for the data.

3.2.2. Model Comparison

We compare models with AIC. Results in Table 6 show FF5-SSAEPD-EGARCH model has smaller AIC value than FF5-Normal model. That is, this new model is better than the one in Fama and French(2015). However, FF5-GARCH model or FF5-SSAEPD-GARCH model seem to be the best because it has the smallest AIC values. That means, models with GARCH is better than the ones with EGARCH. Also, since AIC values are the same for both FF5-GARCH and FF5-SSAEPD-GARCH, we conclude that, to capture fat-tailedness, GARCH equation is better than SSAEPD, which is consistent with what we found out in our previous researches.

4. Conclusions and Future Extensions

In this paper, US sector of services is studied. A new Fama-French 5-factor model (denoted as FF5-SSAEPD-EGARCH) is empirically tested. This new model uses the non-normal error term of SSAEPD of Zhu and Zinde-Walsh

Table 5. P-values of KS test.

Note: The null hypothesis of KS test is ${H}_{0}$ : Data follow a specified distribution. We set the significance level of all tests at 5%. If the P-value of KS test is bigger than 5%, then we do not reject the null hypothesis. Otherwise, we reject the null hypothesis. For example, We apply KS test for the FF5-SSAEPD-EGARCH model residuals with the null hypothesis of ${H}_{0}$ : FF5-SSEAPD-EGARCH model residuals are distributed as $\text{SSAEPD}\left(\stackrel{^}{\alpha}\mathrm{,}{\stackrel{^}{p}}_{1}\mathrm{,}{\stackrel{^}{p}}_{2}\right)$ . For Global five factors, its P-value is 0.07, which is bigger than 0.05. That means, under 5% significance level, we cannot reject the null hypothesis and conclude that the residuals from FF5-SSEAPD-EGARCH model follow SSAEPD.

Table 6. AIC values.

Note: North Ame.=North American.

(2009) and EGARCH-type volatility of Nelson (1991) to extend the 5 factor model of Fama and French (2015). The return of services industry and 3 types of 5 factors (US five factors, North American five factors, Global five factors) from French’s Data Library are analyzed and compared. Sample period is from Jul. 1990 to Feb. 2017. Likelihood Ratio test (LR) is used for parameter restriction test, Kolmogorov-Smirnov test (KS) for residual check and AIC for model comparison. Maximum Likelihood Estimation method (MLE) is used to estimate models via MatLab.

Empirical results show: 1) With EGARCH-typed volatilities and non-normal errors, the Fama-French 5 factors are still alive; 2) The new model fits the data better than Fama-French (2015)’s 5-factor model; 3) Models with GARCH-typed volatility are a little bit better than the ones with EGARCH-typed volatility; 4) To capture fat-tailedness, GARCH equation is better than SSAEPD; 5) Using SSAEPD, model can capture stronger left fat-tailedness.

Future extensions will include but not limited to the following: First, we can construct a new index for services industry; Second, other sectors can be analyzed; Last, different factors can be considered.

NOTES

^{1}
$HML{O}_{t}$ is the sum of the intercept and the residual from the regression of
$HM{L}_{t}$ on
${R}_{mt}-{R}_{ft}\mathrm{,}SM{B}_{t}\mathrm{,}RM{W}_{t}\mathrm{,}CM{A}_{t}$ . The reason we use
$HML{O}_{t}$ is that Fama and French (2015) show
$HM{L}_{t}$ (the high book-to-market ratio minus low book-to-market ratio) is redundant in following 5-factor model.

${R}_{t}-{R}_{ft}=$

${\beta}_{0}+{\beta}_{1}\ast \left({R}_{mt}-{R}_{ft}\right)+{\beta}_{2}\ast SM{B}_{t}+{\beta}_{3}\ast HM{L}_{t}+{\beta}_{4}\ast RM{W}_{t}+{\beta}_{5}\ast CM{A}_{t}+{u}_{t},$ (3)

${u}_{t}~\text{Normal}\left(\mu \mathrm{,}{\sigma}^{2}\right)\mathrm{,}t=\mathrm{1,2,}\cdots \mathrm{,}T\mathrm{.}$ (4)

^{2}A convenient reparametrization of Equation (10) is obtained by rescaling, where

^{3}In this case,
$\mu $ and
$\sigma $ are not the notations for the population mean and the population standard deviation.

^{4}A convenient reparametrization of Equation (14) is obtained by rescaling, where

${f}_{\text{SAEPD}}\left(x\right)=\{\begin{array}{l}\mathrm{exp}\left(-\frac{1}{{p}_{1}}{\left|\frac{x}{2\alpha K\left({p}_{1}\right)}\right|}^{{p}_{1}}\right),\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{if}\text{\hspace{0.17em}}x\le 0,\hfill \\ \mathrm{exp}\left(-\frac{1}{{p}_{2}}{\left|\frac{x}{2\left(1-\alpha \right)K\left({p}_{2}\right)}\right|}^{{p}_{2}}\right),\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{if}\text{\hspace{0.17em}}x>0.\hfill \end{array}$(13)

^{5}Global 5 factors are constructed by 23 countries, there are Australia, Austria, Belgium, Canada, Switzerland, Germany, Denmark, Spain, Finland, France, Great Britain, Greece, Hong Kong , Ireland, Italy, Japan, Netherlands, Norway, New Zealand, Portugal, Sweden, Singapore, United States.

^{6}Data source is mba.tuck.dartmouth.edu/pages/faculty/ken.french/data_library.html.

^{7}Excess returns are got by portfolio returns minus the risk free rate.

^{8}LR formula is from Neyman and Pearson (1993) [22] . The equation is:

$\text{LR}=-2\mathrm{ln}\left(\text{likelihoodfornull}\right)+2\mathrm{ln}\left(\text{likelihoodforalternative}\right)$ .

^{9}The null hypothesis of the joint significance test is
${H}_{0}:{\beta}_{1}={\beta}_{2}={\beta}_{3}={\beta}_{4}={\beta}_{5}=0$ .

^{10}See column T4.

^{11}The null hypothesis
${H}_{0}$ : FF5-Normal residuals are distributed as
$\text{Normal}\left(\stackrel{^}{\mu}\mathrm{,}{\stackrel{^}{\sigma}}^{2}\right)$ .

References

[1] Fama, E.F. and French, K.R. (2015a) A Five-Factor Asset Pricing Model. Journal of Financial Economics, 116, 1-22.

[2] Sharpe, W.F. (1964) Capital Asset Prices: A Theory of Market Equilibrium under Conditions of Risk. Journal of Finance, 19, 425-442.

[3] Lintner, J. (1965) The Valuation of Risk Assets and the Selection of Risky Investments in Stock Portfolios and Capital Budgets. Review of Economics and Statistics, 47, 13-37.

[4] Carhart, M.M. (1997) On Persistence in Mutual Fund Performance. Journal of Finance, 52, 57-82.

[5] Chan, H.W. and Faff, R.W. (2005) Asset Pricing and the Illiquidity Premium. Financial Review, 40, 429-458.

[6] Gregory, C., Hagmann, M. and Linton, O. (2012) Efficient Semiparametric Estimation of the Fama-French Model and Extensions. Journal of the Econometric, 80, 713-754.

[7] Hou, K., Xue, C. and Zhang, L. (2015) Digesting Anomalies: An Investment Approach. Review of Financial Studies, 28, 650-705.
https://doi.org/10.1093/rfs/hhu068

[8] Hou, K, Xue, C. and Zhang, L. (2014) A Comparison of New Factor Models. Social Science Electronic Publishing, Rochester, NY.

[9] Singh, S. and Yadav, S.S. (2015) Indian Stock Market and the Asset Pricing Models. Procedia Economics & Finance, 30, 294-304.
https://doi.org/10.1016/S2212-5671(15)01297-6

[10] Chiah, M., Chai, D., Zhong, A. and Li, S. (2016) A Better Model? An Empirical Investigation of the Fama-French Five-Factor Model in Australia. International Review of Finance, 16, 595-638.

[11] Fama, E.F. and French, K.R. (2017) International Tests of a Five-Factor Asset Pricing Model. Journal of Financial Economics, 123, 441-463.
https://doi.org/10.1016/j.jfineco.2016.11.004

[12] Guo, B., Zhang, W., Zhang, Y.J. and Zhang, H. (2017) The Five-Factor Asset Pricing Model Tests for the Chinese Stock Market. Pacific-Basin Finance Journal, 43, 84-106.

https://doi.org/10.1016/j.pacfin.2017.02.001

[13] Li, L., Zhu, Q.Y. and Yang, Y.J. (2017) An Extension of FAMA and French 5-Factor Model Based on the SSAEPD Errors and the EGARCH-Type Volatilities. Asian Academic Research Association, Journal of Multidisciplinary (AARJMD), 4, 184-203.

[14] Zhu, D. and Zinde-Walsh, V. (2009) Properties and Estimation of Asymmetric Exponential Power Distribution. Journal of Econometrics, 148, 86-99.

https://doi.org/10.1016/j.jeconom.2008.09.038

[15] Nelson, D.B. (1991) Conditional Heteroskedasticity in Asset Returns: A New Approach. Econometrica, 59, 347-370. https://doi.org/10.2307/2938260

[16] Fama, E.F. and French. K.R. (1992) The Cross-Section of Expected Stock Returns. Journal of Finance, 47, 427-465. https://doi.org/10.1111/j.1540-6261.1992.tb04398.x

[17] Griffin, J.M. (2002) Are the Fama and French Factors Global or Country Specific? Review of Financial Studies, 15, 783-803. https://doi.org/10.1093/rfs/15.3.783

[18] Chan, H.W. and Faff, R.W. (2003) An Investigation into the Role of Liquidity in Asset Pricing: Australian Evidence. Pacific-Basin Finance Journal, 11, 555-572.
https://doi.org/10.1016/S0927-538X(03)00003-9

[19] Fama, E.F. and French, K.R. (2012) Size, Value, and Momentum in International Stock Returns. Journal of Financial Economics, 105, 457-472.
https://doi.org/10.1016/j.jfineco.2012.05.011

[20] Chai, D., Faff, R. and Gharghori, P. (2013) Liquidity in Asset Pricing: New Australian Evidence Using Low-Frequency Data. Australian Journal of Management, 38, 375-400.

https://doi.org/10.1177/0312896213489143

[21] Fama, E.F. and French, K.R. (2013) A Four-Factor Model for the Size, Value, and Profitability Patterns in Stock Returns. SSRN Electronic Journal.

[22] Neyman, J. and Pearson, E.S. (1993) On the Problem of the Most Efficient Tests of Statistical Hypotheses. Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences, 231, 289-337.
https://doi.org/10.1098/rsta.1933.0009