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

Error in MixSIAR (using single isotope): inconsistent of raw data and data presented in plot #352

Open
Waiting-William opened this issue Aug 7, 2023 · 0 comments

Comments

@Waiting-William
Copy link

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)

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

1 participant