1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 172 173 174 175 176 177 178 179 180 181 182 183 184 185 186 187 188 189 190 191 192
|
Software and datasets to support 'Modern Applied Statistics with S',
fourth edition, by W. N. Venables and B. D. Ripley.
Springer, 2002, ISBN 0-387-95457-0.
This file documents software changes since the third edition.
Version 6.1 (Apr 1999) release for the third edition. Support for S-PLUS
5.1 and 2000, added more na-handling.
Version 6.2 (Aug 2000) tidy up documentation. Would be 6.1pl010, but that
causes R problems.
(Oct 2000) support for S-PLUS 6.0 for Unix.
(July 2001) support for S-PLUS 6.0 for Windows.
Version 6.3 (Dec 2001) add material destined for fourth edition.
Version 7.0 (Mar 2002) pre-release for the fourth edition. S-PLUS 6.x
and R only.
Version 7.1 (Mar 2003) add namespaces, eliminate data() calls
Version 7.2 (Apr 2004) for R 1.9.0.
Library MASS
------------
- eqscplot has new arguments ratio and uin.
- stepAIC will not drop strata terms in coxph or survreg models.
- profile.glm will report inadequate supplied glm fits, not just fail.
- new method confint.lm.
- fractions/rational allow missing values.
- mvrnorm has an 'empirical' argument.
- predict.lda and predict.qda try harder to avoid exponential underflow.
- new function fitdistr for ML estimation of univariate distributions.
- new function glmmPQL to use lme to fit GLMMs by PQL
- truehist allows rule for nbins to be specified as a character string.
- parcoord function.
- new datasets bacteria, epil, nlschools, SP500
- polr allows control argment for optim, reports lack of convergence.
- stepAIC works again if formula has an offset (R had changed).
- biplot.correspondence now shows the origin as a cross.
- polr was not preserving contrasts to put in the fit object.
- vcov methods for lme, gls, coxph and survReg.
- Added 'tol' argument to isoMDS.
- stepAIC now allow 'direction=both' starting from a full model.
- glm.nb allows R-style 'start' argument.
- truehist passes ... on to both plot.default() and rect().
- isoMDS now uses the C interface to optim.
- addterm, dropterm, stepAIC now work with lme and gls fits.
- huber checks for MAD equal to zero.
- glmmPQL now loads nlme if not already loaded.
- glmmPQL handles list 'random' arguments (7.0-11).
- The MASS datasets no longer require data(foo) to load them. (7.0-11)
- mvrnorm uses eigen(EISPACK=TRUE) for back-compatibility (7.0-11, R 1.7.0)
- print.summary.polr could lose dimnames for 1 coefficient.
- remove heart as survival in R now has it.
- confint.{lm,glm} didn't handle specifying parm in all cases.
- confint and confint.lm have been migrated to base in R.
- addterm.default, dropterm.default and stepAIC work better inside functions.
- glm.nb now sets AIC in the object, and has a logLik() method.
- truehist now accepts a 'ylab' argument.
- negative.binomial and neg.bin no longer generate objects with
package:MASS in their environment.
- stepAIC now drops (if allowed) 0-df terms sequentially from the right.
- lda(CV=TRUE) now works for rank-deficient fits.
- predict methods for lda, polr now check newdata types.
- model.frame.lda/polr now look for the environment of the original formula.
- polr has a new `model' argument defaulting to TRUE.
- fitdistr supports the trivial case of a Normal distribution.
- sammon and isoMDS now allow missing values in the dissimilarity matrix, and
isoMDS allows Minkowski distances in the configuration space.
- cov.trob works better if wts are supplied, and may converge a little faster
in any case.
- The ch11.R script now uses mclust not mclust1998.
- The default xlab for boxcox() is now greek lambda.
- glmmPQL now handles offset terms.
- add predict.rlm method to correct predict.lm in the case se.fit=TRUE.
- weighted rlm fits are handled better, and default to "inv.var".
- logtrans works without specifying 'data'.
- predict() method for glmmPQL.
- polr() has an option for probit or proportional hazard fits.
- neg.bin() and negative.binomial() had an error in the aic() formula.
- The ch05.R script now includes the code for Figure 5.8.
- Datasets austres, fdeaths, lh, mdeaths, nottem and rock
are now visible in the 'datasets' package of R 2.0.0 and so have
been removed here.
- Script ch07.R now gives details using the gam() function in package gam as
well as that in package mgcv.
- rlm's fitted component is now always unweighted.
- theta.{md,ml,mm} now have one help file with examples.
- polr() has a new method "cauchit" suggested by Roger Koenker.
(Requires R >= 2.1.0)
- polr() now works with transformed intercepts, and usually converges
better (contributed by David Firth).
- polr() handles a rank-deficient model matrix.
- polr() now returns the method used, and uses it for predictions.
- anova() method for polr (contributed by John Fox).
- predict.glmmPQL was not using the na.action in the object as intended.
- The default methods for addterm and dropterm and anova.polr now check
for changes in the number of cases in use caused e.g. by na.action=na.omit.
- Added vcov() method for rlm fits.
- eqscplot() accepts reversed values for xlim and ylim.
- Script ch10.R uses se.contrast to calculate se's missing from model.tables.
- profile() and confint() methods for polr().
- glm.convert() was not setting the `offset' component that R's glm objects
have.
- sammon() now checks for duplicates in the initial configuration.
- isoMDS() and sammon() work around dropping of names.dist in 2.1.0
- lda() now gives an explicit error message if all group means are the same.
- fitdistr() now has a logLik() method, chooses the optim() method if not
supplied, handles the log-normal by closed-form and no longer attempts to
handle the uniform.
- glm.nb() now accepts 'mustart'.
- glm.nb() now supports weights: they used to be ignored when estimating
theta.
- fitdistr() now supports geometric and Poisson distributions, and
uses closed-form results for the exponential.
- lm.ridge, lqs and rlm allow offset() terms.
- the 'prior' argument of predict.qda is now operational.
- script ch12.R now has b1() adapted for R's contour().
- anova.polr() quoted model dfs, not residual dfs.
- stepAIC() applied to a polr fit now gets the correct rdf's in the
anova table.
- lm.gls() now returns fitted values and residuals on the original
coordinates (not the uncorrelated ones).
- parcoord() now allows missing values and has a new argument
'var.label' to label the variable axes. (Contributed by Fabian Scheipl.)
- rlm() has a 'lqs.control' argument passed to lqs() where used for
initialization.
- rlm() could fail with some psi functions (e.g. psi.hampel) if 'init' was
given as a numeric vector.
- rlm() handles weighted fits slightly differently, in particular trying
to give the same scale estimate if wt.method="case" as if replicating the
cases.
- confint.nls copes with plinear models in R (now profile.nls does).
- The wrappers lmsreg() etc have been adapted to work in the MASS namespace.
- qda() accepts formulae containing backquoted non-syntactic names.
- polr() gives an explicit error message if 'start' is misspecified.
- glmmPQL() evaluates the formulae for 'fixed' and 'random', which may
help if they are given as variables and not values.
- There are anova() and logLik() methods for class "glmmPQL" to stop misuse.
- profile.polr() now works for a single-coefficient model.
- The print and print.summary methods for polr and rlm make use of
naprint() to print a message e.g. about deleted observations.
- Class "ridgelm" now has a coef() method, and works for n < p.
Library class
-------------
- fixed a bug in knn(use.all = FALSE)
- added SOM, batchSOM and ancillary functions.
- Bugfix in SOM (7.0-5)
- knn* and lvq* explicitly check for missing values in the inputs.
- knn1 was not handling ties correctly in some cases (but knn was).
- added check for user error in inputs in knn.cv.
- SOM segfaulted for 1D inputs.
Library nnet
------------
- no copying of datasets even in R.
- model.frame method for multinom (even in R).
- nnet now uses the C interface to optim.
- nnet.Hess has been renamed nnetHess.
- vcov.multinom now computes the Hessian analytically (thanks to David Firth).
- predict methods for multinom, nnet now check newdata types
- model.frame.multinom now looks for the environment of the original formula
- multinom has a new `model' argument defaulting to TRUE.
- the multinom methods for add1, dropterm and anova now check
for changes in the number of cases in use caused e.g. by na.action=na.omit.
- added confint() method for multinom.
- added logLik() method for multinom.
Library spatial
---------------
- Methods for class "surf.ls" contributed by Roger Bivand.
- surf.ls and surf.gls now check for x,y,z cols in a data frame given for x.
- PP datasets moved from data to ppdata
- no longer duplicate topo dataset from MASS
- Strauss and SSI are now interruptible
- the extractAIC method for "trls" objects was using incorrect edf
|