similar to: query

I ran across this problem when playing with ccf(). Its function call is >function (x, y, lag.max = NULL, type = c("correlation", "covariance"), plot = TRUE, na.action = na.fail, ...) Internally, ccf() calls plot(), which digs up plot.acf() whose default style is type='h' . I wanted to pass the argument type='l' to the plotting routine, but of

Hi all, I am getting different results from ccf and cor, Here is a simple example: set.seed(100) N <- 100 x1 <- sample(N) x2 <- x1 + rnorm(N,0,5) ccf(x1,x2)$acf[ccf(x1,x2)$lag == -1] cor(x1[-N], x2[-1]) Results: > ccf(x1,x2)$acf[ccf(x1,x2)$lag == -1] [1] -0.128027 > cor(x1[-N], x2[-1]) [1] -0.1301427 Thanks, Tal ----------------Contact

Dear All, I am studying some process measurement time series in R and trying to identify time delays using cross-correlation function ccf. The results have however been bit confusing. I found a couple of years old message about this issue but unfortunately wasn't able to find it again for a reference. For example, an obvious time shift is observed between the measurements y1 and y2 when the

Hi. It's my understanding that a cross-correlation function of vectors x and y at lag zero is equivalent to their correlation (or covariance, depending on how the ccf is defined). If this is true, could somebody please explain why I get an inconsistent result between cov() and ccf(type = "covariance"), but a consistent result between cor() and ccf(type = "correlation")? Or

Dear all, Can anybody tell me how to export a 3d figure made with the plot3d function? I'm careless about whether it's still interactive or not in another format, as long I can get it out of R. Thanks! Alejandro Gonz?lez Departamento de Biodiversidad y Conservaci?n Real Jard?n Bot?nico Consejo Superior de Investigaciones Cient?ficas Claudio Moyano, 1 28014 Madrid, Spain Tel +0034

Hello All, I'm having trouble with the ccf() function. I am trying to do cross-correlation between two time-series, but I keep getting an error message I don't know what to do with. This what I type and the error message I get: > ccf(ts(mod[,1]),ts(mod[,2]),na.action='na.exclude',type='cor') Error in "colnames<-"(`*tmp*`, value = c("ts(mod[,

Hi I'm computing the correlation between two time-series x_t and y_t-1 (time-series lagged using the lag(y,-1) function) using the cor() function and the returned value is different from the value of ccf() function at the same lag. Any ideas why this is so? Thanks in advance for any hints. Mihnea [[alternative HTML version deleted]]

Dear all, given two numeric vectors x and y, the ACF(x) at lag k is cor(x(t),x(t+k)) while the CCF(x,y) at lag k is cor(x(t),y(t-k)). See below for a simple example. > set.seed(1) > x <- rnorm(10) > y <- rnorm(10) > x [1] -0.6264538 0.1836433 -0.8356286 1.5952808 0.3295078 -0.8204684 0.4874291 0.7383247 0.5757814 -0.3053884 > y [1] 1.51178117 0.38984324

Dear R gurus, I would like to use the ccf function on two matrices that are each 196000 x 12. Ideally, I want to be able to go row by row for the two matrices using apply for the ccf function and get one 196000 X 1 array output. The apply function though wants only one array, no? Basically, is there a way to use apply when there are two arrays in order to do something like correlation on a row

On 11/28/2006 11:50 AM, A.I. McLeod wrote: > Hi Duncan, Hi Ian. > > ccf(x,y) does not explain whether c(k)=cov(x(t),x(t+k)) or d(k)=cov(x(t),x(t-k)) is calculated. The following example demonstrates > that the c(k) definition is used: > ccf(c(-1,1,rep(0,8)),c(1,rep(0,9))) > However S-Plus acf uses the d(k) definition in their acf function. I don't think our code looks

Hi, May be this helps: #Creating some dummy data.? set.seed(24) lst1<-lapply(1:8,function(x) ts(sample(1:25,20,replace=TRUE))) set.seed(49) lst2<-lapply(1:8,function(x) ts(sample(1:45,20,replace=TRUE))) ??Find_Max_CCF() #No vignettes or demos or help files found with alias or concept or #title matching ?Find_Max_CCF? using regular expression matching. Found a function with the same

Hello everybody, I am sorry if my questions are too simple or not easily understandable. I’m not a native English speaker and this is my first analysis using this function. I have a problem with a cross correlation function and I would like to understand how I have to perform it in R. I have yearly data of an independent variable (x) from 1982 to 2010, and I also have yearly data of a variable

Hi all, I have two time series data (say x and y). I am interested to calculate the correlation between them and its confidence interval (or to test no correlation). Function cor.test(x,y) does the test of no correlation. But this test probably is wrong because of autocorrelated data. ccf() calculates the correlation between two series data. But it does not

Hello, I am a very new R user and not a statistician so please excuse any over explanation, I'm just trying to be as clear as possible. I have performed a cross correlation of two time series (my columns) in a single data setusing: ccf(ts(A[rows,columnX]),(A[rows,columnY]), lag=NULL, type="correlation",plot=F) I?am able to get the results (for example): Autocorrelations of

I can't understand the results from cross-correlation function ccf() even though it should be simple. Here's my short example: ********* a<-rnorm(5);b<-rnorm(5) a;b [1] 1.4429135 0.8470067 1.2263730 -1.8159190 -0.6997260 [1] -0.4227674 0.8602645 -0.6810602 -1.4858726 -0.7008563 cc<-ccf(a,b,lag.max=4,type="correlation") cc Autocorrelations of series 'X',

I am new to R and new to time series modeling. I have a set of variables (var1, var2, var3, var4, var5) for which I have historical yearly data. I am trying to use this data to produce a prediction of var1, 3 years into the future. I have a few basic questions: 1) I am able to read in my data, and convert it to a time series format using 'ts.' data_ts <- ts(data, start = 1988, end =

Dear R-users, I have been using R and its core-packages with great satisfaction now for many years, and have recently started using the "ccf" function (part of the "stats" package version 2.16.0), about which I have a question. The "ccf"-algorithm for calculating the cross-correlation between two time series always calculates the mean and standard deviation per time

Hi, I have been running the ccf() function to find cross-correlations of time series across various lags. When I give the option of plot=TRUE, I get a plot that gives me 95% confidence interval cut-offs (based on sample covariances) for my cross-correlations at each lag. This gives me a sense of whether my cross-correlations are statistically significant or not. However, I am unable to get R to

I am trying to run a cross-correlation using the "ccf()" function. When I select plot = TRUE in the ccf() I get a graph which has ACF on the y-axis, which would suggest that these y-values are the auto-correlation values. How should I adjust the code to produce a plot that provides the cross-correlation values? Here is my code: w002dat <-

Dear all, I am trying to avoid a for loop here and wonder if the following is possible: I have a data.frame with 6 columns and i want to get a cross-correlogram (by using ccf) . Obivously ccf only accepts two columns at once and then returms a list. In fact, with a for loop i?d do the following for (i in 1:6) { x[[i]]=ccf(mydf[,i],mydf[,6]) } Is there any chance to the same with