Scarlett
WAIC is not computed by the ds()
function. There is a function
AIC
() in the Distance package that can compute Bayes Information Criterion (BIC) by setting the argument
k
to be the natural log of the number of detections.
Example here if modelobj
is the object returned by ds()
and
nval
is the number of detections.
BIC <- AIC(modelobj, k = log(nval))$AIC
As I understand WAIC, it requires a posterior distribution of parameters, hence is Bayesian in nature. The Distance package does not carry out Bayesian estimation.