MixSIAR
MixSIAR copied to clipboard
using discrim data with raw source data
Hi Brian,
How do you run the model using raw source data and experimentally derived TDF values for each prey type? I have only been able to get the model to run by loading a file of source means/SD/n for each prey type and a file with the corresponding TDF values for each prey type. Based on the errors I've included below, I'm guessing this is because the row names in the mean source data file and the TDF file match. However, based on what I've read in the MixSIAR manual and your response to issue #102, if I want the model to account for the covariance between C and N isotopes, I would need to use raw source data in the model. When I try to run the model after loading in raw source data and a discrimination file containing TDFs for each prey type, I get the following error:
Error in jags.model(model.file, data = data, inits = init.values, n.chains = n.chains, :
RUNTIME ERROR:
Non-conforming parameters in function inprod
Additionally, when I try to use plot_data
using the same files I get this error:
Error in plot_data(filename = "isospace_plot_raw", plot_save_pdf = TRUE, :
*** Error: Source names do not match in source and discr
data files. Please check your source and discr data file row names.
So, is there a way for me to run the model using raw source data and experimentally derived TDFs for each prey type?
Thanks so much! Hannah