compute_diet_matrix: Computes the posterior diet matrix

In Irstea/escroc: EscrocR: a R package implementing the models ESCROC and EscroPath

Computes the posterior diet matrix

Description

Gives the quantiles of posterior diet

Usage

compute_diet_matrix(mydata, myres, quant = c(0.025, 0.5, 0.975))

Arguments

mydata	a list as returned by function prepare_data
myres	a fit as returned by fit_escroc
quant	a vector of quantiles

Value

a 3D array of dim length(quant) x nb species x nb species

Examples

#importing data
data(signature_data)
data(prior_diet_matrix)
data(LOQ)
data(prior_signature_data)
prior_delta <- data.frame(tracer=c("X15N","X13C"),mean=c(3,0),sd=c(1,1))

#check that everything is ok
mydata <- prepare_data(prior_diet_matrix,signature_data,
LOQ,prior_signature_data,prior_delta)

#build the model
mymodel <- building_model(mydata)

#fit the model
myresults <- fit_escroc(mydata, mymodel)

#diet
compute_diet_matrix(mydata, myresults)

Irstea/escroc documentation built on June 18, 2022, 11:22 p.m.