next up previous [pdf]

Next: TWO-STAGE LINEAR LEAST SQUARES Up: PEF ESTIMATION WITH MISSING Previous: Internal boundaries to multidimensional

Finding the prediction-error filter

The first stage of the least-squares estimation is computing the prediction-error filter. The second stage will be using it to find the missing data. The input data space contains a mixture of known data values and missing unknown ones. For the first stage of finding the filter, we generally have many more fitting equations than we need so we can proceed by ignoring the fitting equations that involve missing data values. We ignore them everywhere that the missing inputs hit the filter.

The codes here do not address the difficulty that maybe too much data is missing so that all weights are zero. To add stabilization we could supplement the data volume with a ``training dataset'' or by a ``prior filter''. With things as they are, if there is not enough data to specify a prediction-error filter, you should encounter the error exit from cgstep() [*].

user/gee/pef.c
void find_pef(int nd       /* data size */, 
	      float* dd    /* data [nd] */, 
	      sf_filter aa /* estimated filter */, 
	      int niter    /* number of iterations */)
/*< find PEF >*/ 
{
    hconest_init( dd, aa);
    sf_solver(hconest_lop, sf_cgstep, aa->nh, nd, aa->flt, dd, 
	      niter, "x0", aa->flt, "end");
    sf_cgstep_close();
}


next up previous [pdf]

Next: TWO-STAGE LINEAR LEAST SQUARES Up: PEF ESTIMATION WITH MISSING Previous: Internal boundaries to multidimensional

2013-07-26