Number of knots | 0 | 1 | 2 | 3 | 4 | 5 |
Placement points | 50% | 33.3% | 25% | 20% | 16.6% | |
66.6% | 50% | 40% | 33.30% | |||
75% | 60% | 50% | ||||
80% | 66.6% | |||||
83.3% |
In this paper, we focused on the Caputo-Hadamard fractional reaction sub-diffusion equations. By using the nonuniform L1 scheme and nonuniform Alikhanov scheme in the temporal domain, we formulated two efficient numerical schemes, where the second order difference method was used in the spatial dimension. Furthermore, we derived the stability and convergence of these proposed schemes. Remarkably, both derived numerical methods exhibited α-robustness, that is, it remained valid when α→1−. Numerical experiments were given to demonstrate the theoretical statements.
Citation: Xingyang Ye, Xiaoyue Liu, Tong Lyu, Chunxiu Liu. α-robust error analysis of two nonuniform schemes for Caputo-Hadamard fractional reaction sub-diffusion problems[J]. Electronic Research Archive, 2025, 33(1): 353-380. doi: 10.3934/era.2025018
[1] | Tahir Afzal, Muhammad Asim Afridi, Muhammad Naveed Jan . Integrating LSTM with Fama-French six factor model for predicting portfolio returns: Evidence from Shenzhen stock market China. Data Science in Finance and Economics, 2025, 5(2): 177-204. doi: 10.3934/DSFE.2025009 |
[2] | Xiaozheng Lin, Meiqing Wang, Choi-Hong Lai . A modification term for Black-Scholes model based on discrepancy calibrated with real market data. Data Science in Finance and Economics, 2021, 1(4): 313-326. doi: 10.3934/DSFE.2021017 |
[3] | Nitesha Dwarika . Asset pricing models in South Africa: A comparative of regression analysis and the Bayesian approach. Data Science in Finance and Economics, 2023, 3(1): 55-75. doi: 10.3934/DSFE.2023004 |
[4] | Yi Chen, Zhehao Huang . Measuring the effects of investor attention on China's stock returns. Data Science in Finance and Economics, 2021, 1(4): 327-344. doi: 10.3934/DSFE.2021018 |
[5] | Moses Khumalo, Hopolang Mashele, Modisane Seitshiro . Quantification of the stock market value at risk by using FIAPARCH, HYGARCH and FIGARCH models. Data Science in Finance and Economics, 2023, 3(4): 380-400. doi: 10.3934/DSFE.2023022 |
[6] | Zimei Huang, Zhenghui Li . What reflects investor sentiment? Empirical evidence from China. Data Science in Finance and Economics, 2021, 1(3): 235-252. doi: 10.3934/DSFE.2021013 |
[7] | Tsumbedzo Mashamba, Modisane Seitshiro, Isaac Takaidza . A comprehensive high pure momentum equity timing framework using the Kalman filter and ARIMA forecasting. Data Science in Finance and Economics, 2024, 4(4): 548-569. doi: 10.3934/DSFE.2024023 |
[8] | Kuo-Shing Chen . Interlinkages between Bitcoin, green financial assets, oil, and emerging stock markets. Data Science in Finance and Economics, 2024, 4(1): 160-187. doi: 10.3934/DSFE.2024006 |
[9] | Qian Shen, Yifan Zhang, Jiale Xiao, Xuhua Dong, Zifei Lin . Research of daily stock closing price prediction for new energy companies in China. Data Science in Finance and Economics, 2023, 3(1): 14-29. doi: 10.3934/DSFE.2023002 |
[10] | Katleho Makatjane, Ntebogang Moroke . Examining stylized facts and trends of FTSE/JSE TOP40: a parametric and Non-Parametric approach. Data Science in Finance and Economics, 2022, 2(3): 294-320. doi: 10.3934/DSFE.2022015 |
In this paper, we focused on the Caputo-Hadamard fractional reaction sub-diffusion equations. By using the nonuniform L1 scheme and nonuniform Alikhanov scheme in the temporal domain, we formulated two efficient numerical schemes, where the second order difference method was used in the spatial dimension. Furthermore, we derived the stability and convergence of these proposed schemes. Remarkably, both derived numerical methods exhibited α-robustness, that is, it remained valid when α→1−. Numerical experiments were given to demonstrate the theoretical statements.
Since the Capital Asset Pricing Model and beta were introduced in modern finance, there has been arguments on whether beta is adequate to explain the complex nature of systematic risk in stock market. The complexity of systematic risk is primarily due to its asymmetry and time-varying nature. The asymmetry of systematic risk, particularly downside risk, begins drawing attention of researchers in recent two decades. Most studies of downside risk follow the classic approach, employing the linear market model to estimate beta. As a result, the time-varying nature of beta is either ignored or being weighted inappropriately. Therefore, adopting appropriate data smoothing technique is crucial to preserve the true time-varying nature of beta. In this paper, two models, the cubic piecewise polynomial function (CPPF) model and the Fourier Flexible Form (FFF) model are employed to model portfolio returns in order to examine the significance of beta, downside beta and upside beta estimates. Both models take flexible approaches, yet are parsimonious, allowing beta estimates to be time-varying with appropriate weight. Innovatively, various numbers of knots and orders are applied on the CPPF model and the FFF model, respectively, to smooth the sample. Also, the Akaike Information Criteria (AIC) (Akaike, 1974) is adopted to determine the most appropriate number of knots and order for the sample. With the AIC, the best fitted estimates of beta, upside beta and downside beta for both models are generated. These estimates are sorted into portfolios to examine the risk-return relationship. Moreover, Fama-Macbeth regressions are performed to discover the significance of the estimates in a cross-sectional context.
Taking the CPPF and the FFF approach is motivated by their flexibility, with both approaches allowing the beta estimates to vary over time with heavier weight on more recent data. The CPPF approach is analogous to cubic spline approach but with no constraints of intercept columns, and the estimates at each point in time are the product of a vector of initial estimates and a piecewise polynomial matrix. For the FFF approach, Sine and Cosine functions are adopted to construct a matrix which creates a non-linear pattern bounded between −1 and 1. The pattern is finally presented on the estimates at each point in time to allow time-variation. The importance of time-varying estimates is that estimates can present the true relationship between variables at each point of time, which allows us to discover the variation of co-movements among variables rather than a single estimate over the whole sample. Compared to the moving window approach, the advantage of the CPPF model and the FFF model is that the whole sample is considered, while the moving window approach is limited to past data and the length of the window used.
We find the estimates of beta, downside beta and upside beta estimates of both models to be highly significant to drive stock returns. The beta estimates positively drive stock returns. The downside and upside beta estimates demonstrate reversed impacts on stock returns, the downside beta has a negative impact on stock returns, while the upside beta, consistent with beta estimates, has a positive impact. This paper is arranged as follows: Section 2 provides literature reviews of both models, followed by Section 3 describes the data. Section 4 explains the econometric models and methods applied. Section 5 and Section 6 provide the empirical results and results of the Fama-Macbeth regressions. Section 7 concludes.
The existing literature of time-varying beta is extensive. There is a long history of literature that has argued that the Capital Asset Pricing Model (CAPM) proposed by Sharpe (1964) and Lintner (1965) is inadequate to explain the risk-return relationship due to the assumption of constant beta. To solve this drawback, relaxing the constant beta assumption and allowing time-varying beta is one possible method. There are a number of approaches to obtain time-varying betas, for instance, in Fama and Macbeth's (1973) study, moving window Ordinary Least Square (OLS) regression is applied to the market factor model to obtain time-varying beta. According to Härdle et al (1988), Härdle (1992), Wand and Jones (1995), Ang and Kristensen (2012), and Li and Yang (2011), various nonparametric approaches which are based on simulation are alternative methods to obtain the time-varying beta. Moreover, the time-varying beta also can be obtained by the multivariate GARCH based model proposed by Engle (2002), Andersen et al. (2002) and Nieto et al. (2011). More recently, Horvath et al. (2020) employ functional data regression to estimate time-varying beta in Chinese stock market, and Chakrbarti and Das (2021) adopt modified multivariate GARCH to capture time-varying risk in Indian and American stock market.
There is relatively fewer literature focusing on downside risk of stock market although research emerge in recent decade. Since Ang et al. (2006) point out that stocks with higher downside risk is compensated with additional risk premium, focus has primarily been lying on how to measure the downside beta. Giglio et al. (2016), Min and Kim (2016) and Li (2021) all adopt unique approaches to estimate downside risk in various financial markets.
Although existing literature are considerably well built on time-varying beta and downside risk separately, there is limited research focusing on time-varying feature of downside risk. Huang (2019) applies linear market model to estimate time-varying downside beta. More recently, Dobrynskaya (2021) investigates downside risk in cryptocurrency market by estimating downside beta at each point in time using high frequency data. The aforementioned studies use linear model to estimate time-varying downside beta, while fail to put appropriate timing weights (heavier weight on more recent data) when smoothing due to a lack of non-linear feature of the models. In this paper, to fill gap in the literature, specifically on both investigating and pricing the time variation and asymmetric nature of systematic risk. Instead of regressing the stock return upon the market portfolio return in a linear fashion, two alternative methods with non-linear feature are proposed to estimate the corresponding time-varying beta coefficients, namely the CPPF regression and the FFF regression, respectively. A brief literature of development of both methods is provided in following sections.
The advantages of the CPPF approach are, firstly, data can be flexibly adjusted without considering the sample size. Secondly, for research with particular focus on data smoothing, it allows time weight to be considered when estimating among various selected knots. Thirdly, apart from the time weight, the nature of the original data is retained and there are no extra functions or patterns to be put into the model.
In order to introduce CPPF approach, a review of cubic spline method is essential. The cubic spline method was originally used in mathematics and engineering (Ferguson, 1963). Mathematically, as a third order piecewise polynomial function, the cubic spline is used to smooth discrete points into a continuous curve. According to Rorres and Anton (1984), a cubic spline can be expressed mathematically in the following form:
S(x)={s1(x) if x1≤x≤x2s2(x) if x2≤x≤x3⋮sn−1(x) if xn−1≤x≤xn | (1) |
where it is assumed that si is the third order polynomial function defined by
si(x)=ai(x−xi)3+bi(x−xi)2+ci(x−xi)+di fori=1,2,3…….n−1. | (2) |
The first and second order derivative of Equation (1) defines the fundamentals of the process. These derivatives are given by
s′i(x)=3ai(x−xi)2+2bi(x−xi)+ci | (3) |
s″i(x)=6ai(x−xi)+2bi fori=1,2,3…….n−1 | (4) |
The piecewise polynomial function has the following properties:
1.The piecewise polynomial function interpolates all data points.
2.The S(x) function is continuous in the interval [x1, xn].
3.The first derivative of the S(x) function is continuous in the interval [x1, xn].
4.The second derivative of the S(x) function is continuous in the interval [x1, xn].
There are a number of studies which employ the cubic spline approach in financial modeling, mainly focusing on estimation of the term structure, autoregressive conditional duration (ACD) models and volatility of high frequency data. Vasicek and Fong (1982) and Jarrow et al. (2004) employ the cubic spline approach in estimating interest rate term structure. Engle and Russell (1998) proposed an ACD model which treats the time between transactions as a stochastic process. Within the ACD model, a daily seasonal factor is modelled by a cubic spline series. Moreover, in Zhang et al. (2001), a threshold autoregressive conditional duration (TACD) model is proposed and shown to be superior to the classic ACD model.1 More recently, Taylor(2004a, b) and adopted cubic splines in their studies in the context of modelling the volatility of high frequency data via the ACD model. Aside from the above studies, there are a few studies that have also adopted cubic splines as a modelling tool, such as Engle and Rivera (1991) who estimated the density factor by using cubic splines in an autoregressive conditional heterosecdasticity (ARCH) context. Yu and Ruppert (2002) introduced the cubic spline approach into the estimation of the single index model. Evans and Speight (2010) employed the cubic spline approach to model intraday exchange rate volatility.
1The cubic spline approach was particularly used to approximate seasonal factors within the model.
According to Eilers and Marx (2004), there are mainly two approaches used in cubic spline regression: the B-spline basis and truncated power functions basis.
For the B-spline basis approach, Eilers and Marx (2004) use equally-spaced knots and spline function B. Mathematically, the B-spline model can be written as
E(y)=μ=Bα | (5) |
and the objective function to be minimized is
QB=|y−Bα|2+λ|Ddα|2 | (6) |
which λ is a non-negative parameter, and Dd is the d-th difference of α, it can be written as
Dd=Δdα | (7) |
and
Δαj=αj−αj−1 | (8) |
Δ2αj=αj−2αj−1+αj−2 | (9) |
and so on for higher orders. So the objective function of QB leads the B-spline model to
(B'B+λD′dDd)ˆα=B'y | (10) |
It can be seen from Equation (10) that when λ = 0, it becomes the classic equation of linear regression.
For the truncated power functions basis, according to Ruppert et al. (2003), for a given asset i, column j and degree p, the truncated power function of F is written as
Fij=(xi−tj)PI(xi>tj) | (11) |
where I(u) is an indicator function, it is 0 when u < 0 and 1 otherwise. The vector t contains the knots, and the knots are placed as quantiles of x. Consequently, the model for E(y) can be written as
E(yi)=∑Pk=0βkxki+∑n−1j=1Fijbj | (12) |
And the objective function to be minimized is given by
QF=|y−βx−Fb|2+κ|b|2 | (13) |
The increasing of κ will increase the smoothness.
Eliers and Marx (2004) point out that both bases allow a mixed model approach, and the B-spline basis can be derived from the truncated power basis. They also show that the truncated power basis has bad numerical properties, and could cause discontinuities in estimation, while the B-spline basis approach has no such issue. However, according to Taylor (2004), the truncated power basis is employed in the spline-based periodical GARCH model on high frequency commodity future return data, which produces excellent smooth estimates. Therefore, in light of Taylor's (2004) study, the truncated power basis is employed in the CPPF approach in this paper, and the detail of the piecewise polynomial matrix used will be introduced in Section 4.2
2The CPPF approach is derived from the cubic spline approach described above. However, we do not refer to CPPF as a spline because we allow for discontinuities at each knot.
The advantages of the FFF approach are: firstly, in the context of normal and high frequency data, the macroeconomic news announcement effect has been filtered by the periodic pattern of the FFF, so there is no need to model the macroeconomic news announcement effect; secondly, the FFF approach creates a smooth pattern for volatility dynamics and changes; thirdly, the FFF approach is based on sound mathematics and the fit of the periodicity of financial data is widely agreed.
The FFF which was first proposed and refined by Gallant (1981, 1982 and 1984). This mathematical function, based on a Fourier series, was initially used to approximate the utility function and derive an appropriate expenditure system for the whole economy. Mathematically, it can be written as
∑Aα=1∑Jj=−Jajαeijk′αx=∑Aα=1{u0α+2∑Jj=1[ujαcos(jk′αx)−υjαsin(jk′αx)]} | (14) |
where
ajα=ujα+iυjα, α=1,2,3...A, j=0,±1,±2...± |
whereas i is defined as the imaginary unit, k is the order of the expansion, aj is the coefficient given by
aj=uj+iυj | (15) |
Recently, the FFF was widely applied in two aspects of economics and finance: estimation of production and banking efficiency, and modeling high frequency volatility periodicity.
In the former aspect, Chung et al. (2001) and Huang and Wang (2001) both applied the FFF in estimating the scale and scope of the Asian banking industry. Huang and Wang (2004) expanded the FFF and applied it to panel data to estimate multiproduct banking efficiency. Featherstone and Cader (2005) employed the FFF in a Bayesian econometrics context to evaluate agricultural production. And Yu et al. (2007) adapted the FFF to estimate agricultural banking efficiency.
Within a volatility context, Andersen and Bollerslev (1998) introduced the FFF into high frequency data volatility modelling. In their study, under GARCH framework, the FFF was used to estimate an intraday periodicity component in order to capture volatility reactions to macroeconomic announcements. The FFF within their study has been simplified as
f(θ,t,n)=μ0+∑Dk=1λk⋅Ik(t,n)+∑Pp=1(δc,p⋅cosp2πNn+δs,p⋅sinp2πNn) | (16) |
where Ik(t, n) is the indicator of event k during time interval n on day t, θ is the parameter vector to be estimated, and μ0, λk, δc, p and δs, p are the fixed coefficients to be estimated (Andersen and Bollerslev, 1998). Moreover, Andersen et al. (2000) applied the FFF in the Japanese stock market, while Bollerslev et al. (2000) employed the FFF in analyzing the US bond market. More recently, Evans and Speight (2010) further adopted the FFF in the foreign exchange market.
Although the cubic spline approach and FFF approach are widely used in the financial literature, the majority of studies use high frequency data in a financial derivatives market, banking industry or foreign exchange market. There is very few study using both approaches to estimate the downside and upside components of risk in stock markets. This paper employs the CPPF and FFF as tools, with various numbers of knots and the AIC used to uncover the best fit of beta, downside beta and upside beta estimates of monthly data with a long span in the US stock market and to improve the goodness of fit of asset pricing models.
As an experiment for CPPF and FFF models, data used in this paper are taken from Center for Research in Security Price (CRSP) database. This paper focuses on the ordinary common stocks listed on the New York Stock Exchange (NYSE), American Stock Exchange (AMEX) and NASDAQ measured on a monthly frequency from January 1960 to December 2010.3 American depositary receipts (ADR), real estate investment trust (REIT), closed-end funds, foreign firms and other securities which do not have a CRSP share code of 10 or 11 are excluded from the sample. Each stock is required to have at least 5 years of consecutive monthly adjusted return observations with at most 5 missing observations. The return of each stock is adjusted for stock splits, mergers and acquisitions, and dividends (dividends are subtracted from stock prices for adjustment), giving 13,557 stocks. The value-weighted return of all listed stocks is taken as a measure of the market portfolio, and the one month Treasury bill rate represents the risk free rate.4 A summary of stocks selected is shown in Table A.1 in Appendix.
3The NASDAQ data are only available from January 1972.
4Using the same criteria Ang et al. (2007) adopted.
All econometric models in this paper are based on the CAPM proposed by Sharpe (1964) shown as Equation (17)
Rit−Rft=αit+βit⋅(RMt−Rft)+ϵit | (17) |
where Rit is the rate of return of stock i at time t, Rft is the risk free rate at time t, αit is the constant at time t, βit is the coefficient to be estimated and represents the co-movement between stock i and the market at time t, RMt is the rate of return of market portfolio at time t, (RMt – Rft) is the excess return of the market portfolio, and εit is the error term of stock i at time t. It is this equation that will be estimated using the CPPF and FFF models.
For convenience, we define
xRit=Rit−Rft | (18) |
and
xRMt=RMt−Rft | (19) |
where xRit is the excess rate of return of stock i at time t, and xRMt is the excess rate of return of the market portfolio at time t.
In this section, the CPPF (with knots) model is described. By using the CPPF model, the excess rate of return on the market portfolio xRM will be divided into different numbers of series depending on the number of knots selected, thus allowing the betas to vary over time.
Deciding the number of knots to use is an interesting tradeoff (Stone, 1986). If a small number of knots are chosen, the estimates will be over-smooth with less variability, and could also be biased. By contrast, if a high number of knots is selected, the bias can be avoided, however, it will also lead to a high variability of estimates in the fit and could result in overfitting. Eilers and Marx (1996) discovered that up to 4 to 5 knots is most appropriate for most applications, therefore, the number of knots selected for the CPPF model will vary from 0 to 5.
Placement of knots follows the quintile method proposed by Stone (1986). In his study, he found that placing knots according to the quintile point with respect to the total number of observations results in less bias than placing knots according to a fixed number of observations. Therefore, the knots are placed at the quintile points as follows:
Number of knots | 0 | 1 | 2 | 3 | 4 | 5 |
Placement points | 50% | 33.3% | 25% | 20% | 16.6% | |
66.6% | 50% | 40% | 33.30% | |||
75% | 60% | 50% | ||||
80% | 66.6% | |||||
83.3% |
The econometric models used here take advantage of the CPPF approach, and apply it to the classic market model. To estimate the beta coefficient of each stock, the model, in matrix terms, can be written as
XRi=ααi+(XRMΘSN)⋅Bi+εεiN=0,1,2,3,4,5 | (20) |
where XRi is a (t × 1) column vector of excess returns of stock i, αi is a t × 1 column vector, XRM is a (t × 1) column vector of excess returns of the market portfolio, SN is a (t × n) piecewise polynomial matrix with N representing the number of knots, ʘ is the element to element multiplication sign which results in (XRM ʘ SN) becoming a (t × n) matrix, 5 Bi is the (n × 1) estimated beta column vector, and εi is the (t × 1) column vector error term.
5ʘ is conventionally used as an element to element multiplication sign when two matrices are in the same rank, we borrow it here for different rank matrices for the sake of simplicity.
Specifically, the piecewise polynomial matrix, SN, varies along with the number of knots selected. When the CPPF has no knots, S0 can be expressed as:
S0=[1011213202222330332334044243⋮⋮⋮⋮t0tt2t3]. | (21) |
Moreover, when one knot is selected, the knot will be placed at the 50% point of observations, with the S0 elements remaining in S1, plus new elements added in with elements valued 0 above the knot, therefore S1 can be written as
S1=[10112130000202222300003033233000040442430000⋮⋮⋮⋮⋮⋮⋮⋮(t2)0(t2)(t2)2(t2)30000⋮⋮⋮⋮10111213⋮⋮⋮⋮20212223⋮⋮⋮⋮30313233⋮⋮⋮⋮40414244⋮⋮⋮⋮⋮⋮⋮⋮t0tt2t3(t2)0(t2)(t2)2(t2)3] | (22) |
The expression for SN with two knots (N = 2), three knots (N = 3), four knots (N = 4) and five knots (N = 5) can be found in Illustration A.1–A.4 in Appendix. It can be seen from the expression for S0 and S1 that as the number of knots increases, the number of columns in SN will increase. More precisely, for every one extra knot placed, the number of columns in SNwill increase by 4, so the dimensions of S0, S1, S2, S3, S4 and S5 will be (t × 4), (t × 8), (t × 12), (t × 16), (t × 20) and (t × 24) respectively.
The OLS regression is applied to each stock to get the vector of beta estimates. Then, beta estimates for each stock at each point in time (BS in vector form) can be calculated as follows:
BS=SN⋅Bi | (23) |
It can be seen from Equation (14) that BS is a product of a (t × n) matrix SN and a (n × 1) estimated beta vector Bi, therefore regardless of the number of knots placed in the function, the rank of BS will always be t × 1. Since the number of knots varies from 0 to 5, there will be 6 possible Bi vectors for each stock corresponding to the number of knots used. In order to find the best fit for each stock, we follow Eilers and Marx's (1996) study, and use the AIC.6 The AIC can be expressed as
6The Schwarz Information Criteria (SIC) can also be used to determine the appropriate number of knots, in this paper, AIC is chosen instead of SIC since the SIC shows less tolerance when the number of parameters in the model is high, according to Eilers and Marx's study (1996).
AIC=−2ln(L)+2k | (24) |
where L is the maximum value of the likelihood function, and k is the number of parameters within the model. There are discontinuities at knots points, while the fitted values are smooth between each knot. These discontinuities are due to the column of ones in the piecewise polynomial matrix, and in these cases, we let the data to decide the appropriate value of the estimates.
In order to calculate the downside and upside beta estimates by using the CPPF model, the same logic is used with Equation (11) modified. Referring to Ang et al. (2006), the downside beta and upside beta are calculated as
β−=cov(xRi,xRM|xRM<¯xRM)var(xRM|xRM<¯xRM) | (25) |
and
β+=cov(xRi,xRM|xRM≥¯xRM)var(xRM|xRM≥¯xRM) | (26) |
where ¯xRM is the average market excess return over the sample period of the stock, and previous notations hold. In light of Ang et al. (2006), dummy variables (vectors) D1i and D2i are created and employed for each stock. D1 and D2 (with time subscript t) can be expressed as
D1i=1andD2i=0if xRM,t<¯xRM | (27) |
and
D1i=0andD2i=1if xRM,t≥¯xRM | (28) |
It can be seen from Equations (18) and (19) that D1i = 1 and D2i = 0 if the market excess return at time t is below the average market excess return, while D1i = 0 and D2i = 1 if the market excess return at time t is above the average market excess return. Then two more variables are created as follows:
D1ixRM=D1i⊙xRM | (29) |
D2ixRM=D2i⊙xRM | (30) |
It can be seen from Equations (20) and (21) that two new variables D1ixRM and D2ixRM are the element to element products of dummy variables of stock i and the corresponding excess market return over the sample period of the stock. For the former, observations are excess market returns if they are below the average excess market return over the sample period, and 0 otherwise. For the latter, observations are excess market returns if they are above the average excess market return over the sample period, and 0 otherwise.
The econometric model used to estimate downside and upside betas for each stock, in matrix form, can be written as
XRi=D1i+D2i+(D1iXRMΘSN)×Bi−+(D2iXRMΘSN)⋅Bi++εεi | (31) |
where N = {0, 1, 2, 3, 4, 5}, Bi- and Bi+ are the (n × 1) estimated downside and upside beta estimate column vectors. Since the number of knots varies from 0 to 5, there will be 6 pairs of Bi- and Bi+vectors. The best downside and upside beta estimates as determined by the AIC for a stock at each point of time, (B−∗S and B+∗S in vector form) can be conducted as follows
B−∗S=SN⋅B−i | (32) |
and
B+∗S=SN⋅B+i | (33) |
As mentioned in the previous paragraph, regardless of the number of knots placed in the function, the dimensions of B−∗S and B+∗S will always be t × 1. Both downside and upside betas can be interpreted in an analogous manner to classic beta regarding to downside and upside market.
In this section, the FFF model is described in detail. In light of Andersen and Bollerslev (1998), Andersen et al. (2000), Bollerslev et al. (2000), and Evans and Speight (2010), the econometric model with the FFF specification employed in this paper is defined as
xRit=αit+∑Pp=1[βcos,p⋅(cosp2πNn⋅xRMt)+βsin,p⋅(sinp2πNn⋅xRMt)]+ϵit | (34) |
where αit is the constant, βcos, p and βsin, p are the coefficients to be estimated for stock i, N is the total number of observations of stock i, n is the order of observations with n = {1, 2, 3…t}, εit is the error term of stock i at time t, and p is the order of the FFF. The order of the FFF can vary from 1 to infinity. However, in order to provide efficient and unbiased estimates, according to previous studies, we chose up to 4.7 In this paper, the order from 1 to 4 is selected to examine and discover the best fit of the estimates.
7See Andersen and Bollerslev (1998), and Evans and Speight (2010) for similar assumptions.
The OLS regression is applied to each stock to get the βcos, p and βsin, p estimates. The AIC will then be used for each regression. Since an order of 1 to 4 is examined, there are 4 AICs for each stock. Taking advantage of the nature of the AIC, the regression that produces the least AIC will indicate the optimal fit. To calculate the best estimate for a stock at each point of time, the AIC supported estimates of βcos, p and βsin, p for each stock are used to get the best estimates at each point of time, specifically,
β∗F=∑Pp=1(βcos,p⋅cosp2πNn+βsin,p⋅sinp2πNn) | (35) |
In order to calculate the downside and upside beta estimates using the above FFF model, the same logic is followed as in the CPPF case. The same variables created in Equations (20) and (21) are created and employed for each stock in the new FFF model and the new market model is given by
xRit=D1i+D2i+P∑p=1[β−cos,p⋅(cosp2πNn⋅D1ixRMt)+β−sin,p⋅(sinp2πNn⋅D1ixRMt)] |
+∑Pp=1[β+cos,p⋅(cosp2πNn⋅D2ixRMt)+β+sin,p⋅(sinp2πNn⋅D2ixRMt)]+ϵit | (36) |
where β−cos,p and β−sin,p are the downside market coefficients to be estimated for stock i, and β+cos,p and β+sin,p are the upside market coefficients to be estimated for stock i, and previous notations hold. As in Equation (22), there is no conventional constant term in the model, rather, full set of dummy variables are instead used. Since the order of the FFF examined varies from 1 to 4, there will be 4 groups of β−cos,p, β−sin,p, β+cos,p and β+sin,p for each stock. For each group an AIC value is calculated and the lowest value indicates the best fit group of β−cos,p, β−sin,p, β+cos,p and β+sin,p.
Furthermore, the best downside and upside beta estimates for each point in time are given by
β−∗F=∑Pp=1(β−cos,p⋅cosp2πNn+β−cos,p⋅sinp2πNn) | (37) |
and
β+∗F=∑Pp=1(β+cos,p⋅cosp2πNn+β+cos,p⋅sinp2πNn) | (38) |
and all previous notations hold. In the next section, the empirical results will be demonstrated and analyzed in detail.
As the method explained in the previous section, the best fits for both the CPPF and FFF models are obtained. In order to illustrate the results in a clearer way, having the best fitted estimates β∗S, β−∗S and β+∗S for the CPPF model and β∗F, β−∗F and β+∗F for the FFF model, the number of stocks with corresponding numbers of knots or orders and the percentage of the whole sample are shown in Table 2 and Table 3, respectively.
Knots | 0 | 1 | 2 | 3 | 4 | 5 | |
β∗S | Number of Stocks | 9409 | 1455 | 861 | 655 | 614 | 563 |
Percentage to Whole sample | 69.40% | 10.73% | 6.35% | 4.83% | 4.53% | 4.15% | |
β−∗S and β+∗S | Number of Stocks | 9399 | 903 | 483 | 416 | 729 | 1627 |
Percentage to Whole sample | 69.33% | 6.66% | 3.56% | 3.07% | 5.38% | 12.00% | |
Note: This table reports the number and percentage of stocks with different knots to construct the best fit estimates of CPPF model. |
Order | 1 | 2 | 3 | 4 | |
β∗F | Number of Stocks | 6204 | 2746 | 2099 | 2508 |
Percentage to Whole sample | 45.76% | 20.26% | 15.48% | 18.50% | |
β−∗F and β+∗F | Number of Stocks | 8377 | 2293 | 1429 | 1458 |
Percentage to Whole sample | 61.79% | 16.91% | 10.54% | 10.75% | |
Note: This table reports the number and percentage of stocks in different orders to construct the best fit estimates of the FFF model. |
For the CPPF model, it can be seen from Table 2 that for 9409 stocks, best estimates are obtained when no knots are used. For other knot values, the number of stocks decreases. Typically, when 5 knots are used, just 563 stocks produced the best estimates. Similar results are obtained when the downside and upside beta estimates are constructed.
For the FFF model, it is clear from Table 3 that to construct β∗F, 6204 stocks have an order of 1. Orders 2, 3 and 4 are generally selected less often. This pattern is even more obvious when constructing β−∗F and β+∗F, 8377 stocks produce the best estimates with order 1, 2293 stocks with order 2, and the number of stocks with order 3 and 4 are 1429 and 1458, respectively.
Furthermore, the relationships among stock returns and corresponding beta, downside beta and upside beta estimates for the CPPF model and the FFF models are examined. In order to uncover the relationship in a cross-sectional fashion, stocks at each point of time are cross-sectionally assigned into five portfolios according to the value of the estimate. Since the beta, upside beta and downside beta estimates for both models are not independent of each other, to distinguish the effects among them, the relative upside beta, denoted by (β+−β) and relative downside beta denoted by (β-−β) are considered. To sort the portfolio, at each point of time, all stocks are sorted into five quintiles according to the value of the beta estimate. Therefore, portfolio 1 contains stocks with the lowest 20% of estimates, portfolio 2 contains stocks with the second lowest 20% of estimates, and accordingly. When stocks are sorted into 5 portfolios at each point of time, 8 the equally weighted average of the estimate for each portfolio and the corresponding average annualized stock returns are calculated. The results of both models are summarized in Table 4 and Table 5, respectively.
Portfolio | Return | βS* | βS−* | βS+* | Portfolio | Return | βS* | βS−* | βS+* | ||||||||
Panel 1 Stocks Sorted by βS* | Panel 2 Stocks Sorted by βS−* | ||||||||||||||||
1 Low | 1.53% | −0.24 | 0.47 | −0.21 | 1 Low | 25.00% | 0.48 | −7.2 | 0.85 | ||||||||
2 | 6.76% | 0.64 | 0.80 | 0.32 | 2 | 11.08% | 0.77 | 0.68 | 0.68 | ||||||||
3 | 8.95% | 0.99 | 1.23 | 0.84 | 3 | 10.53% | 1.03 | 1.1 | 0.83 | ||||||||
4 | 11.59% | 1.41 | 1.54 | 1.25 | 4 | 9.73% | 1.33 | 1.61 | 0.99 | ||||||||
5 High | 24.84% | 2.68 | 2.18 | 2.4 | 5 High | −2.67% | 1.86 | 10.01 | 1.23 | ||||||||
High − Low | 23.30% | 2.92 | 1.72 | 2.62 | High − Low | −27.67% | 1.38 | 17.21 | 0.39 | ||||||||
Panel 3 Stocks Sorted by βS+* | Panel 4 Stocks Sorted by (βS−*−βS*) | ||||||||||||||||
1 Low | −11.93% | 0.42 | 1.32 | −3.87 | 1 Low | 35.87% | 1.64 | −6.6 | 2.03 | ||||||||
2 | 6.67% | 0.71 | 0.91 | 0.34 | 2 | 14.74% | 1.11 | 0.88 | 1.09 | ||||||||
3 | 9.85% | 0.96 | 1.21 | 0.81 | 3 | 9.72% | 0.96 | 1.05 | 0.72 | ||||||||
4 | 13.43% | 1.28 | 1.35 | 1.37 | 4 | 5.72% | 0.95 | 1.38 | 0.6 | ||||||||
5 High | 35.58% | 2.11 | 1.42 | 5.94 | 5 High | −12.35% | 0.83 | 9.49 | 0.16 | ||||||||
High − Low | 47.51% | 1.69 | 0.1 | 9.81 | High − Low | −48.22% | −0.81 | 16.09 | −1.87 | ||||||||
Panel 5 Stocks Sorted by (βS+*−βS*) | Panel 6 Stocks Sorted by (βS−*−βS+*) | ||||||||||||||||
1 Low | −6.16% | 1.35 | 1.88 | −3.38 | 1 Low | 36.79% | 1.43 | −6.32 | 5.09 | ||||||||
2 | 7.24% | 1.03 | 1.37 | 0.52 | 2 | 14.15% | 1.05 | 0.89 | 1.2 | ||||||||
3 | 9.94% | 0.96 | 1.16 | 0.84 | 3 | 10.07% | 0.96 | 1.05 | 0.83 | ||||||||
4 | 13.12% | 1 | 1.1 | 1.21 | 4 | 6.43% | 0.98 | 1.35 | 0.53 | ||||||||
5 High | 29.51% | 1.14 | 0.71 | 5.4 | 5 High | −13.76% | 1.05 | 9.23 | −3.06 | ||||||||
High − Low | 35.67% | −0.2 | −1.18 | 8.78 | High − Low | −50.55% | −0.38 | 15.55 | −8.15 | ||||||||
Note: This table presents the relationship between excess stock returns and factor loading of the CPPF model. The column labeled "return" reports the average stock returns over a one month T-bill rate. "High-Low" reports the difference between portfolio 5 and portfolio 1. Notably, in a perfect market, the average value of βS* is assumed to be 1. |
Portfolio | Return | βF* | βF−* | βF+* | Portfolio | Return | βF* | βF−* | βF+* | ||||||||
Panel 1 Stocks Sorted by βF* | Panel 2 Stocks Sorted by βF−* | ||||||||||||||||
1 Low | 3.54% | −1.19 | −0.74 | −0.87 | 1 Low | 21.88% | −0.67 | −1.32 | −0.31 | ||||||||
2 | 7.51% | −0.37 | −0.26 | −0.28 | 2 | 13.09% | −0.24 | −0.38 | −0.12 | ||||||||
3 | 9.11% | −0.01 | −0.01 | −0.03 | 3 | 10.46% | −0.01 | −0.01 | −0.01 | ||||||||
4 | 11.24% | 0.36 | 0.24 | 0.23 | 4 | 8.22% | 0.23 | 0.37 | 0.11 | ||||||||
5 High | 22.26% | 1.22 | 0.76 | 0.95 | 5 High | −0.03% | 0.69 | 1.33 | 0.33 | ||||||||
High − Low | 18.72% | 2.41 | 1.49 | 1.82 | High − Low | −21.91% | 1.36 | 2.65 | 0.64 | ||||||||
Panel 3 Stocks Sorted by βF+* | Panel 4 Stocks Sorted by (βF−*−βF*) | ||||||||||||||||
1 Low | −7.72% | −0.81 | −0.29 | −1.28 | 1 Low | 33.18% | 0.52 | −0.71 | 0.65 | ||||||||
2 | 5.02% | −0.27 | −0.13 | −0.36 | 2 | 16.29% | 0.11 | −0.18 | 0.19 | ||||||||
3 | 9.44% | −0.02 | −0.01 | −0.02 | 3 | 10.34% | −0.02 | −0.02 | −0.01 | ||||||||
4 | 14.58% | 0.24 | 0.13 | 0.33 | 4 | 4.66% | −0.14 | 0.17 | −0.21 | ||||||||
5 High | 32.33% | 0.86 | 0.3 | 1.33 | 5 High | −10.84% | −0.48 | 0.74 | −0.61 | ||||||||
High − Low | 40.05% | 1.67 | 0.59 | 2.61 | High − Low | −44.02% | −1 | 1.45 | −1.27 | ||||||||
Panel 5 Stocks Sorted by (βF+*−βF*) | Panel 6 Stocks Sorted by (βF−*−βF+*) | ||||||||||||||||
1 Low | −4.27% | 0.36 | 0.57 | −0.63 | 1 Low | 36.12% | 0.17 | −0.81 | 0.85 | ||||||||
2 | 5.48% | 0.12 | 0.21 | −0.14 | 2 | 16.44% | −0.02 | −0.24 | 0.17 | ||||||||
3 | 9.74% | 0 | −0.01 | −0.01 | 3 | 10.15% | −0.02 | −0.01 | −0.02 | ||||||||
4 | 14.99% | −0.14 | −0.23 | 0.12 | 4 | 4.34% | −0.02 | 0.22 | −0.2 | ||||||||
5 High | 27.69% | −0.34 | −0.54 | 0.66 | 5 High | −13.41% | −0.12 | 0.85 | −0.81 | ||||||||
High − Low | 31.95% | −0.7 | −1.1 | 1.29 | High − Low | −49.53% | −0.3 | 1.66 | −1.66 | ||||||||
Note: This table presents the relationship between excess stock returns and factor loading of the FFF model. The column labeled "return" reports the average stock returns over one month T-bill rate. "High-Low" reports the difference between portfolio 5 and portfolio 1. Notably, in a perfect market, the average value of βF* depends on the average value of the intercept. |
8Since monthly data are used in this paper, and the whole sample is from January 1960 to December 2010, there are 612 time points.
It can be seen from the CPPF model results in Table 4 that when sorting by β∗S, portfolio 1 has an average β∗S of −0.24 while on the other hand, portfolio 5 has an average β∗S of 2.68. Consistent with the literature, the average annualized realized rates of return of each portfolio show an ascending order as the average β∗S increases, portfolio 1 yields a return of 1.53% while portfolio 5 shows a return of 24.84%. The average β−∗S and β+∗S values of each portfolio follow the same trend as β∗S, an average β−∗S is 0.47 in portfolio 1 and increases to 2.18 in portfolio 5. Similarly, the average β+∗S is −0.21 in portfolio 1 and increases to 2.4 in portfolio 5.
Interestingly, a different pattern in returns is demonstrated when stocks are sorted by β−∗S. It is clear that average returns demonstrate a reversed trend while average β∗S shows the same ascending trend from portfolio 1 to portfolio 5 along with the increase of β−∗S. β−∗S is −7.2 in portfolio 1 with an average β∗S of 0.48 and an average return of 25%, while in portfolio 5, β−∗S grows to 10.01 with an average β∗S increasing to 1.86 and the average return drops to −2.67%. The pattern of β+∗S is generally increasing but with a subtle variation in that, it drops from 0.85 to 0.62 from portfolio 1 to portfolio 2, and then keeps growing to portfolio 5 ending up with a value of 1.23.
Notably, although β∗S and β+∗S still have increasing trends in this panel, the difference between values for portfolio 1 and portfolio 5 (1.38 and 0.39, respectively) are narrower than the ones in Panel 1 (2.92 and 2.62 respectively).
When stocks are sorted by β+∗S, a similar pattern appears to those in Panel 1. It can be seen from Panel 3 that, β+∗S is −3.87 in portfolio 1 with an average β∗S of 0.42 and an average return of −11.93%, and in portfolio 5, β+∗S grows to 5.94 with average β∗S increasing to 2.11 and average returns increasing to 35.58%. The pattern of β+∗S is also generally increasing but with a sudden drop from 1.32 to 0.91 between portfolio 1 to portfolio 2, and then keeps increasing to portfolio 5 and ends up with a value of 1.42. Compared to Panel 1, the spread of β∗S, and β−∗S between portfolio 5 and portfolio 1 is less, however the spread of β+∗S and average return is much higher (9.81 and 47.51% in Panel 3, while 2.62 and 23.3% in Panel 1). It is clear that β+∗S has the same positive impact on portfolio returns as β∗S.
In order to examine how β−∗S is driving the return not considering the impact of β∗S, a new estimate (β−∗S − β∗S) is employed in the analysis. Using this estimate to sort portfolios could discover the unique property of β−∗S after controlling for β∗S. When stocks are sorted by (β−∗S − β∗S), an unfamiliar pattern appears in Panel 4. From portfolio 1 to portfolio 5, all average returns, β∗S and β+∗S are in descending order while only β−∗S increases from −6.6 to 9.49. Although in Panel 1, Panel 2 and Panel 3, β−∗S are also in ascending order, the spread of average returns in Panel 4 is the highest and reaches −48.22%. So controlling for β∗S, it can be seen that β−∗S shows a negative relationship with portfolio returns and β∗S.
As with (β−∗S − β∗S), (β+∗S − β∗S) is employed to uncover the unique property of β+∗S after controlling for β∗S. It can be seen from Panel 5 that from portfolio 1 to portfolio 5, both average return and β+∗S are in ascending order, starting at −6.16% and −3.38, increasing to 29.51% and 5.4, respectively. β−∗S exhibits a descending trend for the first time within these panels. It drops from −1.88 to 0.71. A U-shaped pattern in β∗S is apparent, it starts at 1.35 in portfolio 1 and drops to 0.96 in portfolio 3, but restores to 1.14 in portfolio 5.
In Panel 6, (β−∗S−β+∗S) is adopted to sort the portfolio in order to control β+∗S from β−∗S and for the sake of precision. In can be seen from Panel 6 that both average returns and β+∗S are in descending orders, starting at 36.79% and 5.09 and dropping to 13.76% and −3.06 respectively, while β−∗S exhibits an ascending trend increasing from -6.32 to 9.23. As in Panel 5, a U-shaped pattern appears in β∗S, starting at 1.43 in portfolio 1 dropping to 0.96 in portfolio 3 and recovering to 1.05 in portfolio 5. Notably, the spread of returns in Panel 6 is the highest among all 6 panels at −50.55%.
Regarding the FFF model, it can be seen from Table 5 that when sorting by β∗F, portfolio 1 has an average β∗F of −1.19 while on the other hand, portfolio 5 shows an average β∗F of 1.22. Again, consistent with the literature, the average annualized rates of return to each portfolio are presented in an ascending order with the β∗F. Portfolio 1 yields a return of 3.54%, while portfolio 5 has a return of 22.26%. Average β−∗F and β+∗F for each portfolio follow the same trend as β∗F, with an average β−∗F of −0.74 in portfolio 1 and 0.76 in portfolio 5. Similarly, the average β+∗F is −0.87 in portfolio 1 and increases to 0.95 in portfolio 5. A different pattern was demonstrated when stocks are sorted by β−∗F. It is clear that average returns exhibit a reversed trend while average β∗F shows the same ascending trend from portfolio 1 to portfolio 5 along with the increase of β−∗F. β−∗F is −1.32 in portfolio 1 with an average β∗F of −0.67 and average return of 21.88%, while in portfolio 5, β−∗F grows to 1.33, average β∗F increases to 0.69 and the average return drops to −0.03%. The pattern of β+∗F is consistent with β−∗F, it starts at −0.31 in portfolio 1, and then keeps growing to portfolio 5 ending up with a value of 0.33. Notably, although β∗F and β+∗F still have an increasing trend in this panel, the difference between values in portfolio 1 and portfolio 5 (1.36 and 0.64, respectively) are narrower than the ones in Panel 1 (2.41 and 1.82, respectively).
When sorting by β+∗F, a similar pattern to that in Panel 1 appears. It can be seen from Panel 3 that β+∗F is −1.28 in portfolio 1 with an average β∗F of −0.81 and an average return of −7.72%, and in portfolio 5, β+∗F grows to 1.33 with an average β∗F increasing to 0.86 and an average return increasing to 32.33%. The pattern of β−∗F is consistent with β+∗F, it starts at −0.29 in portfolio 1, and then keeps growing to portfolio 5 ending up with a value of 0.3. Compared to Panel 1, the spread of β∗F, and β−∗F between portfolio 5 and portfolio 1 are less. However, the spread of β+∗F and average returns is much higher (2.61 and 40.05% in Panel 3, and 1.82 and 18.72% in Panel 1, respectively).
As with the CPPF model, we consider (β−∗F − β∗F) in the analysis. Using this estimate to sort portfolios could uncover further properties of β−∗F after controlling for β∗F. When stocks are sorted by (β−∗F−β∗F), an unfamiliar pattern appears in Panel 4. From portfolio 1 to portfolio 5, average returns, β∗F and β+∗F all decrease only with β−∗F increasing from −0.71 to 0.74. As in Panel 1, Panel 2 and Panel 3, β−∗F are also in ascending order, the spread of average returns in Panel 4 is highest and reaches −44.02%. Referring back to when stocks are sorted by (β−∗S−β∗S), a similar pattern appears. Therefore, it can be concluded that when controlling for the effect of beta, the relative downside beta estimates of both models have a negative relationship with portfolio returns, which can be interpreted as stocks tend to suffer a loss if they have large downside betas.
When (β+∗F−β∗F) is employed (to uncover the unique property of β+∗F after controlling for β∗F), it can be seen from Panel 5 that from portfolio 1 to portfolio 5, both average returns and β+∗F increase, starting at −6.16% and −3.38% and increasing to 29.51% and 5.4, respectively, while both β∗F and β−∗F decrease from 0.36 and 0.57 to −0.34 and −0.54, respectively.
Also in Panel 6, (β−∗F − β+∗F) is adopted to sort the portfolio to control β+∗F from β−∗F and for the sake of precision. It can be seen from this panel that the same pattern as in Panel 4 appears. From portfolio 1 to portfolio 5, average returns, β∗F and β+∗F all decrease while only β−∗F increases from −0.81 to 0.85. Notably, the spread of returns in Panel 6 is the highest among all 6 panels at −49.53%.
To sum up, the results of the CPPF and the FFF models, β∗S and β∗F as classic risk estimates, still have a clear impact stock returns. Specifically, when stocks are sorted by β∗S and β∗F, average returns follow exactly the same increasing trend with the β∗S and β∗F presented, even the portfolio return is inversely related to the market return. More importantly, it can be seen from these panels that β−∗S, β+∗S, β−∗F and β+∗F do have an impact on stock returns. When stocks are sorted by β−∗S, (β−∗S−β∗S), β−∗F and (β−∗F−β∗F), clearly, downside related estimates have a negative relationship with the realized returns, this is even more obvious in Panel 6 when stocks are sorted by (β−∗S−β+∗S) and (β−∗F−β+∗F), while when stocks are sorted by β+∗S, (β+∗S − β∗S), β+∗F and (β+∗F − β∗F), positive relationships appear between upside related estimates and realized returns. Moreover, the classic estimates β∗S and β∗F appear to have a similar impact as the upside related estimates β+∗S, (β+∗S−β∗S), β+∗F and (β+∗F − β∗F). To rationalize that, when downside beta is calculated, the return of the market portfolio is below the average, and very likely to be negative. The stock expected excess return is the product of beta and excess returns to the market portfolio, so when stocks are sorted by downside beta into portfolios, the larger the downside beta, the lower the return. In addition to that, the panic on the falling market of investors' could also be a reason for aggravating the negative returns.
In this section, in order to uncover the premia of β∗S, β−∗S, β+∗S, β∗F, β−∗F and β+∗F on stock returns from a cross-sectional point of view, a series of Fama-Macbeth regressions are performed which employ various combinations of the above estimates as independent variables.
For the purposes of comparison and to demonstrate the importance of placing an appropriate number of knots for the CPPF model and choosing an appropriate order for the FFF model, additional variables βS0, β−S0, β+S0, βF1, β−F1 and β+F1 are introduced. For βS0, β−S0 and β+S0, they are the beta estimate, downside beta estimate and upside beta estimate, respectively, for each stock at each point in time estimated with the CPPF model without placing a knot. However for βF1, β−F1 and β+F1, they are the beta estimate, downside beta estimate and upside beta estimate, respectively, for each stock at each point in time estimated with the FFF model with order 1.
The correlation coefficient matrices for groups of variables are presented in Table 6 and Table 7. It can be seen from Table 6 that βF1 is highly correlated with both β−F1 and β+F1 (the correlation coefficients are 0.7331 and 0.8093, respectively). The other pairs of variables are correlated with each other to some extent, but not as highly as the two mentioned pairs (above 0.5), for instance, βS0 and β+F1 exhibits the highest correlation with a coefficient of 0.462 after the two peak values. Therefore, in the following Fama-Macbeth regression, βF1 will not appear in the same regression with β−F1 or β+F1, and the other variables will form different combinations of independent variables.
βS0 | β−S0 | β+S0 | βF1 | β−F1 | β+F1 | |
βS0 | 1.0000 | |||||
β−S0 | 0.3363 | 1.0000 | ||||
β+S0 | 0.4326 | 0.0365 | 1.0000 | |||
βF1 | 0.4397 | 0.1345 | 0.1770 | 1.0000 | ||
β−F1 | 0.3209 | 0.2611 | 0.0458 | 0.7331 | 1.0000 | |
β+F1 | 0.4620 | 0.0336 | 0.2737 | 0.8093 | 0.2886 | 1.0000 |
Note: This table reports the correlation coefficients between factor loadings of both the CPPF model with zero knots and the FFF model in order one. To avoid unnecessary repetition, only the lower triangle of the matrix is shown. |
β∗S | β−∗S | β+∗S | β∗F | β−∗F | β+∗F | |
β∗S | 1.0000 | |||||
β−∗S | 0.0130 | 1.0000 | ||||
β+∗S | 0.0237 | -0.0001 | 1.0000 | |||
β∗F | 0.3154 | 0.0029 | 0.0085 | 1.0000 | ||
β−∗F | 0.2609 | 0.0082 | 0.0025 | 0.5239 | 1.0000 | |
β+∗F | 0.3099 | -0.0018 | 0.0098 | 0.6689 | 0.1889 | 1.0000 |
Note: This table reports the correlation coefficients between factor loadings of both the CPPF model with appropriate number of knots and the FFF model in appropriate order, to avoid unnecessary repetition, only the lower triangle of the matrix is shown. |
It is clear from Table 7 that, as in Table 6, high correlations appear between the FFF based estimates, with β∗F highly correlated with β−∗F and β+∗F (with correlation coefficients of 0.5239 and 0.6689, respectively). The remaining variables exhibit a weaker correlation with each other. Thus in the following Fama-Macbeth regression, β∗F will not appear in the same regression with β−∗F or β+∗F, and the other variables will form different combinations of independent variables. Notably, in Table 8, for both the CPPF and the FFF models, the upside beta estimates β+∗S and β+∗F are negatively correlated with the downside beta estimates of the CPPF model β−∗S with correlation coefficients of −0.0001 and −0.0018, respectively.
1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 | 11 | |
Intercept |
0.00155 [1.14] |
0.00743***
[3.82] |
0.00154 [1.21] |
0.00852***
[3.48] |
0.00907***
[3.90] |
0.00212 [1.57] |
0.00705***
[3.88] |
0.00213*
[1.70] |
0.00345***
[2.89] |
0.00801***
[4.06] |
0.00251**
[2.24] |
βS0 | 0.00663*** [3.25] |
0.00824***
[3.68] |
0.00614***
[3.02] |
0.00768***
[3.46] |
0.00525***
[2.73] |
0.00626***
[2.98] |
|||||
β−S0 | 0.0000637 [0.08] |
−0.000721 [−1.28] | 0.00187***
[2.99] |
−0.000678 [−1.20] | −0.0119*** [−13.04] | −0.000103 [−0.15] | 0.00147***
[2.95] |
||||
β+S0 | 0.000937*
[1.90] |
−0.000977*** [−2.13] | −0.000884** [−2.21] | −0.000983** [−2.14] | 0.0115***
[14.41] |
0.000520 [1.14] |
−0.00205*** [−4.72] | ||||
βF1 | 0.00542***
[3.35] |
0.00132**
[2.50] |
0.00141***
[2.66] |
0.00509***
[3.77] |
|||||||
β−F1 | −0.0113*** [−9.43] | −0.0126*** [−12.50] | −0.0131*** [−15.27] | ||||||||
β+F1 | 0.0148***
[10.83] |
0.0156***
[12.60] |
0.0127***
[17.27] |
||||||||
Number of obs | 2398103 | 2398103 | 2398103 | 2398103 | 2398103 | 2398103 | 2398103 | 2398103 | 2398103 | 2398103 | 2398103 |
Average R2 | 0.055 | 0.060 | 0.098 | 0.017 | 0.029 | 0.056 | 0.077 | 0.100 | 0.065 | 0.069 | 0.106 |
Note: This table reports the results of the Fama-Macbeth regression of factor loadings without knot and order. The t-statistics in the square brackets are calculated by using Newey-West (1987) heteroscedastic robust standard error with 12 lags. * denotes significance at the 10% level, ** denotes significance at the 5% level and ***denotes significance at the 1% level. |
This phenomenon potentially shows that the downside beta estimates do have an opposite impact on stock returns compared to the upside beta estimates which complies with the conclusion made in the previous sections.
After deciding on the possible combinations of variables, the Fama-Macbeth regressions are performed on both groups of variables which are demonstrated in Table 8 and Table 9. Since the data are at monthly frequency from January 1960 to December 2010, there are 612 cross-sectional time points and 2,398,103 observations of each regression. Newey-West (1987) heterosecdasticity robust standard errors with 12 lags are employed to calculate the t-statistics and the adjusted R2 values obtained from the cross-sectional regressions are provided.
1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 | 11 | |
Intercept |
0.00186 [1.52] |
0.00898***
[7.05] |
0.00538***
[5.04] |
0.00826***
[3.56] |
0.00943***
[4.42] |
0.00279**
[2.19] |
0.00965***
[7.54] |
0.00618***
[5.76] |
0.00369***
[3.38] |
0.00947***
[7.14] |
0.00623***
[6.23] |
β∗S | 0.00610***
[3.42] |
0.00448***
[3.82] |
0.00481***
[3.00] |
0.00383***
[3.54] |
0.00490***
[3.20] |
0.00384***
[3.42] |
|||||
β−∗S | −0.0120***
[−9.82] |
−0.00948***
[−11.19] |
−0.0106***
[−10.25] |
−0.00944***
[−11.19] |
−0.00994***
[−12.94] |
−0.00850***
[−11.30] |
|||||
β+∗S | 0.0134***
[14.47] |
0.00911***
[14.41] |
0.0117***
[15.25] |
0.00890***
[14.63] |
0.0107***
[15.74] |
0.00823***
[14.44] |
|||||
β∗F | 0.00673***
[4.16] |
0.00334***
[5.67] |
0.00238***
[4.73] |
0.00342***
[3.91] |
|||||||
β−∗F | −0.0113***
[−9.13] |
−0.00419***
[−8.21] |
−0.0111***
[−10.53] |
−0.00485***
[−13.74] |
|||||||
β+∗F | 0.0155***
[13.00] |
0.00647***
[11.34] |
0.0121***
[15.64] |
0.00586***
[13.51] |
|||||||
Number of obs |
2398103 | 2398103 | 2398103 | 2398103 | 2398103 | 2398103 | 2398103 | 2398103 | 2398103 | 2398103 | 2398103 |
Adjusted R2 | 0.092 | 0.13 | 0.147 | 0.039 | 0.055 | 0.096 | 0.138 | 0.150 | 0.11 | 0.137 | 0.153 |
Note: This table reports the results of the Fama-Macbeth regression of factor loadings with appropriate knots and orders on stock excess returns. The t-statistics in the square brackets are calculated by using the Newey-West (1987) heteroscedastic robust standard error with 12 lags. * denotes significance at the 10% level, ** denotes significance at the 5% level and ***denotes significance at the 1% level. |
As βF1 is excluded from the regression when β−F1 or β+F1 are employed, therefore, there are 11 possible combinations among βS0, β−S0, β+S0, βF1, β−F1 and β+F1 as independent variables. It can be seen from Table 8 that regression 1, 2 and 3 examine the impact of estimates of the CPPF model without a knot. Generally, these three regressions exhibit poor fit with the intercept in regression 1 and coefficients of β−S0 and β+S0 in regression 2 not significant at the 5% significance level and with adjusted R2 values of 0.06 and 0.055, respectively. In regression 3, the coefficient on β−S0 and the intercept are not significant even at the 10% significance level with an adjusted R2 value of 0.098. Regression 4 and 5 examine the impact of estimates of the FFF model with order 1. Estimates of regression 4 and 5 are all significant at the 1% significance level, however both regressions present low adjusted R2 values of 0.017 and 0.029, respectively.
The remaining regressions 6 to 11 employ variables from both the CPPF model and the FFF model to examine the impact of these variables on stock returns. It is clear from Table 8 that with the exception of the coefficient on β−S0 and coefficients on β−S0 and β+S0 in regression 10, remaining estimates are all significant at the 5% level. Among the regressions in Table 8, regression 11 shows the highest adjusted R2 value at 0.106 and it also contains the most variables. Notably, consistent with the literature, the estimated coefficients of βS0 and βF1 are always positive among regressions, and illustrate that the beta estimates for both models have a positive impact on stock returns. Moreover, the estimated coefficients on β−F1 and β+F1 are always significant at the 1% significance level, and their signs are constantly negative and positive, and show that the downside and upside risk estimates of the FFF model have negative and positive impacts on stock returns, respectively. However, the significance and sign of the estimated coefficients on β−S0 and β+S0 vary across regressions in Table 9, therefore, it is difficult to provide a definitive conclusion.
Table 9 shows the 11 possible combinations among β∗S, β−∗S, β+∗S, β∗F, β−∗F and β+∗F. It can be concluded from Table 9 that, similar to Table 8, regressions 1 to 3 examine the impact of estimates of the CPPF model with appropriate numbers of knots according to the AIC. Moreover, regressions 4 and 5 examine the impact of best estimates of the FFF model (according to the AIC). The remaining regressions 6 to 11 employ variables from both the best CPPF and FFF models to examine the impact of these variables on stock returns. Unlike Table 8, all estimated coefficients except the intercept in regression 1 are significant at the 1% significance level. These best fit estimates are all highly significant, and the beta coefficients are consistently positive, which is consistent with the classic literature (as in Table 8). Moreover, the estimated coefficients of downside and upside beta estimates show negative and positive signs, respectively, over all regressions, which is consistent with the conclusions made regarding Table 4 and Table 5. Furthermore, showing better fit than classic CAPM, among all 11 regressions in Table 9, regression 11 exhibits the highest adjusted R2 value at 0.153.9
9The outcome of CAMP is not shown due to limited space.
However, since β∗F is excluded from regressions employed β−∗F or β+∗F, checking the alterative regression (regression 8) shows that employing β∗F instead of β−∗F and β+∗F, generates the second highest adjusted R2 among all regressions at 0.15. Therefore, it can be concluded from Table 9 that the variables produce a much higher adjusted R2 value than the variables used in Table 8, thus indicating that placing appropriate numbers of knots in the CPPF model and selecting the appropriate order in the FFF model produces better beta estimates. Moreover, all the variables in regression 8 and regression 11 show significant effects on excess stock returns, while the regression that employs β−∗F and β+∗F outperforming the one that employs β∗F.
It can be concluded from this paper that the beta, upside beta and downside beta estimates produced by the CPPF model and the FFF model do have a significant impact on cross-sectional stock returns. The beta estimates, whose role has been doubted in the literature for several decades, are significant in driving the stock returns for both models. Moreover, the downside and upside beta estimates of both models demonstrate reversed impacts on stock returns. The reason for that is when downside beta is calculated, the return of the market portfolio is below the average, and likely to be negative. The expected excess stock return is the product of beta and excess returns to the market portfolio, so when stocks are sorted by downside beta into portfolios, the larger downside beta, the lower the return, and vice versa. The former ones show negative impacts on stock returns, while the latter ones, consistent with the beta estimates, have positive effects (both are significant). For stocks with negative downside beta, they are inversely related with downside risk and more desirable in a downside market, therefore positive returns are rewarded. Moreover, placing the appropriate number of knots in the CPPF model and selecting the correct order of the FFF model are crucial procedures to generate the best fit estimates according to the AIC. It has been shown in this paper that estimates with the appropriate number of knots (or order) deliver more significant impacts on stock returns within the cross-sectional return regressions with respect to those based on non-optimal knots or orders. Furthermore, in order to avoid potential multicollinearity, beta estimates based on the FFF model β∗F can be treated as an alternative variable of downside and upside beta estimates (β−∗F and β+∗F, respectively). However, employing β−∗F and β+∗F in the regression produces higher adjusted R2 values than employing β∗F.
The findings of this paper fill the gap in literature, specifically on both investigating and pricing the time variation and asymmetric nature of systematic risk. The methods and models proposed in this paper embed advanced mathematical techniques of data smoothing and widen the current options of asset pricing models. The application of proposed models provides high degree of explanatory power to capture and price risk in stock market. Nonetheless, as an experiment, there are certain limitations of this paper, e.g. the number of knots and orders are determined relatively arbitrarily, and the sample could be extended.
The authors of this paper are indebted to all colleagues who provided valuable comments.
The authors declare no conflict of interest.
[1] | K. Diethelm, The Analysis of Fractional Differential Equations, Springer-Verlag, 2010. https://doi.org/10.1007/978-3-642-14574-2 |
[2] | A. Kilbas, H Srivastava, J. Trujillo, Theory and Applications of Fractional Differential Equations, Elsevier, 2006. https://doi.org/10.3182/20060719-3-PT-4902.00008 |
[3] | K. Miller, B. Ross, An Introduction to the Fractional Calculus and Fractional Differential Equations, Wiley, 1993. |
[4] | I. Podlubny, Fractional Differential Equations, Acad. Press, 1999. |
[5] |
R. Garra, F. Mainardi, G. Spada, A generalization of the Lomnitz logarithmic creep law via Hadamard fractional calculus, Chaos Solitons Fractals, 102 (2017), 333–338. https://doi.org/10.1016/j.chaos.2017.03.032 doi: 10.1016/j.chaos.2017.03.032
![]() |
[6] |
C. Lomnitz, Application of the logarithmic creep law to stress wave attenuation in the solid earth, J. Geophys. Res., 67 (1962), 365–368. https://doi.org/10.1029/JZ067i001p00365 doi: 10.1029/JZ067i001p00365
![]() |
[7] |
S. Denisov, H. Kantz, Continuous-time random walk theory of superslow diffusion, Europhys. Lett., 92 (2010), 30001. https://doi.org/10.1209/0295-5075/92/30001 doi: 10.1209/0295-5075/92/30001
![]() |
[8] |
J. Dr¨aGer, J. Klafter, Strong anomaly in diffusion generated by iterated maps, Phys. Rev. Lett., 84 (2000), 5998–6001. https://doi.org/10.1103/PhysRevLett.84.5998 doi: 10.1103/PhysRevLett.84.5998
![]() |
[9] |
F. Iglˊoi, L. Turban, H. Rieger, Anomalous diffusion in aperiodic environments, Phys. Rev. E, 59 (1999), 1465. https://doi.org/10.1103/PhysRevE.59.1465 doi: 10.1103/PhysRevE.59.1465
![]() |
[10] |
L. Sanders, M. Lomholt, L. Lizana, K. Fogelmark, R. Metzler, T. Ambjornsson, Severe slowing-down and universality of the dynamics in disordered interacting many-body systems: Ageing and ultraslow diffusion, New J. Phys, 16 (2014), 113050. https://doi.org/10.1088/1367-2630/16/11/113050 doi: 10.1088/1367-2630/16/11/113050
![]() |
[11] |
T. Sandev, A. Iomin, H. Kantz, R. Metzler, A. Chechkin, Comb model with slow and ultraslow diffusion, Math. Model. Nat. Phenom., 11 (2016), 18–33. https://doi.org/10.1051/mmnp/201611302 doi: 10.1051/mmnp/201611302
![]() |
[12] |
C. Li, D. Li, J. Wang, L1/LDG method for the generalized time-fractional Burgers equation, Commun. Appl. Math. Comput., 5 (2023), 1299–1322. https://doi.org/10.1016/j.matcom.2021.03.005 doi: 10.1016/j.matcom.2021.03.005
![]() |
[13] |
E. Fan, C. Li, Z. Li, Numerical approaches to Caputo-Hadamard fractional derivatives with applications to long-term integration of fractional differential systems, Commun. Nonlinear Sci. Numer. Simul., 106 (2022), 106096. https://doi.org/10.1016/j.cnsns.2021.106096 doi: 10.1016/j.cnsns.2021.106096
![]() |
[14] |
T. Zhao, C. Li, D Li, Efficient spectral collocation method for fractional differential equation with Caputo-Hadamard derivative, Fract. Calc. Appl. Anal., 26 (2023), 2903–2927. https://doi.org/10.1007/s13540-023-00216-6 doi: 10.1007/s13540-023-00216-6
![]() |
[15] |
X. Ye, J. Cao, C. Xu, A high order scheme for fractional differential equations with the Caputo-Hadamard derivative, J. Comput. Math., 434 (2025), 615–640. https://doi.org/10.4208/jcm.2312-m2023-0098 doi: 10.4208/jcm.2312-m2023-0098
![]() |
[16] |
H. Liao, D Li, J Zhang, Sharp error estimate of the nonuniform L1 formula for linear reaction-subdiffusion equations, SIAM J. Numer. Anal., 56 (2018), 1112–1133. https://doi.org/10.1137/17M1131829 doi: 10.1137/17M1131829
![]() |
[17] |
M. Stynesy, E. O'riordanz, Z. Graciax, Error analysis of a finite difference method on graded meshes for a time-fractional diffusion equation, SIAM J. Numer. Anal., 55 (2017), 1057–1079. https://doi.org/10.1137/16M1082329 doi: 10.1137/16M1082329
![]() |
[18] |
H. Liao, W Mclean, J. Zhang, A second-order scheme with nonuniform time steps for a linear reaction-sudiffusion problem, Commun. Comput. Phys., 30 (2021), 567–601. https://doi.org/10.4208/cicp.OA-2020-0124 doi: 10.4208/cicp.OA-2020-0124
![]() |
[19] |
H. Liao, W Mclean, J. Zhang, A discrete Gronwall inequality with applications to numerical schemes for subdiffusion problems, SIAM J. Numer. Anal., 57 (2019), 218–237. https://doi.org/10.1137/16M1175742 doi: 10.1137/16M1175742
![]() |
[20] |
B. Zhou, X. Chen, D. Li, Nonuniform Alikhanov linearized Galerkin finite element methods for nonlinear time-fractional parabolic equations, J. Sci. Comput., 85 (2020), 39. https://doi.org/10.1007/s10915-020-01350-6 doi: 10.1007/s10915-020-01350-6
![]() |
[21] |
X. Yang, H. Zhang, The uniform l1 long-time behavior of time discretization for time-fractional partial differential equations with nonsmooth data, Appl. Math. Lett., 124 (2022), 107644. https://doi.org/10.1016/j.aml.2021.107644 doi: 10.1016/j.aml.2021.107644
![]() |
[22] |
D. Hou, M. Azaiez, C. Xu, M¨untz spectral method for two-dimensional space-fractional convection-diffusion equation, Commun. Comput. Phys., 26 (2019), 1415–1443. https://doi.org/10.4208/cicp.2019.js60.04 doi: 10.4208/cicp.2019.js60.04
![]() |
[23] |
X. Yang, W. Qiu, H. Chen, H. Zhang, Second-order BDF ADI Galerkin finite element method for the evolutionary equation with a nonlocal term in three-dimensional space, Appl. Numer. Math., 172 (2022), 497–513. https://doi.org/10.1016/j.apnum.2021.11.004 doi: 10.1016/j.apnum.2021.11.004
![]() |
[24] |
Z. Sun, L. Ling, A kernel-based meshless conservative Galerkin method for solving Hamiltonian wave equations, SIAM J. Sci. Comput., 44 (2022), A2789–A2807. https://doi.org/10.1137/21M1436981 doi: 10.1137/21M1436981
![]() |
[25] |
H. Zhang, Y. Liu, X. Yang, An efficient ADI difference scheme for the nonlocal evolution problem in three-dimensional space, J. Appl. Math. Comput., 69 (2023), 651–674. https://doi.org/10.1007/s12190-022-01760-9 doi: 10.1007/s12190-022-01760-9
![]() |
[26] |
X. Yang, Z. Zhang, Superconvergence analysis of a robust orthogonal Gauss collocation method for 2D fourth-order subdiffusion equations, J. Sci. Comput., 100 (2024), 62. https://doi.org/10.1007/s10915-024-02616-z doi: 10.1007/s10915-024-02616-z
![]() |
[27] |
Z. Sun, L. Ling, A high-order meshless linearly implicit energy-preserving method for nonlinear wave equations on {R}iemannian manifolds, SIAM J. Sci. Comput., 46 (2024), A3779–A3802. https://doi.org/10.1137/24M1654245 doi: 10.1137/24M1654245
![]() |
[28] |
C. Li, Z. Li, Z. Wang, Mathematical analysis and the local discontinuous Galerkin method for Caputo-Hadamard fractional partial differential equation, J. Sci. Comput., 85 (2020), 41. https://doi.org/10.1007/s10915-020-01353-3 doi: 10.1007/s10915-020-01353-3
![]() |
[29] |
Z. Wang, C. Ou, S. Vong, A second-order scheme with nonuniform time grids for Caputo-Hadamard fractional sub-diffusion equations, J. Comput. Appl. Math., 414 (2022), 114448. https://doi.org/10.1016/j.cam.2022.114448 doi: 10.1016/j.cam.2022.114448
![]() |
[30] | J. Zhang, Z. Zhang, C. Zhao, α-robust error estimates of general non-uniform time-step numerical schemes for reaction-subdiffusion problems, preprint, arXiv: 2305.07383. https://doi.org/10.48550/arXiv.2305.07383 |
[31] |
Z. Yang, X. Zheng, H. Wang, Well-posedness and regularity of Caputo-Hadamard fractional stochastic differential equations, Z. Angew. Math. Phys., 72 (2021), 141. https://doi.org/10.1007/s00033-021-01566-y doi: 10.1007/s00033-021-01566-y
![]() |
1. | Fangzhou Huang, Jiao Song, Nick J. Taylor, The impact of business conditions and commodity market on US stock returns: An asset pricing modelling experiment, 2022, 6, 2573-0134, 433, 10.3934/QFE.2022019 |
Number of knots | 0 | 1 | 2 | 3 | 4 | 5 |
Placement points | 50% | 33.3% | 25% | 20% | 16.6% | |
66.6% | 50% | 40% | 33.30% | |||
75% | 60% | 50% | ||||
80% | 66.6% | |||||
83.3% |
Knots | 0 | 1 | 2 | 3 | 4 | 5 | |
β∗S | Number of Stocks | 9409 | 1455 | 861 | 655 | 614 | 563 |
Percentage to Whole sample | 69.40% | 10.73% | 6.35% | 4.83% | 4.53% | 4.15% | |
β−∗S and β+∗S | Number of Stocks | 9399 | 903 | 483 | 416 | 729 | 1627 |
Percentage to Whole sample | 69.33% | 6.66% | 3.56% | 3.07% | 5.38% | 12.00% | |
Note: This table reports the number and percentage of stocks with different knots to construct the best fit estimates of CPPF model. |
Order | 1 | 2 | 3 | 4 | |
β∗F | Number of Stocks | 6204 | 2746 | 2099 | 2508 |
Percentage to Whole sample | 45.76% | 20.26% | 15.48% | 18.50% | |
β−∗F and β+∗F | Number of Stocks | 8377 | 2293 | 1429 | 1458 |
Percentage to Whole sample | 61.79% | 16.91% | 10.54% | 10.75% | |
Note: This table reports the number and percentage of stocks in different orders to construct the best fit estimates of the FFF model. |
Portfolio | Return | βS* | βS−* | βS+* | Portfolio | Return | βS* | βS−* | βS+* | ||||||||
Panel 1 Stocks Sorted by βS* | Panel 2 Stocks Sorted by βS−* | ||||||||||||||||
1 Low | 1.53% | −0.24 | 0.47 | −0.21 | 1 Low | 25.00% | 0.48 | −7.2 | 0.85 | ||||||||
2 | 6.76% | 0.64 | 0.80 | 0.32 | 2 | 11.08% | 0.77 | 0.68 | 0.68 | ||||||||
3 | 8.95% | 0.99 | 1.23 | 0.84 | 3 | 10.53% | 1.03 | 1.1 | 0.83 | ||||||||
4 | 11.59% | 1.41 | 1.54 | 1.25 | 4 | 9.73% | 1.33 | 1.61 | 0.99 | ||||||||
5 High | 24.84% | 2.68 | 2.18 | 2.4 | 5 High | −2.67% | 1.86 | 10.01 | 1.23 | ||||||||
High − Low | 23.30% | 2.92 | 1.72 | 2.62 | High − Low | −27.67% | 1.38 | 17.21 | 0.39 | ||||||||
Panel 3 Stocks Sorted by βS+* | Panel 4 Stocks Sorted by (βS−*−βS*) | ||||||||||||||||
1 Low | −11.93% | 0.42 | 1.32 | −3.87 | 1 Low | 35.87% | 1.64 | −6.6 | 2.03 | ||||||||
2 | 6.67% | 0.71 | 0.91 | 0.34 | 2 | 14.74% | 1.11 | 0.88 | 1.09 | ||||||||
3 | 9.85% | 0.96 | 1.21 | 0.81 | 3 | 9.72% | 0.96 | 1.05 | 0.72 | ||||||||
4 | 13.43% | 1.28 | 1.35 | 1.37 | 4 | 5.72% | 0.95 | 1.38 | 0.6 | ||||||||
5 High | 35.58% | 2.11 | 1.42 | 5.94 | 5 High | −12.35% | 0.83 | 9.49 | 0.16 | ||||||||
High − Low | 47.51% | 1.69 | 0.1 | 9.81 | High − Low | −48.22% | −0.81 | 16.09 | −1.87 | ||||||||
Panel 5 Stocks Sorted by (βS+*−βS*) | Panel 6 Stocks Sorted by (βS−*−βS+*) | ||||||||||||||||
1 Low | −6.16% | 1.35 | 1.88 | −3.38 | 1 Low | 36.79% | 1.43 | −6.32 | 5.09 | ||||||||
2 | 7.24% | 1.03 | 1.37 | 0.52 | 2 | 14.15% | 1.05 | 0.89 | 1.2 | ||||||||
3 | 9.94% | 0.96 | 1.16 | 0.84 | 3 | 10.07% | 0.96 | 1.05 | 0.83 | ||||||||
4 | 13.12% | 1 | 1.1 | 1.21 | 4 | 6.43% | 0.98 | 1.35 | 0.53 | ||||||||
5 High | 29.51% | 1.14 | 0.71 | 5.4 | 5 High | −13.76% | 1.05 | 9.23 | −3.06 | ||||||||
High − Low | 35.67% | −0.2 | −1.18 | 8.78 | High − Low | −50.55% | −0.38 | 15.55 | −8.15 | ||||||||
Note: This table presents the relationship between excess stock returns and factor loading of the CPPF model. The column labeled "return" reports the average stock returns over a one month T-bill rate. "High-Low" reports the difference between portfolio 5 and portfolio 1. Notably, in a perfect market, the average value of βS* is assumed to be 1. |
Portfolio | Return | βF* | βF−* | βF+* | Portfolio | Return | βF* | βF−* | βF+* | ||||||||
Panel 1 Stocks Sorted by βF* | Panel 2 Stocks Sorted by βF−* | ||||||||||||||||
1 Low | 3.54% | −1.19 | −0.74 | −0.87 | 1 Low | 21.88% | −0.67 | −1.32 | −0.31 | ||||||||
2 | 7.51% | −0.37 | −0.26 | −0.28 | 2 | 13.09% | −0.24 | −0.38 | −0.12 | ||||||||
3 | 9.11% | −0.01 | −0.01 | −0.03 | 3 | 10.46% | −0.01 | −0.01 | −0.01 | ||||||||
4 | 11.24% | 0.36 | 0.24 | 0.23 | 4 | 8.22% | 0.23 | 0.37 | 0.11 | ||||||||
5 High | 22.26% | 1.22 | 0.76 | 0.95 | 5 High | −0.03% | 0.69 | 1.33 | 0.33 | ||||||||
High − Low | 18.72% | 2.41 | 1.49 | 1.82 | High − Low | −21.91% | 1.36 | 2.65 | 0.64 | ||||||||
Panel 3 Stocks Sorted by βF+* | Panel 4 Stocks Sorted by (βF−*−βF*) | ||||||||||||||||
1 Low | −7.72% | −0.81 | −0.29 | −1.28 | 1 Low | 33.18% | 0.52 | −0.71 | 0.65 | ||||||||
2 | 5.02% | −0.27 | −0.13 | −0.36 | 2 | 16.29% | 0.11 | −0.18 | 0.19 | ||||||||
3 | 9.44% | −0.02 | −0.01 | −0.02 | 3 | 10.34% | −0.02 | −0.02 | −0.01 | ||||||||
4 | 14.58% | 0.24 | 0.13 | 0.33 | 4 | 4.66% | −0.14 | 0.17 | −0.21 | ||||||||
5 High | 32.33% | 0.86 | 0.3 | 1.33 | 5 High | −10.84% | −0.48 | 0.74 | −0.61 | ||||||||
High − Low | 40.05% | 1.67 | 0.59 | 2.61 | High − Low | −44.02% | −1 | 1.45 | −1.27 | ||||||||
Panel 5 Stocks Sorted by (βF+*−βF*) | Panel 6 Stocks Sorted by (βF−*−βF+*) | ||||||||||||||||
1 Low | −4.27% | 0.36 | 0.57 | −0.63 | 1 Low | 36.12% | 0.17 | −0.81 | 0.85 | ||||||||
2 | 5.48% | 0.12 | 0.21 | −0.14 | 2 | 16.44% | −0.02 | −0.24 | 0.17 | ||||||||
3 | 9.74% | 0 | −0.01 | −0.01 | 3 | 10.15% | −0.02 | −0.01 | −0.02 | ||||||||
4 | 14.99% | −0.14 | −0.23 | 0.12 | 4 | 4.34% | −0.02 | 0.22 | −0.2 | ||||||||
5 High | 27.69% | −0.34 | −0.54 | 0.66 | 5 High | −13.41% | −0.12 | 0.85 | −0.81 | ||||||||
High − Low | 31.95% | −0.7 | −1.1 | 1.29 | High − Low | −49.53% | −0.3 | 1.66 | −1.66 | ||||||||
Note: This table presents the relationship between excess stock returns and factor loading of the FFF model. The column labeled "return" reports the average stock returns over one month T-bill rate. "High-Low" reports the difference between portfolio 5 and portfolio 1. Notably, in a perfect market, the average value of βF* depends on the average value of the intercept. |
βS0 | β−S0 | β+S0 | βF1 | β−F1 | β+F1 | |
βS0 | 1.0000 | |||||
β−S0 | 0.3363 | 1.0000 | ||||
β+S0 | 0.4326 | 0.0365 | 1.0000 | |||
βF1 | 0.4397 | 0.1345 | 0.1770 | 1.0000 | ||
β−F1 | 0.3209 | 0.2611 | 0.0458 | 0.7331 | 1.0000 | |
β+F1 | 0.4620 | 0.0336 | 0.2737 | 0.8093 | 0.2886 | 1.0000 |
Note: This table reports the correlation coefficients between factor loadings of both the CPPF model with zero knots and the FFF model in order one. To avoid unnecessary repetition, only the lower triangle of the matrix is shown. |
β∗S | β−∗S | β+∗S | β∗F | β−∗F | β+∗F | |
β∗S | 1.0000 | |||||
β−∗S | 0.0130 | 1.0000 | ||||
β+∗S | 0.0237 | -0.0001 | 1.0000 | |||
β∗F | 0.3154 | 0.0029 | 0.0085 | 1.0000 | ||
β−∗F | 0.2609 | 0.0082 | 0.0025 | 0.5239 | 1.0000 | |
β+∗F | 0.3099 | -0.0018 | 0.0098 | 0.6689 | 0.1889 | 1.0000 |
Note: This table reports the correlation coefficients between factor loadings of both the CPPF model with appropriate number of knots and the FFF model in appropriate order, to avoid unnecessary repetition, only the lower triangle of the matrix is shown. |
1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 | 11 | |
Intercept |
0.00155 [1.14] |
0.00743***
[3.82] |
0.00154 [1.21] |
0.00852***
[3.48] |
0.00907***
[3.90] |
0.00212 [1.57] |
0.00705***
[3.88] |
0.00213*
[1.70] |
0.00345***
[2.89] |
0.00801***
[4.06] |
0.00251**
[2.24] |
βS0 | 0.00663*** [3.25] |
0.00824***
[3.68] |
0.00614***
[3.02] |
0.00768***
[3.46] |
0.00525***
[2.73] |
0.00626***
[2.98] |
|||||
β−S0 | 0.0000637 [0.08] |
−0.000721 [−1.28] | 0.00187***
[2.99] |
−0.000678 [−1.20] | −0.0119*** [−13.04] | −0.000103 [−0.15] | 0.00147***
[2.95] |
||||
β+S0 | 0.000937*
[1.90] |
−0.000977*** [−2.13] | −0.000884** [−2.21] | −0.000983** [−2.14] | 0.0115***
[14.41] |
0.000520 [1.14] |
−0.00205*** [−4.72] | ||||
βF1 | 0.00542***
[3.35] |
0.00132**
[2.50] |
0.00141***
[2.66] |
0.00509***
[3.77] |
|||||||
β−F1 | −0.0113*** [−9.43] | −0.0126*** [−12.50] | −0.0131*** [−15.27] | ||||||||
β+F1 | 0.0148***
[10.83] |
0.0156***
[12.60] |
0.0127***
[17.27] |
||||||||
Number of obs | 2398103 | 2398103 | 2398103 | 2398103 | 2398103 | 2398103 | 2398103 | 2398103 | 2398103 | 2398103 | 2398103 |
Average R2 | 0.055 | 0.060 | 0.098 | 0.017 | 0.029 | 0.056 | 0.077 | 0.100 | 0.065 | 0.069 | 0.106 |
Note: This table reports the results of the Fama-Macbeth regression of factor loadings without knot and order. The t-statistics in the square brackets are calculated by using Newey-West (1987) heteroscedastic robust standard error with 12 lags. * denotes significance at the 10% level, ** denotes significance at the 5% level and ***denotes significance at the 1% level. |
1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 | 11 | |
Intercept |
0.00186 [1.52] |
0.00898***
[7.05] |
0.00538***
[5.04] |
0.00826***
[3.56] |
0.00943***
[4.42] |
0.00279**
[2.19] |
0.00965***
[7.54] |
0.00618***
[5.76] |
0.00369***
[3.38] |
0.00947***
[7.14] |
0.00623***
[6.23] |
β∗S | 0.00610***
[3.42] |
0.00448***
[3.82] |
0.00481***
[3.00] |
0.00383***
[3.54] |
0.00490***
[3.20] |
0.00384***
[3.42] |
|||||
β−∗S | −0.0120***
[−9.82] |
−0.00948***
[−11.19] |
−0.0106***
[−10.25] |
−0.00944***
[−11.19] |
−0.00994***
[−12.94] |
−0.00850***
[−11.30] |
|||||
β+∗S | 0.0134***
[14.47] |
0.00911***
[14.41] |
0.0117***
[15.25] |
0.00890***
[14.63] |
0.0107***
[15.74] |
0.00823***
[14.44] |
|||||
β∗F | 0.00673***
[4.16] |
0.00334***
[5.67] |
0.00238***
[4.73] |
0.00342***
[3.91] |
|||||||
β−∗F | −0.0113***
[−9.13] |
−0.00419***
[−8.21] |
−0.0111***
[−10.53] |
−0.00485***
[−13.74] |
|||||||
β+∗F | 0.0155***
[13.00] |
0.00647***
[11.34] |
0.0121***
[15.64] |
0.00586***
[13.51] |
|||||||
Number of obs |
2398103 | 2398103 | 2398103 | 2398103 | 2398103 | 2398103 | 2398103 | 2398103 | 2398103 | 2398103 | 2398103 |
Adjusted R2 | 0.092 | 0.13 | 0.147 | 0.039 | 0.055 | 0.096 | 0.138 | 0.150 | 0.11 | 0.137 | 0.153 |
Note: This table reports the results of the Fama-Macbeth regression of factor loadings with appropriate knots and orders on stock excess returns. The t-statistics in the square brackets are calculated by using the Newey-West (1987) heteroscedastic robust standard error with 12 lags. * denotes significance at the 10% level, ** denotes significance at the 5% level and ***denotes significance at the 1% level. |
Number of knots | 0 | 1 | 2 | 3 | 4 | 5 |
Placement points | 50% | 33.3% | 25% | 20% | 16.6% | |
66.6% | 50% | 40% | 33.30% | |||
75% | 60% | 50% | ||||
80% | 66.6% | |||||
83.3% |
Knots | 0 | 1 | 2 | 3 | 4 | 5 | |
β∗S | Number of Stocks | 9409 | 1455 | 861 | 655 | 614 | 563 |
Percentage to Whole sample | 69.40% | 10.73% | 6.35% | 4.83% | 4.53% | 4.15% | |
β−∗S and β+∗S | Number of Stocks | 9399 | 903 | 483 | 416 | 729 | 1627 |
Percentage to Whole sample | 69.33% | 6.66% | 3.56% | 3.07% | 5.38% | 12.00% | |
Note: This table reports the number and percentage of stocks with different knots to construct the best fit estimates of CPPF model. |
Order | 1 | 2 | 3 | 4 | |
β∗F | Number of Stocks | 6204 | 2746 | 2099 | 2508 |
Percentage to Whole sample | 45.76% | 20.26% | 15.48% | 18.50% | |
β−∗F and β+∗F | Number of Stocks | 8377 | 2293 | 1429 | 1458 |
Percentage to Whole sample | 61.79% | 16.91% | 10.54% | 10.75% | |
Note: This table reports the number and percentage of stocks in different orders to construct the best fit estimates of the FFF model. |
Portfolio | Return | βS* | βS−* | βS+* | Portfolio | Return | βS* | βS−* | βS+* | ||||||||
Panel 1 Stocks Sorted by βS* | Panel 2 Stocks Sorted by βS−* | ||||||||||||||||
1 Low | 1.53% | −0.24 | 0.47 | −0.21 | 1 Low | 25.00% | 0.48 | −7.2 | 0.85 | ||||||||
2 | 6.76% | 0.64 | 0.80 | 0.32 | 2 | 11.08% | 0.77 | 0.68 | 0.68 | ||||||||
3 | 8.95% | 0.99 | 1.23 | 0.84 | 3 | 10.53% | 1.03 | 1.1 | 0.83 | ||||||||
4 | 11.59% | 1.41 | 1.54 | 1.25 | 4 | 9.73% | 1.33 | 1.61 | 0.99 | ||||||||
5 High | 24.84% | 2.68 | 2.18 | 2.4 | 5 High | −2.67% | 1.86 | 10.01 | 1.23 | ||||||||
High − Low | 23.30% | 2.92 | 1.72 | 2.62 | High − Low | −27.67% | 1.38 | 17.21 | 0.39 | ||||||||
Panel 3 Stocks Sorted by βS+* | Panel 4 Stocks Sorted by (βS−*−βS*) | ||||||||||||||||
1 Low | −11.93% | 0.42 | 1.32 | −3.87 | 1 Low | 35.87% | 1.64 | −6.6 | 2.03 | ||||||||
2 | 6.67% | 0.71 | 0.91 | 0.34 | 2 | 14.74% | 1.11 | 0.88 | 1.09 | ||||||||
3 | 9.85% | 0.96 | 1.21 | 0.81 | 3 | 9.72% | 0.96 | 1.05 | 0.72 | ||||||||
4 | 13.43% | 1.28 | 1.35 | 1.37 | 4 | 5.72% | 0.95 | 1.38 | 0.6 | ||||||||
5 High | 35.58% | 2.11 | 1.42 | 5.94 | 5 High | −12.35% | 0.83 | 9.49 | 0.16 | ||||||||
High − Low | 47.51% | 1.69 | 0.1 | 9.81 | High − Low | −48.22% | −0.81 | 16.09 | −1.87 | ||||||||
Panel 5 Stocks Sorted by (βS+*−βS*) | Panel 6 Stocks Sorted by (βS−*−βS+*) | ||||||||||||||||
1 Low | −6.16% | 1.35 | 1.88 | −3.38 | 1 Low | 36.79% | 1.43 | −6.32 | 5.09 | ||||||||
2 | 7.24% | 1.03 | 1.37 | 0.52 | 2 | 14.15% | 1.05 | 0.89 | 1.2 | ||||||||
3 | 9.94% | 0.96 | 1.16 | 0.84 | 3 | 10.07% | 0.96 | 1.05 | 0.83 | ||||||||
4 | 13.12% | 1 | 1.1 | 1.21 | 4 | 6.43% | 0.98 | 1.35 | 0.53 | ||||||||
5 High | 29.51% | 1.14 | 0.71 | 5.4 | 5 High | −13.76% | 1.05 | 9.23 | −3.06 | ||||||||
High − Low | 35.67% | −0.2 | −1.18 | 8.78 | High − Low | −50.55% | −0.38 | 15.55 | −8.15 | ||||||||
Note: This table presents the relationship between excess stock returns and factor loading of the CPPF model. The column labeled "return" reports the average stock returns over a one month T-bill rate. "High-Low" reports the difference between portfolio 5 and portfolio 1. Notably, in a perfect market, the average value of βS* is assumed to be 1. |
Portfolio | Return | βF* | βF−* | βF+* | Portfolio | Return | βF* | βF−* | βF+* | ||||||||
Panel 1 Stocks Sorted by βF* | Panel 2 Stocks Sorted by βF−* | ||||||||||||||||
1 Low | 3.54% | −1.19 | −0.74 | −0.87 | 1 Low | 21.88% | −0.67 | −1.32 | −0.31 | ||||||||
2 | 7.51% | −0.37 | −0.26 | −0.28 | 2 | 13.09% | −0.24 | −0.38 | −0.12 | ||||||||
3 | 9.11% | −0.01 | −0.01 | −0.03 | 3 | 10.46% | −0.01 | −0.01 | −0.01 | ||||||||
4 | 11.24% | 0.36 | 0.24 | 0.23 | 4 | 8.22% | 0.23 | 0.37 | 0.11 | ||||||||
5 High | 22.26% | 1.22 | 0.76 | 0.95 | 5 High | −0.03% | 0.69 | 1.33 | 0.33 | ||||||||
High − Low | 18.72% | 2.41 | 1.49 | 1.82 | High − Low | −21.91% | 1.36 | 2.65 | 0.64 | ||||||||
Panel 3 Stocks Sorted by βF+* | Panel 4 Stocks Sorted by (βF−*−βF*) | ||||||||||||||||
1 Low | −7.72% | −0.81 | −0.29 | −1.28 | 1 Low | 33.18% | 0.52 | −0.71 | 0.65 | ||||||||
2 | 5.02% | −0.27 | −0.13 | −0.36 | 2 | 16.29% | 0.11 | −0.18 | 0.19 | ||||||||
3 | 9.44% | −0.02 | −0.01 | −0.02 | 3 | 10.34% | −0.02 | −0.02 | −0.01 | ||||||||
4 | 14.58% | 0.24 | 0.13 | 0.33 | 4 | 4.66% | −0.14 | 0.17 | −0.21 | ||||||||
5 High | 32.33% | 0.86 | 0.3 | 1.33 | 5 High | −10.84% | −0.48 | 0.74 | −0.61 | ||||||||
High − Low | 40.05% | 1.67 | 0.59 | 2.61 | High − Low | −44.02% | −1 | 1.45 | −1.27 | ||||||||
Panel 5 Stocks Sorted by (βF+*−βF*) | Panel 6 Stocks Sorted by (βF−*−βF+*) | ||||||||||||||||
1 Low | −4.27% | 0.36 | 0.57 | −0.63 | 1 Low | 36.12% | 0.17 | −0.81 | 0.85 | ||||||||
2 | 5.48% | 0.12 | 0.21 | −0.14 | 2 | 16.44% | −0.02 | −0.24 | 0.17 | ||||||||
3 | 9.74% | 0 | −0.01 | −0.01 | 3 | 10.15% | −0.02 | −0.01 | −0.02 | ||||||||
4 | 14.99% | −0.14 | −0.23 | 0.12 | 4 | 4.34% | −0.02 | 0.22 | −0.2 | ||||||||
5 High | 27.69% | −0.34 | −0.54 | 0.66 | 5 High | −13.41% | −0.12 | 0.85 | −0.81 | ||||||||
High − Low | 31.95% | −0.7 | −1.1 | 1.29 | High − Low | −49.53% | −0.3 | 1.66 | −1.66 | ||||||||
Note: This table presents the relationship between excess stock returns and factor loading of the FFF model. The column labeled "return" reports the average stock returns over one month T-bill rate. "High-Low" reports the difference between portfolio 5 and portfolio 1. Notably, in a perfect market, the average value of βF* depends on the average value of the intercept. |
βS0 | β−S0 | β+S0 | βF1 | β−F1 | β+F1 | |
βS0 | 1.0000 | |||||
β−S0 | 0.3363 | 1.0000 | ||||
β+S0 | 0.4326 | 0.0365 | 1.0000 | |||
βF1 | 0.4397 | 0.1345 | 0.1770 | 1.0000 | ||
β−F1 | 0.3209 | 0.2611 | 0.0458 | 0.7331 | 1.0000 | |
β+F1 | 0.4620 | 0.0336 | 0.2737 | 0.8093 | 0.2886 | 1.0000 |
Note: This table reports the correlation coefficients between factor loadings of both the CPPF model with zero knots and the FFF model in order one. To avoid unnecessary repetition, only the lower triangle of the matrix is shown. |
β∗S | β−∗S | β+∗S | β∗F | β−∗F | β+∗F | |
β∗S | 1.0000 | |||||
β−∗S | 0.0130 | 1.0000 | ||||
β+∗S | 0.0237 | -0.0001 | 1.0000 | |||
β∗F | 0.3154 | 0.0029 | 0.0085 | 1.0000 | ||
β−∗F | 0.2609 | 0.0082 | 0.0025 | 0.5239 | 1.0000 | |
β+∗F | 0.3099 | -0.0018 | 0.0098 | 0.6689 | 0.1889 | 1.0000 |
Note: This table reports the correlation coefficients between factor loadings of both the CPPF model with appropriate number of knots and the FFF model in appropriate order, to avoid unnecessary repetition, only the lower triangle of the matrix is shown. |
1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 | 11 | |
Intercept |
0.00155 [1.14] |
0.00743***
[3.82] |
0.00154 [1.21] |
0.00852***
[3.48] |
0.00907***
[3.90] |
0.00212 [1.57] |
0.00705***
[3.88] |
0.00213*
[1.70] |
0.00345***
[2.89] |
0.00801***
[4.06] |
0.00251**
[2.24] |
βS0 | 0.00663*** [3.25] |
0.00824***
[3.68] |
0.00614***
[3.02] |
0.00768***
[3.46] |
0.00525***
[2.73] |
0.00626***
[2.98] |
|||||
β−S0 | 0.0000637 [0.08] |
−0.000721 [−1.28] | 0.00187***
[2.99] |
−0.000678 [−1.20] | −0.0119*** [−13.04] | −0.000103 [−0.15] | 0.00147***
[2.95] |
||||
β+S0 | 0.000937*
[1.90] |
−0.000977*** [−2.13] | −0.000884** [−2.21] | −0.000983** [−2.14] | 0.0115***
[14.41] |
0.000520 [1.14] |
−0.00205*** [−4.72] | ||||
βF1 | 0.00542***
[3.35] |
0.00132**
[2.50] |
0.00141***
[2.66] |
0.00509***
[3.77] |
|||||||
β−F1 | −0.0113*** [−9.43] | −0.0126*** [−12.50] | −0.0131*** [−15.27] | ||||||||
β+F1 | 0.0148***
[10.83] |
0.0156***
[12.60] |
0.0127***
[17.27] |
||||||||
Number of obs | 2398103 | 2398103 | 2398103 | 2398103 | 2398103 | 2398103 | 2398103 | 2398103 | 2398103 | 2398103 | 2398103 |
Average R2 | 0.055 | 0.060 | 0.098 | 0.017 | 0.029 | 0.056 | 0.077 | 0.100 | 0.065 | 0.069 | 0.106 |
Note: This table reports the results of the Fama-Macbeth regression of factor loadings without knot and order. The t-statistics in the square brackets are calculated by using Newey-West (1987) heteroscedastic robust standard error with 12 lags. * denotes significance at the 10% level, ** denotes significance at the 5% level and ***denotes significance at the 1% level. |
1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 | 11 | |
Intercept |
0.00186 [1.52] |
0.00898***
[7.05] |
0.00538***
[5.04] |
0.00826***
[3.56] |
0.00943***
[4.42] |
0.00279**
[2.19] |
0.00965***
[7.54] |
0.00618***
[5.76] |
0.00369***
[3.38] |
0.00947***
[7.14] |
0.00623***
[6.23] |
β∗S | 0.00610***
[3.42] |
0.00448***
[3.82] |
0.00481***
[3.00] |
0.00383***
[3.54] |
0.00490***
[3.20] |
0.00384***
[3.42] |
|||||
β−∗S | −0.0120***
[−9.82] |
−0.00948***
[−11.19] |
−0.0106***
[−10.25] |
−0.00944***
[−11.19] |
−0.00994***
[−12.94] |
−0.00850***
[−11.30] |
|||||
β+∗S | 0.0134***
[14.47] |
0.00911***
[14.41] |
0.0117***
[15.25] |
0.00890***
[14.63] |
0.0107***
[15.74] |
0.00823***
[14.44] |
|||||
β∗F | 0.00673***
[4.16] |
0.00334***
[5.67] |
0.00238***
[4.73] |
0.00342***
[3.91] |
|||||||
β−∗F | −0.0113***
[−9.13] |
−0.00419***
[−8.21] |
−0.0111***
[−10.53] |
−0.00485***
[−13.74] |
|||||||
β+∗F | 0.0155***
[13.00] |
0.00647***
[11.34] |
0.0121***
[15.64] |
0.00586***
[13.51] |
|||||||
Number of obs |
2398103 | 2398103 | 2398103 | 2398103 | 2398103 | 2398103 | 2398103 | 2398103 | 2398103 | 2398103 | 2398103 |
Adjusted R2 | 0.092 | 0.13 | 0.147 | 0.039 | 0.055 | 0.096 | 0.138 | 0.150 | 0.11 | 0.137 | 0.153 |
Note: This table reports the results of the Fama-Macbeth regression of factor loadings with appropriate knots and orders on stock excess returns. The t-statistics in the square brackets are calculated by using the Newey-West (1987) heteroscedastic robust standard error with 12 lags. * denotes significance at the 10% level, ** denotes significance at the 5% level and ***denotes significance at the 1% level. |