similar to: Correlation structure in AR(2) and ARMA(p,q)

Hello Following some standard textbooks on ARMA(1,1)-GARCH(1,1) (e.g. Ruey Tsay's Analysis of Financial Time Series), I try to write an R program to estimate the key parameters of an ARMA(1,1)-GARCH(1,1) model for Intel's stock returns. For some random reason, I cannot decipher what is wrong with my R program. The R package fGarch already gives me the answer, but my customized function

Thank you so much all for your invaluable inputs. On Sat, Jun 4, 2011 at 3:36 AM, Patrick Burns <patrick at burns-stat.com> wrote: > A common thing to do is the Ljung-Box > test on the residuals. ?For garch it > would be the residuals squared. > > Actually for garch it should be the > rank of the squared residuals -- see >

Dear all R user, Please forgive me if my problem is too simple. Actually my problem is basically Statistical rather directly R related. Suppose I have return series ret with mean zero. And I want to fit a Garch(1,1) on this. my is r[t] = h[i]*z[t] h[t] = w + alpha*r[t-1]^2 + beta*h[t-1] I want to estimate the three parameters here; the R syntax is as follows: #

Hi all, I have some residuals from regression, and i suspect they have correlations in them... I am willing to cast the correlation into a ARMA(p, q) framework, what's the best way to identify the most suitable p, and q, and fit ARMA(p, q) model and then correct for the correlations in regression? I know there are functions in R, I have used them before, but I just want to see if I can do

I have the following time series model for prediction purposes *Loss_t = b1* Loss_(t-1) + b2*GDP_t + b3*W_(t-1)* where W_t is the usual white noise variable. So this is similar to ARMA(1,1) except that it also contains an extra predictor, GDP at time t. I have only 20 observations on each variable except GDP for which I know till 100 values. And most importantly,I have also calculated

Dear Madams/Sirs, Hello. I am using the gls function to specify an arma correlation during estimation in my model. The parameter values which I am sending the corARMA function are from a previous fit using arima. I have had some success with the method, however in other cases I get the following error from gls: "All parameters must be less than 1 in absolute value". None of

Hi I am trying to fit a time series data.It gives a AR(2) model using the ar function and ARMA(1,1) model using autoarmafit function in timsac package.How do I know which is the correct underlying model? pls help -- View this message in context: http://r.789695.n4.nabble.com/AR-vs-ARMA-model-tp4639015.html Sent from the R help mailing list archive at Nabble.com.

Hello - Here's what I'm trying to do. I want to fit a time series y with ARMA(1,1) + GARCH(1,1), there are also an exogeneous variable x which I wish to include, so the whole equation looks like: y_t - \phi y_{t-1} = \sigma_t \epsilon_t + \theta \sigma_{t-1} \epsilon_{t-1} + c x_t where \epsilon_t are i.i.d. random variables \sigma_t^2 = omega + \alpha \sigma_{t-1}^2 + \beta

Hi, I need a function to convert an ARMA process to an "infinite" AR process. I know that in the Stats package exist one to convert an ARMA to MA, but i can't find one function to convert to AR. There exists one? THANKS

Dear All, Let: dp: depth of the?river tp: temperature with respect to depth We can?have a simple scatter plot, between depth as y-axis and temperature as x-axis,?by using a?plot function as shown?below. #####################? dp <- c(1,4,3,2,5,7,9,8,9,2) tp <- 1:10 plot(tp,dp, type= 'l') ##################### ? Could someone advice me on?the way to?plot the same pair of

Dear All, Let a, b and c are?three matrices with same no. of column but different no. of row. a <- matrix(1, 1, 2) b <- matrix(2, 2, 2) c <- matrix(3, 3, 2) Could someone help me to combine these matrices together as a single matrix? Thank you Fir

Hello, I would like to compute ARMA and AR using arima-function in R. My question is: If I have Null=zero values in my data, what should I do? Remove ? or doesn't matter for ARIMA-models and I can estimate my coefficients including zero values in data in arima-function in R ? What is the better way? How to manage the data for ARIMA estimation? Thank you. [[alternative HTML version

Dear All, I was looking for the color index in?image function,?such as from topo.colors(n) and etc. but still never found it. For instance, from the help menu. ########################################### # Volcano data visualized as matrix. Need to transpose and flip # matrix horizontally. image(t(volcano)[ncol(volcano):1,]) # A prettier display of the volcano x <- 10*(1:nrow(volcano)) y

Thanks a lot to everybody that helped me out with this. Conclusions: (1) In order to edit arima in R: >fix(arima) or alternatively: >arima<-edit(arima) (2) This is not contained in the "Introduction to R" manual. (3) A "productive" fix of arima is attached (arma coefficients printed out and error catched so that it doesn't halt parent loops to search for

Dear All, I have a series of data in which the first column consist of a combination of date and time, for instance 17 April 2008 at 4.01pm, such data is recorded as: 4/17/2008 16:01 I'd like to seperate it into four different columns which consist of Day, Month,Year and Time, respectively. Could someone please advice me on this mater? Thank you, Fir

Dear All, I've been?searching for?appropriate codes to compute the rate of change and the curvature?of ?nonparametric regression model whish was denoted by a smooth function?but?unfortunately?don't manage to?do?it. I presume that such characteristics from a smooth curve can be determined by the first and second derivative operators. The following are the example of fitting a

Dear All, Could anyone?give me a hand?to suggest few packages in R to running Kalman prediction and filtration ? Thanks Fir

Dear R users, I need to fit an ARMA model. As far as I've seen, EACF (extended ACF) is not available in R. 1. Let's say I fit a series of ARMA models in a loop. Given the code/output included below, how do I pull 'Model' and 'Fit' (AIC) from each summary() so that I can combine them into an array/data frame to be sorted by AIC? 2. Apart from EACF, are you aware perhaps

Hello, I am running an ARMA model to run forecast for changes in S&P 500 prices. My ARMA calculations look as follows armacal <- arma( spdata, order = c(0,4), lag = list(ma = c(1,2,4)) ) Output: Call: arma(x = spdata, order = c(0, 4), lag = list(ma = c(1, 2, 4)) ) Coefficient(s): ma1 ma2 ma4 intercept -0.073868 0.058020 -0.081292 0.007082 All's

Hello all, I would like how to modelized a time serie with AR-ARCH process. It can be used arma and garch functions in tseries package for build ar process or a garch process, but how can it be modelized a ar-garch model ? Thanks [[alternative HTML version deleted]]