5

我对使用单调样条很感兴趣,但是当 R 尝试使用它时出现错误。我正在使用 R 2.12.0,并且方法“monoH.FC”表示自 2.8.0 起就支持它

可重现的示例(更复杂的 (x,y) 关系的结果相同)

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

我试过的

?spline返回:

...
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"’.
...

spline函数本身表明不支持“monoH.FC”方法:

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

问题

如何将 method = 'monoH.FC' 与样条曲线一起使用?

4

1 回答 1

9

使用splinefun;它支持method=monoH.FC.

中的最后一个示例?spline向您展示了如何做到这一点。

## 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)
于 2011-01-04T17:15:15.453 回答