Re: [R] ccf and covariance

From: Prof Brian Ripley <ripley_at_stats.ox.ac.uk>
Date: Wed, 23 Apr 2008 19:31:24 +0100 (BST)

On Wed, 23 Apr 2008, Bob Farmer wrote:

> 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).

The ratio of your values is

> MASS::fractions(282568.5/259021)

[1] 12/11

? Do you recognize it?

There is an explanation in MASS4, p. 390, for example.

> 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 have I misunderstood what is a cross-correlation?
> (unfortunately, I can't seem to get a look at the ccf code, since I
> think it's buried in some C function outside of the main environment)

It is in the R sources, not 'buried' at all - that is what 'Open Source' means. You can browse them at https://svn.r-project.org/R/trunk, or download them for study.

> Thanks very much.
> --Bob Farmer
> PhD candidate, Dalhousie University
> Halifax, NS, Canada
>
> Example:
> d1<-data.frame(matrix(ldeaths, nrow = 6, byrow = T))
> seventy_4<-as.numeric(d1[1,])
> seventy_5<-as.numeric(d1[2,])
>
> ccf(x=seventy_4, y=seventy_5,
> plot = F, lag.max = 0, type = "covariance"
> )
> cov(seventy_4, seventy_5) #inconsistent
>
> ccf(x=seventy_4, y=seventy_5,
> plot = F, lag.max = 0, type = "correlation"
> )
> cor(seventy_4, seventy_5) #consistent
>
> ______________________________________________
> R-help_at_r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.
>

-- 
Brian D. Ripley,                  ripley_at_stats.ox.ac.uk
Professor of Applied Statistics,  http://www.stats.ox.ac.uk/~ripley/
University of Oxford,             Tel:  +44 1865 272861 (self)
1 South Parks Road,                     +44 1865 272866 (PA)
Oxford OX1 3TG, UK                Fax:  +44 1865 272595

______________________________________________
R-help_at_r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
Received on Wed 23 Apr 2008 - 19:32:00 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Wed 23 Apr 2008 - 21:30:40 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.

list of date sections of archive