统计软件和r语言.ppt

demo if (dev.cur() = 1) get(getOption(device))() opar - par(ask = interactive() (.Device %in% c(X11, GTK, gnome, windows, Macintosh))) x - rnorm(50) opar - c(opar, par(bg = white)) plot(x, ann = FALSE, type = “n”)#出现“Hit Return to see next plot:” `ann If set to `FALSE, high-level plotting functions do not annotate the plots they produce with axis and overall titles. The default is to do annotation. demo abline(h = 0, col = gray(0.9)) lines(x, col = green4, lty = dotted) points(x, bg = limegreen, pch = 21) title(main = Simple Use of Color In a Plot, xlab = Just a Whisper of a Label, col.main = blue, col.lab = gray(0.8), cex.main = 1.2, cex.lab = 1, font.main = 4, font.lab = 3) demo par(bg = gray) pie(rep(1, 24), col = rainbow(24), radius = 0.9)#Hit Return to see next plot: title(main = A Sample Color Wheel, cex.main = 1.4, font.main = 3) title(xlab = (Use this as a test of monitor linearity), cex.lab = 0.8, font.lab = 3) pie.sales - c(0.12, 0.3, 0.26, 0.16, 0.04, 0.12) names(pie.sales) - c(Blueberry, Cherry, Apple, Boston Cream, Other, Vanilla Cream) pie(pie.sales, col = c(purple, violetred1, green3, cornsilk, cyan, white)) demo pie.sales - c(0.12, 0.3, 0.26, 0.16, 0.04, 0.12) names(pie.sales) - c(Blueberry, Cherry, Apple, Boston Cream, Other, Vanilla Cream) pie(pie.sales, col = c(purple, violetred1, green3, cornsilk, cyan, white))#Hit Return to see next plot: title(main = January Pie Sales, cex.main = 1.8, font.main = 1) title(xlab = (Dont try this at home kids), cex.lab = 0.8, font.lab = 3) demo par(bg = cornsilk) n - 10 g - gl(n, 100, n * 100) x - rnorm(n * 100) + sqrt(codes(g)) boxplot(split(x, g), col = lavender, notch = TRUE)#Hit Return to see next plot: title(main = Notched Boxplots, xlab = Group, font.main = 4, font.lab = 1) demo par(bg = white) n - 100 x - c(0, cumsum(rnorm(n))) y - c(0, cumsum(rnorm(n))) xx - c(0:n, n:0) yy - c(x, rev(y)) plot(xx, yy, type = n, xlab = Time, ylab = Distance)#Hit Return to see next plot: polygon(xx, yy, col = gray) t