BIC_bayesbr {bayesbr} | R Documentation |
A function that receives data from the estimated model, uses the information from the loglik, the number of observations of the model and the number of estimated parameters and returns the BIC, an estimator for the quality of the estimation of a model.
BIC_bayesbr(x)
x |
an object of the class bayesbr, containing the list returned from the |
Proposed by Stone (1979) the BIC (Bayesian Information Criterion) measures the quality of the adjustment made by the model, when comparing adjusted models with the same data, the smaller the BIC the better the adjustment.
The BIC theory requires that the log-likelihood has been maximized, but as we are in the context of Bayesian statistics, the log-likelihood as explained in the logLik.bayesbr
is made with the average of the a priori distribution for each theta and applying this value in the formula to calculate the loglik.
The BIC is calculated by
BIC = log(n)*k - 2 * L ,
where n
is the number of observations of the model variables, k
is the number of covariates used in the model, and L is the average of the loglik chain returned by the function logLik.bayesbr
.
A number corresponding to the BIC (Bayesian Information Criterion) of the estimated model.
Schwarz, G. (1978). Estimating the dimension of a model. The annals of statistics, 6(2), 461-464.
bayesbr
, AIC_bayesbr
, DIC_bayesbr
data("CarTask",package = "bayesbr") car_bayesbr <- bayesbr(probability ~ NFCCscale + task, data = CarTask, iter =100) bic = BIC_bayesbr(car_bayesbr)