0.3-1 ongoing
o bootCoef(): deleted.  No longer needed as bootCase() from car was modified to
    meet the same needs.
o confint.bootCase(), htest.bootCase(), hist.bootCase(), plot.bootCase():  Added,
    but were all originally *.bootCoef.
o compIntercepts: changed a cat()-based warning to warning().  Set the digits for
    printing the common covariate to be the same as digits=.
o Summarize(): slightly modified two warnings about factors on the RHS.


0.3-0 8-Nov-12
o Moved several functions to FSA that are used quite often for fisheries
    analyses.  Ultimately, I want to remove the dependency on NCStats.

o binCI(): moved to FSA.
o cltSim(): changed color of bars in both the population and sampling distribution.
o confint.nlsBoot(): moved to FSA.
o fact2num(): moved to FSA.
o gofCI(): removed dependence on binCI() (which is now in FSA) by using binconf()
    from Hmisc directly.  Added a type="goodman" that will perform the method
    from Goodman (1965) (thanks to Paul Rabie for the discussion and code!).
o htest.nlsBoot(): moved to FSA.
o hist.formula(): modified.  Slightly changed mgp() in par().
o hyperCI(): moved to FSA.
o ks2d(): moved to FSA.
o ks2dp(): moved to FSA.
o ksTest(): moved to FSA.
o lagratio(): moved to FSA.
o lsmean(), and related internals: moved to FSA.
o plotH(): moved to FSA.
o poiCI(): moved to FSA.
o popSizesPlot(): moved to FSA.
o pos2adj(): moved to FSA.
o rcumsum(): moved to FSA.


0.2-8 21Jun12
o Switched to testing under R version 2.15.0 (32-bit).
o Changed license specification from "GPL version 2 or newer" to "GPL (>= 2)" to
    avoid warning on check.


o bootCoef(): Added.  Basically a wrapper to bootCase() in the car package with
    the ability to call some S3 methods (see below).
o confint.bootCoef(): Added.  This is essentially the (very) old ci.bc().
o faqNC(): Added.
o fitPlot.glm(): Added.
o fitPlot.logreg(): Added.
o hist.bootCoef(): Added.
o htest.bootCoef(): Added.  This is essentially the (very) old ho.bc().
o levenesTest(): Added.  A convenience wrapper to leveneTest().
o logregPlot(): Removed.  Functionality is now in fitPlot.logreg().
o plot.bootCoef(): Added.
o plotBinResp(): Added.  Used in fitPlot.logreg() and as stand-alone.
o poiCI.rd: Updated details and references in the help page based on information
    provided by Jerry Lewis about the historical use of the method described in
    Ulm.


0.2-7 2Mar12
o added an importFrom for multcomp for glhtSig() below.

o .onLoad(): Modified.  Moved the startup message into packageStartupMessage()
    in hopes of eliminating the warning when checking the package.
o addMargins(): Modified.  Corrected error when dealing with 1-dimenstional
    tables and handling call from percTable() correctly with no margin=.
o addSigLetters(): Modified.  Added code to deal with a bug that appeared if
    the x-axis variable when adding to a two-way ANOVA interaction plot was an
    ordered, rather than regular, factor.
o chisqPostHoc(): Added.
o distrib(): Modified.  Added a digits= argument for printed answers.
o glhtSig(): Added.
o percTable(): Modified.  Added addMargins= argument.  Corrected problem with
    2-dimensional table and call to addMargins with margin=NULL argument.
o sumTable(): Added.
o transChooser(): Modified.  Added alpha= argument and set outlier test p-values
    that are greater than 1 to 1.
o typeoflm(): Modified.  Fixed a call to the deprecated lgrep().


0.2-6 1Oct11
o Switched to compiling under R version 2.13.1 (32-bit).
o Removed utils from importFrom because it is a base package and due to the
    removal of functions below.
o Removed tools from importFrom due to the removal of functions below.
o Removed suggested packages from description.

o findUpdates(): Removed.  Moved to miscOgle package. 
o ksTest(): Added.
o swvCode(): Removed.  Moved to miscOgle package.
o swvCounts(): Removed.  Moved to miscOgle package.
o swvPvalue(): Removed.  Moved to miscOgle package.
o swvANOVA(): Removed.  Moved to miscOgle package.
o swvGLHT(): Removed.  Moved to miscOgle package.
o swvHtest(): Removed.  Moved to miscOgle package.
o swvREG(): Removed.  Moved to miscOgle package.
o swv2ad(): Removed.  Moved to miscOgle package.
o swv2PDF(): Removed.  Moved to miscOgle package.
o updateNCStats(): Removed.


