Logo Questions Linux Laravel Mysql Ubuntu Git Menu
 

Problem with spline method = 'monoH.FC''

Tags:

function

r

spline

I am interested in using the monotone spline, but I get an error when R tries to use it. I am using R 2.12.0, and the method 'monoH.FC' says that it has been supported since 2.8.0

Reproducible example (same result for more complicated (x,y) relationships)

x<-1:2
y<-1:2
spline(x,y,method="monoH.FC")
    Error in spline(x, y, method = "monoH.FC") : invalid interpolation method

What I have tried

?spline returns:

...
Usage:
...
        spline(x, y = NULL, n = 3*length(x), method = "fmm",
        xmin = min(x), xmax = max(x), xout, ties = mean)
...
Arguments:
  method: specifies the type of spline to be used.  Possible values are
      ‘"fmm"’, ‘"natural"’, ‘"periodic"’ and ‘"monoH.FC"’.
...

But the spline function itself indicates that the 'monoH.FC' method is not supported:

...
method <- pmatch(method, c("periodic", "natural", "fmm"))
if (is.na(method)) 
    stop("invalid interpolation method")
...

Question

How can I use method = 'monoH.FC' with spline?

like image 974
David LeBauer Avatar asked Jan 04 '11 17:01

David LeBauer


1 Answers

Use splinefun; it supports method=monoH.FC.

The last example in ?spline shows you how to do it.

## An example of  monotone  interpolation
n <- 20
set.seed(11)
x. <- sort(runif(n)) ; y. <- cumsum(abs(rnorm(n)))
plot(x.,y.)
curve(splinefun(x.,y.)(x),                add=TRUE, col=2, n=1001)
curve(splinefun(x.,y., method="mono")(x), add=TRUE, col=3, n=1001)
legend("topleft", paste("splinefun( \"", c("fmm", "monoH.CS"), "\" )", sep=''),
         col=2:3, lty=1)
like image 192
Joshua Ulrich Avatar answered Sep 30 '22 12:09

Joshua Ulrich