Changes to FSA — Estimating Abundance
Want to share your content on R-bloggers? click here if you have a blog, or here if you don't.
I mentioned previously, that I have been updating the Mark-Recapture vignettes. That has morphed into a document that is an update of the Mark-Recapture Closed and Open vignettes and the depletion/removal vignettes and associated FSA functions. Some of the changes are cosmetic, others are more important. I attempt to highlight the changes below but you should also see the NEWS file for the FSA package. The new vignette is here.
FSA changes
capHistConvert()
— major changes, should consider this to be a new function.capHistSum()
— added anm.array
object to the output, changed column and row labels to the methodB objects, and added aplot()
object to make the diagnostic plots identified in the Rcapture package. Old code should still work.depletion()
— Changedtype=
tomethod=
and addedDeLury
as an option tomethod=
(and leftDelury
). Changedricker.mod=
toRicker.mod=
. Old code that explicitly usedtype=
will be broken (need to change tomethod=
).mrClosed()
— Changedincl.inputs=
toverbose=,
changedtype=
tomethod=
. Added the ability to construct a CI for the overall PE when multiple groups are used in a Petersen family method (thus, added aincl.all=
toconfint()
). Changed default forincl.all=
fromFALSE
toTRUE
. Old code that explicitly usedtype=
will be broken (need to change tomethod=
).mrOpen()
— Changedci.type=
totype=
andphi.type=
tophi.full=
. Removedtype=
fromsummary()
and added averbose=
which will print only the estimates ifFALSE
or both observables and estimates ifTRUE
. Added averbose=
toconfint()
to control whether the message about the type of confidence interval is printed or not. Added another function,jolly()
, that is the same asmrOpen()
. Old code should work unlessci.type=
andphi.type=
were used inmrOpen()
ortype=
was used insummary()
.removal()
— Added options to perform Moran (1951) and Schnute (1983) removal methods. Completely modified the code so that the examples withapply()
andlapply()
would also provide confidence intervals. Changed so that if catch is a one row or one column matrix or data.frame then it will be converted to a vector, which allows one to extract one row from a data.frame to send toremoval()
without having to useas.numeric()
. Old code should still work.
New Abundance Vignette
Again, this is now a combination of the original M-R closed, M-R open, and depletion/removal events. All sections are like the original except for the changes in FSA from above and the theoretical development is greatly reduced.
The M-R closed section now includes an illustration to estimate N using the Petersen methods if the data is grouped (i.e., separate estimates by length groups) and for fitting the Otis et al. (1978) log-linear models to multiple M-R data with possible capture probability heterogeneities (using the Rcapture package).
The M-R open section now includes a section on fitting log-linear models (using Rcapture package).
The depletion/removal now includes is a section on using the method of Schnute (1983) to identify a statistically different capture probability for the first sample.
Filed under: Administration, Fisheries Science, R Tagged: Abundance, Depletion, FSA, Jolly-Seber, Mark-Recapture, R, Removal
R-bloggers.com offers daily e-mail updates about R news and tutorials about learning R and many other topics. Click here if you're looking to post or find an R/data-science job.
Want to share your content on R-bloggers? click here if you have a blog, or here if you don't.