0.2-5 19Aug11
o Modified description file to show my e-mail address.
o Added an ImportFrom from tools for texi2dvi used in swv2pdf().

o ABCens90: Added this data file.
o binCI(): Removed reference to an unused type in the .Rd file.
o c.region(): Modified.  Changed all arguments (except cex.ans=) so that they
    did not have defaults.  Also added an argument, lbl.col, to allow the axis
    label to be a different color then the shaded area.  Changed the show.axis=
    argument to lbl.axis=.  Re-ordered the arguments.
o distrib(): Modified.  Changed label for normal distribution so that the sd=
    value is rounded to three decimal places.  This makes for a "cleaner" plot
    when the provided sd value is really a computed standard error.  Added a 
    shade.trans= argument to allow the shaded area to have a transparency value;
    this will allow the plots to be a bit softer to the eye.  Added a lbl.col=
    argument for the axis label so that it is not printed in the transparent
    color if transparency is used.  Modified the call to c.region() per the 
    changes above.
o expandTable(): Modified.  Added methodology to deal with a "contingency table"
    with only one row or one column.  Also fixed a bug that implied that three
    dimensional tables could be used.
o geomean(): Added.
o geosd(): Added.
o gofCI(): Added.
o meanMedian(): Modified.  Added an argument that allows one to plot user-derived
    data rather than just the data simulated from a beta distribution.  Also
    streamlined the code by adding two internal functions -- meanMedianPlots()
    and meanMedianBeta().
o percTable(): Added.
o plot.htest(): Modified.  Modified the call to c.region() per the changes above.
o plotH(): Modified.  Modified to plot variable names on axes by default.  Fixed
    error in how x- and y-axes were labeled if a formula was not used.
o powerSim(): Modified.  Modified the call to c.region() per the changes above.
o rSquared(): Added.
o swvCode():  Modified.  Modified to allow file= argument to include the .rnw
    extension and still output a file with a .r extension.
o swv2PDF(): Added.


0.2-4 15Jun11
o Switched to compiling under R version 2.13.0.

o accuracyPrecision():  Modified.  Added several arguments that allowed more control
    over the plottings, changed the default number of points plotted, changed the
    default pch used for plotting the mean, renamed the internal functions (to
    camel code), removed the make.labels() internal function (incorporated into
    the makeTargets() internal function), and changed the default coloring for the
    rings on the targets.  Functionality is unchanged.  Hopefully it is a bit more
    general.
o ABDLakes: Added this data file.
o Frogs: Added this data file.
o hist.formula(): Modified.  Added the ylab= argument with a default to "Frequency".
    Added the right= argument with a default to FALSE; this is opposite of the
    default used in base hist().  Added warn.unused=FALSE arguments to the two
    calls to hist() with plot=FALSE so that the warnings would be suppressed.
o lsmean(): Added.  This is the exact lsmean.default() and lsmean.lme() functions
    from Yandell's pda package, including three dependent functions which were
    kept as internal functions in NCStats.  The pda package does not seem to work
    with R >2.13.0.  To make this function available to those that use it I put
    it in NCStats.
o SquareLakePopn: Added this data file (from FSAdata).
o srsdf(): Added.
o updateNCStats(): Modified.  Added the tools package to the packages that should
    not be detached when updating.  Also made sure that the package.dependencies
    function is from the tools package.  This eliminated a warning when checking
    the packages.


0.2-3 18Apr11
o Added importFrom utils for installed.packages() and available.packages().

o compIntercepts(): Modified.  Changed catches at front of function to note if
    a model with an interaction was sent and, if so, that the interaction term
    would be ignored.  Also added an if that will allow the RHS of the model
    formula to have the covariate and factor in either order.
o findUpdates(): Added.
o lgrep(): Modified.  Changed to be a pass-through to grepl() and added a message
    telling the user to use grepl() instead.  Basically I just found grepl() and
    no longer need lgrep().
o pos2adj(): Added.  Allows ability to place point labels according to directional
    coordinates (e.g., "NE").
