display "BIC = " -2 * e(ll) + ln(e(N)) * (e(df_m)+1)