arrangeParam.Gibbs.2D: Used in the loop of a Gibbs sampler, setting parameters not...

In forestgeo/ctfs: Tools for the Analysis of Forest Dynamics

Description

Used in the loop of a Gibbs sampler, setting parameters not yet tested (j and above) to previous value (i-1), and other parameters (<j) to new value (i).

But unlike arrangeParam.Gibbs, on which this is based, this is for a matrix, not a vector. Also unlike arrangeParam.Gibbs, this accepts only two versions of the parameter matrix, the previous and the next. This is crucial for memory: otherwise, allparam would be a 3D array of large size.

The scalar j refers to one row, and the scalar k a column. The function assumes that Gibbs updates occur by row, so every element on rows < j is set to the next parameter, plus all those in row j up to column k.

Usage

arrangeParam.Gibbs.2D(j, k, allparam.prev, allparam.next)

forestgeo/ctfs documentation built on May 3, 2019, 6:44 p.m.