o residPlot(): Modified.  Corrected an error related to labeling points if they
    were an outlier.  Basically needed to realize that the studentized residuals
    are a named vactor.
o Subset(): Modified.  Added a catch to send a warning message if the resultant
    data frame has no rows in it.
o swvCode(): Modified.  Added the moreItems= argument to allow the user to more
    easily add more items to be removed shen Stangling.  Basically the items in
    moreItems get appended to the items in itemsToRemove.  In addition, changed
    the ItemsToRemove= argument to itemsToRemove (note capitalization).
o swvHtest(): Added.
o updateNCStats(): Added.
o z.test():  Added.  Basically same as z.test() from TeachingDemos package.
    Modified for my use at Northland.
    

0.2-2 3Mar11
o moved to compiling under 2.12.1 (32 bit)
o changed dependency to >2.11.1
o added a dependency to tcltk package (see transChoooser below).

o addMargins(): Added.  Largely the same as addmargins() except that I changed
    the margin= argument such that margin=1 gives a marginal value at the end of
    the rows and margin=2 gives a marginal value at the end of the columns.  This
    change makes addMargins() fit more closely (in my opinion) with the use of
    margin= in prop.table().
o fact2num(): Added.  Changes a factor to numeric values.  Useful for example
    if age or year is a factor but you need it to be numeric.
o fitPlot.nls(): Modified.  Deleted the col= argument and added col.pt= and 
    col.mdl= arguments to allow for different colors for the points and the
    fitted line.
o plotH(): Modified.  Added catches for determining if Y variable is not quantitative
    and if more than one X (or RHS) variable is sent (i.e., function only works
    for one explanatory variable).  Created separate internal functions for whether
    the X variable is quantitative (the old plotH()) or categorical (simply a 
    pass through to barplot).
o predictPlot(): Added.  Simply a pass-through to (i.e., does the exact same thing
    as) predictionPlot().  Simply for convenience.
o residualPlot(): Added.  Simply a pass through to (i.e., does the exact same thing
    as) residPlot().  Sipmly for convenience.
o Summarize(): Modified.  Removed dimnames attribute on the level labels in the
    resulting dataframe when a quantitative variable is summarized by a categorical
    variable.  Mostly for aesthetic reasons.  Added a digits= argument to be more
    compatible with other R functions.  Retained the numdigs= argument for
    backward compatability.
o swvCode(): Modified.  Added two more items to ItemsToRemove= argument to handle
    changes to Stangle() output.  Added a blanks= argument that can be used to
    remove all or just "extra" blank lines in the Stangle() output.  Also added
    the annotate=FALSE argument to set as default for sending to Stangle().  This
    forces Stangle() to not include code block comments.
o transChooser(): Modified.  Added require(tcltk) to attempt to ensure that this
    package is loaded.  Corrected an error which caused the Levene's test
    results not to be printed above the residual plot when show.stats=TRUE for
    the ONEWAY and TWOWAY functions.  Corrected an error which caused overprinting
    of the outlier test p-value if outlierTest() returned more than one value 
    (only the minimum will be now be printed).  Created internal functions for
    placing the p-value labels above each graph (reduced redundancy).  Removed
    several calls to round() inside of a formatC() call.
    

0.2-1 31Jan11
o removed dependency on nortest package (see adTest() below).

o adTest(): brought, with name change, ad.test() from nortest package into NCStats.
    I changed to 'camel-case' naming convention and then removed dependency on 
    the nortest package.
o Summarize(): moved a bunch of code into internal functions.  Modified the code
    for a quantitative variable by a factor variable so that the result prints
    as a data.frame (rather than a matrix) that includes one or two columns that
    contain the levels of the factor variables.  This provides a nicer looking
    output and an output that could be subsetted if needed.  Also modified the
    code so that multiple "values" of a factor variable can be excluded.  The
    help code examples were modified to illustrate the changes.
o transChooser(): updated ad.test() references to adTest().
    

0.2-0 1-Dec-10
o added a work-around .Rd files for the three functions imported from the car
    package namespace and exported as part of the NCStats namespace.  This seems
    to have corrected the warning produced during Rcmd check.

o distrib(): modified.  Added the se= argument in the hopes that this might be
    somewhat easier for students that are learning to distinguish between
    sampling and population distributions.
