2.4-22 (27 Nov 2023) ------ fixed double from int in printf for Kurt, and updated citation format 2.4-21 (06 Jan 2023) ------ sprintf -> snprintf 2.4-20 (19 Oct 2022) ------ providing missing headers 2.4-19 (17 Aug 2022) ------ removed & and other latex Rd fixes 2.4-18 (02 Apr 2022) ------ updated JSS DOI and now USE_FC_LEN for BLAS 2.4-17 (11 Sep 2020) ------ Fixed clang-11 issue with linalg.h and CXX headers pointed out by Brian 2.4-16 (05 Sep 2020) ------ put all btgp and tgp.design examples in donttest 2.4-15 (02 Sep 2020) ------ updated for modern R requirements implemented drop=FALSE in tgp.partition for Jan simplified BLAS and Lapack headers in linalg.c/.h for monomvn now using R_PosInf and R_NegInf 2.4-14 (06 Feb 2016) ------ forgot to do compact vignettes 2.4-13 (06 Feb 2016) ------ removed Rinterface.h; not sure where it came from 2.4-12 (05 Feb 2016) ------ new web page on bobby.gramacy.com, and title case moved R.h and Rmath.h outside of rhelp.h to be compatible with new C++ headers inside of newest R removed beta functionality comment from old vignette 2.4-11 (28 November 2014) ------ removed akima from package 2.4-10 (27 November 2014) ------ added mean0.range1 to NAMESPACE for use in the plgp package cleaned up some warnings in new pedantic gcc compile moved akima and maptree to Imports 2.4-9 (01 April 2013) ----- Changed printed tables to 15 digits to comply with new read.table checks removed zzz.R and moved LICENSE to License global exp2d and data call made fancy to pass checks 2.4-8 (13 Oct 2012) ----- fixed bimprov/improv bug overlooked in NULL issue updates removeed bessel_k.c and assocated functions; now using built-in bessel_k_ex function which didn't exist before 2.4-7 (18 April 2012) ----- fixed NULL issue in lh_sample 2.4-6 (14 April 2012) ----- Added check that E divides T-B Corrected bug in slice grid causing overflow Corrected NULL problems pointed out by BR 2.4-5 (30 January 2012) ----- fixes to vigneetes after NAMESPACE change, and removing direct references to stdout and stderr 2.4-4 (23 January 2012) ----- Added NAMESPACE file 2.4-3 (18 December 2011) ----- Checked for assert problems with NDEBUG and updated R version (2.14) and email addresses 2.4-2 (30 March 2011) ----- added dp.sim to tgp.default.params so that MVN proposal covariance matrix can be specified by the user in GP-SIM models 2.4-1 (23 Feb 2011) ----- Changed my_r_process_events to always ProcessEvents and check for interrupts (every second) Made new sobol_indices function to clean up tgp::Sens, and using shape to drive categorical calculations rather than bmax minor edit to tgp2 doc to correct fried.bool description, and edit to fried.bool to fix fried.bool(1) problem 2.4 (20 Sept 2010) --- added sim (single index model) rank 1 correlation structure overlooked fixing nugget in treed split proposals; now fixed 2.3-4 (19 July 2010) ----- allowed the entire $tree component from tgp.default.params to be written over by the tree argument to b* functions fixed a bug in the linear initialization of correlation functions (for non-LLM models) 2.3-3 (23 April 2010) ----- fixed another bug pointed out by Ripley 2.3-2 (21 April 2010) ----- replace isinf by R_FINITE as suggested by Brian Ripley 2.3-1 (8 Feb 2010) ----- Added reference to second JSS paper 2.3 (23 Dec 2009) --- final version of the tgp2 vignette, accepted at JSS removed mrtgp documentation 2.2-4 (20 Nov 2009) ----- added a check for a full rank design matrix in the pre-processing of inputs corrected some syntax errors in .Rd files caught by new checks on CRAN changed sens warning to check for nn > 0 instead of null XX which was causing a warning to be printed every time regardless of whether or not there was a violation 2.2-3 (22 Jul 2009) ----- added nug.p=0 option to cause the nugget to stay fixed at the params$gd[1] value. bcart, btlm, and blm now automatically use nug.p=0 since the nugget is not identified in these models 2.2-2 (27 May 2009) ----- caught a signifigant typo/bug in Xsplit pre-processing in R which was causing the bounding rectangle to be computed improperly removed combinat dependency 2.2-1 (20 Feb 2009) ----- added Tree::Clear() to rejected prunes in Tree::Prune() tgp.default.params and tgp.check.params now allow the min partition size to depend on basemax rather than d the bounding rectangle and the valid tree-split locations are now governed by an Xsplit automatically generated in the tgp R function, and gathered from the tgp object in predict.tgp 2.2 (16 Jan 2009) --- added lower bound for s2_g0, and added DEFAULT CPP macros to help prevent assertions which would have backed out of bad parameterizations anyways. This more liberal approach essentially prevents assertions in the no-noise regression case now allow basemax != d for constant mean functions added categorical LHS for sensitivity analyses finished draft of second vignette, to be submitted to JSS 2.1-6 (25 Nov 2008) ----- caught important bug in b0-prior pre-processing in R that was causing the emperical Bayes prior to be used when the b0 hierarchical one is requested 2.1-5 (24 Oct 2008) ----- added check for >= d+1 unique rows of X in Tree::Singular changed how d an nug are printed in the progress meter for mrtgp added par(mfrow=c(1,1)) to default 1-d mrtgp plots 2.1-4 (17 Oct 2008) ----- made it possible to cap the number of input locations which are ranked by the expected improvement statistic by allowing the improv argument to b* functions to take a 2-vector input. Also no longer ranks points with an improvement statistic of zero 2.1-3 (7 Oct 2008) ----- added check for >1 unique Z in Tree::Singular allow (undocumented) specification of minpart through the ellipses (...) argument to the b* functions m0r1=TRUE is now the default for all b* functions, the Rnw documentation has been updated (but the pdfs have not been regenerated in the current version) 2.1-2 (17 April 2008) ----- introduced params$Wi and made it so that params${b,Wi} set both the starting an prior values for the corresponding parameters. added bprior="bmznot" so that tau2 can be fixed at its starting value 2.1-1 (4 April 2008) ----- separaing tcounts and cum_tcounts and resetting tcounts to the average of cum_tcounts at beginning of each round fixed bug in returning of tcounts back to R by encuring that its->UpdateProbs() is always called after the last of R rounds whenever its->Numit() > 1 return ess information at each inverse temperature back to R for storage in the tgp-class object 2.1 (22 March 2008) --- moved CHANGES to ChangeLog for the new CRAN packages page cleaned up tgp vignette into one .Rnw file including several .iRnw files added linear=FALSE to akima interp commands stopped printing both separate ess and essN in IT method added new itemps.barplot function for visualizing the number of visits to each inverse temperature in the ST-MCMC chain now tcounts (observation counts) for IT are accumulated across repeats (R>=2) -- not sure if this is a good idea now recording tree acceptance rates (Grow, Prune Change and Swap accepts over total) propogated splitmin to up to the R interface, and added a basemax paramter also specifiable in the R inferface -- the parameters allow control of which X colums can be partitioned upon (by the tree) and which appear in the base (GP) model added heights="map" argument to tgp.trees for plotting only the tree corresponding to the maximum a' posteriori tree wrote two sections in the new vignette (tgp2): one on splitmin and basemax (i.e., categorical inputs) and one on importance tempering 2.0-4 (23 Jan 2008) ----- made a change suggested by Andrea Spano to make the partition function more robust when new X values (on a different range from the original data) are provided 2.0-3 (27 Dec 2007) ----- commented the moving average code (for sens) and made it more efficient 2.0-2 (12 Dec 2007) ----- allow mode in the sens function to be be within the bounds of rect -- i.e., stop checking for positive mode fixed dyslexic error in man page for sens describing total effect indices fixed missing braces in some .Rd files shortened to LGPL in DESCRIPTION 2.0-1 (1 Oct 2007) ----- Memory savings: preds->ZZ and preds->Zp only allocated when pred.n=TRUE, preds->w and preds->itemp only allocated when ST or IT or IS corrected default.itemps() example barplot and changed k.min to 0.1 to agree with IT paper fixed zcov bug in tgp.postprocess 2.0 (21 Sept 2007) --- consolidated R functions from individual files into those grouped by task made gridlen argument a 2-vector so the grid size can be specified in both x and y new economy argument to tree and Base duplicators so that the entire set of covariance matrices is not copied unless necessary Matt Taddy officially joins the list of authors, adding multi-resolution tgp and sensitivity analysis importance tempering, with importance sampling and simulated tempering support Latin hypercube sampling with respect to a Beta distribution moved most of the reading of parameters out of the Tgp constructor and into Tgp::Init so that tgp_cleanup() can safely delete the Tgp module even it it wasn't successfully initialized ordered multiple improv statistics predictive covariances now calculated instead of variances 1.2-7 (19 Aug 2007) ----- tgp.partition 1-d data handling made more robust removed assertion on w that tends to fail for unknown reasons 1.2-6 (17 June 2007) ----- added number of leaves to traces fixed CART bug in Gp::predict that was already handled in v1.3, thanks to Taddy, but was bad in the CRAN version 1.2-5 (29 May 2007) ----- allowing btgp to do a linburn $parts is a now a matrix, and $linburn now {TRUE,FALSE} $tree and $hier removed from tgp-class output, and $trace correctly removed when input trace=FALSE got rid of some hanging code fragments in documentation of btgp and dopt.gp 1.2-4 (07 May 2007) ----- fixed a leak when krige=FALSE appearing in Tree::Predict fixed sub-vignette cutting for preds and traces corrected documentation of nug.p in tgp.default.params() added verbosity and iter arguments to dopt and thus dopt.gp and tgp.design 1.2-3 (09 Apr 2007) ----- updates to the vignette for the JSS publication check.matrix modified to allow one-column matrix Z arguments which have column-name information tgp.trees calculation of rows and columns corrected for situations when there is a height 1 tree 1.2-2 (19 Dec 2006) ----- slice.interp bugs fixed (is.null(NULL) and plotting of axis labels) and tick marks added added R argument to predict.tgp and changed the krige argument to be called MAP added new krige argument to b* and predict.tgp to allow the user to control whether Z*k* kriging samples are to be obtained (possibly saving on memory) bug involving the order of trees read from files in the tgp.get.trees function is fixed; as a result the $trees output field in "tgp"-class objects is listed by height and so may contain NULL entries documented plot.tgp and predict.tgp with the method markup 1.2-1 (03 Dec 2006) ----- Added "cart-" option to bprior in order to allow a CART style mean rather than a LM mean 1.2 (21 Nov 2006) --- Now allowing BTE[2]-BTE[1] % BTE[3] != 0 Added predict.tgp and necessary functionality to save and load in MAP tgp model for later prediction Corrected par calculation in tgp.trees for non-default which argument when there is an unprinted height 1 tree Traces not recorded for linarea, parts, and ego, when they are off or fixed to unchanging settings Fixed EGO calculations; now calling them "improv", made them thread-safe, and added improv samples to traces Also now renamed ds2x to Ds2x, and made Ds2x (i.e., ALC) calculations thread-safe Ellipses to b* function allow other params arguments to be passed to replace defaults in tgp.default.params() Added traces of hierarchical parameters, and all predictive summaries (including kriging means and vars) Initial implementation of Importance Tempering idea Calculating variances of predictive distribution, and saving kriging means and variances Fixed numPrunable (P+1) calculation in grow attempts EGO calculations made thread-safe 1.1-11 (20 Sept 2006) ------ Bad compilation flags caused bad bug in log_HierPrior to go unnoticed until assertion failure on Windows 1.1-10 (19 Sept 2006) ------ Fixed big compute_ego() bug Tweaks to RNG seeds so that they give same sequence across platforms (tested on OSX/PB & Linux/Xeons) Allow X=NULL in dopt.gp() Removal of constants LOG_2_PI covered by R internal M_LOG_SQRT_2PI Added mvnpdf_log_dup and wishpdf_log functions Corrected full log-pdf function for tree prior, and added log-pdf of W and B0 1.1-9 (03 Sept 2006) ----- Many comments added in predict[_linear].c and exp_sep.cc Fixed ds2x (predict_linear) memory leak Fixed Matern bug -- now using correct pwr=1 distance computation; added comments to gen_covar.c Compilation warnings addressed in update of compat- ibility with C/C++ adaptive sampling (as) code Now getting BestPartitions from Tgp instead of at the end of Model::Rounds 1.1-8 (29 Aug 2006) ----- Added trace capability for Matern corr family Fixed bug that num_consumed is not updated on final NULL read from tlist Using new log_bessel_k() function to bypass allocs of floor(nu)+1 double each time bessel_k() was called 1.1-7 (23 Aug 2006) ----- Added the ability to get traces of all parameters, except the hierarchical ones Fixed bug to undo.m0r1 so that quantile differences don't have the mean subtracted out twice Changes to allow "..." arguments to mapT, and plotting of 2-d slices of input locations Minor fix to do params->Print() within model->Print(). Begin addition of mr_tgp, by Matt Taddy. When these changes are complete, we should be moving to version 2.0 1.1-6 (26 Jun 2006) ----- Added Latin Hypercube design function Added "mapT" function for stand-alone plotting of MAP partitions via projection and slice. Fixed bug in "tgp.plot.slice" in the process Also now calling "tgp.get.partitions" by the shorter name "partition" 1.1-5 (27 May 2006) ----- Fixed memory leak with RNG state in predict_consumer Solved -DPARALLEL adaptive sampling mystery exp2d.rand now returns more information about responses including ZZ, Ztrue, and ZZtrue 1.1-4 (22 May 2006) ----- Instructions for compiling with Pthreads. Got rid of duplicate printing of d[a,b] in tgp header. Added verb arguments to specify how verbose the R-console printing should be be. 1.1-3: (16 Apr 2006) ------ Necessary changes in order to get read_ctrlfile to work for the command-line (non-R) interface to tgp, and other minor changes. Noticed bug for bcart and b0tau options due to bad T- matrix initialization -- fixed. 1.1-2: (04 Mar 2006) ------ 1. Beta version of Matern correlation function working (abiet a bit shakily) 2. Tgp class and cleanup function now implemented in order to let tgp be interrupted, and then clean up to avoid memory leaks "on.exit()" 1.1-1: (26 Feb 2006) ------ Minor edits to documentation. 1.1: (22 Feb 2006) ---- 1. Modularization overhall of C/C++ code for passing to Matt Taddy, who will be responsible for such future additions as Matern correlation functions, and multi- resolution GP implementations. 2. Plotting enhancements, including 1-d projections and (now default) loess interpolation for 2-d surfaces and projections, plotting of ALC & EGO statistics, and ability to control plot layout to show mean surface, error surface, or both. 3. URL/Email now to UCSC/AMS address. 1.0-2: (Dec 2005) ------ 1. Wrote generic random deviate function which can use unif_rand(), randomkit, and erand48. 2. Fixed up some documentation bugs and typos. 3. system command in model.cc no longer calls "rm" which doesn't work on Windows. 1.0-1: (Dec 2005) ------ 1. No longer using erand48; using unif_rand() from R for Windows compatibility.