The value and degrees of freedom are the engine's (the tulpa_fit
method); the "nobs" attribute is nobs.tobs_fit() on the same fit, so
stats::BIC() and anything else reading the attribute counts the
observations nobs() reports.
Usage
# S3 method for class 'tobs_fit'
logLik(object, ...)