Code to plot continuous effect + random effect of individual (PeerJ Fig 6)

[brianstock]

Hi Brian,

I've been working through the MixSIAR code you posted on GitHub that deals with the alligator data from Nifong et al 2015 and would like to apply a similar approach for my fish data. In particular, I would like to replicate Figure 6 from your 2018 PeerJ paper. The alligator script seems to work well for my data but, unless I missed something, I didn't see the necessary code for Figure 6. Would it be possible for you to provide the code for this? If not, do you have any suggestions on how to proceed?

Sorry if this is answered somewhere on GitHub. I am new to both it and MixSIAR.

[brianstock]

Sorry, this should have been available already, it's not answered on Github. I thought I'd included all the code for the alligator case study, but I checked and it only runs the 8 models in Table 1 (not the Length + Individual model), then makes Table 1, Figure 5 and Figure 7 (not 6 or 8).

I've uploaded code to run the Length + Individual model and reproduce Figures 6 and 8. If you re-install MixSIAR from GitHub instead of CRAN, you'll find a new script in the example_scripts folder, mixsiar_script_alligator_length_ind.R:

devtools::install_github("brianstock/MixSIAR", build_vignettes = TRUE)
library(MixSIAR)
find.package("MixSIAR") # location of MixSIAR install, example_scripts is in here

Alternatively, view and download from GitHub here.