R 等张/递增平滑函数回归

我想拟合一个严格递增的光滑函数。可按如下方式生成数据

set.seed(100)
x  <- sort( runif(100) )
y  <- sort( runif(100)^2 )
y[sample.int(100,10)]  <- runif(10)
plot(x, y)
loessMod  <- loess("y~x", data=data.frame(x,y), span=3)
lines(x, predict(loessMod, x), col="red")
min( diff(loessMod$fitted) )

set.seed（100）
x