close Warning: Can't synchronize with repository "(default)" (/var/svn/tolp does not appear to be a Subversion repository.). Look in the Trac log for more information.

Opened 16 years ago

Closed 16 years ago

Last modified 16 years ago

#688 closed enhancement (fixed)

New optional method marginalDraw in non linear filters

Reported by: Víctor de Buen Remiro Owned by: Víctor de Buen Remiro
Priority: highest Milestone: BSR hierarchy and priors
Component: Math Version: 1.1.7
Severity: blocker Keywords: BSR, non lineal filter
Cc:

Description (last modified by Víctor de Buen Remiro)

It's needed a new optional method in non linear filters with next default behaviour

  Matrix marginalDraw(Real unused)
  {
    SetCol(For(1,Card(_.colNames),Real(Real i) { ? })
  };

This method will be called just before to call draw method of corresponding internal non linear block in order to select which parameters will be drawn. It will return a column matrix with marginal sampling in selected rows and ? in other ones that will be sampled by non linear block in the standard way.

You can use it to apply marginal simulation to some selected parameters. For example, in order to apply a mixed-discrete distribution with 10% of probability of zero for all parameters:

  Matrix marginalDraw(Real unused)
  {
    SetCol(For(1,Card(_.colNames),Real(Real i) 
    { 
      If(Rand(0,1)<.10, 0, ?)
    }))
  };

Change History (6)

comment:1 Changed 16 years ago by Víctor de Buen Remiro

Resolution: fixed
Status: newclosed

(In [962]) Fixes #688

comment:2 Changed 16 years ago by Víctor de Buen Remiro

(In [963]) Fixes #688

comment:3 Changed 16 years ago by Víctor de Buen Remiro

Description: modified (diff)
Summary: New optional method selectDrawIndex in non linear filtersNew optional method marginalDraw in non linear filters

comment:4 Changed 16 years ago by Víctor de Buen Remiro

(In [964]) Fixes #688

comment:5 Changed 16 years ago by Víctor de Buen Remiro

(In [965]) Fixes #688

comment:6 Changed 16 years ago by Víctor de Buen Remiro

(In [968]) Fixes #688

Note: See TracTickets for help on using tickets.