# NOT RUN {
ourModel <- ces(rnorm(100,0,1),h=10)
pointLik(ourModel)
# Bias correction
pointLik(ourModel) - nParam(ourModel)
# Bias correction in AIC style
2*(nParam(ourModel)/nobs(ourModel) - pointLik(ourModel))
# BIC calculation based on pointLik
log(nobs(ourModel))*nParam(ourModel) - 2*sum(pointLik(ourModel))
# }
Run the code above in your browser using DataLab