https://github.com/cran/pracma
Raw File
Tip revision: d50cda4811eb3cdb8d2ce9e01ddb5e07d4a8c24f authored by HwB on 11 September 2013, 00:00:00 UTC
version 1.5.5
Tip revision: d50cda4
pchip.R
##
##  p c h i p . R  Test suite
##


pchip <- pracma::pchip

x <- c(1, 2, 3, 4, 5, 6)
y <- c(16, 18, 21, 17, 15, 12)

xs <- c(1.5, 2.5, 3.5, 4.5, 5.5)
ys <- pchip(x, y, xs)
# ys <- interp1(x, y, xs, method="cubic")       # the same
# 16.88750 19.80000 19.33333 15.96667 13.63750

yml <- c(16.887499999999999, 19.800000000000001, 19.333333333333332,
         15.966666666666667, 13.637499999999999)

all.equal(ys, yml, tolerance = 1e-7)
back to top