brianstock / MixSIAR

A framework for Bayesian mixing models in R:
http://brianstock.github.io/MixSIAR/
90 stars 75 forks source link

Another TDF question - dealing with variability & uncertainty #280

Closed amchiapella closed 3 years ago

amchiapella commented 3 years ago

The more I get to know my data, run test models, and explore reviews of TDFs, the more apparent the incredible importance of proper TDFs in mixing models becomes. The range of possible TDFs is quite large for my consumer group (fish), and no lab studies exist for my species, or even a similar species (let alone their natural diets). How do folks reconcile this? Just take the mean of all relevant literature TDFs I can find and have a large SD in my discrimination file? Try to use C:N ratios of my prey to somehow adjust lit TDFs accordingly? Even with pretty robust prior diet information (hundred of stomach samples over the course of 6 months) and good model diagnostics, my MixSiar results just really don't make sense, and it's clear to me this is due to the TDFs. Any insight welcome! Thanks

amchiapella commented 3 years ago

This is one of those fun instances where I figured out a good answer (for my case at least) within minutes of posting this issue... maybe putting it down on paper helped me think through my options better. Posting the solution here in case it's helpful for others. Out of my three fish species that each had three size classes, exactly one of those size classes of one of those species had a specialized diet (ie, only ate one thing, as indicated by stomach contents). In this case, small alewife exclusively ate zooplankton. So I calculated the TDF between zooplankton and small alewife, and was even able to account for tissue turnover by taking the mean d13c and d15n of zooplankton for the first three months of the study and subtracting those from the mean d13c and d15n of small alewife from the latter 3 months. I decided these TDFs were more accurate to extrapolate to my other fish species (which were a similar size & fed on roughly similar prey) than some random lab study (or compilation of lab TDFs). In case it's helpful: TDFs were 2.47 (SD 1.15) for d13C and 2.3 (SD 1.08) for d15N. System is a large temperate oligotrophic lake.

ChrisHarrod commented 3 years ago

Hi Butting in here. What fish is it and what tissues and foods are you using?

Cheers Chris

From: amchiapella @.> Sent: 05 May 2021 11:31 To: brianstock/MixSIAR @.> Cc: Subscribed @.***> Subject: [brianstock/MixSIAR] Another TDF question - dealing with variability & uncertainty (#280)

The more I get to know my data, run test models, and explore reviews of TDFs, the more apparent the incredible importance of proper TDFs in mixing models becomes. The range of possible TDFs is quite large for my consumer group (fish), and no lab studies exist for my species, or even a similar species (let alone their natural diets). How do folks reconcile this? Just take the mean of all relevant literature TDFs I can find and have a large SD in my discrimination file? Try to use C:N ratios of my prey to somehow adjust lit TDFs accordingly? Even with pretty robust prior diet information (hundred of stomach samples over the course of 6 months) and good model diagnostics, my MixSiar results just really don't make sense, and it's clear to me this is due to the TDFs. Any insight welcome! Thanks

— You are receiving this because you are subscribed to this thread. Reply to this email directly, view it on GitHubhttps://github.com/brianstock/MixSIAR/issues/280, or unsubscribehttps://github.com/notifications/unsubscribe-auth/ACNZSNE42WN7LNRHGQSEXETTMFQDXANCNFSM44FGICKA.

amchiapella commented 3 years ago

Hi Chris,

I ended up posting a comment to my original post then closing it because I resolved it on my own (had one specialized consumer that I realized I could calculate a TDF from). But my consumers are alewife, rainbow smelt, and slimy sculpin. Sources are zooplankton, mysis, and benthic invertebrates. Also trying to run some models with mysis (detritus, benthic inverts, zoops, and phytoplankton as sources) but will most likely have to use a lit TDF for those.

On Wed, May 5, 2021 at 4:42 PM ChrisHarrod @.***> wrote:

Hi Butting in here. What fish is it and what tissues and foods are you using?

Cheers Chris

From: amchiapella @.> Sent: 05 May 2021 11:31 To: brianstock/MixSIAR @.> Cc: Subscribed @.***> Subject: [brianstock/MixSIAR] Another TDF question - dealing with variability & uncertainty (#280)

The more I get to know my data, run test models, and explore reviews of TDFs, the more apparent the incredible importance of proper TDFs in mixing models becomes. The range of possible TDFs is quite large for my consumer group (fish), and no lab studies exist for my species, or even a similar species (let alone their natural diets). How do folks reconcile this? Just take the mean of all relevant literature TDFs I can find and have a large SD in my discrimination file? Try to use C:N ratios of my prey to somehow adjust lit TDFs accordingly? Even with pretty robust prior diet information (hundred of stomach samples over the course of 6 months) and good model diagnostics, my MixSiar results just really don't make sense, and it's clear to me this is due to the TDFs. Any insight welcome! Thanks

— You are receiving this because you are subscribed to this thread. Reply to this email directly, view it on GitHub< https://github.com/brianstock/MixSIAR/issues/280>, or unsubscribe< https://github.com/notifications/unsubscribe-auth/ACNZSNE42WN7LNRHGQSEXETTMFQDXANCNFSM44FGICKA>.

— You are receiving this because you modified the open/close state. Reply to this email directly, view it on GitHub https://github.com/brianstock/MixSIAR/issues/280#issuecomment-832993050, or unsubscribe https://github.com/notifications/unsubscribe-auth/AJZNYKHXWPN6KZRWLN7RW3LTMGURZANCNFSM44FGICKA .

--

Ariana Chiapella, PhD

Postdoctoral Associate

Rubenstein Ecosystem Science Laboratory

University of Vermont

http://arianachiapella.weebly.com/