MixSIAR
MixSIAR copied to clipboard
Error in MixSIAR (using single isotope): inconsistent of raw data and data presented in plot
Hello, I am trying to calculate endmember contributions using a single isotope case in MixSIAR. Then, I follow the code from this example and apply my data to do the calculations. Howerve, I found that the data in the plot is not consistent with my original data, and the results of each run are different from run to run. I tried many ways but failed. Do you have any suggestion to solve this problem?
library(MixSIAR)
Load mix data
mix.filename <- system.file("extdata", "soil_consumer.csv", package = "MixSIAR") mix <- load_mix_data(filename=mix.filename, iso_names=c("d13C"), factors=NULL, fac_random=NULL, fac_nested=NULL, cont_effects=NULL)
Load source data
source.filename <- system.file("extdata", "soil_sources.csv", package = "MixSIAR") source <- load_source_data(filename=source.filename, source_factors=NULL, conc_dep=FALSE, data_type="raw", mix)
Load discrimination/TDF data
discr.filename <- system.file("extdata", "soil_discrimination.csv", package = "MixSIAR") discr <- load_discr_data(filename=discr.filename, mix)
Make isospace plot
plot_data(filename="isospace_plot", plot_save_pdf=TRUE, plot_save_png=FALSE, mix,source,discr)