У меня есть данные, которые строго увеличиваются, и я хотел бы соответствовать сглаживающему сплайну, который также монотонно увеличивается сsmooth.spline() Функция, если возможно, из-за простоты использования этой функции.

Например, мои данные могут быть эффективно воспроизведены на следующем примере:

testx <- 1:100
testy <- abs(rnorm(length(testx)))^3
testy <- cumsum(testy)
plot(testx,testy)
sspl <- smooth.spline(testx,testy)
lines(sspl,col="blue")

который не обязательно увеличивается везде. Какие-либо предложения?