data(ple4)
hs = srrTMB(as.FLSR(ple4,model=segreg),spr0=spr0y(ple4),lplim=0.07,uplim=0.5)
plotsrs(hs)
jitter = srrjitter(hs)
#> Warning: NaNs produced
plotsrs(jitter$groups)
plotsrs(jitter$best) # Best
# Relax lower bound
plotsrs(FLSRs(init=hs,best=jitter$best))