3.3缺失值处理
R中缺失值以NA表示，判断数据是否存在缺失值的函数有两个，最基本的函数是is.na()它可以应用于向量、数据框等多种对象，返回逻辑值。
> attach(data)
The following objects are masked fromdata (pos = 3):
city, price, salary
> data$salary=replace(salary,salary>5,NA)
> is.na(salary)
[1] FALSEFALSE TRUE FALSE TRUE TRUE TRUE FALSE FALSE FALSE FALSEFALSE
> sum(is.na(salary))
[1] 4
另一个判断缺失值的函数是complete.cases()，它同样返回逻辑值向量，但值与is.na()的相反:缺失值为FALSE，正常数据为TRUE，利用它来选取无缺失数据的行非常方便。
> complete.cases(data$salary)
[1] TRUE TRUE FALSE TRUE FALSE FALSEFALSE TRUE TRUE TRUE TRUE TRUE
3.3.2判断缺失模式
存在缺失数据时，需要进一步判断数据的缺失模式，判断是否是随机的，然后才能确定处理的方法。
程序包mice，利用链式方程进行多元插补，可以处理混合变量类型的数据缺失，自动产生填补变量的预测变量，是处理缺失值的重要工具。
> library(mice)
> data$price=replace(price,price>5,NA)
> md.pattern(data)
price salary city
5 1 1 0 1
3 0 1 0 2
4 1 0 0 2
3 4 12 19
输出结果中的“1”表示没有缺失数据，“0”表示存在缺失数据。第1列第1行的“5”表示有5个样本是完整的，下面的“3”表示有3个样本缺少了salary这一变量的值，第1列最后一个数字“4”表示有4条记录在salary和price上都有缺失。最后一行表示各个变量缺失的样本数合计。
程序包VIM提供了在R中探索数据缺失情况的新工具，实现缺失模式的可视化
> library(VIM)
> aggr(data)
第一个图由小条形的长度显示各变量缺失数据比例
第二个图显示了综合的缺失模式，可以与md.pattern()生成的结果对照观察，其中浅色方框表示完整数据，深色框表示缺失值。底部的颜色框高度反映了相应组合的频率。
3.3.3处理缺失数据
(1)删除缺失样本
过滤掉缺失样本是最简单的方式，其前提是缺失数据的比例较少，而且缺失数据是随机出现的，这样删除缺失数据后对分析结果影响不大。 R可以使用complete.cases()指令选取完整的记录，有缺失值的行则删去不要。
> data1=data[complete.cases(data$salary),]
> dim(data1)
[1] 8 3
或
> data2=data[!is.na(salary),]
> dim(data2)
[1] 8 3对于有多个变量缺失的数据，如果想直接删除所有的缺失值，可以通过na.omit()函数来完成，
> data3=na.omit(data)
> dim(data3)
[1] 5 3
(2)替换缺失值
> data[is.na(data)]=mean(salary[!is.na(salary)])
(3)多重插补法
多重插补(Multiple Imputation)是用于填补复杂数据缺失值的一种方法，该方法通过变量间关系来预测缺失数据，利用蒙特卡罗随机模拟方法生成多个完整数据集，再对这些数据集分别进行分析，最后对这些分析结果进行汇总处理。FSC是基于链式方程的插补方法，因此也称为MICE (Multiple Imputation by Chained Equations )。它与其他多重插补算法的本质区别是，它在进行插补时不必考虑被插补变量和协变量的联合分布，而是利用单个变量的条件分布逐一进行插补。在R语言中通过程序包mice中的函数mice()可以实现该方法，它随机模拟多个完整数据集并存入imp，再对imp进行线性回归，最后用pool函数对回归结果进行汇总。
3.4数据整理
3.4.1数据合并
（1）函数cbind(),rbind()
> a=c("hk",12,10)
> data1=rbind(data,a)
> data1
cityprice salary
………
12 qa 6 5
13 hk 12 10
(2)构造data.frame
对数据“整容”最简单的思路是把数据向量化，再按要求用向量构建其他类型的对象。一些结构相似的对象，如向量(数值型、字符型、逻辑型)、因子、数值矩阵、列表或其他数据框等，可以被合并为一个数据框。
> weight=c(150,135,210,140)
> height=c(65,61,70,65)
> gender=c("F","F","M","F")
> stu=data.frame(weight,height,gender)
> stu
weightheight gender
1 150 65 F
2 135 61 F
3 210 70 M
4 140 65 F
合并时，变量名称就白动变成了新数据框的列名，也可以用names()重新给其赋值。
> row.names(stu)=c("Alice","Bob","Cal","David")
> stu
weightheight gender
Alice 150 65 F
Bob 135 61 F
Cal 210 70 M
David 140 65 F
(3)函数merge()
在R中合并两个数据集可以通过专门的函数merge()来实现。merge通过相同的列或行名来识别，合并两个数据框或列表，其调用格式如下:
merge(x, y, by = intersect(names(x),names(y)),by.x = by, by.y = by, all = FALSE, all.x = all, all.y = all,sort =TRUE, suffixes = c(".x",".y"),
incomparables = NULL, ...)
x,y 要合并的数据集
by指定合并的依据(相同的行或列)
by.x by.y分别为第一个数据框和第二个数据框要连接的列名
all, all.x, all.y逻辑值，默认为FALSE。
> index=list("city"=data$city,"index"=1:12)
> index
$city
[1]"bj" "sh" "gz" "ab" "cd""as" "ac" "fa" "ff" "ee""er" "qa"
$index
[1] 1 2 3 4 5 6 7 8 9 10 11 12
> data.index=merge(data,index,by="city")
> data.index
cityprice salary index
1 ab 6 5 4
2 ac 5 NA 7
3 as 3 NA 6
4 bj 1 2 1
5 cd 1 NA 5
6 ee 3 4 10
7 er 5 3 11
8 fa 6 1 8
9 ff 1 2 9
10 gz 5 NA 3
11 qa 6 5 12
12 sh 3 4 2
3.4.2选取数据的子集
在R中，选取数据子集用中括号[]
> data[data$salary>6]
3.4.3数据排序
R中的排序函数sort()只能对向量进行简单的排序，对含有多变量的数据集，需要用order指令来完成，其调用格式如下:
order(..., na.last = TRUE, decreasing =FALSE)
> order.price=order(data$price)
或
> sort.list(data$price)
指令order返回向量排序后各数字的原始位置，与之非常相关的指令是秩(rank )，它返回每个数字在整个向量中的秩，可以简单地理解为各个数字的大小顺序。
> rank(data$price)
3.5长宽格式的转换.
>t (data)
3.5.1揉数据函数
R中有两个揉数据函数stack()和unstack|()，用于数据长格式和宽格式之间的转换.
stack()把一个数据框转换成两列:一列为数据，另一列为数据对应的列名称。
unstack()是stack的逆过程，被转换的对象包含两列，它把数据列按照因子列的不同水平重新排列，分离为不同的列。
3.5.2揉数据的最佳伴侣
程序包reshape2是reshape的重写版，是专门用于数据集形状转换的，一般用户常使用melt(), acast()和dcast()，它们却可以把数据“揉成各种形状。
melt本身的意思是溶解、分解，其作用在一个数据集上其实就是拆分数据，它的对象一可以是数组(array )、数据框或列表。
> library(reshape2)
> data(airquality)
> str(airquality)
'data.frame': 153obs. of 6 variables:
$Ozone : int 41 36 12 18 NA 28 23 19 8 NA ...
$ Solar.R:int 190 118 149 313 NA NA 299 99 19 194...
$Wind : num 7.4 8 12.6 11.5 14.3 14.9 8.6 13.8 20.1 8.6...
$Temp : int 67 72 74 62 56 66 65 59 61 69 ...
$Month : int 5 5 5 5 5 5 5 5 5 5 ...
$ Day : int 1 2 3 4 5 6 7 8 9 10 ...
> longdata=melt(airquality,id.vars=c("Ozone",'Month',"Month","Day"),measure.vars=2:4)
> str(longdata)
'data.frame': 459obs. of 6 variables:
$Ozone : int 41 36 12 18 NA 28 23 19 8 NA ...
$Month : int 5 5 5 5 5 5 5 5 5 5 ...
$Month : int 5 5 5 5 5 5 5 5 5 5 ...
$ Day : int 1 2 3 4 5 6 7 8 9 10 ...
$variable: Factor w/ 3 levels "Solar.R","Wind",..: 1 1 1 1 11 1 1 1 1 ...
$value : num 190 118 149 313 NA NA 299 99 19 194 ...
利用ggplot2在一个图形中多维度地展示value值
> library(ggplot2)
> p=ggplot(data=longdata,aes(x=Ozone,y=value,color=factor (Month)))
> p+geom_point(shape=20,size=4)+facet_wrap(~variable,scales="free_y")+geom_smooth(aes(group=1),fill="gray80")
和stack()一样，melt()也有对应的函数用来还原数据:acast()用于数组，dcast()用于数据框，其中的参数formula是一个公式，左边的每个变量都会成为新数据集中的一列，右边的变量是因子，其每个水平行在新数据集中成为一列，从而把长格式数据转换为短格式。
