作者：包寒吳霜

中科院心理所碩士在讀(名字/人格/社會/文化心理學(xué))

知乎：https://www.zhihu.com/people/psychbruce

社會科學(xué)研究經(jīng)常會遇到“超多變量”的情況——多量表、多維度、多題項，以及復(fù)雜的正反計分題……如何更高效地計算量表總分？如何更簡潔地進(jìn)行反向計分？

本文將為大家分享如何使用R語言(data.table包 + 自編函數(shù))優(yōu)雅地計算多變量。

當(dāng)我們的數(shù)據(jù)中存在成百上千個變量時，不僅變量管理存在一定的難度，而且變量計算也會變得比較復(fù)雜。如果使用R來處理，有哪些現(xiàn)存的方法呢？

• dplyr包的mutate函數(shù)(需再次賦值給data)

1data=mutate(data,?X=x1+x2+x3,?Y=y1+y2+y3)

• data.table包的“:=”函數(shù)(原地更新，無需賦值)

1#?單變量計算2data[,?X:=x1+x2+x3]3#?多變量同時計算4data[,?":="(X=x1+x2+x3,?Y=y1+y2+y3)]

變量少還好說，然而，一旦遇到幾十個甚至幾百個變量，我們都希望利用更簡便的方法來計算總分或平均分，而不是一個個敲變量名。

首先，直接使用sum(...)或mean(...)是行不通的，大家可以自行嘗試(算出來的其實(shí)是每一列的總分而不是每一行的總分)。實(shí)際上，我們需要用mapply把相關(guān)函數(shù)施加于各個變量上，這樣才能分別對每一行計算總分：

1data[,?":="(Xsum=mapply(sum,?x1,?x2,?x3))]

這只是一個初步的解決思路，但依然尚未解決多變量的問題——我們并不想一個個敲變量名。

于是乎，針對大家經(jīng)常遇到的一些變量計算需求，我編制了相應(yīng)的R函數(shù)，既可以用于普通的data.frame(mutate)，也可以用于data.table(:=)。

下面以計算平均值為例，介紹自編的MEAN函數(shù)的基本用法：

1data[,?":="(Xmean1=MEAN(data,?"x",?1:50),2????????????Xmean2=MEAN(data,?vars=c("x1",?"x2",?"x3")),3????????????Xmean3=MEAN(data,?varrange="x1:x50"),4????????????Xmean4=MEAN(data,?"x",?1:50,?rev=41:50,?likert=1:7))]

MEAN的第一個參數(shù)為原來的數(shù)據(jù)(data.frame或data.table)，后面的若干參數(shù)用來定義變量范圍和反向計分題：

• 對于有規(guī)律的變量名，如x1-x50，只需要定義var="x"和item=1:50即可(此為推薦用法，并且var和item參數(shù)名可省略)

• 或者可以通過vars具體列出參與計算的變量(字符串向量)

• 另外還可以通過varrange="x1:x50"或varrange=c("x1","x50")的形式定義變量的起止范圍
  【提示：此時的1和50不代表數(shù)字上的連續(xù)范圍，也就是說，varrange定義的起止范圍是數(shù)據(jù)中變量的原始位置，如果原始變量的后綴數(shù)字是亂序排列的，則推薦使用var和item定義】

• 如果涉及反向計分題，只需要通過rev參數(shù)定義哪些題目反向計分(可以是單個數(shù)字、單個字符串、數(shù)字向量、字符串向量)，并通過likert參數(shù)定義題目是幾點(diǎn)量表(例如1-7點(diǎn)的量表，設(shè)置likert=1:7或者likert=c(1, 7)均可)

MEAN函數(shù)的源代碼如下(復(fù)制并運(yùn)行即可使用)：

