3

如何按物种规范化 Sepal.Length 的值？

    iris
    Sepal.Length Sepal.Width Petal.Length Petal.Width    Species
1            5.1         3.5          1.4         0.2     setosa
...

# i have to divide by 
tapply(iris$Sepal.Length, iris$Species, max)
    setosa versicolor  virginica 
       5.8        7.0        7.9

换句话说，我想将所有值Species=="setosa"除以 5.8，依此类推，最后我希望在 Sepal.Length 列中有一个标准化值为 0..1 的数据框。

最后它应该返回

    iris
    Sepal.Length Sepal.Width Petal.Length Petal.Width    Species
1      0.8793103         3.5          1.4         0.2     setosa
...
4

4 回答 4

7

显然有很多方法可以做到这一点。我最喜欢ave()（参见 DWin 的答案）或data.table包的语法：

library(data.table)
dt <- data.table(iris)
dt[, Sepal.Length:=(Sepal.Length)/max(Sepal.Length), by="Species"]
dt
#      Sepal.Length Sepal.Width Petal.Length Petal.Width   Species
#   1:    0.8793103         3.5          1.4         0.2    setosa
#   2:    0.8448276         3.0          1.4         0.2    setosa
#   3:    0.8103448         3.2          1.3         0.2    setosa
#   4:    0.7931034         3.1          1.5         0.2    setosa
#   5:    0.8620690         3.6          1.4         0.2    setosa
# 146:    0.8481013         3.0          5.2         2.3 virginica
# 147:    0.7974684         2.5          5.0         1.9 virginica
# 149:    0.7848101         3.4          5.4         2.3 virginica
# 150:    0.7468354         3.0          5.1         1.8 virginica

df <- data.frame(dt) ## It's possible (but not necessary) to coerce back to
                     ## a plain old data.frame
于 2012-10-25T20:19:01.817 回答
5

我严格解释您除以最大值的愿望。

一种选择：

aggregate(iris$Sepal.Length,list(iris$Species),FUN = function(x) x/max(x))

另一个，使用ddplyfrom plyr（并一次缩放所有列：

ddply(iris,.(Species),colwise(function(x){x / max(x)}))

还有一个更像@Dwinave示例的变体，它使其他列保持不变，但使用ddply：

ddply(iris,.(Species),transform,Sepal.Length = Sepal.Length / max(Sepal.Length))
于 2012-10-25T20:20:20.227 回答
3 
  iris$ratio_to_max <- ave( iris$Sepal.Length, list(iris$Species), 
                                                     FUN= function(x) x/max(x))
#-------------
> str(iris)
'data.frame':   150 obs. of  6 variables:
 $ Sepal.Length: num  5.1 4.9 4.7 4.6 5 5.4 4.6 5 4.4 4.9 ...
 $ Sepal.Width : num  3.5 3 3.2 3.1 3.6 3.9 3.4 3.4 2.9 3.1 ...
 $ Petal.Length: num  1.4 1.4 1.3 1.5 1.4 1.7 1.4 1.5 1.4 1.5 ...
 $ Petal.Width : num  0.2 0.2 0.2 0.2 0.2 0.4 0.3 0.2 0.2 0.1 ...
 $ Species     : Factor w/ 3 levels "setosa","versicolor",..: 1 1 1 1 1 1 1 1 1 1 ...
 $ ratio_to_max: num  0.879 0.845 0.81 0.793 0.862 ...

如果你想更换Sepal.Length柱子，你可以这样做，但我通常会避免这种破坏性的做法，直到我真的确定我得到了我想要的。（即使这样我也感到内疚。）如果你希望它在单独的列表“数据包”中并丢弃原来的“Sepal.Length”列，你可以使用split：

 spl.iris <- split(iris[-1], iris$Species)
 str(spl.iris)
于 2012-10-25T20:27:08.227 回答
0

我确信有更好的 plyr 或数据表，甚至是基本方法：

L1 <- lapply(split(iris[, -5], iris$Species), function(x) apply(x, 2, scale))
L2 <- lapply(seq_along(L1), function(i) {
    data.frame(SPecies=names(L1)[i], L1[[i]])
})
do.call(rbind, L2)
于 2012-10-25T20:15:26.507 回答