similar to: Optimization and garch

Hello R-list-members, I'm trying to model ARMA(0,2) & GARCH(1,1) process using the code below, but according to my textbook, the estimated parameters are wrong. The MA-parameters should be negative. (I've got the same problem using garchFit()). Can anyone tell me what I'm doing wrong? And how can I calculate the hessian matrix? Many thanks, Desislava Kavrakova Code:

Hello, I'm quite new to R but tried to learn as much as possible in the last few months. My problem is that I would like to estimate the model of Leon et al. (2005). I have shortly summarised the most important equations in the following pdf file: http://hannes.fedorapeople.org/leon2005.pdf My main question is now how could I introduce these two additional terms into the Likelihood

I'm trying to estimate the parameters for GARCH(1,1) process. Here's my code: loglikelihood <-function(theta) { h=((r[1]-theta[1])^2) p=0 for (t in 2:length(r)) { h=c(h,theta[2]+theta[3]*((r[t-1]-theta[1])^2)+theta[4]*h[t-1]) p=c(p,dnorm(r[t],theta[1],sqrt(h[t]),log=TRUE)) } -sum(p) } Then I use nlminb to minimize the function loglikelihood: nlminb(

Dear R-users, I am wondering if anyone can tell me how Garch coefficients are calculated in R. What is algorithm for that? If anyone give me a detail desceiption I will be very grateful. Thanks and Regards, stat700004 thanks in advance --------------------------------- [[alternative HTML version deleted]]

Hello, i want to estimate a complex GARCH-model (see below). http://r.789695.n4.nabble.com/file/n4112396/GJR_Garch.png W stands for the Day of the Week Dummies. r stands for returns of stock market indices. I stands for the GJR-term. I need some help with three problems: 1.) implementation of the GJR-term in the variance equation 2.) compute robust covariance matrix

I tried to use LME (on a fairly large dataset, so I am not including it), and I got this error message: Error in lme.formula(formula(paste(c(toString(TargetName), "as.factor(nodeInd)"), : nlminb problem, convergence error code = 1 message = false convergence (8) Is there any way to get more information or to get the potentially wrong estimates from LME? (Also, the page in the

Hello everyone, I'm trying to estimate the parameters of the returns series attached using the GARCH code below, but I get the following error message: Error in solve.default(Hessian) : system is computationally singular: reciprocal condition number = 0 Error in diag(solve(Hessian)) : error in evaluating the argument 'x' in selecting a method for function 'diag' Can

Hello, I have a rather complex problem... I will have to explain everything in detail because I cannot solve it by myself...i just ran out of ideas. So here is what I want to do: I take quotes of two indices - S&P500 and DJ. And my first aim is to estimate coefficients of the DCC-GARCH model for them. This is how I do it: library(tseries) p1 = get.hist.quote(instrument =

Hi All, I,m trying to do a GARCH simulation in R 2.3.0 release in Windows XP. I've seen garchsim function but that is for garch (1,1) and ?garch gives an example for ARCH simulation. Can anyone help me how can i extend the help shown in ?garch to GARCH simulation? Please help me in this regard. Thanks, Sumanta Basak.

I'm using R 2.1.1 and just successfully installed packages tseries, fseries. I try to run example http://www.maths.lth.se/help/R/.R/library/tseries/html/garch.html But it shows > x.arch <- garch(x, order = c(0,2)) # Fit ARCH(2) Error: couldn't find function "garch" Then I run command > help.search("garch") it shows the R information.

Dear R users, I am trying to use tseries' garch function in order to determine the volatility of a return series generated by quantmod. Here is the code that I am using: > library(quantmod) > getSymbols("AAPL") convert daily closing prices into continuous log returns > dret<-dailyReturn(AAPL,type='log') check to see that the autocorrelations decay >

Dear R users, Could you please help me out. I am in trouble as I am unable to model graphs to explain the GARCH (1 , 1) model, the Hill estimator (of alpha), and the Pareto estimator. I just got introduce to R. I am working on a paper which must be worked from R. You look at the difficulty I had from the text below. [1] "DAX" "DAX_CAC" "DAX_CAC40"

I need to change the code of Garch to the FCGARCH (a non-linear multi-regime GARCH). I don't know nothing about R. I'd like to know how can I get the code of the garch in order to change it and make the fit for the FC-GARCH. Any non-linear code will be helpfull because if doesn't help in the programming it helps in getting familiar with R. Thank you Renato -- PhD Student Renato

Dear r users, Few days ago I posted the same topic but unable to receive any suggestion. So I am asking this same question. I was trying to fit a garch(1,1) model to my dataset. But while executing I got a warning message "NaNs produced in: sqrt(pred$e)". And got the estimated sd's along with five "NA", but as per my best knowledge I should get only one

I have done the usual estimation of GARCH models, applied to my historical dataset (commodities futures) with a maximum likelihood function and selected the best model on the basis of information criteria such as Akaike and Bayes. Can somebody explain me please the calibration scheme for a GARCH model? I was not able to find a paper, dealing with exactly this algorithm for my case. I only

tseries_0.3-0 at CRAN now contains the following new features: NelPlo Nelson-Plosser Macroeconomic Time Series garch Fit GARCH Models to Time Series get.hist.quote Download Historical Finance Data jarque.bera.test Jarque-Bera Test na.remove NA Handling Routines for Time Series garch contains a GARCH estimation routine together

tseries_0.3-0 at CRAN now contains the following new features: NelPlo Nelson-Plosser Macroeconomic Time Series garch Fit GARCH Models to Time Series get.hist.quote Download Historical Finance Data jarque.bera.test Jarque-Bera Test na.remove NA Handling Routines for Time Series garch contains a GARCH estimation routine together

Dear All, I have loaded package(tseries), but when I run predict.garch(...) R tells me could not find function "predict.garch", however ?predict.garch shows me something. I am confused about this. How can I forecast garch volatility? I have tried: predict(...,n.ahead=...),give me fitted value predict(...,n),give me NA,NA

Dear all R-users, I have a GARCH related query. Suppose I fit a GARCH(1,1) model on a dataframe dat >garch1 = garch(dat) >summary(garch1) Call: garch(x = dat) Model: GARCH(1,1) Residuals: Min 1Q Median 3Q Max -4.7278 -0.3240 0.0000 0.3107 12.3981 Coefficient(s): Estimate Std. Error t value Pr(>|t|) a0 1.212e-04 2.053e-06 59.05 <2e-16 *** a1

Dear all R-users, I have a GARCH related query. Suppose I fit a GARCH(1,1) model on a dataframe dat >garch1 = garch(dat) >summary(garch1) Call: garch(x = dat) Model: GARCH(1,1) Residuals: Min 1Q Median 3Q Max -4.7278 -0.3240 0.0000 0.3107 12.3981 Coefficient(s): Estimate Std. Error t value Pr(>|t|) a0 1.212e-04 2.053e-06 59.05 <2e-16 *** a1