Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

species-specific group size estimation #9

Closed
amandalbradford opened this issue Jul 26, 2023 · 1 comment
Closed

species-specific group size estimation #9

amandalbradford opened this issue Jul 26, 2023 · 1 comment

Comments

@amandalbradford
Copy link

Aloha @smwoodman - After finding some species-specific group sizes (GsSpBestx) being higher than the total group sizes (GsSchoolBest), @mchill808 and I realized that swfscDAS is averaging observer best estimates for GsSchoolBest as expected. However, for computing GsSpBestx for mixed-species sightings, instead of 1) averaging the observer best estimates, then 2) averaging the available percentages, and finally 3) multiplying those average values for each species (as we have always done for the group sizes that feed into abundance estimation), we found that swfscDAS is 1) multiplying each observer's GsSchoolBest by their GsSpBestx and then 2) averaging those values together for the final GsSpBestx, even though the included GsSchoolBest and SpPercx represent the averaged values that we think should be multiplied instead.

If all observers provide percentages along with their estimates, then the two approaches above are more consistent (but still differ). However, the discrepancy becomes greater with fewer estimates associated with percentages, with some species-specific group sizes even ending up higher than the total group size estimate. I know you provide the complete format to allow for differences in group size estimation, but we wondered if you would consider revising the current species-specific group size estimation approach in swfscDAS? Thank you!

@smwoodman
Copy link
Owner

In short, I don't remember a specific reason for the current approach, and I am happy to update das_sight as you propose. I agree with y'all that sum(GsSpBestx) should equal GsSchoolBest, and same for the Highs and Lows, and imagine that I just didn't notice this potential discrepancy because the files that I was using for testing had sightings where all observers provided percentages.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants