R语言学习笔记——BaseGraphics

OStack程序员社区-中国程序员成长平台 › 门户 › 编程› R语言›R语言教程

原作者: [db:作者] 来自: [db:来源] 收藏邀请

做exploratory data annalysis的作业，差点被虐死了，R从头开始，边做边学，最后搞到一点多才弄完，还有一个图怎么画都不对，最后发现是数据读取的时候有问题。

用来画图的数据来自：http://archive.ics.uci.edu/ml/datasets/Individual+household+electric+power+consumption

数据属性如下：

Attribute Information:

1.date: Date in format dd/mm/yyyy 
2.time: time in format hh:mm:ss 
3.global_active_power: household global minute-averaged active power (in kilowatt) 
4.global_reactive_power: household global minute-averaged reactive power (in kilowatt) 
5.voltage: minute-averaged voltage (in volt) 
6.global_intensity: household global minute-averaged current intensity (in ampere) 
7.sub_metering_1: energy sub-metering No. 1 (in watt-hour of active energy). It corresponds to the kitchen, containing mainly a dishwasher, an oven and a microwave (hot plates are not electric but gas powered). 
8.sub_metering_2: energy sub-metering No. 2 (in watt-hour of active energy). It corresponds to the laundry room, containing a washing-machine, a tumble-drier, a refrigerator and a light. 
9.sub_metering_3: energy sub-metering No. 3 (in watt-hour of active energy). It corresponds to an electric water-heater and an air-conditioner.

1.读取数据：

在R中，要从txt读取数据，需要使用read.table函数，它的几个常用的参数列表如下：

参数	含义
file	文件的名字，可以用文件的绝对路径
header	逻辑值，为1表示文件第一行包含各个变量的名称，如果不包含，可以用colnames参数指定
sep	各个变量之间的分隔符，可能是空格，分号，逗号等等，根据文件自定义就可以了
col.names	是一个向量，指定每一列的名称
colClasses	是一个类向量，指定每一列的类别
na.strings	表明数据中什么样的值被定义为NA

本例中用如下代码读入数据：

x<-read.table("D:/coursera/data analysis/household_power_consumption.txt",sep=";", header=T, colClasses = c('character', 'character', 'numeric',
                                    'numeric', 'numeric', 'numeric','numeric', 'numeric', 'numeric'),na.strings='?')

2. 合并年月日和时分秒，并把它们转换成Date类型存放到新增加的一列DateTime中

上述代码中是将年月日和时分秒按照character类型读入的，接下来的运算中使用Date类型会比较方便。

主要有两个步骤：1.利用paste函数拼接年月日和时分秒；2.将拼接后的串转换成Date型。

代码如下:

x$DateTime <- strptime(paste(x$Date, x$Time),"%d/%m/%Y %H:%M:%S")

这里，strptime把粘贴Date和Time得到的串转换成Date格式，其中第二个参数format非常重要，原来的数据形如16/12/2006 17:24:00，所以format中的"%d/%m/%Y %H:%M:%S"要和数据一一对应。比如这里用了Y而不是y，原因在于Y表示带世纪的年份，而y只能表示不带实际的年份，它的取值范围只有00~99。还有数据中用于分隔的“/”和":"都不能省略，一旦格式和数据不同，就会发现转换后的数据都变成NA。strptime中format常用的几个参数列出如下：

%d	day（01~31）
%m	month（01~12）
%Y	Year with century(0~9999)
%y	Year without century(00~99)
%H	hour（00~23）
%M	minutes（00~59）
%S	second(00~61)

(ps.我也不知道为什么%S可以到61，帮助文档里面是这么写的=。=)

转换后DateTime值形如：2006-12-16 17:24:00

3.抽取日期范围在“2007-2-1”和“2007-2-2”之间（包括边界）的数据放入s中：

s <- subset(x,as.Date(DateTime) >= as.Date("2007-02-01")&as.Date(DateTime) <= as.Date("2007-02-02"))

这里主要用到as.Date()函数，它把一个串转换为Date格式。

4.直方图

这里主要以global_active_power频度直方图为例，代码如下：

hist(s$Global_active_power,freq=TRUE,col="red",xlab="Global Active Power(kilowatts)",ylab="Frequency",xaxt="n",yaxt="n",main="Global Active Power")
axis(side=2,at=seq(0,1200,200))
axis(side=1,at=NULL)

hist函数有多个参数，可以在帮助文档中查看，常用的几个列举如下：

pch	图例样式，默认为空心的小圆圈
col	颜色，用整数定义，可以用colors()函数查看所有的颜色（查了一下，居然有657种=。=）
xlab	x轴标签
ylab	y轴标签
main	图像名称

上述的参数中还设置了xaxt="n",yaxt="n",是为了之后用axis(side=2,at=seq(0,1200,200))和axis(side=1,at=NULL)分别标注x和y轴的刻度。其中的seq(0,1200,200)表示产生一个0~1200，以200为公差的等差序列。

5.把图像存为png格式

