R语言函数作图 - 开源中国社区
当前访客身份：游客 [
当前位置：
发布于 日 15时,
R语言函数作图
代码片段(1)
1.&[代码]函数作图&&&&
# Goal: To make a panel of pictures.
par(mfrow=c(3,2))
# 3 rows, 2 columns.
# Now the next 6 pictures will be placed on these 6 regions. :-)
# Let me take some pains on the 1st
plot(density(runif(100)), lwd=2)
text(x=0, y=0.2, "100 uniforms")
# Showing you how to place text at will
abline(h=0, v=0)
# All these statements effect the 1st plot.
x=seq(0.01,1,0.01)
par(col="blue")
# default colour to blue.
plot(x, sin(x), type="l")
lines(x, cos(x), type="l", col="red")
plot(x, exp(x), type="l", col="green")
lines(x, log(x), type="l", col="orange")
plot(x, tan(x), type="l", lwd=3, col="yellow")
plot(x, exp(-x), lwd=2)
lines(x, exp(x), col="green", lwd=3)
plot(x, sin(x*x), type="l")
lines(x, sin(1/x), col="pink")
开源中国-程序员在线工具：
相关的代码(25)
5回/11053阅
1回/5582阅
33回/3927阅
0回/4037阅
0回/3583阅
1回/3277阅
1回/2626阅
0回/2386阅
0回/2383阅
7回/2337阅
介绍下怎么和java，php等结合的方法吧，用在网页上应该不错
2楼：刘学炜 发表于
引用来自“jingshishengxu”的评论介绍下怎么和java，php等结合的方法吧，用在网页上应该不错估计他也没看懂这个程序的，只是代码分享而已。
开源从代码分享开始
DTC2的其它代码couldn't connect to host4780人阅读
R图形（2）
常用统计作图函数汇总
hist() 直方图
stem() 茎叶图
boxplot() 箱线图（盒形图）
coplot() 协同图
qqnorm() 正态qq图
qqplot() 两总体qq图
1. 高级低级图形函数的常用选项
高、低级图形函数概述
高级图形函数可以迅速简便地绘制常见类型的图形，但是，某些情况下你可能希望绘制一些有特殊要求的图形。比如，你希望坐标轴按照自己的设计绘制，在已有的图上增加另一组数据，在图中加入一行文本注释，绘出多个曲线代表的数据的标签，等等。
低级图形函数让你在已有的图的基础上进行添加。
低级图形函数一般需要指定位置信息，其中的坐标指的是所谓用户坐标，即前面的高级图形函数所建立的坐标系中的坐标。坐标可以用两个向量x和y给出，也可以由一个两列的矩阵给出。如果交互作图可以用下面介绍的locator()函数来交互地从图形中直接输入坐标位置。
最常用高级函数的是plot()函数。比如，
- plot(x,y)（其中x，y是向量）对两个变量画散点图。
- 用plot(z) （其中z是一个定义了x变量和y变量的列表，或者一个两列的矩阵）也可以达到同样目的。
如果x是一个时间序列对象（时间序列对象用ts()函数生成），plot(x)绘制时间序列曲线图。
- 如果x是一个普通向量，则绘制x的值对其下标的散点图。
- 如果x是复数向量则绘制虚部对实部的散点图。
- 如果f是一个因子，则plot(f)绘制f的条形图（每个因子水平的个数）。
- 如果f是因子，y是同长度的数值向量，则plot(f,y)对f的每一因子水平绘制y中相应数值的盒形图。
- 如果d是一个数值型数据框，则plot(d)对d的每两个变量之间作图（散点图等）
load("F:/R/cl.RData")
plot(cl$Sex)
高级绘图函数常用选项
低级绘图函数常用选项
低级图形函数一般需要指定位置信息，其中的坐标指的是所谓用户坐标，即前面的高级图形函数所建立的坐标系中的坐标。坐标可以用两个向量x和y给出，也可以由一个两列的矩阵给出。如果交互作图可以用下面介绍的locator()函数来交互地从图形中直接输入坐标位置。
图形参数详解
指定用于绘制散点的符号。绘制的点往往略高于或低于指定的坐标位置，仅pch=“.”无这个问题。
指定画线用的线型。缺省值lty=1是实线。从2开始是各种虚线。
指定线粗细，以标准线粗细为单位。这个参数影响数据曲线的线宽以及坐标轴的线宽。
指定颜色，可应用于绘点、线、文本、填充区域、图象。颜色值也可以用象”red”，”blue” 这样的颜色名指定。
用来指定字体的整数。一般font=1是正体，2是 黑体，3是 斜体，4是 黑斜体。
font.axisfont.labfont.mainfont.sub
分别用来指定坐标刻度、坐标轴标签、标题、小标题所用的字体。
指定文本相对于给定坐标的对齐方式。取0表示左对齐，取1表示右对齐，取0.5表示居中。此参数的值实际代表的是出现在给定坐标左边的文本的比例，所以adj=-0.1的效果是文本出现在给定坐标位置的右边并空出相当于文本10%长度的距离。
指定字符放大倍数。
坐标轴与坐标刻度
许多高级图形带有坐标轴，还可以先不画坐标轴然后用axis()单独加。
函数box() 用来画坐标区域四周的框线。
坐标轴包括三个部件：轴线（用lty可以控制线型），刻度线，刻度标签。它们可以用如下的图形参数来控制：
lab=c(5, 7, 12)
第一个数为x轴希望画几个刻度线，第二个数为y轴希望画几个刻度线，这两个数是建议性的；第三个数是坐标刻度标签的宽度为多少个字符，包括小数点，这个数太小会使刻度标签四舍五入成一样的值。
坐标刻度标签的方向。0表示总是平行于坐标轴，1表示总是水平，2表示总是垂直于坐标轴。
mgp=c(3,1,0)
坐标轴各部件的位置。第一个元素为坐标轴位置到坐标轴标签的距离，以文本行高为单位。第二个元素为坐标轴位置到坐标刻度标签的距离。第三个元素为坐标轴位置到实际画的坐标轴的距离，通常是0。
坐标轴刻度线长度，单位是绘图区域大小，值为占绘图区域的比例。tck小于0.5时x轴和y 轴的刻度线将统一到相同的长度。取1时即画格子线。取负值时刻度线画在绘图区域的外面。
xaxs=”s”yaxs=”d”
控制x轴和y轴的画轴方法。取值为”s”（即standard）或”e”（即extended）的时候数据范围控制在最小刻度和最大刻度之间。取”e”时如果有数据点十分靠近边缘轴的范围会略微扩大。这种画轴方式有时会在轴的一边留下太大的空白。取值为“i”（即internal）或“r”（此为缺省）使得刻度线都落在数据范围内部，而“r”方式所留的边空较小。取值设为”d”时会锁定此坐标轴，后续的图形都使用与它完全相同的坐标轴，这在要生成一系列可比较的图形的时候是有用的。要解除锁定需要把这个图形参数设为其它值。
一页多图用mfrow参数或mfcol参数规定，如：
par(mfrow=c(3,2))
表示同一页有三行两列共六个图，而且次序为按行填放。类似地，
par(mfcol=c(3,2))
规定相同的窗格结构，但是次序为按列填放，即先填满第一列的三个再填第二列。
par(mfrow=c(1,1))即可。
缺省时无外边空。为了规定外边空大小，可以用omi参数或oma参数。omi参数使用英寸为单位，oma参数以文本行高为单位，两个参数均为四个元素的向量，分别给出下、左、上、右方的边空大小。如：
par(oma=c(2,0,3,0))
函数mtext用来在外边空加文字标注。其用法为
mtext(text, side = 3, line = 0, outer = FALSE)
在多图环境中还可以用mfg参数来直接跳到某一个窗格，比如
par(mfg=c(2,2,3,2))
表示在三行两列的多图环境中直接跳到第二行第二列位置。mfg参数的后两个表示多图环境的行、列数，前两个表示要跳到的位置。
可以不使用多图环境而直接在页面中的任意位置产生一个窗格来绘图，参数为fig，如：
par(fig=c(4,9,1,4)/10)
此参数为一个向量，分别给出窗格的左、右、下、上边缘的位置，取值为占全页面的比例，比如上面的例子在页面的右下方开一个窗格作图。
S中一个单独的图由绘图区域（绘图的点、线等画在这个区域中）和包围绘图区域的边空组成，边空中可以包含坐标轴标签、坐标轴刻度标签、标题、小标题等，绘图区域一般被坐标轴包围。
边空的大小由mai参数或mar参数控制，它们都是四个元素的向量，分别规定下方、左方、上方、右方的边空大小，其中mai取值的单位是英寸，而mar的取值单位是文本行高度。例如：
par(mai=c(1, 0.5, 0.5, 0))
par(mar=c(4, 2, 2, 1))
这两个图形参数不是独立的，设定一个会影响另一个。S缺省的图形边空常常太大，以至于有时图形窗口较小时边空占了整个图形的很大一部分。
通常我们可以取消右边空，并且在不用标题时可以大大缩小上边空。例如下例可以生成十分紧凑的图形：
oldpar &- par(mar=c(2,2,1,0.2))
在一个页面上画多个图时边空自动减半，但我们往往还需要进一步减小边空才能使多个图有意义。
交互图形函数
函数locator(n, type)运行时会停下来等待用户在图中点击，然后返回图形中鼠标点击的位置的坐标。等待点击时用鼠标中键点击可以选择停止等待，立即返回。参数n指定点击多少次后自动停止，缺省为500次；参数type如果使用则可指定绘点类型，与plot()函数中的type 参数用法相同，在鼠标点击处绘点（线、垂线，等等）。locator()的返回值是一个列表，有两个变量（元素）x和y，分别保存点击位置的横坐标和纵坐标
例如，为了在已经绘制的曲线图中找一个空地方标上一行文本，只要使用如下程序：
text(locator(1), "Normal density", adj=0)
text()函数的adj参数用一个数字表示文本串相对于给定的坐标的画法，adj=0表示给定坐标为文本串左侧的坐标，adj=1表示给定坐标为文本串右侧的坐标，adj=0.5表示给定坐标为文本串中间的坐标。
函数identify(x, y, labels)在运行时也会停下来等待用户点击，直到按了鼠标中键，任何返回用户在图形中用鼠标点击的点的序号，点击时对点击的点加标签。参数x和y给出要识别的各个点的坐标。labels参数指定点击某个点时要在旁边绘制的文本标签，缺省时标出此点的序号，如果只需要返回值而不想画任何标记则可以在调用此函数时加一个plot=F参数。
注意identify()与locator()不同，locator()返回图中任意点击位置的坐标，而identify() 只返回离点击位置最近的点的序号。
例如，我们在向量x和y中有若干个点的坐标，运行如下程序：
attach(cl)
plot(Height, Weight) ;
identify(Height,Weight)
这时显示转移到图形窗口，进入等待状态，用户可以点击图中特别的点，该点的序号就会在旁边标出。为了结束，只要单击鼠标中键或单击右键并选择停止。返回结果为你点击的各个点的序号：
警告: 已经找到了最近的点
警告: 没有0.25英尺的点
2 10 15 16
2. hist直方图
hist(x)作向量x的直方图。缺省时自动确定分组，也可以用nclass=参数指定分组个数，或者用breaks=参数指定一个分组点向量。如果指定了prob=T则纵轴显示密度估计。
3. stem茎叶图
stem(x, scale = 1, width = 80, atom = 1e-08)
Arguments:
x: a numeric vector.
scale: This controls the plot length.
width: The desired width of plot.
atom: a tolerance.
References:
Becker, R. A., Chambers, J. M. and Wilks, A. R. (1988) _The New S
Language_. Wadsworth & Brooks/Cole.
stem(islands)
stem(log10(islands))
stem(cl$Weight, scale = 1, width = 80, atom = 1e-08)
The decimal point is 1 digit(s) to the right of the |
8 | 3445508
10 | 0332233
4. boxplot箱线图
箱线图的图形中的参数意义如下：
例子: 作了两个图叠加,一个为subset= supp == “VC“,另一个为subset= supp == ”OJ“.数据ToothGrowth有两个数量变量(len,dose)和一个属性变量(supp); 对每个supp,dose各有三个值(各形成三个盒形图). len ~ dose意味着box为len所做.
data(ToothGrowth)
boxplot(len ~ dose, data = ToothGrowth, boxwex = 0.25, at = 1:3 - 0.2, subset= supp == "VC", col="yellow", main="Guinea Pigs' Tooth Growth", xlab="Vitamin C dose mg", ylab="tooth length", ylim=c(0,35))
boxplot(len ~ dose, data = ToothGrowth, add = TRUE, boxwex = 0.25, at = 1:3 + 0.2, subset= supp == "OJ", col="orange")
legend(2, 9, c("Ascorbic acid", "Orange juice"), fill = c("yellow", "orange"))
boxplot(cl$Height~cl$Sex,data=cl,col = "lightgray");
5. coplot() 协同图
attach(cl)
coplot(Weight ~ Height | Sex)
如果z是一个数值型变量，则coplot() 先对z的取值分组，然后对z的每一组取值分别绘图。
甚至可以用如coplot(y~x | x1+x2)表示对x1和x2的每一水平组合绘图。
coplot()和pairs()函数缺省绘制散点图，但可以用一个panel= 参数指定其它的低级绘图函数，如lines，panel.smooth等。
tsplot(x)绘制时间序列曲线图。多个参数时tsplot(x1, x2, ...)表示绘制多条曲线，自动统一曲线取值范围。如果参数非时间序列对象则以下标1，2，3等为横坐标绘图。
6. qqnorm() 正态qq图与qqplot() 两总体qq图
- qqnorm(x), qqline(x), qqplot(x,y)作分位数－分位数图。
- qqnorm(x)对向量x作正态概率（纵轴为次序统计量值，横轴为对应该次序统计量的标准正态分布分位数值）。
- qqline(x) 除作qqnorm(x)图之外还画一条拟合曲线。
qqplot(x,y)把x和y的次序统计量分别画在x轴和y 轴以比较两个变量的分布。
x &- rnorm(50)
y &- rexp(60)
qqplot(x,y)#可以判断是不是同一个类型的分布
par(mfrow=c(2,2))#准备画2?2的4个图
plot(compression, distance,main= "Hooke's Law") #只有标题的图
plot(compression, distance,main= "Hooke's Law", xlab= "x",ylab= "y") #标题+x,y标记
identify(compression,distance) #标出点号码
plot(compression, distance,main="Hooke's Law") #只有标题的图
text(46,120, "f=1/2*k*s")#在指定位写入文字
plot(compression, distance,main="Hooke's Law") #只有标题的图
text(locator(2), "I am here!") #在点击的两个位置写入文字
Par(cex=0.7,mex=0.7)
plot(log(body),log(brain))
text(x=log(body), y=log(brain),labels=row.names(Animals), adj=1.5)
plot(log(body),log(brain))
identify(log(body),log(brain),row.names(Animals))
3 颜色画图函数
view.colours = function(){
plot(1, 1, xlim=c(0,14), ylim=c(0,3), type="n", axes=F, xlab="",ylab="")
text(1:6, rep(2.5,6), paste(1:6), col=palette()[1:6], cex=2.5)
text(10, 2.5, "Default palette", adj=0)
rainchars = c("R","O","Y","G","B","I","V")
text(1:7, rep(1.5,7), rainchars, col=rainbow(7), cex=2.5)
text(10, 1.5, "rainbow(7)", adj=0)
cmtxt = substring("cm.colors", 1:9,1:9)
text(1:9, rep(0.5,9), cmtxt, col=cm.colors(9), cex=3)
text(10, 0.5, "cm.colors(9)", adj=0)
view.colours()
高级低级函数画图
x &- seq(-pi, pi, len = 65)
plot(x, sin(x), type = "l", ylim = c(-1.2, 1.8), col = 3, lty = 2)
points(x, cos(x), pch = 3, col = 4)
lines(x, tan(x), type = "b", lty = 1, pch = 4, col = 6)
title("legend(..., lty = c(2, -1, 1), pch = c(-1,3,4), merge = TRUE)", cex.main = 1.1)
# 注意图例书写方式
legend(-1, 1.9, c("sin", "cos", "tan"), col = c(3,4,6), lty = c(2, -1, 1), pch = c(-1, 3, 4), merge = TRUE, bg='gray90')
matplot画图
matplot ：用矩阵的列画多线图
matplot(x, y, type = "p", lty = 1:5, lwd = 1, lend = par("lend"),
pch = NULL,
col = 1:6, cex = NULL, bg = NA,
xlab = NULL, ylab = NULL, xlim = NULL, ylim = NULL,
..., add = FALSE, verbose = getOption("verbose"))
require(grDevices)
matplot((-4:5)^2, main = "Quadratic")
sines &- outer(1:20, 1:4, function(x, y) sin(x / 20 * pi * y))
matplot(sines, pch = 1:4, type = "o", col = rainbow(ncol(sines)))
matplot(sines, type = "b", pch = 21:23, col = 2:5, bg = 2:5,
main = "matplot(...., pch = 21:23, bg = 2:5)")
5 时间序列图
=runif(100)
x=ts(x,start=1960+(3/12),frequency=12)#构造时间序列frequency=12即按个月份构造
options(digits=2)
0.444 0.849 0.459 0.659 0.277 0.268 0.958 0.926 0.595
1961 0.250 0.769 0.111 0.944 0.593 0.459 0.670 0.331 0.045 0.694 0.473 0.227
1962 0.037 0.181 0.562 0.455 0.419 0.968 0.288 0.187 0.143 0.548 0.437 0.280
1963 0.149 0.623 0.937 0.108 0.755 0.077 0.963 0.781 0.660 0.458 0.787 0.176
1964 0.997 0.650 0.045 0.849 0.018 0.744 0.148 0.964 0.844 0.660 0.080 0.303
1965 0.954 0.217 0.909 0.658 0.991 0.797 0.020 0.363 0.183 0.332 0.793 0.676
1966 0.459 0.270 0.813 0.818 0.138 0.780 0.679 0.029 0.553 0.895 0.548 0.480
1967 0.898 0.212 0.299 0.223 0.637 0.167 0.102 0.424 0.438 0.058 0.048 0.474
1968 0.254 0.381 0.532 0.089 0.825 0.337 0.382
plot(x,main="Time series")
参考知识库
* 以上用户言论只代表其个人观点，不代表CSDN网站的观点或立场
访问：70171次
积分：2737
积分：2737
排名：第10256名
原创：128篇
转载：81篇
评论：16条
(2)(8)(10)(1)(2)(6)(7)(13)(11)(31)(23)(12)(30)(36)(21)下次自动登录
现在的位置:
& 综合 & 正文
R语言的graphics功能（画图）
找了好久，终于找到了画各种统计图的画图工具，那就是R，本文主要介绍R中的各种图形绘制方法。R中有非常强大的绘图功能，本文也不能一一列举，如果想要具体了解，可在R中输入demo(graphic)或者demo(persp)查看。绘图函数需要将图输出到绘图设备上，绘图设备的形式有很多种，可以是窗口，可以是写入一个jpeg、bmp、tiff等图片文件格式，具体可用help(device)查看。
打开多个绘图设备
前面已经提到，绘图设备有多重，如果是R中的窗口，可以用函数x11()来打开窗口，虽然x11是在linux下的叫法，但是在windows中也可以用该函数，当然也可以用windows打开；另外如果要打开不同的图形文件设备，可以用jpeg()/bmp()、tiff()等函数打开不同类型的文件设备。dev.list()可返还当前打开的所有的图形设备。dev.cur()返回当前的图形设备，还有dev.set(number)用于激活相应的图形设备,dev.next(),dev.prev()用于切换。dev.off()用于关闭当前设备，这个函数可以指定设备号来关闭，比如dev.off(1);用graphic.off()可以关闭所有的设备。
2. Graphic的分割
下面介绍的几个函数可以实现graphic的分割功能，但是这些都不是兼容的（即不能同时使用），
用x11()等函数打开图形设备之后，就可以使用相应的screen函数了。
split.screen(c(2,2))该函数用于将screen分成2行2列，标识符按行分别从1到2*2标记，参数向量用于指定每个screen的标识符，注意不能分割太多，会造成图画不下的情况。
split.screen(c(1,2), screen=2)可以将子screen再次进行分割screen这个参数用于指定需要再次分割的标识符。
加下去就是画图了
screen(2)用于选中标识符为2的子屏幕，以后的画图函数都在2这个子屏幕中生效，如果2这个子屏幕中已经被画过了，它会清空已经画过的。如果不想清空 只需 screen(2, FALSE)即可
erase.screen(n=) 用于清空指定标识符上的图形
close.screen(n, all.screens =FALSE) 关闭屏幕，相当于清空了标识符，因此再也不能定位到相应的屏幕了
相比screen函数，还有一个Layout，比screen更加灵活（理论上可以设置出任意的布局），也更加方便。
layout(mat, widths, heights):mat用于划分整个布局，如果mat中两块数字相同则布局显示时候会将那两块显示在一起。Widths和heights都是一个向量，长度分别和列数和行数相同，用来表示每一块列或者行的相对长度。如果要用绝对的，可用lcm(5)这个函数，返回5cm
layout.show(n)用来显示子窗口的个数。
在画图的时候，Layout中默认是bycolonm的，即画图时候，先画a[1,1].再a[1,2],再a[2,1]a[2,2]，如果要a[1,1]之后是a[1,2]只要在matrix这个参数中指定byrow=TRUE即可。
3. 绘图函数
高级绘图函数：
& mat &- matrix(c(1:4), 2, 2)
& layout(mat)
& plot(1:10)
& text(5, 5, labels = "plot(1:10)")
& plot(1:10, cos(1:10))
& text(5, 5, labels = "plot(1:10, cos(1:10))")
& sunflowerplot(1:10, 1:10)
& pie(c(0.1, 0.2, 0.3, 0.4))
& layout(mat)
& boxplot(1:10)
& stripchart(1:10)
& x = 1:10
& y = 3:12
& z = x + y
& coplot(x~y|z)
好多函数啊！！！脖子酸了！！就不一一列举了，因为暂时还用不到，以后等用到了再贴进来吧！！
低级绘图函数：
4. 绘图参数
绘图参数可以作为图形函数的参数来指定相应的属性，这里来介绍一个另外的方法，par();
Par这个函数有如下参数，可以修改图形中的不同的属性：
绘图参数和低级作图函数使我们可以进一步改善图形。前面我们已经看到，一些绘图参数不允许作为plot这样的函数的自变量。我们可以用par()来修改这些参数，这样就必须输入多行的命令。在改变绘图参数时，预先保存它们的初始值以便以后恢复十分有用，如下：
opar &- par()
par(bg="lightyellow", col.axis="blue", mar=c(4, 4, 2.5, 0.25))
plot(x, y, xlab="Ten random values", ylab="Ten other values",
xlim=c(-2, 2), ylim=c(-2, 2), pch=22, col="red", bg="yellow",
bty="l", tcl=-.25, las=1, cex=1.5)
title("How to customize a plot with R (bis)", font.main=3, adj=1)
R的graphics就介绍到这里，待实际使用的时候再慢慢去领悟。
&&&&推荐文章:
【上篇】【下篇】[R语言绘图]气泡图symbols
绘制气泡图主要使用函数symbols(x,y,circle=r).其中x、y是坐标轴，r是每个点的半径。
x&-rnorm(10)
y&-rnorm(10)
r&-abs(rnorm(10))
symbols(x,y,circle = r,
bg=rainbow(10))
此外，我们还可以给每个气泡加上文字。
attach(mtcars)
r&-sqrt(disp/pi)
symbols(wt,mpg,circle=r,
inches=0.3,
bg=lightblue)
text(wt,mpg,row.names(mtcars),
(window.slotbydup=window.slotbydup || []).push({
id: '2467140',
container: s,
size: '1000,90',
display: 'inlay-fix'
(window.slotbydup=window.slotbydup || []).push({
id: '2467141',
container: s,
size: '1000,90',
display: 'inlay-fix'
(window.slotbydup=window.slotbydup || []).push({
id: '2467142',
container: s,
size: '1000,90',
display: 'inlay-fix'
(window.slotbydup=window.slotbydup || []).push({
id: '2467143',
container: s,
size: '1000,90',
display: 'inlay-fix'
(window.slotbydup=window.slotbydup || []).push({
id: '2467148',
container: s,
size: '1000,90',
display: 'inlay-fix'