Version 0.14
-------------------------------------------------------------------------

FEATURES

* new release on CRAN.

* performs Bayesian non-parametric modeling on a (rectangular) set of functional data observations.

* The collection of functions may be modeled under Gaussian process (GP) or intrinsic Gaussian Markov

* random field (iGMRF) prior formulations.

* The covariance and precision parameters of the GP and iGMRF formulations, respectively, are placed under

* a Dirichlet process (DP) prior to allow the data to discover dependence among the estimated functions

* where co-clusters functions are drawn from distributions sharing the same covariance and precision parameters.

* the GP prior formulation is invoked with gpdpgrow()

* any number of additive covariance terms may be specified with gpdpgrow().

* for example, if there are 4 terms, then the input variable, gp_cov = c("rq","se","sn","sn")

* if the covariance functions for the 4 terms are structured as (rational quadratic, squared exponential,

* seasonal, seasonal), respectively.  The input variable, sn_order = c(3,12), sets the order for each seasonality

* term; in this case, 3 months and 12 months (assuming the data time scale is denoted by month).

* the iGMRF prior is invoked with gmrfdpgrow(), also allowing any number of additive precision terms

* the input variable, q_type = c("tr","sn","sn"), denotes "tr' = trend, and "sn" = seasonality terms.

* input, q_order = c(2,3,12) denotes the order for the associated term; for example, the second term

* is specified as seasonal of order = 3 (e.g. months).


CHANGES

08/10/2014
-----------
* version 0.1 launched on CRAN.

12/09/2014
----------
* replaced srand() with arma_rng::set_seed_random() to initialize random seed for RcppArmadillo.
* removed use of arma::sp_mat (sparse matrix) formulation for normalized CAR adjacency matrix from
  gmrfdpgrow() due to memory initialization issue in RcppArmadillo.

02/20/2015
----------
* removed use of arma_rng::set_seed_random() to avoid warning about setting R seed from C++. 

10/16/2015
----------
* Minor changes in anticipation of ggplot2 1.1.0

07/27/2016
----------
* return pop_plot and samp_plot ggplot2 objects in gen_informative_sample() that plot generated synthetic functions in the population 
and sample, respectively.
* employ uniform(0,1) starting values in the MCMC sampler when co-sampling the GP functions, bb, in gpdpgrow() to prevent numerical instability
producing NaNs in auxclusterstep() for sampling cluster assignments under Debian Linux.

04/05/2017
----------
* added new option to gmrfdpgrow() to allow input of N x R predictor matrix, ksi, that is used to update the prior probability of cluster assignments,
s.   This dependent product partition model is enabled by placing a distribution on the ksi, though, we don't believe they're random to esablish a coherence
function multiplicative input to the prior distribution on cluster assignment.
* Note that we had to switch the sampling algorithm for cluster assignment from a conjugate multinomial to the non=conjugate Neal's algorithm 8 since the 
inclusion of ksi produces a non-conjugate mixture.
* added new option to gmrfdpgrow() to support count data response if user inputs non-NULL N x T offset matrix, E.

06/23/2021
----------
* Update imports to remove package mvtnorm and use spam::rmvnorm.