o hist.formula(): modified.  Added the ability to create a histogram of a 
    variable not separated by groups with code such as ~x, x~1, or x~0.  This
    allows one to use the data= argument for a single histogram.  Added a
    main= argument for the situation where a single histogram is plotted and set
    the default to "" to avoid the annoying main plot label.  Added a byrow=
    argument so that the user can choose the order in which the plots are populated.
    Changed ylmts= to ymax= argument so as not to confuse ylmts= with ylim=.
    Changed main.pre= to pre.main= to avoid confusion with partial matching with
    main= argument.  Changed cols= argument to ncol= argument (and rows= to nrow=)
    to avoid confusion with partial matching of col= that may be sent in the
    "dots" argument (and to make consistent with notation in matrix()).  Modified
    code such that if pre.main=NULL then no main title will be printed on the
    histograms.  Removed the breaks= argument (will work as part of dots argument).
    Removed the code that "caught" if the formula= argument was not a formula --
    this was redundant because this is an S3 generic for formula types.
o .onLoad(): modified.  Changed to include version number of loaded version.
o plotH(): added.
o residPlot(): modified.  Added code to plot multiple points returned by 
    outlierTest().  Add the use of thigmophobe.labels to place outlier text marks
    with a minimum of overlapping.  Added three internal functions -- 
    addOutlierTestResults(), addLoessLine(), and getMainTitle() to reduce and
    streamline code among the different generics.
o spois(): modified.  Fixed small bug (extra comma).
o summarize.formula(): modified.  Added the ability to create summaries of a 
    variable not separated by groups with code such as ~x, x~1, or x~0.  This
    allows one to use the data= argument for a single variable.  Modified default
    exclude= value to better handle these situations.
        

0.1-5 23-Aug-10
o renamed functions to camel-type case naming.

o accuracyPrecision(): renamed from accuacyPrecision().
o addSigLetters(): renamed from add.sig.letters().
o binCI(): renamed from bin.ci().
o chooseColors(): renamed from choose.colors().
o ciSim(): renamed from ci.sim().
o cltSim(): renamed from clt.sim(). 
o compIntercepts(): renamed from comp.intercepts().
o compSlopes(): renamed from comp.slopes().
o detachAll(): renamed from detach.all().
o diagPlot(): renamed from diag.plot().  Updated for fitPlot() and 
    residualPlot() changes.
o expandTable(): renamed from expand.table().
o extract.data(): deleted.  Capabilities are now found in the digitize package.
o fitPlot(): renamed from fit.plot().
o fit.plot(): deprecated.
o hoCoef(): renamed from ho.coef().
o hyperCI(): renamed from hyper.ci().
o legendHelp(): renamed from legend.help().  This is an internal function.
o logregPlot(): renamed from logreg.plot().
o meanMedian(): renamed from meanVSmedian().
o mrnorm(): renamed from mult.rnorm().
o paletteChoices(): renamed from palette.choice().  This is an internal function.
o poiCI(): renamed from poi.ci().
o powerSim(): renamed from powerSim().
o popSizesPlot(): renamed from pop.sizes.plot().
o predictionPlot(): renamed from prediction.plot().
o rcumsum(): added from FSA package.
o residPlot(): renamed from residual.plot().  Did not rename as residualPlot()
    because of a conflict with residualPlots() from car package.  Updated for
    fitPlot() change.
o residual.plot(): deprecated.
o sdCalc(): renamed from sd.calc().
o transChooser(): renamed from trans.chooser().  Updated for fitPlot() and 
    residualPlot() changes.
o trans.chooser(): deprecated.

0.1-4 22-Aug-10
o moved to compiling under 2.11.1.
o updated namespace to handle name changes from car package.

o ci.bc(): Deleted.  No longer use boot.case(), use nlsBoot() instead.
o confint.nlsBoot(): Added.
o diag.plot(): Modified to handle name changes from car package.  Also had to
    modify to handle different output from outlier.test.
o ho.bc(): Deleted.  No longer use boot.case(), use nlsBoot() instead.
o htest.nlsBoot(): Added.
o ks2d(): Added along with generics.
o ks2dp(): Added along with generics.
o lagratio(): Added.
o pop.sizes.plot(): Added.
o residual.plot(): Modified to handle name changes from car package.  Also had to
    modify to handle different output from outlier.test.
