add imputation to `rolling_mean`/`rolling_sd`

[dsweber2]

This is a minor variation on either of those techniques where any actually missing data is imputed based on the neighborhood average. Might be better filed under a separate function called impute that can also get used as a way of dealing with latency, though that is frequently extrapolation, rather than interpolation