2

我想基于两列拆分数据帧，但我希望输出是数据帧的二维矩阵，而不是数据帧的平面列表。我可以实现我想要使用的东西by()，subset但我被告知（我认为是 Ripley）应该避免subset在包开发中使用。是否有split保留暗名的优雅替代方案（可能使用 ）？

# sample data
df <- data.frame(x=rnorm(20), y=rnorm(20), v1=rep(letters[1:5],each=4), v2=rep(LETTERS[6:9]))

# what I did previously
submat <- by(df, list(df$v1,df$v2), subset)
dim(submat) # 5 x 4
dimnames(submat) # "a" "b" "c" "d" "e" ; "F" "G" "H" "I"
4

1 回答 1

2

要获得您所要求的数据帧矩阵，请tapply与返回特定数据帧子集但行名称与因子级别匹配的函数一起使用。

> dfmat <- with(df, tapply(1:NROW(df), list(v1,v2), function(idx) df[idx,] ) )
> dfmat[1,1]  # items that are in a single dataframe accessed via matrix indexing
[[1]]
           x         y v1 v2
1 -0.5604756 -1.067824  a  F

> dfmat
  F      G      H      I     
a List,4 List,4 List,4 List,4
b List,4 List,4 List,4 List,4
c List,4 List,4 List,4 List,4
d List,4 List,4 List,4 List,4
e List,4 List,4 List,4 List,4

将列表作为条目的矩阵被print-ed 以仅显示对象类型和条目数（在本例中为列）。注意每个条目是一个包含一个项目的列表，这样dataframe属性保持不变，但需要“向下钻取”才能获得宝藏：编辑：添加了dfmat的属性：

>  attributes(dfmat)
$dim
[1] 5 4

$dimnames
$dimnames[[1]]
[1] "a" "b" "c" "d" "e"

$dimnames[[2]]
[1] "F" "G" "H" "I"    
#------------
> attributes( dfmat[1,1])
NULL
#------------
> attributes( dfmat[1,1][[1]])
$names
[1] "x"  "y"  "v1" "v2"

$row.names
[1] 1

$class
[1] "data.frame"
于 2012-11-29T19:52:48.603 回答