在分析数据时我们经常会遇到将变量值转换成其他的值的情况(如:将连续变量转成分类变量)这时就需要我们对原有数据进行重新编码。本文将介绍R软件中常用的三种重编吗方法:
1、使用逻辑判断式编码。
2、使用cut函数编码。
3、使用car程序包的recode函数。
(一)使用逻辑判断式
(1)现假设我们需要将下面的连续型变量x按照10与20分成三个组,新的分组名称为1、2、3:
> x2=1*(x<=10)+2*(x>10&x<=20)+3*(x>20) > x2 [1] 1 2 3 2 3 3 3 3 1 3 3 2 1 2 3 3 3 2 3 3
将上述变量的数字编码改为字符编码
> labels=c("A","B","C") > x3=labels[x2] > x3 [1] "A" "B" "C" "B" "C" "C" "C" "C" "A" "C" "C" "B" "A" "B" "C" "C" "C" "B" "C" "C"
假设如下将以下范例月收入数据分成“低收入”,“中等收入”,“高收入”三个组:
> income<-c(130065,82961,133076,123028,108945,173466,17477) > income [1] 130065 82961 133076 123028 108945 173466 17477 > newcodes=c("低收入","中等收入","高收入") Error: unexpected input in "newcodes=c("低收入"? > newcodes=c("低收入","中等收入","高收入") > index=1*(income<20000)+2*(income>=20000&income<=60000)+3*(income>60000) > income=newcodes[index] > income [1] "高收入" "高收入" "高收入" "高收入" "高收入" "高收入" "低收入"
(2)使用ifelse函数
基本语法:ifelse(逻辑判断式,TRUE-表达式,FALSE-表达式)
编码成两个分组:
> x [1] 4 12 50 18 50 22 23 46 8 46 36 18 10 14 35 48 23 17 29 30 > (x2=ifelse(x<=30,1,2)) [1] 1 1 2 1 2 1 1 2 1 2 2 1 1 1 2 2 1 1 1 1 > (x3=ifelse(x<=30,"A","B")) [1] "A" "A" "B" "A" "B" "A" "A" "B" "A" "B" "B" "A" "A" "A" "B" "B" "A" "A" "A" "A"
搭配%int%运算符,将"A",“C"重编码为"Group1”,“B”,“D"重编码为"Group2”:
> y [1] "B" "A" "C" "C" "B" "A" "D" "B" "C" "D" > (y2=ifelse(y %in% c("A","C"),"Group1","Group2")) [1] "Group2" "Group1" "Group1" "Group1" "Group2" "Group1" "Group2" "Group2" "Group1" "Group2"
当编码成三个或者三个以上的组时需要多次使用ifelse 函数:
将x按照10与20两个分割点分成1、2、3三组:
> x [1] 4 12 50 18 50 22 23 46 8 46 36 18 10 14 35 48 23 17 29 30 > (x2=ifelse(x<=10,1,ifelse(x<20,2,3))) [1] 1 2 3 2 3 3 3 3 1 3 3 2 1 2 3 3 3 2 3 3
将“A”“E”编码为1,"C"编码为2,“B”“D”编码为3:
> y [1] "B" "A" "C" "C" "B" "A" "D" "B" "C" "D" > y2=ifelse(y%in%c("A","E"),1,ifelse(y=="C",2,3)) > y2 [1] 3 1 2 2 3 1 3 3 2 3
(二)使用cut 函数
cut函数可以根据我们设置的分割点(breaks)将数据重编码,将一个数值向量变量转换为分组形态的factors变量。
基本语法:
cut(x,breaks,labels,include.lowest=F,right=T)
其中
x为数值向量
breaks为分割点信息。若breaks为向量,则根据向量中的数字进行分割。若breaks为大于1正整数k,则将
x分成均等的k组。
labels为分割后各组的名称,若为null,则输出数字向量,否则输出factor变量。
include.lowest=FALSE表示分割时不含各区间端点的最小值。
right=T表示各区间为左端open,右端closed的区间
使用cut函数将x向量依照0、10、20,max(x)分成3组
> x [1] 4 12 50 18 50 22 23 46 8 46 36 18 10 14 35 48 23 17 29 30 > x2=cut(x,breaks = c(0,10,20,max(x)),labels = c(1,2,3)) > x2 [1] 1 2 3 2 3 3 3 3 1 3 3 2 1 2 3 3 3 2 3 3 Levels: 1 2 3 > as.vector(x2) [1] "1" "2" "3" "2" "3" "3" "3" "3" "1" "3" "3" "2" "1" "2" "3" "3" "3" "2" "3" "3"
若没有给定labels参数,cut函数自动按照分割点生成分组名称:
x3=cut(x,breaks = c(0,10,20,max(x))) > x3 [1] (0,10] (10,20] (20,50] (10,20] (20,50] (20,50] (20,50] (20,50] (0,10] (20,50] (20,50] (10,20] (0,10] [14] (10,20] (20,50] (20,50] (20,50] (10,20] (20,50] (20,50] Levels: (0,10] (10,20] (20,50]
现在我们模拟产生10个N(60,10)的随机成绩,并且使用cut函数的breaks选项将其分成5个组:
> score=round(rnorm(10,60,10)) > score [1] 39 65 60 69 58 69 70 62 61 75 > score.cut=cut(score,breaks=5) > score.cut [1] (39,46.2] (60.6,67.8] (53.4,60.6] (67.8,75] (53.4,60.6] (67.8,75] (67.8,75] (60.6,67.8] [9] (60.6,67.8] (67.8,75] Levels: (39,46.2] (46.2,53.4] (53.4,60.6] (60.6,67.8] (67.8,75]
由以上结果可知,cut()函数默认输出一个factor变量,并且自动将五个分组命名为“ (39,46.2]”…“ (67.8,75]”。如果cut()的选项labels=FALSE,则输出的结果是数字编码的一半向量变量:
> score.cut=cut(score,breaks=5,labels = F) > score.cut [1] 1 4 3 5 3 5 5 4 4 5 > score.cut=as.factor(score.cut) > score.cut [1] 1 4 3 5 3 5 5 4 4 5 Levels: 1 3 4 5
(三)使用car程序包中的recode函数
car程序包的recode函数可以将数值或者字符向量、factor变量重新编码。
基本语法:recode(x,recodes,as.factor.result,levels)
其中:
x为数值向量,字符向量或者factor 变量。
recode为设定重新编码规则的字符串。
as.factor.result为是否输出factor变量。若是则为TRUE,不是为FALSE。
levels为排序向量。指定新的编码分组的顺序(默认是按照分组名称排序)。
recodes参数编码规则的写法
recodes参数的值是一个字符串,字符串里面是以分号分隔的编码规则:
recodes=“规则1;规则2…”
每一个编码规则的格式为旧码列表=新码,“旧码列表”部分可用lo代表旧码的最小值(low)、hi代表旧码的最大值(high)撰写规则如下:
(1)旧码=新码 旧码只有单一数值。例如:“0=NA”表示将0改为NA。
(2)旧码向量=新码 多个旧码改为一个新码。例如:“c(7,8,9)=‘high'”,将7,8,9改为high。
(3)start:end=新码 有序数字改码。例如:“lo:19=‘C'”。
(4)else=新码 所有其他情况。例如:“else=NA”。
程序范例:
> library(carData) > library(car) > x [1] 1 2 3 1 2 3 1 2 3 > recode(x,"c(1,2)='A';else='B'") [1] "A" "A" "B" "A" "A" "B" "A" "A" "B"
将成绩0~40分之间的分数编码为1,41-60分之间为2,61-80分为3,81以上为4,其他情况为NA
> score [1] 75 70 66 65 55 69 75 69 82 83 > recode(score,"lo:40=1;41:60=2;61:80=3;81:hi=4;else=NA") [1] 3 3 3 3 2 3 3 3 4 4
上例改用‘A'‘B'‘C''‘D'
> recode(score,"lo:40='A';41:60='B';61:80='C';81:hi='D';else=NA") [1] "C" "C" "C" "C" "B" "C" "C" "C" "D" "D"
补充:R语言字符串处理时的编码的一些问题问题(Windows环境)
R在windows中文系统时默认使用GB编码字符,而许多函数是只支持UTF-8的,例如
不作任何处理读入数据时,数据会以系统默认编码读入,tolower()函数工作失败
如果有干预的读入数据,即指定字符串的编码,则不会出错
数据库读写中文出现乱码的原理也是相似的。
以上为个人经验,希望能给大家一个参考,也希望大家多多支持脚本之家。如有错误或未考虑完全的地方,望不吝赐教。