sign_changes: Sign Changes

In dewittpe/cpr: Control Polygon Reduction

Sign Changes

Description

Count the number of times the first, or second, derivative of a spline changes sign.

Usage

sign_changes(
  object,
  lower = min(object$bknots),
  upper = max(object$bknots),
  n = 1000,
  derivative = 1L,
  ...
)

Arguments

object	a cpr_cp object
lower	the lower limit of the integral
upper	the upper limit of the integral
n	number of values to assess the derivative between lower and upper.
derivative	integer value denoted first or second derivative
...	pass through

Value

the number of times the sign of the first or second derivative changes within the specified interval.

See Also

wiggle

Examples

xvec <- seq(0, 6, length = 500)

# Define the basis matrix
bmat1 <- bsplines(x = xvec, iknots = c(1, 1.5, 2.3, 4, 4.5))
bmat2 <- bsplines(x = xvec)

# Define the control vertices ordinates
theta1 <- c(1, 0, 3.5, 4.2, 3.7, -0.5, -0.7, 2, 1.5)
theta2 <- c(1, 3.4, -2, 1.7)

# build the two control polygons
cp1 <- cp(bmat1, theta1)
cp2 <- cp(bmat2, theta2)
plot(cp1, cp2, show_cp = FALSE, show_spline = TRUE)

sign_changes(cp1)
sign_changes(cp2)

dewittpe/cpr documentation built on Feb. 16, 2024, 1:11 p.m.