1MEAN=function(data,?var=NULL,?items=NULL, 2??????????????vars=NULL, 3??????????????varrange=NULL, 4??????????????rev=NULL,?likert=NULL, 5??????????????na.rm=TRUE)?{ 6??Mean=function(...)?mean(c(...),?na.rm=na.rm) 7??if(!is.null(varrange))?{ 8????dn=names(data) 9????if(length(varrange)==1)?varrange=strsplit(varrange,?":")[[1]]10????varMin=varrange[1]11????varMax=varrange[length(varrange)]12????vars=dn[which(dn==varMin):which(dn==varMax)]13??}?else?{14????if(is.null(vars))?vars=paste0(var,?items)15??}16??if(is.character(rev))?rev=which(vars?%in%?rev)17??vars=paste(deparse(substitute(data)),?vars,?sep="$")18??pre=rep("",?length(vars))19??pre[rev]=ifelse(is.null(likert),?"",?paste0(min(likert)+max(likert),?"-"))20??varlist=paste0(pre,?vars)21??eval(parse(text=paste0("mapply(Mean,?",?paste(varlist,?collapse=",?"),?")")))22}

除了MEAN之外，還有其他的幾個自編函數(shù)，用法與MEAN類似，包括SUM(計算總分)、COUNT(統(tǒng)計某個值在多個變量中的出現(xiàn)次數(shù))、CONSEC(統(tǒng)計多個變量中連續(xù)相同數(shù)字出現(xiàn)最多的個數(shù))：

