broman (version 0.86)

runningratio2: Running ratio with adaptive window

Description

Calculates a running ratio; a ratio sum(top)/sum(bottom) in a sliding window, but rather than a fixed-width window, use nearest positions to give some target denominator

Usage

runningratio2(pos, numerator, denominator, at = NULL, window_denom = 100)

Value

A vector with the same length as the input at (or pos, if at is NULL), containing the running ratio.

Arguments

pos

Positions for the values.

numerator

Values for numerator in ratio.

denominator

Values for denominator in ratio.

at

Positions at which running ratio is calculated. If NULL, pos is used.

window_denom

Target denominator for window for calculating ratio

Author

Karl W Broman

See Also

runningmean(), runningratio()

Examples

x <- 1:1000
y <- runif(1000, 1, 5)
z <- runif(1000, 1, 5)
plot(x, runningratio2(x, y, z, window_denom=10), type="l", lwd=2)
lines(x, runningratio2(x, y, z, window_denom=50), lwd=2, col="blue")
lines(x, runningratio2(x, y, z, window_denom=100), lwd=2, col="red")