o trans.chooser(): Modified to handle name changes from car package.  Also had to
    modify to handle different output from outlier.test.


0.1-3 5Jun10
o hist.formula(): added the ability to have two factors on the RHS of the formula
    to allow users to separate the quantitative response variable according to
    two factors.  Also modified the main.pre= argument.  Modified the code so
    that if the product of the number of rows and number of columns supplied by
    the user is less than the total number of histograms to be plotted then
    multiple pages of histograms will be produced (each requiring the user to
    click on it to go to the next one).
o identify.formula(): added this function so that identify() could be used with
    a formula which is easier to do following using a formula in plot.
o power.sim(): added a label for beta, set minimum for n slider at 2.
o Summarize.formula(): added the ability to have two factors on the RHS of the
    formula when working with a quantitative response variable.  Also removed an
    extraneous (from debugging) print() in the "categorical" portion.
o swvCode(): added this function to remove some "unwanted" lines from files
    created with Stangle().  The function currently removes lines with the
    words "Stangle", "SweaveHooks", and "swvCode".  The function will have to be
    modified to allow removal of other words.


0.1-2 17Dec09
o moved to compiling under 2.10.1.

o NCStats-package(): updated.
o bin.ci(): updated to better handle a first argument that is a vector.


0.1-1 15Apr09
o added a namespace
o removed dependencies and changed to imports ...
o   multcomp removed completely (only used in an example).
o   nortest and sciplot left as dependencies (they do not have a namespace).
o   added a dependency to plotrix (does not have a namespace) because I needed
      thigmophobe() in highlight.

o .FirstLib(): removed (changed to .onLoad() because of namespace).
o comp.slopes(): added the order.slopes= argument to allow the user to choose whether the
    slopes portion of the output should be ordered from smallest to largest or not.
    Added digits= argument to control printing of significant digits.  Removed attachment
    of data frame in the example in the .rd file.
o highlight(): changed to a S3 method with default and formula methods.  This allowed
    use of the formula with the data= argument.  Also included a new argument that
    allowed labeling the points with values in another variable rather than just the
    number.  Allowed a pos= argument with the default being thigmophobe() positions.
o print.CompInt(): added digits= argument to control printing of significant digits.
o sd.calc(), print.sdcalc(): added.
o Summary(): removed for good (was deprecated).
o Subset.rd: removed drop.levels() example.
o plot.htest.rd: removed z.test() example.
o removed BOD, Consumption, CrabClaw, LakeTroutEggs, MooseBrowse, NAKATPase,
    Opossums, PeakDischarge, Tadpoles, and TomatoeNematode.


0.0-13  15Jan09
o moved to compiling under 2.8.0

o add.sig.letters(): minor corrections, slight modification of examples, and added note to .Rd.
o attached(): added notes to .Rd.
o detach.all(): added notes to .Rd.
o diag.plot(): minor corrections, updated examples, and added note to .Rd.
o distrib(): added notes to .Rd.
o expland.table(): changed use of x$Freq[i] to x[i,"Freq"] in hopes of removing the warning about "no visible binding."
    This did not fix the warning because of the use of "Freq" in the next line.
o extract.data(): added notes to .Rd.
o fit.plot(): minor corrections and added note to .Rd.
o hist.formula(): added notes to .Rd.
o lgrep(): minor corrections and changes to .Rd.
o mult.rnorm(): minor corrections, updated examples (deleted tapply() and used Summarize()), and added a note in .Rd.
o print.anova(): added as a replacement function.
o residual.plot(): minor corrections, updated examples, and added note to .Rd.  Also corrected x-axis label (changed
    from 'fitted values' to 'treatment groups') for residual plots using boxplots for one-way and two-way ANOVAs.
o Subset(): added an error message for sending non-data.frames.  Added notes and correction minor errors in .Rd.
o Summarize(): added notes to .Rd.
o Summarize.default(): removed exclude= argument and moifided code around valid percents
    for the categorical data.  This was required because of a change in the use of exclude=
    in table() for 2.8.0.  Modified code to catch more errors.
o Summarize.formula(): added data= argument.  Modified code to catch more errors and to
    provide better tables for categorical data.
o swvANOVA(): added.
o swvGLHT(): added.
o swvPvalue(): modified so that it may (default) include a 'p=' prefix.
o swvREG(): added.
o trans.chooser(): minor corrections, updated examples, and added note to .Rd.


