

Good Morning,
I investigate the dayof –theweek of 11 indices in 11 countries using EGARCH model and rolling sample method with one sample interval case of 500 days. The main purpose of the empirical analysis is to determine whether there exist calendar effects or not by observing whether the t values exceed the quantile of the standard normal distribution z = 1.96 for confidence level of 5%. I want to use the same method as Zhang, Lai, Lin (2017), “the dayoftheweek effects of stocks markets in different countries”, Finance Research Letters 20.
For the EGARCH model, I used Package “Rugarch”, and everything is good :
spec2b<ugarchspec(variance.model=list(model="eGARCH",garchOrder=c(1,1), external.regressors=L1), mean.model=list(armaOrder=c(0,0),archm=TRUE,archpow=1,external.regressors=L2),
distribution.model = "std")
fit2b<ugarchfit(spec=spec2b,solver="hybrid", data=RENT)
show(fit2b)
But for the rolling sample method, with EGARCH(1,1), I thought using “ugarchroll” or “ugarchdistribution ”, to calculate the student t for each rolling window , but I do not get student t :
dist < ugarchdistribution(fit2b, n.sim = 500, n.start = 1, m.sim = 100, solver = "solnp")
show(dist)
or
roll1=ugarchroll(spec1b,data = RENT,n.start = NULL,window.size = 500,solver = "solnp" )
show(roll1)
Could you orient me to a solution ?
Thank you very much
Sandrine BOULERNE
Maître de Conférences  Habilitation à Diriger des Recherches (HDR)
IAE de Tours  Faculté de Droit, Économie et Sciences Sociales
50 av. Jean Portalis, BP 0607, 37206 TOURS Cedex 3
Téléphone : 02 47 36 10 47
www.iae.univtours.fr
[[alternative HTML version deleted]]
_______________________________________________
[hidden email] mailing list
https://stat.ethz.ch/mailman/listinfo/rsigfinance Subscriberposting only. If you want to post, subscribe first.
 Also note that this is not the rhelp list where general R questions should go.