1SUM=function(data,?var=NULL,?items=NULL, 2?????????????vars=NULL, 3?????????????varrange=NULL, 4?????????????rev=NULL,?likert=NULL, 5?????????????na.rm=TRUE)?{ 6??Sum=function(...)?sum(...,?na.rm=na.rm) 7??if(!is.null(varrange))?{ 8????dn=names(data) 9????if(length(varrange)==1)?varrange=strsplit(varrange,?":")[[1]]10????varMin=varrange[1]11????varMax=varrange[length(varrange)]12????vars=dn[which(dn==varMin):which(dn==varMax)]13??}?else?{14????if(is.null(vars))?vars=paste0(var,?items)15??}16??if(is.character(rev))?rev=which(vars?%in%?rev)17??vars=paste(deparse(substitute(data)),?vars,?sep="$")18??pre=rep("",?length(vars))19??pre[rev]=ifelse(is.null(likert),?"",?paste0(min(likert)+max(likert),?"-"))20??varlist=paste0(pre,?vars)21??eval(parse(text=paste0("mapply(Sum,?",?paste(varlist,?collapse=",?"),?")")))22}232425COUNT=function(data,?var=NULL,?items=NULL,26???????????????vars=NULL,27???????????????varrange=NULL,28???????????????value=NA)?{29??Count=function(...)?sum(c(...),?na.rm=TRUE)30??if(!is.null(varrange))?{31????dn=names(data)32????if(length(varrange)==1)?varrange=strsplit(varrange,?":")[[1]]33????varMin=varrange[1]34????varMax=varrange[length(varrange)]35????vars=dn[which(dn==varMin):which(dn==varMax)]36??}?else?{37????if(is.null(vars))?vars=paste0(var,?items)38??}39??vars=paste(deparse(substitute(data)),?vars,?sep="$")40??if(is.na(value))?{41????varlist=paste0("is.na(",?vars,?")")42??}?else?{43????varlist=paste0(vars,?"==",?value)44??}45??eval(parse(text=paste0("mapply(Count,?",?paste(varlist,?collapse=",?"),?")")))46}474849CONSEC=function(data,?var=NULL,?items=NULL,50????????????????vars=NULL,51????????????????varrange=NULL,52????????????????values=0:9)?{53??Conseq=function(string,?number=values)?{54????#?Consecutive?Identical?Digits55????require(stringr)56????pattern=paste(paste0(number,?"{2,}"),?collapse="|")57????ifelse(grepl(pattern,?string),?max(nchar(str_extract_all(string=string,?pattern=pattern,?simplify=TRUE))),?0)58??}59??if(!is.null(varrange))?{60????dn=names(data)61????if(length(varrange)==1)?varrange=strsplit(varrange,?":")[[1]]62????varMin=varrange[1]63????varMax=varrange[length(varrange)]64????vars=dn[which(dn==varMin):which(dn==varMax)]65??}?else?{66????if(is.null(vars))?vars=paste0(var,?items)67??}68??vars=paste(deparse(substitute(data)),?vars,?sep="$")69??varlist=vars70??eval(parse(text=paste0("mapply(Conseq,?paste0(",?paste(varlist,?collapse=",?"),?"))")))71}

最后，我們以一個具體的例子來綜合演示上述函數(shù)的用法。在下面的這個data.table中，我特意將x2和x4的位置對調(diào)了，大家可以自行體會不同參數(shù)之間的微妙差別(如var、vars、varrange)

1d=data.table(x1=1:5,?x4=c(2,2,5,4,5),?x3=c(3,2,NA,NA,5),?x2=c(4,4,NA,2,5),?x5=c(5,4,1,4,5)) 2#????x1?x4?x3?x2?x5 3#?1:??1??2??3??4??5 4#?2:??2??2??2??4??4 5#?3:??3??5?NA?NA??1 6#?4:??4??4?NA??2??4 7#?5:??5??5??5??5??5 8 9d[,":="(n.na=COUNT(d,?"x",?1:5,?value=NA),10????????n.2=COUNT(d,?"x",?1:5,?value=2),11????????sum=SUM(d,?"x",?1:5),12????????mean1=MEAN(d,?"x",?1:5),13????????mean2=MEAN(d,?vars=c("x1",?"x4")),14????????mean3=MEAN(d,?varrange="x1:x2",?rev="x2",?likert=1:5),15????????cons1=CONSEC(d,?"x",?1:5),16????????cons2=CONSEC(d,?varrange="x1:x5"))]17#????x1?x4?x3?x2?x5?n.na?n.2?sum?mean1?mean2?mean3?cons1?cons218#?1:??1??2??3??4??5????0???1??15???3.0???1.5?????2?????0?????019#?2:??2??2??2??4??4????0???3??14???2.8???2.0?????2?????2?????320#?3:??3??5?NA?NA??1????2???0???9???3.0???4.0?????4?????0?????021#?4:??4??4?NA??2??4????1???1??14???3.5???4.0?????4?????2?????222#?5:??5??5??5??5??5????0???0??25???5.0???5.0?????4?????5?????5

Tips：三種定義變量的方式選擇其一即可，分別適用于不同需求——

• var和item會將變量按照item的數(shù)字順序重新排列(在函數(shù)內(nèi)部其實(shí)是paste了var和item)

• vars則以該參數(shù)實(shí)際定義的順序為準(zhǔn)

• varrange定義的是起止位置的變量(因此在本例中“x1:x2”實(shí)則對應(yīng)了“x1, x4, x3, x2”四個變量)

?

更多實(shí)用函數(shù)，請訪問作者的GitHub：

https://github.com/psychbruce/stats/blob/master/BruceFunctions.R

公眾號后臺回復(fù)關(guān)鍵字即可學(xué)習(xí)

回復(fù)?爬蟲???????? ? ?爬蟲三大案例實(shí)戰(zhàn)??回復(fù)?Python? ? ? ?1小時破冰入門

回復(fù)?數(shù)據(jù)挖掘?? ? R語言入門及數(shù)據(jù)挖掘回復(fù)?人工智能?? ? 三個月入門人工智能回復(fù)?數(shù)據(jù)分析師??數(shù)據(jù)分析師成長之路?回復(fù)?機(jī)器學(xué)習(xí)????? 機(jī)器學(xué)習(xí)的商業(yè)應(yīng)用回復(fù)?數(shù)據(jù)科學(xué)????? 數(shù)據(jù)科學(xué)實(shí)戰(zhàn)回復(fù)?常用算法????? 常用數(shù)據(jù)挖掘算法

萬水千山總是情，點(diǎn)個 “好看” 行不行！！！

?

總結(jié)

以上是生活随笔為你收集整理的R 回归 虚拟变量na_如何优雅地计算多变量的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯，歡迎將生活随笔推薦給好友。