similar to: Inputing data from multiple files as time series objects

Hello, Happy New Year!... I am encountering a problem trying to work on the data that I load in R. I have loaded to R a series of stock data using (csv files are named e.g. IBM.R) length.R <- length(list.files(".", pattern=".R")) # the number of files with one #column in the directory "./" ending to ".R" for (i in 1:length.R) {

Hello, good morning or evening!... After studying some of the examples at S-poetry Document, I tried to implement some of the concepts in my R script, that intensively uses looping constructs. However I did not manage any improvement. My main problem is that I have a list of a lot of data e.g.: > xs [[1]] [1]........................[1000] [[2]] [1]........................[840] ... [[50]]

Hi to all, In other statistical software, such as Eviews, it is possible to regress a model with the Least Squares method, assuming that the residuals follow an AR(q) process. For example the resulting regression is something like y = 1.2154 + 0.2215 x + 0.251 AR(1) How is it possible to do the same in R? Thank you very much in advance, Constantine Tsardounis http://www.costis.name

Dear Sir, First of all Happy Holidays!,... I am writing to you because I am a bit confused about ARCH estimation. Is there a way to find what garch() exactly does, without the need of reading the source code (because I cannot understand it)? In Eviews (the results at the end) I am getting different results than in R (for those that have the program I do: Quick -> Estimage Equation ->

-- jim knoble | jmknoble at pobox.com > On Apr 30, 2017, at 19:40, Darren Tucker <dtucker at zip.com.au> wrote: > >> On Mon, May 1, 2017 at 12:09 PM, Damien Miller <djm at mindrot.org> wrote: >> >> Hi, >> >> I just deleted SSHv1 support in OpenBSD and portable OpenSSH. There's >> probably a little dead code still to be expunged, but

Hi, I just deleted SSHv1 support in OpenBSD and portable OpenSSH. There's probably a little dead code still to be expunged, but all user-visible functionality and the bulk of the supporting infrastructure is gone. Sic transit gloria mundi. -d

Hello to everyone! I am working on my final year project about multivariate time series. There are three variables in the multivariate time series model. I have a few questions: 1. I used acf and pacf plot and find my variables are nonstationary. But in adf.test() and pp.test(), the data are stationary. why? 2.I use VAR to get a model. y is the matrix of data set and I have made a once

Hi, I''m running into some trouble when trying to access some variables which actually are a class parameter. The following works OK: node "common" { $var=''foo'' } node "gratia-base" inherits common { class { ''gratia'': dbhost => ''server.mysite.net'', } } node "node1"

Hello! I do not understand what is meant by: "aic" and "bic" follow a top-down strategy based on the Akaike's and Schwarz's information criteria in the datails to the ADF.test function. What does a "top-down strategy" mean? Probably the respective criterion is minimized and the mode vector contains the lag orders at which the criterion attains it

hi, in 1.9.1, the return value from ARMAacf(pacf=TRUE) is not named by lags, contrary to ?ARMAacf. the simple fix is to move names(Acf) <- down after if(pacf), with an appropriate starting lag as pacf=TRUE appears to start at lag 1 (whereas pacf=FALSE starts at lag 0). for consistency, one could argue to append 1 for lag 0 for pacf=TRUE (or start pacf=F at lag 1). however, given the

pacf in devel seems by default to return a different number of lags than 1.7.1 for $pacf. I don't see any mention of this in the NEWS file, or any change in the documentation, so I suspect it is and error, though it may be an undocumented improvement. (Newbie question: How is the simplest way to display a function like pacf.default that is not exported from a namespace?) Paul

How can I make sure the residual signal, after subtracting the trend extracted through some technique, is actually trend-free ? I would greatly appreciate any suggestion about some Stationarity tests. I'd like to make sure I have got the difference between ACF and PACF right. In the following I am citing some definitions. I would appreciate your thoughts. ACF(k) estimates the correlation

I am getting the following: *** Weave Errors *** Error in driver$runcode(drobj, chunk, chunkopts) : Error in eval(expr, envir, enclos) : couldn't find function "pacf.mts" *** Source Errors *** Error in eval(expr, envir, enclos) : couldn't find function "pacf.mts" make[1]: *** [checkVignettes] Error 1 I don't really understand the new namespace mechanism,

Dear list, according to the documentation of acf{ts} "the partial correlation coefficient is estimated by fitting autoregressive models of successively higher orders up to lag.max. " However, R seems to return the Yule-Walker estimates of the PACF by default. You can check this using c(1:10) as the series: the YW estimates are 0.7000000 and -0.1527035 for lags 1 and 2 . If the PACF

Hi! I want to plot three graphs (residuals, ACF and PACF of a model). Ideally I would use a c(2,2) disposition where the residuals plot would start at position 1,1 and span to position 1,2. Then I would plot the ACF in position 2,1 and the PACF in position 2,2. Maybe is clearer like this: -------------------------- | | | residuals | |

When I use the command for PACF, lags 5, 10, 15, and 20 are labeled. I would like to label lag 1. I would greatly appreciate if someone could tell me how to do this. Below is the command that I am using: pacf(data$R1,main="Series R1 Residuals") [[alternative HTML version deleted]]

R: Is there an easy way to get the acf and pacf for an irregular times series? That is, the acf and pacf with lag lengths that are in units of time, not observation number. Thanks, Jason Higbee Research Associate Federal Reserve Bank of St. Louis The views expressed in this email are the author's and not necessarily those of the Federal Reserve Bank of St. Louis or the Federal Reserve

Dear all, I found the following behaviour with pacf() in the multivariate case, set.seed(10) x <- rnorm(1000,sd=10000) y <- rnorm(1000,sd=1) pacf(ts(cbind(x,y)),plot=FALSE,lag.max=10) Partial autocorrelations of series 'cbind(x, y)', by lag , , x x y 0.047 ( 1) 0.000 ( -1) 0.011 ( 2) 0.000 ( -2) 0.005 ( 3) 0.000 ( -3) 0.013 ( 4)

Estimados amigos Estoy dibujando las funciones acf y pacf de una variable de una serie "zoo": > ls.str(pat="T0.5") T0.5 : 'zoo' series from 2017-11-08 23:00:00 to 2017-11-15 06:59:00   Data: num [1:9120, 1:3] 55 49.8 51 50.1 36.5 ...   Index:  POSIXct[1:9120], format: "2017-11-08 23:00:00" "2017-11-08 23:01:00" "2017-11-08

Dear all, following the thread http://tolstoy.newcastle.edu.au/R/e2/devel/07/09/4338.html regarding the bug in the partial autocorrelation function for multivariate time series. I have prepared a web page with patches and relevant information. http://www2.stat.unibo.it/giannerini/R/pacf.htm Please do not hesitate to contact me for further clarifications regards Simone --