StatLib---Software and extensions for the S (Splus) language

To contribute software to this archive, please obtain a copy of the submission instructions: send mail to and ask send submissions from S The software is not warranteed in any way. Unless otherwise stated the software is exactly what is provided by the submitters. All questions and comments should be directed to the submitter.

Some of the entries are shar archives. If you don't know how to deal with a shar archive, send the message

send shar from general for instructions.

The S archive currently contains:

You might also wish to look at the S DOS/WINDOWS collection.
The instructions for how to submit software to statlib. [1/May/89] (1673 bytes)
Topical index of the many fine software routines for S and S-PLUS which are available from StatLib. Maintained by Charles Roosen, [28/Aug/95][28/Feb/96] (5k)
Piece together an array from several vectors and/or arrays. This is a multidimensional generalization of the cbind and rbind functions. Written by Tony Plate ( and Richard Heiberger ( [22/Oct/96] (33k)
Fits the very general accelerated failure model with log-F error to possibly censored data. (Barry W. Brown) [2/Oct/95] (440k)
TEX source to accompany accflf. (Barry W. Brown) [2/Oct/95] (25k)
Perform ACE transformation detection on Generalized Linear Models. Extension of ace.logit. Simon Byers ( [13/Dec/96] [24/Mar/97](5k)
Calculate and plot the ACE transformations for logistic regression (adapted from linear regression). Adrian E. Raftery, [23/Jan/95] [13/Dec/96] (4k)
An S function that calls Genz's "adapt" suboutine for adaptive quadrature (integration) in from 2 to 20 dimensions. Adapted by Mike Meyer ( [17/Oct/89](31709 bytes)
Additive hazards regression alternative or supplement to the Cox model. It gives plots that are shows the effect of covariates on survival. A gzipped tar file. Submitted by Harald Fekjaer ( [14/Oct/97] (74k) See also the Windows version
Obsolete, see survival2 below.
S functions for classical and resistant analysis of averages and the analysis of variance using the sweep operator. Created by Colin Goodall ( [6/July/89][24/March/92](70342 bytes)
S-function to compute robust estimates of partial autocorrelation for AR processes (upto order 6). This function can also be used to obtain posterior probabilities for the competing AR models. Nhu Le ( [1/Mar/95] (99 k)
A shar archive of three separate sets of programs for 1D, 2D and 3D density estimation. The last, ashn, is combined with a visualization program. David Scott, ( [14/Apr/93] (84kbytes)
Provides asymptotic sample size, power, and detectable effect calculations for a variety of generalized linear models and survival models. Shar file. (Barry W. Brown) [2/Oct/95] (197k)
TEX dvi file for paper on asypow. Will not e-mail well. (Barry W. Brown) [2/Oct/95] (47k)
S-plus software to perform Bayesian estimation for an autologistic model with covariates. A shar archive. Submitted by Jennifer Hoeting ( [3/Mar/99][9/Sep/99][5/Jul/00] (155k)
A shar arhive of new S code for avas (additivity and variance stabilization). It is used to estimate transformations for regression. Contributed by Rob Tibshirani ( [7/Oct/89](51150bytes) {Revised, Oct 7. 89}
Balloon is a package of Splus functions for creating skew-t plots and winds profile plots for rawinsonde data. Submitted by Eric Gilleland ( [9/Nov/99] (196 kbytes)
A shar file containing all S functions and datasets from "Statistics for Long-Memory Processes" by Jan Beran are provided along with a LaTeX file listing the data. Contributed by Brandon Whitcher ( [24/Mar/97] (124k)
Bayesian model selection and accounting for model uncertainty in linear regression models, using the BIC approximation. Contributed by "Adrian E. Raftery" ( [17/Oct/94] [23/Jan/95] [13/Dec/96] [14/Oct/97] [2/Jun/98] (8 kbytes)
An S-PLUS function to implement Bayesian Model Averaging for Generalized Linear Models using BIC. Contributed by Chris T. Volinsky ( [2/Jun/98] (19k)
Logistic regression with Bayesian model selection and accounting for model uncertainty, using the BIC approximation Contributed by "Adrian E. Raftery" ( [29/Nov/94] [28/Feb/96] [2/Jun/98] (15 kbytes)
S-Plus function to implement Bayesian Model Averaging for Cox Proportional Hazards Models, using the BIC Approximation. Includes a summary function. Contributed by Chris T. Volinsky ( [1/Dec/95] [28/Feb/96] [13/Dec/96] [2/Jun/98] (17 kbytes)
Plots a bivariate normal density. It is a dumped file and it can be read into Splus (or R for that matter) by simply typing restore("biv.norm.dmp"). Submitted by Veith Tiemann ( [9/Nov/99] (1053 bytes)
Functions for nonparametric bivariate survival curve analysis. Requires dynamic loading. Contributed by Ron Pruitt ( [19/Nov/91] [27/July/92] [30/Sept/92] (112764 bytes).
Box, Line, and Point Plot. A flexible one-dimensional distribution plot(Am.Stat. 51:353-358,1997). Shar file. (J. Jack Lee) [2/Oct/95][3/Jun/98] (53k)
Bayesian model averaging (BMA) to account for model uncertainty in linear regression models. Contributed by Jennifer Hoeting ( [3/Apr/95] [17/Sep/98] (27kbytes)
Replaced by bootstrap.funs
Replaced by bootstrap.funs
S language software (bootstrap, cross-validation, jackknife), data and errata for the book "An Introduction to the Bootstrap" by B. Efron and R. Tibshirani, 1993, Chapman and Hall. Please send comments or questions to R. Tibshirani ( [15/June/93][18/June/93][25/Aug/93][17/Mar/94][25/Jul/94] [1/Nov/94] (116 kbytes) There is also a precompiled DOS/Windows version
Boxplot Matrix. Extension of the standard oneway boxplot to cross-classified data with multiple observations per cell. A shar file which includes implementations in S and XLISPSTAT and the LaTeX paper. Richard M. Heiberger ( [19/Oct/92](138485 bytes)
A fortran program running through Splus which performs additive, robust nonparametric regression with automatic data transformation using a Bayesian Monte Carlo Markov Chain as the method of estimation. Based on Smith and Kohn, "Nonparametric regression using Bayesian variable selection", (to appear in the Journal of Econometrics.) Submitted by Mike Smith ( [28/Aug/95] [2/Oct/95] [25/Apr/96](305 Kbytes)
A shar archive of FORTRAN and Splus functions to compute shrinkage estimators with bridge penalty (\sum |\beta|^\gamma) for linear models. It uses shooting method for Lasso and modified Newton-Raphson method for general bridge. Submitted by Wenjiang FU ( [23/May/00] (20441bytes)
Box-percentile plots are similiar to boxplots, except box-percentile plots supply more information about the univariate distributions. ( [11/May/92] (14808 bytes)
The famous Monte-Carlo simulation is performed. It is a dumped file and it can be read into Splus (or R for that matter) by simply typing restore("buffon.dmp"). Submitted by Veith Tiemann ( [9/Nov/99] (3189 bytes)
A list of known bugs in S, compiled by Rick Becker. ( [10/Mar/90] (32022 bytes)
Function to compute proportions, means, and other statistics by levels of one or more classification variables. Part of the Harrell/hmisc collection.
S-PLUS functions for calibration and assay development using the four parameter logistic model. The functions are particularly useful for immunoassay (RIA) and enzyme linked immunosorbent assay (ELISA) applications. (Michael A. OConnell) [17/Nov/93] (111k)
An addon to the calibration functions of O'Connell, Belanger and Haaland which allows reuse of the same standard curve to create numerous calibration intervals. Submitted by Jim Robison-Cox ( [28/Feb/96] (14k)
A set of simple summaries of S syntax and expressions. Created by Barry Brown ( [5/Nov/90] (8361 bytes)
Library of functions for classification. Includes robust and predictive linear and quadratic discrimination, k-nearest neighbours and allied methods. Replaces classif, and uses only C. Part of the Modapplstat collection.
Replaced by class
S routines to write out data into common lisp structures, and common lisp routines to write data in a form that S can understand. A shar archive. ( [24/April/90] (22289 bytes)
These functions calculate the likelihood of a time series using the results of the Kalman Filter (distribution of the state given the past) and the results of the Reverse Kalman Filter (distribution of the state given the future data). Submitted by Ross Taplin ( [13/Jul/98] (1940 bytes)
A display for comparing one state (or individual) to many others. Display is described in an upcoming American Statistician. Submitted by Russell Almond ( [1/Dec/99] (20109 bytes)
Two S(plus) functions and help files to explore the use of the cone plot, a graphical object for exploring the geometric structure of high dimensional data. Brian Dawkins ( [16/Jan/92] (14689 bytes).
A collection of simple functions to generate and manipulate confounded and fractional factorial designs. Contr. by Bill Venables ( [8/Aug/94] (23 k)
translate S matrices into UNIX format sas, bmdp or blss data sets and vice versa. Also includes a function for writing an S matrix into source for xlisp. Contr. by Phil Spector ( [21/May/91] (41346 bytes)
S functions for "A note on the graphical analysis of survival data" (Cox, 1979). A shar archive. Submitted by George Tomlinson ( [22/Feb/99] (16768 bytes)
This is the shar file of GAM Cox prop hazards model, implemented in S/Splus (with most stuff in C/ratfor for efficiency). Splus driver comes with support of formula language prevalent in most other fitting software in S (lm, glm, gam etc). Rafal Kustra, rafal@utstat.utoronto. [23/Jan/95] (62k)
This is the shar file of Varying coefficient Cox prop hazards model, implemented in S/Splus (with most stuff in C/ratfor for efficiency). Splus driver comes with support of formula language prevalent in most other fitting software in S (lm, glm, gam etc). Rafal Kustra, rafal@utstat.utoronto. [23/Jan/95] (69k)
2 Functions: cpower, Compute power of two-sample Cox or logrank survival test; ciapower, Compute power of interaction test for exponential survival; Part of the Harrell/hmisc collection.
Software to perform subdistribution regression analysis with competing risks, using methods described in Fine and Gray (JASA, 1999). In shar format. Submitted by Bob Gray ( [24/May/99](44k)
Function to make a vector of character strings from a list of valid S names. Used most frequently for selecting subsets of variables from a data frame or from a SAS dataset. Part of the Harrell/hmisc collection.
Functions for construction of objects of class "ctree", which inherits from "tree". This group of functions isolates the data management capabilities of S "tree" functions. An application defining subsets of a cohort database is provided with functions to display the cohort progression plot. Contributed by Richard M. Heiberger ( [6/Oct/93]
An enhanced version of cut. Submitted by Frank Harrell Part of the Harrell/hmisc collection.
A shar file describing strategies for customizing S, with some functions and scripts used at the Columbia Univ. Psychology installation. Contributed by David H. Krantz ( [25/Oct/90] (116284 bytes)
compute the NPMLE of distribution function from doubly censored data. Contributed by Mai Zhou ( [1/Jun/95] (15kbytes)
A dyadic version of apply() and a generalized outer product function, used to construct distance matrices, cross product matrices, etc., under general conditions. Contributed by Colin Goodall ( [14/Apr/93] (23kbytes)
Functions that define and manipulate S "dataset" objects. A dataset is a matrix whose columns (variables) may be of different data types. Though motivated by a need to interface to SAS, they are useful in any data analysis. ( (Made obsolete by data.frame). [7/may/90][22/May/90] (63433 bytes)
Revised and new functions and methods for handling data frames in S. These allow new classes of objects to be included as variables, improve consistency and efficiency, and make better use of classes and methods. Paper (latex) and detailed documentation included. ( [9/Feb/94] (43k)
A shar archive containing the postscript version ( and the tex files to produce it for an AT&T Bell Labs reseearch report on data managment in S. ( [30/June/92] (260430 bytes)
A shar archive of S functions for dealing with dates. The most useful of them accepts a vector of input dates in any of the forms 8/30/53, 30Aug53, 30 August 1953,..., August 30 53, or any mixture of these. Submitted by Terry Therneau ( [2/April/91] [12/June/91][8/July/91][23/Apr/92] [17/Aug/92][14/Apr/93] (28 kbytes)
A library of functions and data has been written by Angelo Canty to accompany the forthcoming Cambridge University Press publication Bootstrap Methods and Their Application by Anthony Davison and David Hinkley. Compressed shar file. (See DH/ Contributed by Angelo Canty ( Mirrored periodically from EPFL Lausanne.
Nonlinear Calibration Design. A documented shar archive of Splus routines for simulating, displaying, evaluating, and comparing designs for nonlinear calibration experiments. Please send comments or questions to Doug Oman ( [1/Jun/95] (215 k)
S function deldir() to compute the Dirichlet (Voronoi) tesselation and Delaunay triangulation of a planar set of data points. (Rolf Turner) and (Don MacQueen) [28/Feb/96] [13/Dec/96] (87k)
A modification to the Aug '91 S function deriv() including an option to generate a hessian function suitable for input to ms(). Submitted by David Smith ( [13/Feb/92] (4694 bytes)
Biostatistical/epidemiologic modeling, testing, estimation, validation, graphics, and prediction by storing enhanced model design attributes in the fit. Mirrored periodically from Virginia Biostat . Submitted by Frank Harrell ( [31/May/94] [19/Sep/94] [2/Oct/95] [1/Dec/95][ 16/Dec/96] (939 kbytes)
PDF documentation for design (see above) Submitted by Frank Harrell ( [19/Sep/94] [13/Dec/96] (450 kbytes) & [1/Jun/95] (551 kbytes)
concise printed statistical description of a vector, matrix, or data frame. Especially useful with sas.get. Part of the Harrell/hmisc collection.
Function to construct approximately D-optimal experimental designs from candidate points. Only simple blocking is allowed but arbitrary treatment structures are accommodated. Bill Venables (, [1/Jun/95] (67k)
DOS program to convert simple S help files to text format for use in S-PLUS for Windows. NOT(!) available by e-mail. Use FTP, gopher, or WWW. (you might need to "cd" to the ../DOS/S area first). Submitted by dmurdoch@mast.QueensU.CA (Duncan Murdoch). [30/Aug/94] (10k)
Functions d/p/q poisson and binomial, ( [4/Nov/91] (4442 bytes)
(note the capital S ...) `D'etection `O'f `M'ultivariate `O'utliers `U'sing `S' `E'nvironment. Contains 13 functions useful in detection od multivariate outliers. ( (Normand Ranger) [6/Oct/93] (79k)
This is a history function for Version 2.0 Splus for DOS. It returns a menu of past commands that can be executed or edited then executed. Jack Turnock ( [30/Sept/92](3813 bytes)
Clustering algorithm based on dynamic altering of hierarchies. Suitable for large data sets. (See dysect in the general collection). (438 k).
An alternative to dyn.load for S on many BSD style UNIX's. Bill Dunlap ( [26/Feb/90] (28416 bytes)
Empirical cumulative distribution plot. (5k) Part of the Harrell/hmisc collection. Frank Harrell ( [17/Sept/93][8/Aug/94]
Routines for drawing ellipses and ellipse-like confidence regions, Version 1.1 (59k) dmurdoch@pickard.mast.QueensU.CA (Duncan Murdoch) [30/Aug/94] [28/Aug/95]
emu.v10 is an Splus object for fitting the limu model, a generalized mixed linear model for correlated responses on an arbitrary scale, with the (restricted) maximum likelihood techniques. The algorithm includes a maximization step and an updating step. The EM algorithm is used at the maximization step. References attached and comments welcome. Submitted by, [1/Dec/95] (35 kbytes)
ESS (Emacs Speaks Statistics) is updates gnuemacs4 and is the penultimate Emacs and XEmacs user interface for programming in and controlling S-PLUS and other statistical packages (including SAS and XLispStat). A link to the primary distribution site at the University of Wisconsin.
The function event.chart() is useful for plotting multiple time-to-event data at the individual level. Many extensions, such as sorting, subsetting, alignment, and multiple point and line formats, are provided. (Am.Stat. 54:60-73,2000) Submitted by J. Jack Lee ( [3/Jun/98][11/Apr/00] (73 Kbytes)
This article contains answers to some of the most frequently asked questions in Snews. Charles Roosen [12/Nov/92][23/Jan/95][28/Aug/95][28/Feb/96] (59k)
This is a TeXinfo version of the FAQ, which makes it useful as online documentation to be used with Emacs (and S/S-PLUS, as shown in the README file). Charles Roosen [23/Mar/93][23/Jan/95][28/Aug/95] (79k)
Superseded by mda.
Device driver for interfacing to the Fig/XFig public domain interactive graphics program. Shar archive. ( [26/Feb/91] (28740 bytes) (For Splus v.2.x)
Device driver for interfacing to the Fig/XFig public domain interactive graphics program. Shar archive. ( [24/Apr/92] [5/May/92] [24/Nov/92](17806 bytes) (For Splus v.3.x)
Find matches between columns of matrices. This function is useful for drawing matched samples from a reference sample using exact matching, matching within specified tolerances, and Euclidean distance matching with specified scale constants. Part of the Harrell/hmisc collection. (11 kbytes)
Numeric to character conversion with flexibility of C's printf. Submitted by Martin Maechler ( [19/Sep/94]
Maximum likelihood estimation of the parameters of a fractionally-differenced ARIMA (p,d,q) model (Haslett and Raftery, Applied Statistics 1989). Fortran source code in a shar archive. See also fracdiff in the general collection of StatLib. Submitted by Chris Fraley, ( [24/Mar/91]. [12/Aug/91] [9/Oct/91] [4/Nov/91] [13/Nov/91] [24/March/92] [15/Sept/92] [22/Feb/93] [11/Mar/94] [4/May/94] (135 kbytes)
A set of functions that implement fractional polynomial models in S-Plus. Submitted by Gareth Ambler ( [20/Oct/99] (29 kbytes)
A set of Splus function to estimate parameters of a positive stable frailty model for correlated time until event data. Described in: Recurrent Injury Event-Time Analysis. Wassell JT, Wojciechowski WC and DD Landen. Statistics in Medicine, 18(23), December, 1999, pp. 3355-3363. Submitted by James Wassell ( [14/Dec/99] (9445 bytes)
A shar archive containing a framemaker device driver for S. The contents are unix "diff"'s from the postscript driver in Spost.c. Submitted by David Lubinsky ( [22/Nov/90] (16475 bytes)
Introduction to S-Plus 3.1 in French. Written for students at the Swiss Federal Institute of Technology in Lausanne. Useful as an introductory document for beginners in S-Plus. PostScript file of a LaTeX document, 51 pages. Contributed by Marcel Baumgartner ( Contact the author for the LaTeX sources. [27/Jul/94] (988 kbytes)
FUNFITS is a comprehensive S-Plus module for fitting functions and nonlinear time series, including multivariate splines, Kriging and neural networks. Contributed by Doug Nychka ( [25/Apr/96] [24/Mar/97][24/Sep/99] (3 kbytes). The actual compressed tar file is available as funfits23.tar.gz in the S collection. Access this file via FTP, or the WWW, but not e-mail. (596k). Older version avaulable at funfits.tar.Z
gmanova version 2.2 is an Splus object for simultaneous inference of multiple response variables on arbitrary scales. Maximum likelihood is used for estimating a generalized multivariate linear model and results for each response variables are combined with either Fisher's method for combining p-values or measures of the Mahalanobis distances for two different purposes. Submitted by [1/Dec/95] (12 kbytes)
Nonlinear Discriminant Analysis via Scaling and ACE. A shar file. Written by Ross Ihaka ( [25/May/90] (47352 bytes)
A shar file containing S functions and a C engine for the Liang/Zeger generalized estimating equation approach to GLMs for dependent data. This revision includes support of model formulae. ( [5/April/90][5/July/90][27/June/91] [18/July/92] [8/Aug/92] [31/May/94][22/Oct/96][3/Jun/98] (218 kbytes) There is also a precompiled DOS/Windows version
gee.v26 is an Splus object for fitting the generalized linear regression model (GEEs) of Liang and Zeger. The syntax is similar to the glm(). This program was originally for comparing with the limu model (see emu.v10). With this program, the covariance matrix can be as large as 180x180, see Chapter 8 of Xie's dissertation for the analysis of Hummingbird learning data. Submitted by [1/Dec/95] (32 kbytes)
This S function implements the algorithm of Thompson and Taylor for the generation of pseudo data points drawn from the estimated density of an observed data set. Created by Barry Brown ( [30/Aug/89](8773 bytes)
Determines the number of iterations required to achieve a specified accuracy level in Markov chain Monte Carlo. An S translation of gibbsit in then general collection. Submitted by Steven Lewis ( [1/Mar/95] [4/Oct/95](7 kbytes)
A shar archive of glm code for new S. The core is a Newton minimizer. The modular structure allows easy user extension to "nonstandard" likelihoods. (A different version of glim). ( [6/Feb/91] (42753 bytes)
Bayesian generalized linear modeling, model comparison and accounting for model uncertainty. An S-PLUS function. Adrian Raftery, ( [17/May/93] [18/May/93] [13/Dec/96] [2/Jun/98] (33 Kbytes)
Fit one of the models described in McCullagh and Nelder, "Generalized Linear Models". Contributed by Terry Therneau ( [27/Sept/89][Updated 20/June/90,19/Nov/90,29/Jan/91] (22505 bytes)
Functions to allow design matrices to be specified in S using Glim or Genstat-type structure formulae. For example, design(a*b+x,length(a)), generates a (suitably labelled) design matrix for a double classification with covariate x. ([31/Mar/90] [16/Jan/92] (20370 bytes)
See ESS above.
A version of the gnuemacs tags program which has been modified to find S function definitions. You can use this to find function defintions or arbitrary strings in multiple related files. Submitted by David Lubinsky (hoqas!hoqas2! [4/Sept/90] (41607 bytes).
Graphical analysis of fractional factorial designs. A shar archive. Contributed by Perry Haaland ( and Michael O'Connell ( [10/Oct/90] (70041 bytes)
Hazard Regression: S functions to estimate the conditional log-hazard function based on possibly censored data and covariates and obtain corresponding densities, hazard rates, probabilities, quantiles and random samples. Includes proportional hazard model as special case. Charles Kooperberg ( shar archive. [26/April/93][17/Aug/93][1/Dec/95] [28/Feb/96] [24/Mar/97] [3/Jun/98] (193 kbytes)
S and C code from `Smoothing Techniques with Implementation in S' by Wolfgang Haerdle, Springer, 1991. A shar archive. Contributed by Brian Ripley ( [30/March/91][7/April/92] [8/Aug/92][21/Feb/94] (124 kbytes) There is also a precompiled DOS/Windows version
The Hmisc library contains many functions useful for data analysis, high-level graphics, utility operations, functions for computing sample size and power, translating SAS datasets into S-PLUS, imputing missing values, advanced table making, variable clustering, and character string manipulation. A shar archive. Mirrored periodically from Virginia Biostat . Contributed by Frank E Harrell Jr ( [13/Dec/96] (567k) See also the DOS version
Modifications of the high breakdown point estimators LMS, LTS, and MM for logistic regression, Poisson regression or Gamma regression with large groups. Submitted by Andreas Christmann ( [10/Sep/1998] (7k)
A shar archive of functions for estimating highest density regions in one and two dimensions. Contributed by Rob Hyndman ([29/Nov/94] [3/Apr/95](16 kbytes)
Hazard Estimation with Flexible Tails: S functions to estimate an unknown hazard function based on possibly censored data and obtain corresponding densities probabilities, quantiles and random samples. Charles Kooperberg ( archive. [26/April/93][17/Aug/93][9/Oct/93] [31/May/94] [1/Dec/95] [24/Mar/97] [3/Jun/98](141 kbytes)
Finite Difference Hessian Approximation from Function Values. An S implementatoin of Dennis and Schnabel's FDHESSF routine. Submitted by Barry Brown ( [30/May/90] (26302 bytes)
A shar file containing a S library to fit continuous-time hidden Markov models to irregularly spaced longitudinal observations. Contributed by Alexandre Bureau ( [24/May/00] (271k)
Function to compute Hoeffding's D statistic for testing independence of two variables, along with asymptotic P-values. Does computations for a matrix of variables, using pairwise deletion of NAs. Part of the Harrell/hmisc collection.
Functions for determining exact p-values and confidence bounds from stratified samples using the hypergeometric distribution. "John P. Wendell" ( [2/Oct/95] (18k)
Class and methods for imputation of missing values Works with the describe function (counts no. of imputed values) and the transcan function (develops imputation models predicting each predictor variable from all other predictors). Part of the Harrell/hmisc collection.
A shar archive containing New S code for using IMSL routines from S. Includes a patch for error handling and print routines, a Makefile, and examples of interface routines and S functions. Created by Ritei Shibata ( [6/May/89](21202 bytes)
A shar archive giving S functions for hazard smoothing with multiple covariates. Contributed by Robert Gray ( [16/Oct/97] (44k)
Simpler IMSL interface to new S than the previous "imsl". The new dynamic loading routine dyn.load2, contributed by Bill Dunlap, is necessary. Necessary subroutines will be directly searched for the IMSL library. A shar archive. ( ) [14/Dec/89] (19718 bytes)
Numerically integrate a function of a single variable. Includes S functions and Fortran and C source code. Created by Barry Brown ( (Made obsolete by integrate2). [Updated 6/Nov/90] (81709 bytes)
Single and multiple dimension numeric rectangular integration. (Change: Add good multiple dimension routine; improve packaging.) Includes S functions and Fortran and C source code. Shar file. (Barry W. Brown) [2/Oct/95] (190k)
S and Splus routines to compute projections, intersections, orthogonal complements, extensions, etc. for glm and lm model matrices. Brian Junker (, [23/Jan/95] (22k)
Compressed tar file with 3 functions for use with interval censored data, to fit nonparametric survival curves, plot them, and do logrank or Wilcoxon type tests. See also windows version. Submitted by Michael P. Fay (m7f@helix.nih.go). [6/Apr/99] (11k)
(Invert Monotone Function). Finds x such that f(x,...)=y for f montone in x and specified y. shar file (Barry W. Brown) [2/Oct/95] (37k)
Routines for output analysis of multiple sequences of iterative simulation (e.g., the Gibbs sampler or the Metropolis algorithm) ( [3/May/92][28/Feb/96] (9 kbytes)
an implementation of Duan and Li's (1991) inverse regression, handling multiple slices. Written entirely in S. Duan and Li, Slicing regression: A link free regression method. Ann. of Statist. (1991) 505-530. Mai Zhou ( [4/May/94][28/Aug/95] (2685 bytes)
Functions for Kruskal's isotonic multidimensional scaling and Shepard plots. Part of the Modapplstat collection.
A shar archive containing S functions to plot the coast-line of Japan and prefecture boundaries; similar to the usa() function. Contributed by Ritei Shibata ( [6/May/89](73845 bytes)
A pair of S functions, named filter and forecast, implementing Kalman filter recursions using the EM algorithm. (Gary Nan-Tie) [17/Nov/93] (9k)
The function kappa.cfint() will compute four different confidence intervals for the interrater agreement kappa in a 2x2 table. (J. Jack Lee), [17/Oct/94] (11k)
Two-dimensional kernel density estimation with help page. Submitted by [18/June/93] (10kbytes)
Two-dimensional gaussian kernel density estimation. Features: scaling, spherizing, automatic bandwidth, 1st derivatives. Submitted by, and Martin Maechler ([29/Nov/94] (16 kbytes)
S-PLUS/Fortran functions for implementing several of the procedures described in the book `Kernel Smoothing' by M.P. Wand and M.C. Jones, London: Chapman and Hall, 1995. Contributed by Matt Wand ( [28/Feb/96] [25/Jan/99] (139k)
Calculates Kent & O'Quigley's measure of dependence for censored data, when dependence between time and covariates is modelled by Cox's model. Function definion and help file. Contributed by Andrej Blejec ( and Janez Stare ( [13/Dec/96] (13k)
Functions for estimation of linear, generalized linear and the proportional hazards models via the "lasso" selection method. submitted by Rob Tibshirani ( [23/Jan/95] (501 kbytes)
A shar archive of C and S-PLUS functions to solve regression problems in which the L1 norm of the estimated parameters is constrained. Currently handles normal linear models and generalised linear models. It uses S-PLUS modelling functions. Submitted by Berwin Turlach ( [4/Nov/98][14/Dec/99] (954 kbytes)
Convert an S matrix to a LaTeX table, with major groupings, decimal alignment, multi-page tables, and control of number of digits to right of decimal point. See also `mat2tex'. Supersceded by print.display. Frank Harrell ( [6/Oct/92][17/Sept/93] (16kbytes)
Two S-PLUS functions to compute influence diagnostics for assessing case influence on the goodness-of-link test in logistic regression. Submitted by Andy H. Lee ( [28/Feb/96] (4k)
LINREG contains 9 functions useful in linear regression for diagnostics, residual plotting, Box and Cox transformations, an Anova table, bootstrap of a linear regression and a sequential ANOVA table. A shar archive. Submitted by Christian Leger ( [14/Sept/90] [23/Oct/91] (31897 bytes)
replaced by nlme
Modifications of the high breakdown point estimators LMS and LTS proposed by P.J. Rousseeuw for logistic regression and Poisson regression with large groups. Written by Andreas Christmann ( [2/Jun/98] (11k)
Superseded by design
S functions to, fully automaticly, estimate an unknown density based on possibly grouped or censored data and obtain corresponding probabilities, quantiles and random samples. C Kooperberg ( [5/Feb/92][6/Oct/92][26/April/93][17/Aug/93][6/Oct/93] [1/Dec/95] [24/Mar/97] (182 kbytes)
Replaced by logspline (above)
Functions for L-spline smoothing using a penalty term defined by an arbitrary linear differential operator L. Contributed by Jim Ramsay James Ramsay ( [14/Oct/97] (640k)
For both Windows and Unix ls() gives unquoted output that is void of any names beginning with a period. ls(<simplesubstring>) (optional quotes) works the same in Windows and Unix, giving all object names that contain the simple substring. John R. Wallace, [3/Apr/95] (3kbytes)
Logspline estimation of a possibly mixed spectral distribution: S functions to estimate the spectral density, line spectrum and spectral distribution based on time series data. C Kooperberg ( [26/July/93][17/Aug/93][9/Oct/93][1/Dec/95] [24/Mar/97][3/Jun/98] (95k kbytes)
See mar.builder in the General Archive
Functions and a database, to be used with the map library, for maps based on additional search keys, such as fips codes and state abbreviations. Contributed by Colin Goodall ( [14/Apr/93] (147 kbytes)
S-Plus program for a two-dimensional representation of a dendrogram. Submitted by Marcel Hendrickx ( [21/Jan/00] (12kbytes)
Model-based cluster and discriminant analysis, including hierarchical clustering and EM for parameterized Gaussian mixture models. Fortran code interfaced to Splus (UNIX and Windows). Mirrored daily from Submitted by Chris Fraley ( [25/Sept/91] [21/Oct/91] [4/Nov/91] [16/Jan/92] [13/Feb/92] [19/Oct/92] [22/Jan/93] [17/Nov/93] [24/Mar/97] [26/Aug/98]
Classification of binary Gaussian + Poisson mixtures in the plane. Fortran with Splus interface. Submitted by Chris Fraley ( [24/Mar/97] [14/Oct/97] (17k)
Markov chain Monte Carlo code for estimating item response theory (IRT) models, including 2- and 3-parameter logistic models, generalized partial credit models, and models with rater effects. A shar archive. Submitted by Richard J. Patz and Brian W. Junker. Please send comments or questions to R. Patz ( [8/Sep/99] (149k). See also the windows version.
Mixture Discriminant Analysis. Tools for fitting mixture models for classification, with subspace selection as in LDA. This archive supersedes fda, and includes revisions of all the functions there. Submitted by Trevor Hastie and Rob Tibshirani ( [22/Mar/96] (151k)
Menu system with generic menubar and popup menus; has methods for Emacs and can be extended to other systems (e.g., Tcl/Tk). A shar file. Submitted by Jeff Hallman (m1jjh00@FRB.GOV). [3/Jun/98] (36 Kbytes)
SPlus function to merge corresponding rows of two data frames sorted over the values of the set of like-named columns. Contributed by Raymond W. FALK ( [26/Oct/94] (11k)
A modification for S of David Gay's routine, mnfb, that finds parameter values to minimize a user supplied function. Created by Barry Brown ( [10/Oct/89](171388 bytes)
A package of Splus functions to compute the fit of a mixture of linear regressions and implement various associated inferential procedures. Submitted by Rolf Turner ( [15/Sep/98][2/Oct/98][20/Oct/98] (55k)
Robust regression estimator (MM-estimator). ("Plain Splus" without any extra C / Fortran code.) Andreas Ruckstuhl [28/Aug/95] (37 Kb)
Software, scripts, and data from the book Modern Applied Statistics with S-Plus by W.N. Venables & B.D. Ripley. There are unix archives and a Windows version. This material is not available via e-mail. Mirrored from Brian Ripley's collection at Oxford. ( Previous version is also available here.
A tabular summary of functions introduced in "Statistical Models in S". LaTeX document, 6 pages. ( [18/March/92] (10151 bytes)
Set up plots for testing model predictions vs. data. Frank E. Ritter (ritter@psychology.nottingham.AC.UK) [4/Jan/1994] (17k)
Performs a monotone smoothing on the vector x, optionally using weights in a parallel vector wt. Created by Barry Brown ( [1/Sept/89](15812 bytes)
Prints contents of motif() or openlook() window to colour postscript file. Contributed by Thomas Lumley ( [14/Oct/97] (7k)
A shar file of S-plus functions for robust regression model selection using stochastic complexity (minimum description length) criterion, robust AIC and robust BIC. Submitted by Guoqi Qian ( [31/Aug/98](33349 bytes).
Library of functions (S functions, Fortran source code) for hierarchical clustering (includes O(n^2) implementations, and 'see-it- as-it-happens' visualization program), partitioning, bond enery algorithm, Sammon mapping, PCA and correspondence analysis; latter includes fuzzy coding function. Contributed by Fionn Murtagh ( [29/Aug/94](205 kbytes)
A shar archive of functions to provide a model class for fitting log-linear multinomial models without using surrogate Poisson models. Also includes an enhanced nnet library (see Modapplstat). Uses C. Part of the Modapplstat collection.
Functions for fitting and analysing Negative Binomial GLMs. Bill Venables, (33k) [30/Aug/94]
New Environment for Statistical Inference. Various S functions for Classical Inference combined with graphical presentations of the results. A tree structured menu allows the user to interactively selection the method. This is an enhanced version (version 3.0) of the old "nesi". ( [1/Nov/89][10/Mar/90](177k bytes)
Splus-Tools for Model Selection and Construction of Bootstrap Confidence, Prediction and Calibration Intervals in Heteroscedastic Nonlinear Regression A tar file (so only available via FTP and WWW) Submitted by Joerg Polzehl ( [24/Mar/97] (266k)
A set of SPlus functions and programs to estimate the parameters of a non-linear regression model over a given set of observations. This file contains instructions for retrieving the code by ftp, gopher or WWW. Bouvier ( [29/Mar/94]
Functions, classes, and methods to fit and analyze linear and nonlinear mixed-effects models. A gzip'ed tar archive. Contributed by Jose' Pinheiro ( and Douglas Bates ([3/Jun/93] [6/Feb/94] [1/Mar/95] [2/Aug/96] [16/Oct/97] [25/Jun/99][2/Aug/99](1.8Mbytes) There is also a precompiled DOS/Windows version
A function to find features in point processes based on nearest neigbour distances. Simon Byers, ( [5/Jun/96][3/Jun/98] (5k)
Model class for feed-forward neural networks with fitting and prediction functions. Superseded by Modapplstat.
Nonlinear regression code. A corrected version of material from Appendix 3.1.1 of Bates and Watts, "Nonlinear Regression Analysis and Its Applications". A shar archive, Doug Bates ( Fixes by Peter Dalgaard ( [25/Oct/90] [22/Nov/90] (5054 bytes)
Function to do text searches in data and function objects in S. Versions for Unix and MS Windows. Contributed by Duncan Murdoch,dmurdoch@pickard.mast.QueensU.CA [23/Jan/95] (10 kbytes)
Oswald is an S-plus library for the analysis of Longitudinal Data. It includes informative dropout modelling and GEE analysis. A compressed tar archive. Contributed by David M Smith ( [13/Dec/96] (390k)
A shar archive containing S functions to plot the coast-line of Australia and state boundaries; similar to the usa() and jpn() functions. Contributed by Bill Venables ( [4/Aug/92] [1/Jun/95](66 kbytes)
This file contains the source codes for the S-Plus functions described in Levesque, A.C. and Eaves, D.M., Decision Modelling for Lot Acceptance, Chapter 9 in Advances in Plant Pathology (Springer-Verlag). pp 241-277. Submitted by David Eaves ( [8/Oct/98][29/Oct/98][30/Oct/98] (15k)
An implementation of the one-dimensional pool adjacent violators algorithm for monotone regression. Submitted by Anthony Gamst ( [10/May/00] (10k)
Splus 3.2 functions which implement the methodology for periodic autoregression modelling discussed in my article: A. Ian McLeod (1994) "Diagnostic Checking of Periodic Autoregression Models with Application", Journal of Time Series Analysis pp.221--233. Queries to [3/Apr/95](57kbytes)
Splus methods for visualizing non-linear functions (e.g., neural networks, Gaussian process models, and bagged trees) through plots of "additive effects". Tar format. Submitted by Tony Plate ( [4/Feb/99] (71kbytes)
Polytomous logistic regression: fit multinomial response models with possibly cross-classified outcomes, also describable as conditional loglinear models. Submitted by Alan Zaslavsky ([3/Jun/98][31/Jan/00] (215 kbytes)
Fortran, Matlab and Splus code for Univariate Partial Least Squares regression as described in Denham(1994), "Implementing partial least squares", (to appear in Statistics and Computing). Submitted by Mike Denham ( [85 kbytes]
printer-menu-window, a shar archive. A revised printer driver for new S. Emulates graphics windows and pop-up menus under Emacs on an ordinary 24 line ASCII terminal. Contributed by Richard M. Heiberger ( [25/May/89](64813 bytes)
A set of S functions for classification and polychotomous regression using adaptively selected polynomial splines for the model building. Charles Kooperberg ( [24/Mar/97][3/Jun/98] (170k)
A set of S functions for classification and (polychotomous) regression using adaptively selected polynomial splines for the model building. Contributed by Charles Kooperberg & Martin O'Connor, ( [15/Oct/97] (372k)
a small collection of S functions to implement a class for univariate polynomial manipulations. Bill Venables ( [28/Feb/96] (50k) There is also a precompiled DOS/Windows version
Functions to plot strings under postscript() with mixed fonts, character sizes, sub- and super-scripts; and display of available postscript fonts. Updated for Splus 3.3. ( [27/April/92] [7/April/93] [17/Sept/93] [28/Aug/95] [13/Dec/96] (48 kbytes)
A shar archive containing S functions for power and sample size analysis for tests of the general linear hypothesis. Contributed by Daniel F. Heitjan ( [15/May/91][27/Nov/91] (15501 bytes).
A shar archive containing S functions which will print any object as a list. Simple but useful for checking the organization of an object. Contributed by Ritei Shibata ( [6/May/89][Updated 23/Aug/91] (1449 bytes)
A shar archive containing S functions which fit a principal curve a matrix of points in arbitrary dimension Trevor Hastie ( [Aug/7/93] (15k)
A shar archive containing S functions which fit a hierarchical Poisson regression model. Poisson Regression Interactive Multilevel Modeling (PRIMM) Cindy L. Christiansen and Carl N. Morris Cindy Christiansen ( [22/Oct/96] (38k)
A set of object-oriented S functions that harness the automatic printing facility in S to convert an S object to another format, e.g. displaying an S object in its own window, converting a data.frame to a LaTeX table, for displaying an S help file in a LaTeX window. Richard M. Heiberger ( and Frank E. Harrell, Jr. ( [4/May/94] (130 kbytes)
S functions for doing projected scores (approx conditional inference) for GLM's with independent strata. A compressed tar archive. Contributed by Richard Waterman ( [24/Mar/97] (13k) See also the Windows version
A shar archive containing S functions which implement penalized spline smoothing. Contributed by Jim Ramsay ( [2/Aug/96] [13/Dec/96] (78k)
S functions for quadratic programming: see Goodall and Brown (1994), In: A Probability, Statistics and Optimization, A Tribute to Peter Whittle (Wiley). Contributed by Colin Goodall ( [1/Mar/95] (8 k)
A shar archive containing FORTRAN77 and S-PLUS functions for solving quadratic programming problems. Contributed by Berwin Turlach ( [14/Oct/97][6/Jul/98][23/Jul/98] (117k)
Routines to compute quantile regressions including analogues of trimmed means for the linear model. Updated version contains new functions implementing new forms of rank tests for linear models based on the dual quantile regression process. ([23/Aug/91] [1/Jun/95](48 kbytes)
A few small functions to find numerical rational approximations using an obvious continued fraction method. Created by Bill Venables ( [7/Aug/89](12607 bytes)
function to compute a Pearson or Spearman correlation matrix based on pairwise deletion of missing data Part of the Harrell/hmisc collection.
A provisional and largely untested interface between S and routines from "Numerical Recipes in C". The interface routine include methods for root-finding, optimization and so forth. A shar archive. ( [22/April/90] (90230 bytes)
A program that implements relative distribution methods (comparison density). Includes non-parametric estimates of the relative density, CDF and relative polarization indices. Submitted by Mark S. Handcock ( [23/Oct/98][7/Sep/99](47k)
Ian Abramson's S-program and manual for likelihood-free repeated measures/growth curve fits; handles imbalance, multiple fixed/floating effects, dependencies. (abramson@osiris.UCSD.EDU) [21/Dec/89] (30231 bytes)
This is a set of four fumctions in SHAR file that computes the (robust) Greatest Deviation Correlation Coefficient that is usefulon normal and erratic data. ( [31/Oct/91] (4032 bytes)
A Robust estimate of the intercept, slopes, and error scale in a simple or multiple linear regression with the Greatest Deviation Correlation Coefficient. Written by Rudy Gideon ( [5/May/92] (4813 bytes)
Functions and datasets used for teaching at Oxford. Used in sguide. A shar archive. Submitted by Brian Ripley ( [7/Aug/92] [17/May/93][21/Feb/94] [19/Oct/94] (71 kbytes)
S-plus functions for the graphical, parametric, non-parametric and likelihood analysis of repeated measures (i.e. longitudinal data). Contributed by Adelchi Azzalini ( and Monica Chiogna ( [26/Jan/94][28/Feb/96] (93k) See also the Windows version
Postcript source of user manual for (552k).[26/Jan/94][28/Feb/96]
Root mean square relative curvature measures for nonlinear regression. Bill Venables (, [9/July/93] [31/May/94] (6 kbytes)
Produces a matrix in which each row is an independent drawing from a specified multinomial distribution. Submitted by Jonathan Rougier ([3/Jun/98](6 kbytes)
Fortran code and S interface for the ROBETH family of robust statistics functions documented in `Algorithms, Routines and S Functions for Robust Statistics' (Wadsworth & Brooks/Cole, 1993) by A. Marazzi. This file contains instructions for retrieving the code by ftp, or WWW. Contributed by Alfio Marazzi ( [16/Mar/94] [24/Mar/97]
Collection of S functions for robust statistics plus a model class for robust linear regression via Huber's M-estimator with valid standard errors for estimates. Also contains some `dirty' datasets. Superseded by Modapplstat.
Row and column means, sums, variances, calls C for speed. Submitted by Tim Hesterberg ( [19/Sep/94] (10k)
A shar archive of functions for recursive partitioning - an alternative to tree(). Terry Therneau & Beth Atkinson ( [14/Oct/97] [3/Jun/98] (1.2M)
A shar archive of functions for recursive partitioning - an alternative to tree(). Update of rpart1, including case weights. Runs on Splus5.1. Terry Therneau & Beth Atkinson ( [8/Jun/00] (1.3M)
Splus function for reduced-rank regression and softly shrunk reduced-rank regression. Submitted by Magne Aldrin ( [19/Apr/99][8/Mar/00] (14k)
Simple interface from S to f77 or C programs which use Fortran I/O. No modifications to either S or IMSL are required! IMSL and some other libraries use Fortran I/O for errors. Partial replacement for dyn.load() when dyn.load is unavailable. ( [27/Nov/89] (7775 bytes)
Modified versions of bs() and ns() that allow safe predictions, especially in the context of the S modelling functions. New predict() methods as well. Trevor Hastie ( [11/Mar/94] (14k)
Calculates optimal sample sizes for a beta-binomial mixture model in a situation where both priors and utilities differ for two parties in a negotiation process, Based on the paper Wolfson, L.J., Kadane, J.B., and Small, M.J. (1995), "Expected Utility as a Policy-Making Tool: An Environmental Health Example", Lara Wolfson ( [1/Jun/95]
A shar file containing S-plus functions, C and fortran codes,which provide curvilinear principal components analysis (a particular nonlinear PCA) via a Smoothing Alternating with Least Squares Algorithm (SALSA). For more details, see Besse and Ferraty, 1995, Ferraty, 1997. Contributed by Frederic Ferraty, ( [14/Oct/97] [3/Jun/98] (170kbytes)
An S-plus program to calculate sample sizes based on highest posterior density intervals for normal means using three different Bayesian approaches. Submitted by Lawrence Joseph ( and Patrick Belisle. [7/Sep/95] (83 kbytes)
A S-plus program to calculate sample sizes based on highest posterior density intervals when comparing two binomial proportions using three different Bayesian approaches. Submitted by Lawrence Joseph (, Roxane du Berger and Patrick Belisle. [7/Sep/95](53kbytes)
an update to the S-Plus function for converting SAS datasets to S data frames or lists. Part of the Harrell/hmisc collection.
A short (and somewhat outdated) comparision of Sas and S. ( Terry Therneau, [20/Dec/93] (11 kbytes)
S functions for approximate Bayesian analyses. ([12/Aug/90] [21/Jan/92] Updated by Alix Gitelman ( [28/Feb/96] (191k)
One-dimension scattergram (generalized rug()). Part of the Harrell/hmisc collection.
A shar archive containting functions to make use of X or Suntools windowing when editing or getting help. The archive contains substitute versions of help and ed, and new functions tool.options and textedit. Contr. by ( [24/May/89](6453 bytes)
sguide.ps1 sguide.ps2
`Introductory Guide to S-Plus'. A beginners' guide to doing statistics in S-Plus. SGuide.ps1 and SGuide.ps2 are PostScript fullsize and reduced 2-on-1 respectively. Archive ripley is also needed. Submitted by Brian Ripley ( [7/Aug/92] [19/Oct/92] [21/Feb/94][19/Oct/94] (sguide.ps1 485 kbytes; .ps2 474kb)
The Nelder-Mead simplex algorithm for function minimization. The method uses no derivatives, is very robust in the face of local minima, but is slow. A shar file. Submitted by Terry Therneau ( [15/Nov/99] (14 kbytes)
S-Plus function for Sliced Inverse Regression based on the algorithm described by K.C. Li (1991, JASA 86, 316-327). Either the number of slices or the elements per slice can be controlled. Contributed by (Thomas Koetter) [28/Aug/95] (4710 bytes)
A lab course for teaching probability and statistics includes a postscript lab manual and some new functions which make S easier for students to use. Contributed by Doug Nychka ( [2/Mar/95] (3 kbytes). The actual compressed tar file is available as slab.tar.Z in the S collection. Access this file via gopher, FTP, or the WWW, but not e-mail.
The library `sm' is associated with the book `Applied Smoothing Techniques for Data Analysis: the Kernel Approach with S-Plus illustrations" by A.W.Bowman and A.Azzalini, 1997, Oxford University Press. A compressed tar file. This is version 2 of the library. Please send comments or questions to A.W.Bowman ( or A.Azzalini ( There is also a MSwindows version available here [2/Jun/98][18/Dec/00] (191k)
A set of S functions used to build a variety of Tukey median based smoothers and digital filters. Contributed by Colin Goodall ( [16/June/89](10664 bytes)
S interface to Fortran functions for universal kriging; K-fn and pseudo-likelihood analyses of spatial point patterns. Contains several spatial datasets. Superseded by Modapplstat.
A shar achive of spline code for new S. It requires that you be able to dynamically load new C code. Contributed by Doug Bates ( [1/May/89],[11/Jan/90]{Revised 24/Jan/90} (15679 bytes)
Instructions for obtaining the LaTeX (and postscript) source, and associated data, for a short course on Splus. The document talks mostly about plain S features and t does not concentrate on features specific to Splus. Very useful as an introductory document. Supersedes snotes. Created by Bill Venables ( [4/Aug/92][24/Mar/97]
pdf file containing a 429 page book on programming in (version 3 of) the S language. Also a shar file of accompanying functions, many of general utility. Submitted by Pat Burns
An interactive spreadsheet interface to S/Splus. Submitted by, (Richard M. Heiberger). [25/Apr/96] (171k)
The function sshc() calculates sample size for historical control trials using the randomized trial formula, the Makuch-Simon method, and the uniform power method. Submitted by J. Jack Lee ( [10/Nov/00] (8k)
Compactly print STRucture of any S object. S function 'str', help, utilities and examples. Martin Maechler ( [11/Mar/94] [29/Nov/94] [24/Mar/97] (35k)
Fits Stukel's generalization of the logistic link function to yes/no data. Shar archive. (Barry W. Brown) [2/Oct/95] (246k)
Sunflower plots, Martin Maechler ( [1/Jun/95] (14 k)
S and C code for Classification And Regression Trees analysis with survival data. A shar archive. Georg Haupt, Ulrich Mansmann, ( [1/Jun/95] (16 kbytes)
Obsolete, see survival2 below.
Functions for survival analysis including cox regression, Kaplan-Meier estimationm Time-dependent and/or multiple event Cox models. A shar archive. ( [17/June/91] (270786 bytes)
Functions for survival analysis. Encompasses all of survival2, + uses the "models and methods" style of 8/91 S, plus several new functions. ( [20/July/92] [7/Aug/92] [14/Apr/93] (446k)
Ancillary documentation for "survival3". A 30 page writeup of comments and examples (latex) plus some smaller notes on specific topics. ( [20/July/92] [7/Aug/92] [14/Apr/93] (153k)
Functions for survival analysis. A major update to survival3. ( [17/Oct/94][29/Nov/94] (596k) There is also a precompiled DOS/Windows version
Ancillary documentation for "survival4", including a writeup of comments and examples (latex) plus some smaller notes on specific topics. ( [17/Oct/94] [29/Nov/94] (227k)
Survival5, including some updates not yet a part of S-Plus (depending on your release). For Splus 3.4. Terry Therneau ([10/Jun/00] (1.7M)
Survival5, including some updates not yet a part of S-Plus (depending on your release). For Splus 5. Terry Therneau ([10/Jun/00] (1.9M)
Tools for converting S help files and S code to latex, and a corresponding latex .sty file. Submitted by Chambers and Hastie ( [23/Mar/93] (15k)
A library of S functions for estimating, converting and analysing multivariate ARMA and state space models. Especially useful for testing statistics and estimation techniques and for examining small sample properties. The most current version of the archive, user guide, and additional information are also available at and by anonymous ftp to in the directory pub/pgilbert. [22/Oct/96] (shar file 1.19M bytes) P.Gilbert ( See also for a Windows version.
A postscript file with a brief user's guide for the archive time.series.[May/95] (file 250K bytes)
A set of S functions and Fortran code for time-series analysis (modelled after the PC program of the same name). Includes documentation in TeX source. A shar archive. (jnewton@stat.tamu.end) [25/March/92] (164731 bytes).
S code to do one-dimensional data analytic wavelet threshold selection (must be used with "wavethresh" S-package by Nason and Silverman). Todd Ogden (, [23/Jan/95] (10k)
an S function that invokes unix() to use sed or tr to translate text strings. Part of the Harrell/hmisc collection.
Simultaneous transformation and imputation using canonical variates. Part of the Harrell/hmisc collection.
A shar file with S functions and test datasets for logistic regression analysis of data collected using two phase stratified sampling. Submitted by Breslow and Chatterjee ( [9/Aug/98] (26858 bytes)
A shar file containing S and C functions for a graphical display to facilitate selecting a transformation of the response in an ANOVA or regression model. Submitted by Pat Burns ( [6/Feb/92] (20771 bytes)
(Transform Plot Axes). Transform axes by a provided function as per built in log. Can thus do probit or logit plots, for example. Shar archive. (Barry W. Brown) [2/Oct/95] (11k)
A shar archive containing two functions to display S data objects in a tree-structured format. list.tree displays ordinary data objects (usually lists) and expr.tree displays expressions and functions. Contr. by ( [24/May/89][7/Aug/91](8270 bytes)
A shar archive of functions to replace prune.tree and predict.tree, correcting errors and enhancing functionality. Both C-based and pure S versions are included. Part of the Modapplstat collection.
An S-Function that fits the general linear growth curve model using the two-stage method when a linear combination of the growth curve coefficients is being modelled. Written by E. Demidenko and T.Stukel ( [3/Jun/98] (9216 bytes)
A shar file containing S class, methods, and help files, a C engine, and postscript documention for Tree Structured Survival Analysis (TSSA) as described in Segal M.R. (1988), Regression trees for censored data, Biometrics, 44:35--47. Contributed by Mark Segal ( and Carrie Wager ( Version 1.16. [22/Oct/96](680K)
This shar archive Cumulative Probabilities (ptukey) and Quantiles (qtukey) of the Studentized Range, or (more generally) for the maximum of m studentized range statistics. Contributed by John Maindonald ( [28/Aug/95] (36 Kb)
A menu-driven interactive tutorial program for Splus. It provides an introductory overview to many Splus functions and utilities. Contributed by Fu-Chuen Chang, ( [1/Aug/95] (29kbytes)
Exact tests for two independent binomial samples based on a wide range of non-classical Nullhypotheses and several families of test statistics. In data.dump format. Submitted by Ulrich Mansmann ( [19/May/99] (67kbytes)
Hierarchical clustering of variables using various similarity measures (Pearson r^2, Spearman rho^2, Hoeffding D). Part of the Harrell/hmisc collection.
List variables and selected attributes for a data frame or matrix. Submitted by Frank Harrell ( [23/July/92] (5kbytes)
This ascii file created by S contains 22 data sets from the book Visualizing Data published by Hobart Press ( and written by William S. Cleveland ( There is also a README object so there are 26 objects in all. The objects were written into the file by data.dump(), and should be restored by data.restore(). For example, if the dump is in the file xxx, the S command data.restore("xxx") reads them into S. The name of each data set is the name used in the book. To find the description of the data set, look under the entry "data, name" in the index. For example, one data set is barley. To find the description of barley, look in the index under the entry "data, barley". The datasets archive of Statlib also has the data in a table format, but, of course, it is far easier to use this dump to get the data into S. (671 kbytes) [12/Nov/93] [2/Oct/95] [22/Oct/97]
This is a collection of S-scripts for producing figures in W.S.Cleveland's "Visualizing Data" book. Included is a README file and 270+ scripts. [17/Oct/94][1/Dec/95] [22/Oct/97] (128k)
This is a collection of tables from W.S.Cleveland's "Visualizing Data" book. [22/Oct/97] (534k)
Enhanced vi and ed functions that allow multiple S objects to be edited and merged. Contributed by Colin Goodall ( [6/Oct/93] (11k)
S and Fortran 90 code for Nonparametric Regression Analysis of Longitudinal Data (Joan G. Staniswalis and J. Jack Lee, JASA, 93:1403-1418, 1998). A shar file. Submitted by J. Jack Lee ([13/Jan/99][14/May/99] (59k)
A shar archive containing S functions in order to perform a wavelet analysis of variance on univariate and bivariate time series. Submitted by Brandon Whitcher ( [13/Aug/98][5/Mar/99][27/Jul/99] (752k)
Code for doing wavelet transforms and thresholding in 1 and 2D in S. ([7/April/93] [Version 2.2, 7/July/93] (180kbytes)
Window/kernel/moving-average smoothing/estimation /nonparametric regression. Central values, quantiles, derivatives. Shar file. Contr. by Barry Brown ( [24/May/91] (153041 bytes)
S code to perform a weighted likelihood bootstrap. Michael Newton ( [17/Dec/92] (3k)
A simple S script for a world map. Submitted by Steve Wofsy, ( [17/Sept/91] (105019 bytes).
Double resolution world map data. Shar archive, Submitted by Steve Wofsy, ( [1/Oct/91] (238558 bytes).
A complete world political map. Much more detailed (and much larger) than either world or world.hires. Available only via FTP. Submitted by Becker and Wilks. [14/May/92]
A shar archive containing an x11 driver for new S. Contributed by Mike Meyer ( [5/May/89](56629 bytes)

Credit where credit is due

If you use an algorithm, dataset, or other information from StatLib, please acknowledge both StatLib and the original contributor of the material.
Last modified: Tue Dec 19 10:24:33 EST 2000 by Pantelis Vlachos