#加载nitrogen samples
mix.filename <- "E:/R2025222/samples.csv"
mix <- load_mix_data(filename=mix.filename,
iso_names=c("d18O","d2H"),factors="Group",fac_random=FALSE,fac_nested=FALSE,
Error in load_mix_data(filename = mix.filename, iso_names = c("d18O", :cont_effects=NULL)
could not find function "load_mix_data"加载source data
source.filename <- "E:/R2025222/sources.csv"
source <- load_source_data(filename=source.filename,source_factors=NULL,conc_dep=FALSE,data_type="means",
Error in load_source_data(filename = source.filename, source_factors = NULL, :mix)
could not find function "load_source_data"加载discrimination/TDF data
discr.filename <- "E:/R2025222/discrimination.csv"
discr <- load_discr_data(filename=discr.filename, mix)
Error in load_discr_data(filename = discr.filename, mix) :
could not find function "load_discr_data"#绘制散点图
plot_data(filename="isospace_plot",plot_save_pdf=TRUE,plot_save_png=TRUE,
Error in plot_data(filename = "isospace_plot", plot_save_pdf = TRUE, plot_save_png = TRUE, :mix,source,discr)
could not find function "plot_data"Calculate standardized convex hull area
if(mix$n.iso==2) calc_area(source=source,mix=mix,discr=discr)
Error: object 'mix' not found#绘制所有组的先验预测贝叶斯图和pdf
plot_prior(alpha.prior=1,source)
Error in plot_prior(alpha.prior = 1, source) :
could not find function "plot_prior"定义模型结构,编写JAGS模型文件
model_filename <- "MixSIAR_model.txt"
resid_err <- FALSE
process_err <- TRUE
write_JAGS_model(model_filename, resid_err, process_err, mix, source)
Error in write_JAGS_model(model_filename, resid_err, process_err, mix, :
could not find function "write_JAGS_model"运行JAGS模型
jags.1 <- run_model(run="short", mix, source, discr, model_filename, alpha.prior=1)
Error in run_model(run = "short", mix, source, discr, model_filename, :
could not find function "run_model"
jags.1 <- run_model(run="normal", mix, source, discr, model_filename, alpha.prior=1)
Error in run_model(run = "normal", mix, source, discr, model_filename, :
could not find function "run_model"过程诊断,总结统计和后验预测贝叶斯图
output_JAGS(jags.1, mix, source)
Error in output_JAGS(jags.1, mix, source) :
could not find function "output_JAGS"
library(MixSIAR)
Error in library(MixSIAR) : 不存在叫‘MixSIAR’这个名字的程辑包#加载nitrogen samples
mix.filename <- "E:/R2025222/samples.csv"
mix <- load_mix_data(filename=mix.filename,iso_names=c("d18O","d2H"),factors="Group",fac_random=FALSE,fac_nested=FALSE,
Error in load_mix_data(filename = mix.filename, iso_names = c("d18O", :cont_effects=NULL)
could not find function "load_mix_data"打印加载的数据信息,可根据需要选择是否保留
print(mix)
Error in print(mix) : object 'mix' not found检查 MixSIAR 包是否已安装,如果未安装则进行安装
if (!require("MixSIAR", quietly = TRUE)) {
- install.packages("MixSIAR"