Have you tried 'coef(roll1)'?
I think the documentation is reasonably clear (?ugarchroll).
It returns an object of class 'uGARCHroll' (which is clickable in the
documentation and will take you to the class help page).
In the class help page you will see that one of the methods on the class
is:
"coef
signature(object = "uGARCHroll"): Extracts the list of coefficients for
each estimated window in the rolling backtest."
Therefore, coef(roll1) will return the list of coefficient matrices
per roll window, from which you can extract the shape parameter of the
standardized Student distribution.
Similarly, please read the documentation for ugarchdistribution.
Alexios
On 9/11/2017 9:03 AM, Sandrine Boulerne wrote:
> Good Morning,
>
> I investigate the dayof –theweek of 11 indices in 11 countries using
> EGARCH model and rolling sample method with one sample interval case of
> 500 days. The main purpose of the empirical analysis is to determine
> whether there exist calendar effects or not by observing whether the /t/
> values exceed the quantile of the standard normal distribution z = 1.96
> for confidence level of 5%. I want to use the same method as Zhang, Lai,
> Lin (2017), “the dayoftheweek effects of stocks markets in different
> countries”, Finance Research Letters 20.
>
> For the EGARCH model, I used Package “Rugarch”, and everything is good:
>
> spec2b<ugarchspec(variance.model=list(model="eGARCH",garchOrder=c(1,1),
> external.regressors=L1),
> mean.model=list(armaOrder=c(0,0),archm=TRUE,archpow=1,external.regressors=L2),
>
> distribution.model = "std")
>
> fit2b<ugarchfit(spec=spec2b,solver="hybrid", data=RENT)
>
> show(fit2b)
>
> But for the rolling sample method, with EGARCH(1,1), I thought using
> “ugarchroll” or “ugarchdistribution*”, ** to calculate the student /t/
> for each rolling window*, *but **I do not get student /t/* :
>
> dist < ugarchdistribution(fit2b, n.sim = 500, n.start = 1, m.sim = 100,
> solver = "solnp")
>
> show(dist)
>
> or
>
> roll1=ugarchroll(spec1b,data = RENT,n.start = NULL,window.size =
> 500,solver = "solnp" )
>
> show(roll1)
>
> Could you orient me to a solution ?
>
> Thank you very much
>
>
>
> Sandrine BOULERNE
> Maître de Conférences  Habilitation à Diriger des Recherches (HDR)
> IAE de Tours  Faculté de Droit, Économie et Sciences Sociales
> 50 av. Jean Portalis, BP 0607, 37206 TOURS Cedex 3
> Téléphone : 02 47 36 10 47
> www.iae.univtours.fr
>
_______________________________________________
[hidden email] mailing list
https://stat.ethz.ch/mailman/listinfo/rsigfinance Subscriberposting only. If you want to post, subscribe first.
 Also note that this is not the rhelp list where general R questions should go.


Dear Alexios
That's wonderful. Thanks to you I got what I was looking for several weeks. you are right, the documentation is very clear, but sometimes we do not see what may seem obvious to an expert eye.
Thanks a lot
Sandrine BOULERNE
Maitre de Conferences  Habilitation à Diriger des Recherches (HDR)
Coresponsable du thème Finance Contrôle Audit du Laboratoire de recherche VALLOREM
IAE de Tours  Faculté de Droit, Économie et Sciences Sociales
50 av. Jean Portalis, BP 0607, 37206 TOURS Cedex 3
Téléphone : 02 47 36 10 47
www.iae.univtours.fr
 Mail original 
De: "alexios" < [hidden email]>
À: "Sandrine Boulerne" < [hidden email]>
Cc: [hidden email]
Envoyé: Mardi 12 Septembre 2017 06:27:33
Objet: Re: To obtain the t student of each rolling window with EGARCH model
Have you tried 'coef(roll1)'?
I think the documentation is reasonably clear (?ugarchroll).
It returns an object of class 'uGARCHroll' (which is clickable in the
documentation and will take you to the class help page).
In the class help page you will see that one of the methods on the class
is:
"coef
signature(object = "uGARCHroll"): Extracts the list of coefficients for
each estimated window in the rolling backtest."
Therefore, coef(roll1) will return the list of coefficient matrices
per roll window, from which you can extract the shape parameter of the
standardized Student distribution.
Similarly, please read the documentation for ugarchdistribution.
Alexios
On 9/11/2017 9:03 AM, Sandrine Boulerne wrote:
> Good Morning,
>
> I investigate the dayof –theweek of 11 indices in 11 countries using
> EGARCH model and rolling sample method with one sample interval case of
> 500 days. The main purpose of the empirical analysis is to determine
> whether there exist calendar effects or not by observing whether the /t/
> values exceed the quantile of the standard normal distribution z = 1.96
> for confidence level of 5%. I want to use the same method as Zhang, Lai,
> Lin (2017), “the dayoftheweek effects of stocks markets in different
> countries”, Finance Research Letters 20.
>
> For the EGARCH model, I used Package “Rugarch”, and everything is good:
>
> spec2b<ugarchspec(variance.model=list(model="eGARCH",garchOrder=c(1,1),
> external.regressors=L1),
> mean.model=list(armaOrder=c(0,0),archm=TRUE,archpow=1,external.regressors=L2),
>
> distribution.model = "std")
>
> fit2b<ugarchfit(spec=spec2b,solver="hybrid", data=RENT)
>
> show(fit2b)
>
> But for the rolling sample method, with EGARCH(1,1), I thought using
> “ugarchroll” or “ugarchdistribution*”, ** to calculate the student /t/
> for each rolling window*, *but **I do not get student /t/* :
>
> dist < ugarchdistribution(fit2b, n.sim = 500, n.start = 1, m.sim = 100,
> solver = "solnp")
>
> show(dist)
>
> or
>
> roll1=ugarchroll(spec1b,data = RENT,n.start = NULL,window.size =
> 500,solver = "solnp" )
>
> show(roll1)
>
> Could you orient me to a solution ?
>
> Thank you very much
>
>
>
> Sandrine BOULERNE
> Maître de Conférences  Habilitation à Diriger des Recherches (HDR)
> IAE de Tours  Faculté de Droit, Économie et Sciences Sociales
> 50 av. Jean Portalis, BP 0607, 37206 TOURS Cedex 3
> Téléphone : 02 47 36 10 47
> www.iae.univtours.fr
>
_______________________________________________
[hidden email] mailing list
https://stat.ethz.ch/mailman/listinfo/rsigfinance Subscriberposting only. If you want to post, subscribe first.
 Also note that this is not the rhelp list where general R questions should go.


Hello Alexios
I have advanced in my research and now I will want to get the student t from each past rolling window (500 days for each window) from 2001 to 2016 and not the student t from each forcasting rolling window.
My question has evolved and is following your response of September 12th.
For the EGARCH model, I used Package “Rugarch” :
spec2b<ugarchspec(variance.model=list(model="eGARCH",garchOrder=c(1,1),external.regressors=L1), mean.model=list(armaOrder=c(0,0),archm=TRUE,archpow=1,external.regressors=L2), distribution.model = "std")
fit2b<ugarchfit(spec=spec2b,solver="hybrid", data=RENT)
show(fit2b)
roll1=ugarchroll(spec2b,data = RENT,n.start = NULL,window.size=500,solver = "solnp" )
coeff(roll1)
coef(roll1) allows me to get the student t from the 500 forecasting windows but not the student t from past windows (500 days per windows) from 2001 to 2016.
Could you orient me to a solution ? I am looking for a solution for several days.
Thank you very much
Sandrine BOULERNE
University of Tours
www.iae.univtours.fr
 Mail original 
De: "alexios" < [hidden email]>
À: "Sandrine Boulerne" < [hidden email]>
Cc: [hidden email]
Envoyé: Mardi 12 Septembre 2017 06:27:33
Objet: Re: To obtain the t student of each rolling window with EGARCH model
Have you tried 'coef(roll1)'?
I think the documentation is reasonably clear (?ugarchroll).
It returns an object of class 'uGARCHroll' (which is clickable in the
documentation and will take you to the class help page).
In the class help page you will see that one of the methods on the class
is:
"coef
signature(object = "uGARCHroll"): Extracts the list of coefficients for
each estimated window in the rolling backtest."
Therefore, coef(roll1) will return the list of coefficient matrices
per roll window, from which you can extract the shape parameter of the
standardized Student distribution.
Similarly, please read the documentation for ugarchdistribution.
Alexios
On 9/11/2017 9:03 AM, Sandrine Boulerne wrote:
> Good Morning,
>
> I investigate the dayof –theweek of 11 indices in 11 countries using
> EGARCH model and rolling sample method with one sample interval case of
> 500 days. The main purpose of the empirical analysis is to determine
> whether there exist calendar effects or not by observing whether the /t/
> values exceed the quantile of the standard normal distribution z = 1.96
> for confidence level of 5%. I want to use the same method as Zhang, Lai,
> Lin (2017), “the dayoftheweek effects of stocks markets in different
> countries”, Finance Research Letters 20.
>
> For the EGARCH model, I used Package “Rugarch”, and everything is good:
>
> spec2b<ugarchspec(variance.model=list(model="eGARCH",garchOrder=c(1,1),
> external.regressors=L1),
> mean.model=list(armaOrder=c(0,0),archm=TRUE,archpow=1,external.regressors=L2),
>
> distribution.model = "std")
>
> fit2b<ugarchfit(spec=spec2b,solver="hybrid", data=RENT)
>
> show(fit2b)
>
> But for the rolling sample method, with EGARCH(1,1), I thought using
> “ugarchroll” or “ugarchdistribution*”, ** to calculate the student /t/
> for each rolling window*, *but **I do not get student /t/* :
>
> dist < ugarchdistribution(fit2b, n.sim = 500, n.start = 1, m.sim = 100,
> solver = "solnp")
>
> show(dist)
>
> or
>
> roll1=ugarchroll(spec1b,data = RENT,n.start = NULL,window.size =
> 500,solver = "solnp" )
>
> show(roll1)
>
> Could you orient me to a solution ?
>
> Thank you very much
>
>
>
> Sandrine BOULERNE
> Maître de Conférences  Habilitation à Diriger des Recherches (HDR)
> IAE de Tours  Faculté de Droit, Économie et Sciences Sociales
> 50 av. Jean Portalis, BP 0607, 37206 TOURS Cedex 3
> Téléphone : 02 47 36 10 47
> www.iae.univtours.fr
>
_______________________________________________
[hidden email] mailing list
https://stat.ethz.ch/mailman/listinfo/rsigfinance Subscriberposting only. If you want to post, subscribe first.
 Also note that this is not the rhelp list where general R questions should go.


Hi Sandrine,
Not sure if this is what you are after:
#####
# reproducible example:
library(rugarch)
library(xts)
data(sp500ret)
spx = xts(as.numeric(sp500ret[,1]), as.Date(rownames(sp500ret)))
spec2b<ugarchspec(variance.model=list(model="eGARCH",garchOrder=c(1,1)),mean.model=list(armaOrder=c(0,0)),
distribution.model = "std")
roll1=ugarchroll(spec2b,data=spx,n.start=NULL,window.size=500,solver="solnp")
cf=coef(roll1)
st=do.call(rbind, lapply(1:length(cf), function(i){
xts(cf[[i]]$coef["shape",1], as.Date(cf[[i]]$index))
}))
st = na.locf(cbind(spx[paste0(index(st[1]),"/",st[length(st)])], st))[,2]
#####
Alexios
On 11/23/2017 12:30 PM, Sandrine Boulerne wrote:
> Hello Alexios
>
> I have advanced in my research and now I will want to get the student t from each past rolling window (500 days for each window) from 2001 to 2016 and not the student t from each forcasting rolling window.
> My question has evolved and is following your response of September 12th.
>
> For the EGARCH model, I used Package “Rugarch” :
>
> spec2b<ugarchspec(variance.model=list(model="eGARCH",garchOrder=c(1,1),external.regressors=L1), mean.model=list(armaOrder=c(0,0),archm=TRUE,archpow=1,external.regressors=L2), distribution.model = "std")
> fit2b<ugarchfit(spec=spec2b,solver="hybrid", data=RENT)
> show(fit2b)
> roll1=ugarchroll(spec2b,data = RENT,n.start = NULL,window.size=500,solver = "solnp" )
> coeff(roll1)
>
> coef(roll1) allows me to get the student t from the 500 forecasting windows but not the student t from past windows (500 days per windows) from 2001 to 2016.
>
> Could you orient me to a solution ? I am looking for a solution for several days.
>
> Thank you very much
>
> Sandrine BOULERNE
> University of Tours
> www.iae.univtours.fr
>
>
>  Mail original 
> De: "alexios" < [hidden email]>
> À: "Sandrine Boulerne" < [hidden email]>
> Cc: [hidden email]
> Envoyé: Mardi 12 Septembre 2017 06:27:33
> Objet: Re: To obtain the t student of each rolling window with EGARCH model
>
> Have you tried 'coef(roll1)'?
>
> I think the documentation is reasonably clear (?ugarchroll).
> It returns an object of class 'uGARCHroll' (which is clickable in the
> documentation and will take you to the class help page).
> In the class help page you will see that one of the methods on the class
> is:
>
> "coef
> signature(object = "uGARCHroll"): Extracts the list of coefficients for
> each estimated window in the rolling backtest."
>
> Therefore, coef(roll1) will return the list of coefficient matrices
> per roll window, from which you can extract the shape parameter of the
> standardized Student distribution.
>
> Similarly, please read the documentation for ugarchdistribution.
>
>
> Alexios
>
> On 9/11/2017 9:03 AM, Sandrine Boulerne wrote:
>> Good Morning,
>>
>> I investigate the dayof –theweek of 11 indices in 11 countries using
>> EGARCH model and rolling sample method with one sample interval case of
>> 500 days. The main purpose of the empirical analysis is to determine
>> whether there exist calendar effects or not by observing whether the /t/
>> values exceed the quantile of the standard normal distribution z = 1.96
>> for confidence level of 5%. I want to use the same method as Zhang, Lai,
>> Lin (2017), “the dayoftheweek effects of stocks markets in different
>> countries”, Finance Research Letters 20.
>>
>> For the EGARCH model, I used Package “Rugarch”, and everything is good:
>>
>> spec2b<ugarchspec(variance.model=list(model="eGARCH",garchOrder=c(1,1),
>> external.regressors=L1),
>> mean.model=list(armaOrder=c(0,0),archm=TRUE,archpow=1,external.regressors=L2),
>>
>> distribution.model = "std")
>>
>> fit2b<ugarchfit(spec=spec2b,solver="hybrid", data=RENT)
>>
>> show(fit2b)
>>
>> But for the rolling sample method, with EGARCH(1,1), I thought using
>> “ugarchroll” or “ugarchdistribution*”, ** to calculate the student /t/
>> for each rolling window*, *but **I do not get student /t/* :
>>
>> dist < ugarchdistribution(fit2b, n.sim = 500, n.start = 1, m.sim = 100,
>> solver = "solnp")
>>
>> show(dist)
>>
>> or
>>
>> roll1=ugarchroll(spec1b,data = RENT,n.start = NULL,window.size =
>> 500,solver = "solnp" )
>>
>> show(roll1)
>>
>> Could you orient me to a solution ?
>>
>> Thank you very much
>>
>>
>>
>> Sandrine BOULERNE
>> Maître de Conférences  Habilitation à Diriger des Recherches (HDR)
>> IAE de Tours  Faculté de Droit, Économie et Sciences Sociales
>> 50 av. Jean Portalis, BP 0607, 37206 TOURS Cedex 3
>> Téléphone : 02 47 36 10 47
>> www.iae.univtours.fr
>>
>
_______________________________________________
[hidden email] mailing list
https://stat.ethz.ch/mailman/listinfo/rsigfinance Subscriberposting only. If you want to post, subscribe first.
 Also note that this is not the rhelp list where general R questions should go.


Hi Alexios,
Thank for your example with sp500ret, but we get only the t values from 500 rolling windows. I will try to explain myself with a simpler example for me.
If the samples include 7 data: {1, 2, 3, 4, 5, 6, 7} and we fix the rolling interval with length of 4, the first sample subinterval is {1, 2, 3, 4}; the second sample subinterval is {2, 3, 4, 5}, etc. The length of the rolling sample interval used in my paper is 500 days. It is assumed that the average trading days in one year is about 250 days. Therefore, 500 days represent the number of trading days in 2 years,
I investigate the dayoftheweek effects of the 10 indices in 10 countries using GARCH model and rolling sample method with one sample interval case of 500 days, from 2001 to 2016.
The main purpose of the empirical analysis in my paper is to determine whether there exist calendar effects or not by observing whether the t values in the graph exceed the quantile of the standard normal distribution, z(α), for a given confidence level α. For example, z(0.005) = 1.96.
So, I would like to get 3750 t values (15 years * 500 days in 2 years), in using :
spec2b<ugarchspec(variance.model=list(model="eGARCH",garchOrder=c(1,1)),mean.model=list(armaOrder=c(0,0)),
distribution.model = "std")
roll1=ugarchroll(spec2b,data=spx,n.start=NULL,window.size=500,solver="solnp") #perhaps ???
cf=coef(roll1) # from 2001 to 2016
Thank you very much for your help
Sandrine
 Mail original 
De: "alexios" < [hidden email]>
À: "Sandrine Boulerne" < [hidden email]>
Cc: [hidden email]
Envoyé: Vendredi 24 Novembre 2017 04:08:12
Objet: Re: To obtain the t student of each rolling window with EGARCH model
Hi Sandrine,
Not sure if this is what you are after:
#####
# reproducible example:
library(rugarch)
library(xts)
data(sp500ret)
spx = xts(as.numeric(sp500ret[,1]), as.Date(rownames(sp500ret)))
spec2b<ugarchspec(variance.model=list(model="eGARCH",garchOrder=c(1,1)),mean.model=list(armaOrder=c(0,0)),
distribution.model = "std")
roll1=ugarchroll(spec2b,data=spx,n.start=NULL,window.size=500,solver="solnp")
cf=coef(roll1)
st=do.call(rbind, lapply(1:length(cf), function(i){
xts(cf[[i]]$coef["shape",1], as.Date(cf[[i]]$index))
}))
st = na.locf(cbind(spx[paste0(index(st[1]),"/",st[length(st)])], st))[,2]
#####
Alexios
On 11/23/2017 12:30 PM, Sandrine Boulerne wrote:
> Hello Alexios
>
> I have advanced in my research and now I will want to get the student t from each past rolling window (500 days for each window) from 2001 to 2016 and not the student t from each forcasting rolling window.
> My question has evolved and is following your response of September 12th.
>
> For the EGARCH model, I used Package “Rugarch” :
>
> spec2b<ugarchspec(variance.model=list(model="eGARCH",garchOrder=c(1,1),external.regressors=L1), mean.model=list(armaOrder=c(0,0),archm=TRUE,archpow=1,external.regressors=L2), distribution.model = "std")
> fit2b<ugarchfit(spec=spec2b,solver="hybrid", data=RENT)
> show(fit2b)
> roll1=ugarchroll(spec2b,data = RENT,n.start = NULL,window.size=500,solver = "solnp" )
> coeff(roll1)
>
> coef(roll1) allows me to get the student t from the 500 forecasting windows but not the student t from past windows (500 days per windows) from 2001 to 2016.
>
> Could you orient me to a solution ? I am looking for a solution for several days.
>
> Thank you very much
>
> Sandrine BOULERNE
> University of Tours
> www.iae.univtours.fr
>
>
>  Mail original 
> De: "alexios" < [hidden email]>
> À: "Sandrine Boulerne" < [hidden email]>
> Cc: [hidden email]
> Envoyé: Mardi 12 Septembre 2017 06:27:33
> Objet: Re: To obtain the t student of each rolling window with EGARCH model
>
> Have you tried 'coef(roll1)'?
>
> I think the documentation is reasonably clear (?ugarchroll).
> It returns an object of class 'uGARCHroll' (which is clickable in the
> documentation and will take you to the class help page).
> In the class help page you will see that one of the methods on the class
> is:
>
> "coef
> signature(object = "uGARCHroll"): Extracts the list of coefficients for
> each estimated window in the rolling backtest."
>
> Therefore, coef(roll1) will return the list of coefficient matrices
> per roll window, from which you can extract the shape parameter of the
> standardized Student distribution.
>
> Similarly, please read the documentation for ugarchdistribution.
>
>
> Alexios
>
> On 9/11/2017 9:03 AM, Sandrine Boulerne wrote:
>> Good Morning,
>>
>> I investigate the dayof –theweek of 11 indices in 11 countries using
>> EGARCH model and rolling sample method with one sample interval case of
>> 500 days. The main purpose of the empirical analysis is to determine
>> whether there exist calendar effects or not by observing whether the /t/
>> values exceed the quantile of the standard normal distribution z = 1.96
>> for confidence level of 5%. I want to use the same method as Zhang, Lai,
>> Lin (2017), “the dayoftheweek effects of stocks markets in different
>> countries”, Finance Research Letters 20.
>>
>> For the EGARCH model, I used Package “Rugarch”, and everything is good:
>>
>> spec2b<ugarchspec(variance.model=list(model="eGARCH",garchOrder=c(1,1),
>> external.regressors=L1),
>> mean.model=list(armaOrder=c(0,0),archm=TRUE,archpow=1,external.regressors=L2),
>>
>> distribution.model = "std")
>>
>> fit2b<ugarchfit(spec=spec2b,solver="hybrid", data=RENT)
>>
>> show(fit2b)
>>
>> But for the rolling sample method, with EGARCH(1,1), I thought using
>> “ugarchroll” or “ugarchdistribution*”, ** to calculate the student /t/
>> for each rolling window*, *but **I do not get student /t/* :
>>
>> dist < ugarchdistribution(fit2b, n.sim = 500, n.start = 1, m.sim = 100,
>> solver = "solnp")
>>
>> show(dist)
>>
>> or
>>
>> roll1=ugarchroll(spec1b,data = RENT,n.start = NULL,window.size =
>> 500,solver = "solnp" )
>>
>> show(roll1)
>>
>> Could you orient me to a solution ?
>>
>> Thank you very much
>>
>>
>>
>> Sandrine BOULERNE
>> Maître de Conférences  Habilitation à Diriger des Recherches (HDR)
>> IAE de Tours  Faculté de Droit, Économie et Sciences Sociales
>> 50 av. Jean Portalis, BP 0607, 37206 TOURS Cedex 3
>> Téléphone : 02 47 36 10 47
>> www.iae.univtours.fr
>>
>
_______________________________________________
[hidden email] mailing list
https://stat.ethz.ch/mailman/listinfo/rsigfinance Subscriberposting only. If you want to post, subscribe first.
 Also note that this is not the rhelp list where general R questions should go.


 Mail original 
De: "Sandrine Boulerne" < [hidden email]>
À: "alexios" < [hidden email]>
Cc: [hidden email]
Envoyé: Vendredi 24 Novembre 2017 10:55:46
Objet: Re: To obtain the t student of each rolling window with EGARCH model
Hi Alexios,
Thank for your example with sp500ret, but we get only the t values from 500 rolling windows. I will try to explain myself with a simpler example for me.
If the samples include 7 data: {1, 2, 3, 4, 5, 6, 7} and we fix the rolling interval with length of 4, the first sample subinterval is {1, 2, 3, 4}; the second sample subinterval is {2, 3, 4, 5}, etc. The length of the rolling sample interval used in my paper is 500 days. It is assumed that the average trading days in one year is about 250 days. Therefore, 500 days represent the number of trading days in 2 years,
I investigate the dayoftheweek effects of the 10 indices in 10 countries using GARCH model and rolling sample method with one sample interval case of 500 days, from 2001 to 2016.
The main purpose of the empirical analysis in my paper is to determine whether there exist calendar effects or not by observing whether the t values in the graph exceed the quantile of the standard normal distribution, z(α), for a given confidence level α. For example, z(0.005) = 1.96.
So, I would like to get 3750 t values (15 years * 500 days in 2 years), in using :
spec2b<ugarchspec(variance.model=list(model="eGARCH",garchOrder=c(1,1)),mean.model=list(armaOrder=c(0,0)),
distribution.model = "std")
roll1=ugarchroll(spec2b,data=spx,n.start=NULL,window.size=500,solver="solnp") #perhaps ???
cf=coef(roll1) # from 2001 to 2016
Thank you very much for your help
Sandrine
 Mail original 
De: "alexios" < [hidden email]>
À: "Sandrine Boulerne" < [hidden email]>
Cc: [hidden email]
Envoyé: Vendredi 24 Novembre 2017 04:08:12
Objet: Re: To obtain the t student of each rolling window with EGARCH model
Hi Sandrine,
Not sure if this is what you are after:
#####
# reproducible example:
library(rugarch)
library(xts)
data(sp500ret)
spx = xts(as.numeric(sp500ret[,1]), as.Date(rownames(sp500ret)))
spec2b<ugarchspec(variance.model=list(model="eGARCH",garchOrder=c(1,1)),mean.model=list(armaOrder=c(0,0)),
distribution.model = "std")
roll1=ugarchroll(spec2b,data=spx,n.start=NULL,window.size=500,solver="solnp")
cf=coef(roll1)
st=do.call(rbind, lapply(1:length(cf), function(i){
xts(cf[[i]]$coef["shape",1], as.Date(cf[[i]]$index))
}))
st = na.locf(cbind(spx[paste0(index(st[1]),"/",st[length(st)])], st))[,2]
#####
Alexios
On 11/23/2017 12:30 PM, Sandrine Boulerne wrote:
> Hello Alexios
>
> I have advanced in my research and now I will want to get the student t from each past rolling window (500 days for each window) from 2001 to 2016 and not the student t from each forcasting rolling window.
> My question has evolved and is following your response of September 12th.
>
> For the EGARCH model, I used Package “Rugarch” :
>
> spec2b<ugarchspec(variance.model=list(model="eGARCH",garchOrder=c(1,1),external.regressors=L1), mean.model=list(armaOrder=c(0,0),archm=TRUE,archpow=1,external.regressors=L2), distribution.model = "std")
> fit2b<ugarchfit(spec=spec2b,solver="hybrid", data=RENT)
> show(fit2b)
> roll1=ugarchroll(spec2b,data = RENT,n.start = NULL,window.size=500,solver = "solnp" )
> coeff(roll1)
>
> coef(roll1) allows me to get the student t from the 500 forecasting windows but not the student t from past windows (500 days per windows) from 2001 to 2016.
>
> Could you orient me to a solution ? I am looking for a solution for several days.
>
> Thank you very much
>
> Sandrine BOULERNE
> University of Tours
> www.iae.univtours.fr
>
>
>  Mail original 
> De: "alexios" < [hidden email]>
> À: "Sandrine Boulerne" < [hidden email]>
> Cc: [hidden email]
> Envoyé: Mardi 12 Septembre 2017 06:27:33
> Objet: Re: To obtain the t student of each rolling window with EGARCH model
>
> Have you tried 'coef(roll1)'?
>
> I think the documentation is reasonably clear (?ugarchroll).
> It returns an object of class 'uGARCHroll' (which is clickable in the
> documentation and will take you to the class help page).
> In the class help page you will see that one of the methods on the class
> is:
>
> "coef
> signature(object = "uGARCHroll"): Extracts the list of coefficients for
> each estimated window in the rolling backtest."
>
> Therefore, coef(roll1) will return the list of coefficient matrices
> per roll window, from which you can extract the shape parameter of the
> standardized Student distribution.
>
> Similarly, please read the documentation for ugarchdistribution.
>
>
> Alexios
>
> On 9/11/2017 9:03 AM, Sandrine Boulerne wrote:
>> Good Morning,
>>
>> I investigate the dayof –theweek of 11 indices in 11 countries using
>> EGARCH model and rolling sample method with one sample interval case of
>> 500 days. The main purpose of the empirical analysis is to determine
>> whether there exist calendar effects or not by observing whether the /t/
>> values exceed the quantile of the standard normal distribution z = 1.96
>> for confidence level of 5%. I want to use the same method as Zhang, Lai,
>> Lin (2017), “the dayoftheweek effects of stocks markets in different
>> countries”, Finance Research Letters 20.
>>
>> For the EGARCH model, I used Package “Rugarch”, and everything is good:
>>
>> spec2b<ugarchspec(variance.model=list(model="eGARCH",garchOrder=c(1,1),
>> external.regressors=L1),
>> mean.model=list(armaOrder=c(0,0),archm=TRUE,archpow=1,external.regressors=L2),
>>
>> distribution.model = "std")
>>
>> fit2b<ugarchfit(spec=spec2b,solver="hybrid", data=RENT)
>>
>> show(fit2b)
>>
>> But for the rolling sample method, with EGARCH(1,1), I thought using
>> “ugarchroll” or “ugarchdistribution*”, ** to calculate the student /t/
>> for each rolling window*, *but **I do not get student /t/* :
>>
>> dist < ugarchdistribution(fit2b, n.sim = 500, n.start = 1, m.sim = 100,
>> solver = "solnp")
>>
>> show(dist)
>>
>> or
>>
>> roll1=ugarchroll(spec1b,data = RENT,n.start = NULL,window.size =
>> 500,solver = "solnp" )
>>
>> show(roll1)
>>
>> Could you orient me to a solution ?
>>
>> Thank you very much
>>
>>
>>
>> Sandrine BOULERNE
>> Maître de Conférences  Habilitation à Diriger des Recherches (HDR)
>> IAE de Tours  Faculté de Droit, Économie et Sciences Sociales
>> 50 av. Jean Portalis, BP 0607, 37206 TOURS Cedex 3
>> Téléphone : 02 47 36 10 47
>> www.iae.univtours.fr
>>
>
_______________________________________________
[hidden email] mailing list
https://stat.ethz.ch/mailman/listinfo/rsigfinance Subscriberposting only. If you want to post, subscribe first.
 Also note that this is not the rhelp list where general R questions should go.

