r减去平均值并除以几个变量的标准差

时间:2017-04-06 21:41:44

标签: r apply lapply sapply

我正在尝试标准化数据框中的某些列,而不是所有列。通过标准化我的意思是,减去平均值并除以标准偏差。我的问题是如何才能对值进行标准化 第1,2,4和6列假设我正在处理此data(mtcars)数据集。

我可以手动完成此操作,但我很想知道是否有一种有效的方法。

2 个答案:

答案 0 :(得分:3)

scale为您完成此操作。所以

df<-mtcars
df[,c(1,2,4,6)]<-scale(df[,c(1,2,4,6)])

将保持其他变量不变。 scale返回mean和sd作为可用于反转过程的属性。

答案 1 :(得分:0)

mt <- mtcars
str(mt)
# 'data.frame': 32 obs. of  11 variables:
#  $ mpg : num  21 21 22.8 21.4 18.7 18.1 14.3 24.4 22.8 19.2 ...
#  $ cyl : num  6 6 4 6 8 6 8 4 4 6 ...
#  $ disp: num  160 160 108 258 360 ...
#  $ hp  : num  110 110 93 110 175 105 245 62 95 123 ...
#  $ drat: num  3.9 3.9 3.85 3.08 3.15 2.76 3.21 3.69 3.92 3.92 ...
#  $ wt  : num  2.62 2.88 2.32 3.21 3.44 ...
#  $ qsec: num  16.5 17 18.6 19.4 17 ...
#  $ vs  : num  0 0 1 1 0 1 0 1 1 1 ...
#  $ am  : num  1 1 1 0 0 0 0 0 0 0 ...
#  $ gear: num  4 4 4 3 3 3 3 4 4 4 ...
#  $ carb: num  4 4 1 1 2 1 4 2 2 4 ...

诀窍是在*apply调用以及重新分配(<-=的左侧)中对其进行子集化。

mysd <- 3 # something important

mt[c(1,2,4,6)] <- lapply(mt[c(1,2,4,6)], `+`, mysd)
str(mt)
# 'data.frame': 32 obs. of  11 variables:
#  $ mpg : num  24 24 25.8 24.4 21.7 21.1 17.3 27.4 25.8 22.2 ...
#  $ cyl : num  9 9 7 9 11 9 11 7 7 9 ...
#  $ disp: num  160 160 108 258 360 ...
#  $ hp  : num  113 113 96 113 178 108 248 65 98 126 ...
#  $ drat: num  3.9 3.9 3.85 3.08 3.15 2.76 3.21 3.69 3.92 3.92 ...
#  $ wt  : num  5.62 5.88 5.32 6.21 6.44 ...
#  $ qsec: num  16.5 17 18.6 19.4 17 ...
#  $ vs  : num  0 0 1 1 0 1 0 1 1 1 ...
#  $ am  : num  1 1 1 0 0 0 0 0 0 0 ...
#  $ gear: num  4 4 4 3 3 3 3 4 4 4 ...
#  $ carb: num  4 4 1 1 2 1 4 2 2 4 ...

请注意,lapply的返回值为list,而不是data.frame。虽然它的行为通常完全相同,但您可以使用as.data.frame(lapply(...))将其包装起来以将其返回到原始类。

对多列进行单个修改的一种流行方法是形成logical向量(可以比整数更安全),例如这个过于简化的示例。使用矢量使随后的重新分配更容易阅读。

vec <- sapply(mt, function(x) min(x)>10)
mt[vec] <- lapply(mt[vec], `+`, mysd)

(如果整数向量包含低于1或高于列数的任何内容,则使用整数变得不那么可预测/健壮。它可以与integer(0)一起使用,所以如果需要,可以随意使用整数。)

这样做的一个很好的副作用是,如果功能是昂贵的&#34; (时间或资源),然后它只在相关列上运行。如果没有选择任何内容,则不执行任何操作。

vec <- sapply(mt, function(x) min(x) > 300)
any(vec)
# [1] FALSE
system.time( mt[vec] <- lapply(mt[vec], function(x) { Sys.sleep(100); x+1; }) )
#    user  system elapsed 
#       0       0       0