1、R Graphics,SCS Short Course March 14, 2008,Archeology,Archeological expedition Basic graphics easy and flexible Lattice (trellis) graphics powerful but less flexible Rgl nice 3d but challenging Tons of new stuff: e.g. install.packages(“heplots”) library(heplots) library(help=heplots) Special purpose

2、: library(car): avp In development: Plot3d.R,Goals for graphics,Two purposes for graphics: Exploratory need something quick that can produce graphs as fast as your mind can grasp them Presentation need flexibility to meet publication requirements,Why Lattice=Trellis graphics,Lattice (trellis) graphi

3、cs: developed at Bell Labs in early 80s Marginality vs Conditional association Motivation was problem of seeing relevant structures in higher dimensional data, the nub of the problem captured by Simpsons Paradox: two dimensional views of data show marginal association But conditional association can

4、 be entirely different Interaction: conditional association can be different for different values of moderator variables,Smoking and Life Expectancy,Source of demographic data (from the CIA): Two variables: - Life Expectancy - Cigarette Consumption Per Capita dl = read.csv(http:/www.math.yorku.ca/ge

5、orges/Data/CigLE.csv) head(dl),Data, head(dl) X Country Continent LE CigCon 1 1 Afghanistan Asia 35.5 98 2 2 Albania Europe 61.4 NA 3 3 Algeria Africa 60.6 1021 4 4 Andorra Europe 72.2 NA 5 5 Angola Africa 33.4 571 6 6 Antigua and Barbuda South America 61.9 NA tail(dl) X Country Continent LE CigCon

6、184 184 Vietnam Asia 61.3 NA 185 185 Yemen Asia 49.3 NA 186 186 Zambia Africa 34.9 408 187 187 Ethiopia Africa NA 87 188 188 Venezuela South America NA 1079 189 189 Zimbabwe Africa NA 399,Regression, table( LE = !is.na( dl$LE) , CigCon = !is.na(dl$CigCon) CigCon LE FALSE TRUE FALSE 0 3 TRUE 83 103 f

7、it = lm( LE CigCon , dl , na.action = na.omit) summary(fit) Call: lm(formula = LE CigCon, data = dl, na.action = na.omit) Residuals: Min 1Q Median 3Q Max -19.4576 -5.8225 0.8188 5.3636 17.8291 Coefficients: Estimate Std. Error t value Pr(|t|) (Intercept) 4.799e+01 1.371e+00 34.995 2e-16 * CigCon 8.5

8、28e-03 9.007e-04 9.468 1.33e-15 * - Signif. codes: 0 * 0.001 * 0.01 * 0.05 . 0.1 1 Residual standard error: 8.158 on 101 degrees of freedom (86 observations deleted due to missingness) Multiple R-Squared: 0.4702, Adjusted R-squared: 0.465 F-statistic: 89.64 on 1 and 101 DF, p-value: 1.333e-15,Interp

9、reting coefficient,Regression coefficient for CigCon is 0.008529 An increase in 1 cigarette per year increases Life Expectancy by 0.008529 years An increase in 1 cigarette per day increases Life Expectancy by 365 x 0.008529 = 3.11 years,All it takes is 3 cigarettes a day to add 10 years to your life

10、.,Why graphics!,Always look at your data The plot function: The Joys of Object-Oriented Programming: plot is a generic function so it does different things to different kinds of objects What does it do to a regression object? plot ( fit ),Basic regression diagnostics,So plot a fitted regression to g

11、et some diagnostic plots But diagnostics not enough to inspire understanding What happens if we plot a data frame? plot( dl ),Basic graphics functions in R,Axis Generic function to add an Axis to a Plot abline Add Straight Lines to a Plot arrows Add Arrows to a Plot assocplot Association Plots axTic

12、ks Compute Axis Tickmark Locations axis Add an Axis to a Plot axis.POSIXct Date and Date-time Plotting Functions barplot Bar Plots box Draw a Box around a Plot boxplot Box Plots bxp Draw Box Plots from Summaries cdplot Conditional Density Plots contour Display Contours coplot Conditioning Plots curv

13、e Draw Function Plots dotchart Cleveland Dot Plots filled.contour Level (Contour) Plots fourfoldplot Fourfold Plots frame Create / Start a New Plot Frame graphics-package The R Graphics Package grid Add Grid to a Plot hist Histograms hist.POSIXt Histogram of a Date or Date-Time Object,Basic graphics

14、 functions in R,identify Identify Points in a Scatter Plot image Display a Color Image layout Specifying Complex Plot Arrangements legend Add Legends to Plots lines Add Connected Line Segments to a Plot locator Graphical Input matplot Plot Columns of Matrices mosaicplot Mosaic Plots mtext Write Text

15、 into the Margins of a Plot pairs Scatterplot Matrices panel.smooth Simple Panel Plot par Set or Query Graphical Parameters persp Perspective Plots pie Pie Charts,Basic graphics functions in R,Generic function plot and its methods: plot Generic X-Y Plotting plot.data.frame Plot Method for Data Frame

16、s plot.default The Default Scatterplot Function plot.design Plot Univariate Effects of a Design or Model plot.factor Plotting Factor Variables plot.formula Formula Notation for Scatterplots plot.histogram Plot Histograms plot.table Plot Methods for table Objects plot.window Set up World Coordinates

17、for Graphics Window plot.xy Basic Internal Plot Function,Basic graphics functions in R,points Add Points to a Plot polygon Polygon Drawing rect Draw One or More Rectangles rug Add a Rug to a Plot screen Creating and Controlling Multiple Screens on a Single Device segments Add Line Segments to a Plot

18、 spineplot Spine Plots and Spinograms stars Star (Spider/Radar) Plots and Segment Diagrams stem Stem-and-Leaf Plots stripchart 1-D Scatter Plots strwidth Plotting Dimensions of Character Strings and Math Expressions sunflowerplot Produce a Sunflower Scatter Plot symbols Draw Symbols (Circles, Squares, Stars, Thermometers