OverviewTo work efficiently with Splus, you need to have a window open in which you can save the commands (those which worked) and the output.
- Use Emacs Speaks Statistics, ESS, with one window for Splus and another for saved input and output. Start emacs
> emacs &Then open a new window (File| Make New Frame or keys Cntrl-X 5b) Now start Splus with commands: Meta-X S. Return when it asks which directory to start from, or you can specify any directory in which you've created a MySwork directory to hold Splus objects. Default is your root directory where it will build the MySwork directory if none exists.
If the above commands don't work, you need to copy a specially tailored .emacs file to your root directory.
gauss> cp ~jimrc/.emacs .emacsFor more details on running within emacs see Emacs Speaks Statistics, ESS. In both cases, I assume you're running an Xsession so that you can mark text by dragging across it with the left mouse button depressed, and paste it by clicking the center mouse button. You can open a graphics window with the (Splus) command:
After you create a plot, you should be able to print it by clicking on the Graph menu, then the Print choice. If you have any problem, check the Options: Printing choices. You should see `lp' in the command box. If not, change it, click Save and Close, then try Graph: Print again.
Sample SessionStartup Splus uses functions for everything, even quit is a function, q(), with no argument in the ()'s. It handles vectors and matrices nicely. It uses <- to assign a value to a variable, but within a function call, = is used to set parameters.
> tires <- read.table("tire.data",col.names=c("hardness","tensile.strength", "abrasion.loss")) > tiresThis reads the file and prints out the hardness dataset. It has 3 variables, tires$hardness, tires$tensile.strength, and tires$abrasion.loss. The read.table function brought it in as a data frame. Now to see how it looks, we need to open a graphics window and do a pairs plot.
> motif() > pairs(hardness)The motif() command only is done once per session. Other choices for graphics devices are available, eg trellis.device(motif) or postscript("filename") to create a postscript file. The graphics.off() command closes the motif window, or completes the postscript file. Next we'll try fitting a regression model using the lm() (linear models) function. The tilde below means ``is modeled using'', then a ``sum'' of explanatory variables is given.
> tires.fit <- lm(abrasion.loss ~ hardness + tensile.strength, data=tires) > summary(tires.fit)We'll finish with a plot of residuals versus fitted values, add a reference line at 0. If you want a paper copy of the plot, click the graph menu option above the picture and the print option below that.
> plot(fitted(tires.fit),tires.fit$resid) > abline(h=0)Exiting After you've looked at the plot, quit q().
- Trellis Graphics Homepage
- For more details on running within emacs see Emacs Speaks Statistics, ESS.
- Splus is quite slow at looping, so simulations are slow.
- Splus homepage.
- Questions and Answers from the S-news list (old).
Other Stat package demos for unix packages:
Back to MSU Math Department Home Page
Last Updated: Tuesday, 14-Jun-2011 16:52:23 MDT