static void populate_x12a_arma_model (MODEL *pmod, const char *path, const DATASET *dset, arma_info *ainfo) { char fname[MAXLEN]; int err; pmod->t1 = ainfo->t1; pmod->t2 = ainfo->t2; pmod->ncoeff = ainfo->nc; pmod->full_n = dset->n; err = gretl_model_allocate_storage(pmod); if (err) { pmod->errcode = err; return; } sprintf(fname, "%s.est", path); err = get_estimates(fname, pmod, ainfo); if (!err) { sprintf(fname, "%s.rsd", path); err = get_uhat(fname, pmod, dset); } if (!err) { sprintf(fname, "%s.lks", path); err = get_ll_stats(fname, pmod); } if (!err) { sprintf(fname, "%s.rts", path); err = get_roots(fname, pmod, ainfo); } #if 0 if (!err) { sprintf(fname, "%s.acm", path); err = get_x12a_vcv(fname, pmod, nc); /* also .rcm */ } #endif if (err) { fprintf(stderr, "problem reading X-12-ARIMA model info\n"); pmod->errcode = err; } else { write_arma_model_stats(pmod, ainfo, dset); } }
static void handle_null_model (MODEL *pmod, arma_info *ainfo) { int full_n = pmod->full_n; pmod->ncoeff = 1; pmod->full_n = 0; pmod->errcode = gretl_model_allocate_storage(pmod); pmod->full_n = full_n; if (!pmod->errcode) { gretl_model_set_int(pmod, "null-model", 1); pmod->coeff[0] = 0.0; pmod->sigma = pmod->sdy; } }