mwCor {astrochron} | R Documentation |

## Calculate moving window correlation coefficient for two stratigraphic series, using a 'dynamic window'

### Description

Calculate moving window correlation coefficient for two stratigraphic series, using a 'dynamic window'. This routine adjusts the number of data points in the window so it has a constant duration in time or space, for use with unevenly sampled data.

### Usage

```
mwCor(dat,cols=NULL,win=NULL,conv=1,cormethod=1,output=T,pl=1,genplot=T,verbose=T)
```

### Arguments

`dat` |
Your data frame containing stratigraphic data; any number of columns (variables) are permitted, but the first column should be a location identifier (e.g., depth, height, time). |

`cols` |
A vector that identifies the two variable columns to be extracted (first column automatically extracted). |

`win` |
Moving window size in units of space or time. |

`conv` |
Convention for window placement: (1) center each window on a stratigraphic level in 'dat' (DEFAULT), (2) start with the smallest location datum in 'dat', (3) start with the largest location datum in 'dat'. For options 2 and 3, the center of the window will not necessarily coincide with a measured stratigraphic level in 'dat', but edges of the data set are better preserved. |

`cormethod` |
Method used for calculation of correlation coefficient (1=Pearson, 2=Spearman, 3=Kendall) |

`output` |
Output results? (T or F) |

`pl` |
(1) Plot results at center of window, or (2) create "string of points plot" as in Sageman and Hollander (1999) |

`genplot` |
Generate summary plots? (T or F) |

`verbose` |
Verbose output? (T or F) |

### References

B.B. Sageman and D.H. Hollander, 1999,
*Cross correlation of paleoecological and geochemical proxies: A holistic approach to the study of past global change*, in E. Barrera and C.C. Johnson, eds., GSA Special Paper 332, p. 365-384.

### Examples

```
# generate example series
ex <- cycles(freqs=c(1/40,1/20),noisevar=.2)
# add second variable
ex[3] <- cycles(freqs=c(1/40,1/20),noisevar=0.2)[2]
# jitter sampling times
ex[1]=ex[1]+rnorm(500,sd=1)
# sort
ex = ex[order(ex[,1],na.last=NA,decreasing=FALSE),]
# run mwCor
mwCor(ex,win=50)
```

*astrochron*version 1.3 Index]