mwStatsGrid {astrochron} | R Documentation |

## 'Dynamic window' moving average, median and variance of stratigraphic series, using evenly spaced spatial/temporal grid

### Description

'Dynamic window' moving average, median and variance of stratigraphic series. 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. The results are output on an evenly spaced spatial/temporal grid (this contrasts with mwStats).

### Usage

```
mwStatsGrid(dat,cols=NULL,win=NULL,step=NULL,start=NULL,end=NULL,output=T,norm=F,
palette=6,ncolors=100,genplot=1,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 variable column to be extracted (first column automatically extracted). |

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

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

`start` |
Starting point for analysis, in units of space or time. |

`end` |
Ending point for analysis, in units of space or time. |

`norm` |
Normalize density estimates to maximum value? (T or F). If false, density estimates are normalized to unit area. |

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

`palette` |
What color palette would you like to use? (1) rainbow, (2) grayscale, (3) blue, (4) red, (5) blue-white-red, (6) viridis |

`ncolors` |
Number of colors to use in plot. |

`genplot` |
Generate summary plots? (0=none, 1=all time series, 2=kernel density estimates for each window, 3=kernel density estimates with median, 4=kernel density estimates with mean) |

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

### Value

A data frame with four columns: Center of window, Average, Median, Variance

### Examples

```
# generate example series from ar1 noise, 5 kyr sampling interval
ex = ar1(npts=1001,dt=5)
# jitter sampling times
ex[1]=ex[1]+rnorm(1001,sd=1)
# sort
ex = ex[order(ex[,1],na.last=NA,decreasing=FALSE),]
# run mwStats
mwStatsGrid(ex,win=100)
```

*astrochron*version 1.3 Index]