jackknife.SR.Rd
結果のプロットもこの関数で行う
jackknife.SR(
resSR,
is.plot = TRUE,
use.p0 = TRUE,
output = FALSE,
filename = "jackknife",
ylim.range = c(0, 1.5),
pch = 19,
cex = 1.1,
...
)
if (FALSE) { # \dontrun{
data(res_vpa)
SRdata <- get.SRdata(res_vpa)
resSR <- fit.SR(SRdata, SR = c("HS","BH","RI")[1],
method = c("L1","L2")[2], AR = 1,
out.AR = TRUE)
res_jackSR <- jackknife.SR(resSR,output = T)
# if calculate steepness
bio_par <- derive_biopar(res_obj=res_vpa,derive_year = 2010)
resSR <- fit.SR(SRdata, SR = c("HS","BH","RI")[1],
method = c("L1","L2")[2], AR = 1,
out.AR = TRUE,bio_par=bio_par)
res_jackSR <- jackknife.SR(resSR,output = T)
} # }