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

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)