Skip to contents
#devtools::load_all()
## Loading required package: ramp.xds
## Loading required package: ramp.func
## Loading required package: ramp.qa
sp <- makepar_F_sin(bottom = 1)
tp <- makepar_F_spline(seq(0, 3650, length.out=5), c(1,0.2,.25,1.5,1))
ts1 <- ts_setup(.3, sp, tp, tnorm = 3650)  
tm = 1:3650
ts <- F_pr(tm, ts1)
plot(tm, ts, type ="l", main = round(100*mean(ts))/100)