R语言初步统计绘图与编程-----ggg要点.ppt

* 空间分析绘图 图22 gstat程序包实例 * R绘图功能 R具备卓越的绘图功能，通过参数设置对图形进行精确控制。绘制的图形能满足出版印刷的要求，可以输出Jpg、tiff、eps、emf、pdf、png等各种格式。 通过与GhostScript软件的结合，可以生成600dpi，1200dpi的等各种分辨率和尺寸的图形。 绘图是通过绘图函数结合相应的选项完成的。 绘图函数包括高级绘图函数和低级绘图函数。 * 高级绘图函数 plot() 绘制散点图等多种图形,根据数据的类，调用相应的函数绘图 hist() 频率直方图 boxplot() 箱线图 stripchart() 点图 barplot() 柱状图 dotplot() 点图 piechart() 饼图 matplot() 数学图形 * 图23 散点图与箱线图 * 低级绘图函数 lines() 添加线 curve() 添加曲线 abline() 添加给定斜率的线 points() 添加点 segments() 折线 arrows() 箭头 axis() 坐标轴 box() 外框 title() 标题 text() 文字 mtext() 图边文字 * 绘图参数 参数用在函数内部，在没有设定值时使用默认值。 font = 字体,lty = 线类型,lwd = 线宽度,pch = 点的类型, xlab = 横坐标,ylab = 纵坐标, xlim = 横坐标范围, ylim = 纵坐标范围, 也可以对整个要绘制图形的各种参数进行设定 参见 par() * 举例：绘图 生成0到2之间的50个随机数，分别命名为x,y x - runif(50,0,2) y - runif(50,0,2) 绘图：将主标题命名为“散点图”， 横轴命名为”横坐标”， 纵轴命名为“纵坐标” plot(x, y, main=散点图, xlab=横坐标, ylab=纵坐标) text(0.6,0.6,text at (0.6,0.6)) abline(h=.6,v=.6) * 图24 绘图举例 * 例：分步绘图 1.打开绘图窗口，不绘制任何对象 plot(x, y, type=n, xlab=, ylab=, axes=F) 2.添加坐标点 points(x,y) 3.添加坐标轴 axis(1); axis(at=seq(0.2,1.8,0.2), side=2) 4.补齐散点图的边框 box() 5.添加标题、副标题、横轴说明、纵轴说明 title(main=Main title, sub=subtitle, xlab=x-label, ylab=y-label) * 图25 分步绘图 1 2 3 4 5 6 * 一页多图 图26 一页多图 par() par(mfrow=c(2,2) ... * 在原有图形上添加元素 举例： x - rnorm(100) # 生成随机数 hist(x,freq=F) # 绘制直方图 curve(dnorm(x),add=T) # 添加曲线 h - hist(x, plot=F) # 绘制直方图 ylim - range(0, h$density, dnorm(0)) #设定纵轴的取值范围 hist(x, freq=F, ylim=ylim) #绘制直方图 curve(dnorm(x),add=T,col=red) #添加曲线 * 图27 在原有直方图上添加曲线 * 练习八：绘图练习 练习1 将Rplots.r中的代码拷贝到R控制台中，查看R绘制的图形。 练习2 对例进行回归分析，并绘制散点图，并为散点图添加回归直线。 plot(impurity~rate) reg-lm(impurity~rate) abline(reg,col=red) summary(reg) * 五 编写函数 * 编程基础 R可以灵活的编写程序，用户自己编写的程序可以直接调用。编程时无需声明变量的类型，这与C,C++等语言不同。 基本格式 函数名 - function(数据,参数1= 默认值,…) { 异常处理; 表达式(循环/判别); return(返回值); } 函数内部也可用#添加注释 * 函数实例 图28 函数实例 data2mat() * 程序流程控制 if if(条件) 表达式 if(条件) 表达式1 else 表达式2 举例 p = 0.03 { if(p=0.05) print(p = 0.05!) else pr