0.0-12  17Nov08
o added a dependency to tcltk package to make the simulations work properly -- fix something that was changed
    in the TeachingDemos package.
o added CITATION file.

o distrib(): added "forward" and "reverse" to the type= argument.  Modified so that the x-axis labels on the
    normal distribution plots correspond to exact integer standard deviation away from the mean.  This change
    is purely for pedagogical reasons.
o hist.formula(): added a ylmts= argument to allow user to set maximum of y-axis separately for each histogram.
    Removed data(iris) from help file (because all examples use the data= argument).
o Subset(): added.
o swvCounts(): added.
o swvPvalue(): added.


0.0-11	15Sep08
o .First.lib: added to print startup message and check for installed packages.
o c.region(): modified by returning a list of values used in the plotting and adding a plot= argument.  These
    two changes allow the user to use c.region to find the plotting values without making the plot.  This was
    needed for the plot.htest() function.
o d.region(): modified as described for c.region().
o distrib(): modified by adding a show.alt= argument which allows user control of whether the message showing
    alternative code for obtaining the same answer is shown or not.  Defaults to FALSE (don't show message).
    Also added a shade.col= argument to control the color of the shaded area.  Added code to work-around a
    problem of constructing plots when the value of interest is one of the values along the x-axis.
o fit.plot.POLY(): added this S3 method to plot polynomial regressions.  It basically calls the SLR method.
o fit.plot.TWOWAY(): corrected error when more then 6 line types were needed (max is 18, though).
o hist.formula(): added this S3 method of allos plotting of multiple histograms at once.  This replaces mhist().
o mhist(): deleted and replaced with hist.formula().
o plot.htest(): added this S3 method to provide a visual of the p-value computed for the one-sample z-test (from
    the TeachingDemos package), one-sample t-test, paired t-test, two-sample t-test, and chi-square tests.
o residual.plot.POLY(): added this S3 method to plot polynomial residuals.  It basically calls the SLR method.
o Summarize(): this replaces the old Summary() and adds the ability to use a formula to get summary statistics
    by levels of a factor.  Slightly modifed help file.
o Summary(): this is now deprecated ... it simply passes its arguments to Summarize.
o typeoflm(): added the ability to detect polynomials.
0 view(): modified by adding which= argument to allow choosing variables/columns to display.


0.0-10	15May08
o Moved to RForge.net

o distrib(): modified, now returns answer regardless of plot= argument.
o extract.data(): modified, added xdec=, ydec=, and filename= arguments.
o lgrep(): added.
o rhead(): removed, replaced with view()
o view(): added, replaced rhead()


0.0-9	22Apr08
o upgraded to R 2.7.0.
o changed dependencies to include sciplot package.

o choose.colors(): added, to allow users to choose among a variety of palettes.  See also palette.choices().
o fit.plot.ONEWAY(): modified, uses lineplot.CI() from the sciplots package to allow the plotting of confidence
    intervals on the one-way fit.plot.  Other related controls were also implemented.
o fit.plot.TWOWAY(): modified, uses lineplot.CI() from the sciplots package to allow the plotting of confidence
    intervals on the two-way fit.plot.  Other related controls were also implemented.
o fit.plot.IVR(): modified, allowed color choices through choose.colors() and placed legend in topright by default.
o palette.choices(): added, aids choose.colors().
o residual.plot(): modified, changed mdl= argument to object= argument to make more compatible with fit.plot().
    Created SLR, IVR, ONEWAY, and TWOWAY methods.  This caused some over coding but clarified the construction of the
    IVR residual plots and makes more consistent with fit.plot().  Along the way, changed the ONEWAY and TWOWAY to have
    bp= argument that defaults to TRUE so that boxplots of the residuals are shown rather than points.  Also modified
    the colors used in the IVR.
o residual.plot.IVR(): added, see also residual.plot().
o residual.plot.ONEWAY(): added, see also residual.plot().
o residual.plot.SLR(): added, see also residual.plot().
o residual.plot.TWOWAY(): added, see also residual.plot().
o trans.chooser.ONEWAY(): modified, added bp= argument because of changes to residual.plot().
o trans.chooser.TWOWAY(): modified, added bp= argument because of changes to residual.plot().


All prior NEWS in ONEWS file. 