Constructs coefficient matrices for low and high limits on the uniform prior distribution for beta.
gjamPriorTemplate(formula, xdata, ydata, lo = NULL, hi = NULL)A list containing two matrices. lo is a Q x S matrix of lower coefficient limits. hi is a Q x S matrix of upper coefficient limits. Unless specied in lo, all values in lo = -Inf. Likewise, unless specied in hi, all values in hiBeta = -Inf.
object of class formula, starting with ~, matches the formula passed to gjam
n x Q observation by predictor data.frame
n x Q observation by response data.frame
list of lower limits
list of upper limits
James S Clark, jimclark@duke.edu
The prior distribution for a coefficient beta[q,s] for predictor q and response s, is dunif(lo[q,s], hi[q,s]). gjamPriorTemplate generates these matrices. The default values are (-Inf, Inf), i.e., all values in lo equal to -Inf and hi equal to Inf. These templates can be modified by changing specific values in lo and/or hi.
Alternatively, desired lower limits can be passed as the list lo, assigned to names in xdata (same limit for all species in ydata), in ydata (same limit for all predictors in xdata), or both, separating names in xdata and ydata by "_". The same convention is used for upper limits in hi.
These matrices are supplied in as list betaPrior, which is included in modelList passed to gjam. See examples and browseVignettes('gjam').
Note that the informative prior slows computation.
Clark, J.S., D. Nemergut, B. Seyednasrollah, P. Turner, and S. Zhang. 2017. Generalized joint attribute modeling for biodiversity analysis: Median-zero, multivariate, multifarious data. Ecological Monographs 87, 34-56.
gjam