dev.copy(png,filename="plot1.png",height=480, width=480,bg="white")
dev.off()

这里会把图像存放在当前路径下，可以用getwd()查看当前路径。

生成的图像如下：

6. 折线图

这里以DateTime-global_active_power折线图为例，有两种画图方式

5.1 直接用plot：

plot(s$DateTime, s$Global_active_power, xlab="n", ylab="Global Active Power(kilowatt)", type="l",lty=1)

5.2 用plot和lines函数：

plot(s$DateTime, s$Global_active_power, xlab="n", ylab="Global Active Power(kilowatt)", pch=NA)
lines(x$DateTime, x$Global_active_power)

先用pch=NA画出一张空白图，然后用lines在上面增加折线。

第二种方法的时间会比第一种的慢。

画出的图形如下：

7. 在一张图中画出多条折线并增加图例

上述的第二种方法只要多加几个lines()函数就可以实现这一功能了。这里以DateTime-sub_metering_*(*=1,2,3)为例：

plot(s$DateTime, s$Sub_metering_1, yaxt="n", ylab="Energy Sub metering",type="l")
lines(s$DateTime, s$Sub_metering_2, col="red")
lines(s$DateTime, s$Sub_metering_3, col="blue")

增加图例用legend()函数：

legend("topright",legend=c("sub_metering_1","sub_metering_2","sub_metering_3"),col=c("black","red","blue"),cex=0.8)

其中legend=c("sub_metering_1","sub_metering_2","sub_metering_3")规定三个图例的名称，col=c("black","red","blue")规定三个图例的颜色，cex=0.8规定图例大小。

整体折线图如下：

8. 在一张图上画多个小图。

使用par函数设置一幅图的整体特征就可以做到这一点。par()函数常用的参数列表如下：

las	取值在{0,1,2,3}中，规定轴标签的位置，比如与坐标轴平行，垂直等等
bg	背景颜色
mar	边缘大小，一般用mar=c(1,2,3,4)制定，1,2,3,4对应的数字分别制定下，左，上，右的边缘大小
oma	图外部的边缘大小，默认为0
mfrow	每行图的数量
mfcol	每列图的数量

上述有一点关于mfrow和mfcol的区别，通过使用mfrow(2,2)和mfcol(2,2)我们都可以得到2*2张子图，二者的区别在于作图的顺序不同,mfrow(2,2)通过左上->右上->左下->右下顺序作图，而mfcol(2,2)通过左上->左下->右上->右下的方式作图，虽然不知道这种区别有什么用，先放在这里了。

如上表所示，如果我们改变mfrow和mfcol的数值就可以在一幅图上画出mfrow*mfcol张小图，代码如下：

par(oma=c(0,0,0,0),mfrow=c(2,2),mar=c(4,4,2,2))

把上述三幅图加上另外一幅图一共四幅图画在一张图上：

 1 #读入数据
 2 x<-read.table("D:/coursera/data analysis/household_power_consumption.txt",sep=";", header=T, colClasses = c('character', 'character', 'numeric','numeric', 'numeric', 'numeric','numeric', 'numeric', 'numeric'),na.strings='?')
 3 
 4 #整个图中画2*2幅小图
 5 par(oma=c(0,0,0,0),mfrow=c(2,2),mar=c(4,4,2,2))
 6 
 7 #处理时间
 8 x$DateTime <- strptime(paste(x$Date, x$Time),"%d/%m/%Y %H:%M:%S")                                                                                                            
 9 s <- subset(x,as.Date(DateTime) >= as.Date("2007-02-01")&as.Date(DateTime) <= as.Date("2007-02-02"))
10 s$Time <- strptime(s$Time,"%Y-%m-%d %H:%M:%S")
11 
12 #第一张图
13 hist(s$Global_active_power,freq=TRUE,col="red",xlab="Global Active Power(kilowatts)",ylab="Frequency",xaxt="n",yaxt="n",main="Global Active Power")
14 axis(side=2,at=seq(0,1200,200))
15 axis(side=1,at=NULL)
16 
17 #第二张图
18 plot(s$DateTime, s$Global_active_power, xlab="n", ylab="Global Active Power(kilowatt)", type="l",lty=1)
19 axis(side=2,at=seq(0,6,2))
20 
21 #第三张图
22 plot(s$DateTime, s$Sub_metering_1, yaxt="n", ylab="Energy Sub metering",type="l")
23 lines(s$DateTime, s$Sub_metering_2, col="red")
24 lines(s$DateTime, s$Sub_metering_3, col="blue")
25 axis(side=2,at=seq(0,30,10))
26 legend("topright",legend=c("sub_metering_1","sub_metering_2","sub_metering_3"),col=c("black","red","blue"),cex=0.8)
27 
28 #第四章图
29 plot(s$DateTime, as.numeric(s$Global_reactive_power),xlab="datetime", ylab="Global reactive power", type="h",lty=1)
30 axis(side=2,at=seq(0.0,0.5,0.1))
31 dev.copy(png,filename="plot4.png",height=480, width=480,bg="white")
32 dev.off()