2

我正在尝试将一列添加到由一个因子归一化值组成的数据框中。

例如：

'data.frame':   261 obs. of  3 variables:
 $ Area   : Factor w/ 29 levels "Antrim","Ards",..: 1 1 1 1 1 1 1 1 1 2 ...
 $ Year   : Factor w/ 9 levels "2002","2003",..: 1 2 3 4 5 6 7 8 9 1 ...
 $ Arrests: int  18 54 47 70 62 85 96 123 99 38 ...

我想添加一列，该列是按区域分组归一化的逮捕值。

我想出的最好的是：

data$Arrests.norm <- unlist(unname(by(data$Arrests,data$Area,function(x){ scale(x)[,1] } )))

该命令处理但数据被加扰，即标准化值与数据帧中的正确区域不匹配。

感谢您的提示。

编辑：只是为了澄清我所说的加扰数据的含义，在我的代码之后对数据框进行子集化，我得到如下输出，其中标准化值显然属于另一个因子组。

      Area Year Arrests Arrests.norm
199 Larne 2002      92 -0.992843957
200 Larne 2003     124 -0.404975825
201 Larne 2004      89 -1.169204397
202 Larne 2005      94 -0.581336264
203 Larne 2006      98 -0.228615385
204 Larne 2007       8  0.006531868
205 Larne 2008      31  0.418039561
206 Larne 2009      25  0.947120880
207 Larne 2010      22  2.005283518
4

2 回答 2

3

跟进您的by尝试：

df <- data.frame(A = factor(rep(c("a", "b"), each = 4)),
                 B = sample(1:4, 8, TRUE))

ll <- by(data = df, df$A, function(x){
  x$B_scale <- scale(x$B)
  x
  }
   )

df2 <- do.call(rbind, ll)
于 2013-10-13T22:32:46.603 回答
2 
data <- transform(data, Arrests.norm = ave(Arrests, Area, FUN = scale))

会成功的。

于 2013-10-13T22:38:18.400 回答