similar to: time series with NA - acf - tsdiag - Ljung-Box

Dear Everyone, The document of `tsdiag? says > These tests are sometimes applied to the residuals from an ARMA(p, q) fit, in which case the references suggest a better approximation to the null-hypothesis distribution is obtained by setting fitdf = p+q, provided of course that lag > fitdf. This implies that we should pass the `fitdf' parameter when applying `Box.test' to

Dear list members, I have 982 quotations of a given stock index and I want to run a Ljung-Box test on these data to test for autocorrelation. Later on I will estimate 8 coefficients. I do not know how many degrees of freedom should I assume in the formula for Ljung-Box test. Could anyone tell me please? Below the formula: Box.test(x, lag = ????, type = c("Ljung-Box"), fitdf = 0)

I fit a simple linear model y = bX to a data set today, and that produced 24 residuals (I have 24 data points, one for each year from 1984-2007). I would like to test the time-independence of the residuals of my model, and I was recommended by my supervisor to use the Ljung-Box test. The Box.test function in R takes 4 arguments:  x a numeric vector or univariate time series. lag the statistic

Hi, I saw that tsdiag function doesn't provide a correct result for Ljung-Box test. I want to use Box.test function for this, but I don't know how to determine lag parameter for this function. For fitdf, as I'm using a SARIMA model (0,1,1)(0,1,1)12, I decided to set it to 2. Can you confirm me the value for fitdf and give me a way to determine the lag value? Thanks Myriam -- View

HI list, I want to know whether tsdiag uses k-(p+q) as the lag in ljung box test. How is it possible to save those values nuncio -- Nuncio.M Research Scientist National Center for Antarctic and Ocean research Head land Sada Vasco da Gamma Goa-403804 [[alternative HTML version deleted]]

Hello, I'm getting the following error message when I try to run 'tsdiag' on what seems to be a valid time series: > tsdiag(small) returns: [Error in tsdiag(small) : no applicable method for "tsdiag"] where small is a little test series where I have isolated this problem (the original has 30-years worth of daily data) When I print it (small), it looks like the

I am using the ts package to fit ARIMA models, and the tsdiag() function to plot diagnostics. In doing so I'm generating an awful lot of diagnostic plots of different models and different data sets all within the same R session. So my question is, is there an option in tsdiag() similar to <main="Title"> that I can use? This would be quite helpful when I print out the plots,

two ARIMA models, both have several bars signicant in ACF and PACF plots of their residuals, but when run Ljung Box tests, both don't show any significant correlations... however, one model has p-value that is larger than the other model, based on the p-values, can I say the model with larger p-values should be better than the model with smaller p-values? [[alternative HTML version

Hi R-Helpers, I'm dealing with the STL procedure and trying to apply the tsdiag and StructTS onto the ts object to analyse the different parameters which need to be set. How can I use the tsStructure & tsdiag to create a seasonal, trend and cycle subseries plot so that I can select & analyse the correct np,ns, nt and nl? The problem is that too much signal goes into the seasonal

Hi all I'm using the Box-Ljung test (from within R) to test if a time-series in independently distributed. 2 questions: 1) p-value returned by Box-Ljung: IF I want to test if the time-series is independant at say 0.05 sig-level (it means that prob of erroneously accepting that the time-series is independent is 0.05 right?) --> then do I consider time-series as "independant"

Hi All R Experts, I met with below mentioned statistics in paper "Stock Index Volatility Forecasting with High Frequency Data" by Eugenie Hol, Siem Jan Koopman http://ideas.repec.org/p/dgr/uvatin/20020068.html I would like to ask that what is "Box-Ljung portmantacau statistic based on N squared autocorrelation" ? Is it same as "Box-Ljung Statistics" of stats

Hi! I know that there is function in fBasics package for univariate Jarque-Bera test and a funtion for univariate Ljung-Box test in stats package. But I am wondering if there is a function somewhere to do the tests for multivariate time series? Thanks, John [[alternative HTML version deleted]]

Hi, I am newbie user of ggplot and would like some assistance in implementing time series plots. I'd like to know how the tsdiag plot can be made in ggplot? Thanks Harsh Singhal Decisions Systems, Mu Sigma Inc.

The "tsdiag" function in the TSA package overrides the "tsdiag" function in the "stats" package. There are a few annoying bugs in the TSA's version of the function so I would like to use the "stats" function but still have access to other TSA functions. I have tried using stats::tsdiag( ) but as long as the TSA package is attached the function from

Hello, I am using the Ljung Box test in R to compute if the resiudals of my fitted model is random or not. I am not sure though what the results mean, I have looked at various sources on the internet and have come up with contrasting explanations (mainly because these info deal with different program languages, like SAS, SPSS, etc). I know that my residuals should appropriate white noise( is

Dear colleagues, I am new to R and statistics so please keep that in mind. I have doubts on the df calculation of Ljung-Box test (Box.test). The function seems to use always the df=lag=m and not df=m-p-q like suggested in Ljung and Box (1978) paper (that is referenced). Do you agree with this? If so, is there an R package function that computes Ljung-Box test with the degrees of

Full_Name: Peter Solymos Version: 2.8.1 OS: Windows Submission from: (NULL) (129.128.141.92) The help page of the Box.test function (stats) states that the Ljung-Box test was published in: Ljung, G. M. and Box, G. E. P. (1978), On a measure of lack of fit in time series models. Biometrika 65, 553--564. The page numbers are incorrect. The correct citation should be as follows: Ljung, G. M.

After an RSiteSeach("Box.test") I found some discussion regarding the degrees of freedom in the computation of the Ljung-Box test using Box.test(), but did not find any posting about the proper degrees of freedom. Box.test() uses "lag=number" as the degrees of freedom. However, I believe the correct degrees of freedom should be "number-p-q" where p and q are

Hi: I am trying to do a spurious regression in R but I can not find the function. Anybody used it before? The problem I have is try to do a regression with several time series. An alternative is to use the GLS function to fit the linear regression with the correlation structure AR(3) for the response (or residual). I hope the residuals after the GLS regression will be independent judged by

I'm analizing the Argentina stock market (merv) I download the data from yahoo library(tseries) Argentina <- get.hist.quote(instrument="^MERV","1996-10-08","2003-11-03", quote="Close") merv <- na.remove(log(Argentina)) I made the Augmented Dickey-Fuller test to analyse if merv have unit root: adf.test(merv,k=13) Dickey-Fuller = -1.4645,