                 Changes in version 0.8.6 (2026-07-26)                  

    o   nlr and wnl5 no longer fail with "cannot change value of locked
	binding" when the user's workspace happens to contain an object
	named e.

    o   nlr recomputes the predictions at the final estimates, so that
	Prediction, Residual, the run test and dx no longer use the
	values left behind by the last hessian evaluation.

    o   nlr no longer stops with an uninformative error when the
	estimated error variance is very small; a non-positive variance
	is now penalised instead of producing NaN.

    o   nlr with Error="NOSG" works; the count of free error variance
	parameters was computed from a reversed sequence. Error="NOSG"
	is now documented.

    o   ExpandDH no longer overwrites the CMT of other dose records
	that fall inside an infusion interval, which sent a bolus into
	the wrong compartment.

    o   ExpandDH copies ID and covariates into the row it inserts at
	the end of an infusion, instead of setting every column to
	zero. The returned table is always sorted by TIME.

    o   hSkew finds the predictors from the model data rather than from
	terms, which failed or silently used same-named objects from
	the workspace.

    o   nComp handles Ka equal to an eigenvalue, which gave NaN.

    o   cmpChi returns 1 rather than 0 for two models of equal fit, and
	for a larger model that fits worse.

    o   pComp returns quietly when there is no rate to draw, and
	detects elimination by compartment number instead of
	compartment count.

    o   SolComp3 is robust to round-off near repeated roots.

    o   pProf draws the likelihood interval labels slightly above the
	curve minimum, so that they no longer overlap the profile
	curve.

    o   numDeriv moved from Depends to Imports, since only hessian is
	used. It is no longer attached to the search path by
	library(wnl); load it yourself if your own code calls it.

    o   Contact field added for bug reports.

                 Changes in version 0.8.5 (2026-03-07)                  

    o   Minor bugs/errata fix and refactored using Claude Code.

                 Changes in version 0.8.4 (2025-09-07)                  

    o   Default k value is adjusted to be more conservative.

                 Changes in version 0.8.3 (2025-03-23)                  

    o   A bug in ExpandDH is fixed.

                 Changes in version 0.8.2 (2025-02-21)                  

    o   logk in nlr adjusted

                 Changes in version 0.8.1 (2024-06-01)                  

    o   g2inv utility function bug with p=0 fixed.

                 Changes in version 0.8.0 (2024-02-22)                  

    o   g2inv utility function is revised to handle non-square matrix.

                 Changes in version 0.7.3 (2023-04-25)                  

    o   Default k value in nlr() is adjusted.

                 Changes in version 0.7.2 (2023-03-26)                  

    o   If hessian matrix is singular, g2inv() is used instead of
	solve().

    o   If k option in nlr() is missing, recommended k for conf.level
	is used.

    o   pProf() improved and exposed

    o   EnvObj() introduced

                 Changes in version 0.7.1 (2022-11-29)                  

    o   pProf() bug fixed for NaNs

                 Changes in version 0.7.0 (2022-10-28)                  

    o   hSkew() function added

    o   nlr function returns 1/8 likelihood interval which is similar
	to 95 percent CI.

    o   pProf() internal function plots likelihood profile and interval
	for the last nlr() run.

                 Changes in version 0.6.4 (2021-11-22)                  

    o   Minor beautification of the script - more spaces

    o   dx made compatible for rmarkdown; it further checks the default
	graphic parameter.

                 Changes in version 0.6.3 (2021-09-06)                  

    o   pComp option Col for compartment color added.

    o   pComp option asp=1 should be provided if you need.

                 Changes in version 0.6.2 (2021-08-07)                  

    o   pComp for plotting compartment model diagram is added.

                 Changes in version 0.6.1 (2020-08-27)                  

    o   qr.solve is used instead of solve.

    o   Documents enhanced.

                 Changes in version 0.6.0 (2020-07-16)                  

    o   nlr function gained Error="S" (scale) option. This enables GLS.

                 Changes in version 0.5.1 (2019-06-10)                  

    o   DAT table changed.

                 Changes in version 0.5.0 (2019-05-10)                  

    o   Analytical solutions of compartmental models added.

                        Changes in version 0.4.2                        

    o   Combined error with fixed additive variance model added.

                 Changes in version 0.4.1 (2018-06-14)                  

    o   Errata in manual corrected.

                 Changes in version 0.4.0 (2018-04-02)                  

    o   dx for simplest diagnostic plot is added.

    o   nlr and wnl5 returns more.

    o   nlr checks arguments more.

    o   Better p-value from cmpChi.

                 Changes in version 0.3.3 (2018-03-02)                  

    o   nlr returns more information.

    o   cmpChi is newly introduced.

                 Changes in version 0.3.2 (2017-12-01)                  

    o   Method argument added.

                 Changes in version 0.3.1 (2017-10-16)                  

    o   Secondary argument added.

                        Changes in version 0.3.0                        

    o   wnl5 added.

                        Changes in version 0.2.0                        

    o   Changed interface.

                        Changes in version 0.1.0                        

    o   First version
