Histogram smoothing for whitestripe

Source: R/smooth_hist.R

Uses a generalized additive model (GAM) to smooth a histogram for whitestripe

smooth_hist(x, y, deg = 4, k = floor(min(250, length(x)/2)),
  method = "REML", ...)

Arguments

x

values of midpoints from hist
y

values of counts from hist
deg

degree of polynomials used
k

Number of knots
method

Method for smoothing for GAM
...

Arguments passed to gam

Value

List of objects: x and y coordinates of histogram, coefficients from GAM, fitted values from GAM, the GAM model, the knots fitted, and degrees of polynomials

See also

gam

Examples

data(t2.voi.hist)
y = t2.voi.hist$counts
x = t2.voi.hist$mids
x = x[!is.na(y)];
y = y[!is.na(y)]
# 30 used for speed of example
s.hist = smooth_hist(x, y, k=30)
plot(t2.voi.hist, border="red")
lines(s.hist)