# NOT RUN {
library(MEDITS)
merge_TATB <- m.TATB(TA,TB,"ARISFOL")
merge_TATC <- m.TATC(TA,TC,"ARISFOL")
GSA <- unique(TA$AREA)
indices <- index.ts(merge_TATB,GSA,"ARISFOL",index = "abundance",
depth_range=c(500,800), sampling = "RSS",plot=FALSE)
freq <- LFD(merge_TATC,indices,sex="m",LC=1,depth_range=c(500,800))
quant(freq,0.95)
# }
Run the code above in your browser using DataLab