4

Rのグループごとに順序統計を計算するにはどうすればよいですか。列に従って結果を集計し、グループごとに1行だけを返します。その行は、いくつかの順序に従って、グループのn番目の要素である必要があります。理想的には、基本関数のみを使用したいと思います。

x <- data.frame(Group=c("A","A", "A", "C", "C"), 
                Name=c("v", "u", "w", "x", "y"), 
                Quantity=c(3,3,4,2,0))
> x
  Group Name Quantity
1     A    v        3
2     A    u        3
3     A    w        4
4     C    x        2
5     C    y        0

数量、名前の順に並べて、n番目に高い値を取得したいと思います。N = 2の場合、これは

  Group Name Quantity
1     A    u        3
5     C    y        0

For N=1
  Group Name Quantity
3     A    w        4
4     C    x        2

次のことを試しましたが、情報のないエラーメッセージが表示されます。

 aggregate.data.frame(x, list(x$Group), function(y){ max(y[,'Quantity'])})
 Error in `[.default`(y, , "Quantity") (from #1) : incorrect number of dimensions"
4

4 に答える 4

2 
x <- 
    data.frame(
        Group = c("A","A", "A", "C", "C", "A", "A") , 
        Name = c("v", "u", "w", "x", "y" ,"v", "u") , 
        Quantity = c(3,3,4,2,0,4,1)
    )

# sort your data to start..
# note that Quantity vs. Group and Name
# are sorted in different directions,
# so the -as.numeric() flips them
x <- 
    x[ 
        order( 
            -as.numeric( x$Group ) , 
            x$Quantity , 
            -as.numeric( x$Name ) , 
            decreasing = TRUE 
        ) , 
    ]
# once your data frame is sorted the way you want your Ns to occur, the rest is easy

# rank your data..  
# just create the numerical order, 
# but within each group..
# (or you could add those ranks directly to the data frame if you like)
ranks <- 
    unlist( 
        tapply( 
            order( x$Group ) , 
            as.numeric( x$Group ) , 
            order 
        ) 
    )

# N = 1
x[ ranks == 1 , ]

# N = 2
x[ ranks == 2 , ]
于 2013-01-20T03:45:25.900 に答える
1

いくつかの集約-マージマジック：

f <- function(x, N) {
  sel <- function(x) {                                   # Choose the N-th highest value from the set, or lowest element if there < N unique elements.  Is there a built-in for this? 
    z <- unique(x)                                       # This assums that you wan the N-th highest unique value.  Simply don't filter by unique if not.
    z[order(z, decreasing=TRUE)][min(N, length(z))]
  }

  xNq <- aggregate(Quantity ~ Group, data=x,   sel)      # Choose the N-th highest quantity within each "Group"
  xNm <- merge(x, xNq)                                   # Add the matching "Name" values
  x <- aggregate(Name ~ Quantity + Group, data=xNm, sel) # Choose the N-th highest Name in each group
  x[c('Group', 'Name', 'Quantity')]                      # Put into original order
}


> f(x, 2)
##   Group Name Quantity
## 1     A    u        3
## 2     C    y        0

> f(x, 1)
##   Group Name Quantity
## 1     A    w        4
## 2     C    x        2
于 2013-01-20T03:42:30.530 に答える
1 
# define ordering function, increasing on Quantity, decreasing on Name
in.order <- function(group) with(group, group[order(Quantity, -rank(Name)), ])

# set desired rank for each Group
N <- 2

# get Nth row by Group, according to in.order
group.rows <- by(x, x$Group, function(group) head(tail(in.order(group), N), 1))

# collapse rows into data.frame
do.call(rbind, group.rows)

#   Group Name Quantity
# A     A    u        3
# C     C    y        0

このエラーが表示aggregate.data.frameされる理由は、この関数が引数FUNに従って各列に適用されby、完全な各サブセットではなく適用されるためですdata.frame(上記のように、by関数の目的はそのためです)。を使用する場合aggregate、 に提供するものはすべて、 sFUNではなく列を受け入れる必要がありますdata.frame。あなたの例では、のyようにベクトルにインデックスを付けようとしているdata.frameため、次元エラーが発生します。

于 2013-01-20T06:16:06.623 に答える
0

一緒に行きました

do.call(rbind, by(x, x$Group, function(x)
      x[order(-x$Quantity, x$Name),][1,]))

他の誰かの提案に従って。他の投稿されたソリューションよりも私の思考プロセスに少し適していることがわかりました（感謝しています）。

于 2013-01-20T22:47:20.527 に答える