1

これは実際のデータの簡単な例で、N>100 の NxN 行列です。今、私が見つけたいのは、各行の上位 3 つです。次に、上位 3 つに最も頻繁に含まれる列と、それが発生する行の下を見つけます。

        A   B   C   D   E   F   G
    A   0   70  5   73  96  46  58
    B   47  0   20  89  75  50  19
    C   42  98  0   30  30  22  76
    D   66  20  18  0   63  18  60
    E   73  0   63  51  0   23  7
    F   79  34  61  56  12  0   99
    G   25  26  41  86  51  30  0

ここで、セルが行の上部 3 にある場合は 1、それ以外の場合は 0 のバイナリ マトリックスを作成して、同じ結果を見つけようとしました。ただし、長方形ではないテーブルを作成するのが最も困難です。

        A   B   C   D   E   F   G
    A   0   1   0   1   1   0   0
    B   0   0   0   1   1   1   0
    C   1   1   0   0   0   0   1
    D   1   0   0   0   1   0   1
    E   1   0   1   1   0   0   0
    F   1   0   1   0   0   0   1
    G   0   0   1   1   1   0   0

最終的な結果は次のようになります。条件 C、D、E、および F の下で、A が上位 3 位に入っています。

    A       D       E       C       G       B       F
    C (42)  A (73)  A (96)  E (63)  C (76)  A (70)  B (50)
    D (66)  B (89)  B (75)  F (61)  D (60)  C (98)  
    E (73)  E (51)  D (63)  G (41)  F (99)      
    F (79)  G (86)  G (51)

正しい方向へのアドバイスやポイントは大歓迎です。

ありがとう！

4

4 に答える 4

1

あなたが示している構造を見ると、かなり混乱します。代わりに、2 番目の行列の位置に値を代入するとどうなるでしょうか。

foo = structure(c(0L, 47L, 42L, 66L, 73L, 79L, 25L, 70L, 0L, 98L, 20L, 
0L, 34L, 26L, 5L, 20L, 0L, 18L, 63L, 61L, 41L, 73L, 89L, 30L, 
0L, 51L, 56L, 86L, 96L, 75L, 30L, 63L, 0L, 12L, 51L, 46L, 50L, 
22L, 18L, 23L, 0L, 30L, 58L, 19L, 76L, 60L, 7L, 99L, 0L), .Dim = c(7L, 
7L), .Dimnames = list(c("A", "B", "C", "D", "E", "F", "G"), c("A", 
"B", "C", "D", "E", "F", "G")))

bar = structure(c(0L, 0L, 1L, 1L, 1L, 1L, 0L, 1L, 0L, 1L, 0L, 0L, 0L, 
0L, 0L, 0L, 0L, 0L, 1L, 1L, 1L, 1L, 1L, 0L, 0L, 1L, 0L, 1L, 1L, 
1L, 0L, 1L, 0L, 0L, 1L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 
1L, 0L, 1L, 0L), .Dim = c(7L, 7L), .Dimnames = list(c("A", "B", 
"C", "D", "E", "F", "G"), c("A", "B", "C", "D", "E", "F", "G"
)))

bar[as.logical(bar)] = foo[as.logical(bar)]

> bar
   A  B  C  D  E  F  G
A  0 70  0 73 96  0  0
B  0  0  0 89 75 50  0
C 42 98  0  0  0  0 76
D 66  0  0  0 63  0 60
E 73  0 63 51  0  0  0
F 79  0 61  0  0  0 99
G  0  0 41 86 51  0  0

あなたが示している構造は として作成できますがlist、明確なデータ表示にはならないと思います。ただし、出力のニーズを自由に拡張または明確にしてください。

于 2013-01-02T22:26:08.873 に答える
1

これが私がやった方法です：

X <- read.table(stdin())
A   B   C   D   E   F   G
A   0   70  5   73  96  46  58
B   47  0   20  89  75  50  19
C   42  98  0   30  30  22  76
D   66  20  18  0   63  18  60
E   73  0   63  51  0   23  7
F   79  34  61  56  12  0   99
G   25  26  41  86  51  30  0



indices <- which(t(apply(X,1,function(x) x >= sort(x,decr=TRUE)[3])), arr.ind=TRUE)
values  <- X[indices]

letters <- matrix(rownames(X)[indices],ncol=2)


data.frame(letters, values)

まず、データをテーブルとして読み込みます。次に、各行の 3 番目に大きい要素に対して各要素をテストします。残りは出力を少し良くするだけです:

   X1 X2 values
1   C  A     42
2   D  A     66
3   E  A     73
4   F  A     79
5   A  B     70
6   C  B     98
7   E  C     63
8   F  C     61
9   G  C     41
10  A  D     73
11  B  D     89
12  E  D     51
13  G  D     86
14  A  E     96
15  B  E     75
16  D  E     63
17  G  E     51
18  B  F     50
19  C  G     76
20  D  G     60
21  F  G     99
于 2013-01-02T23:24:11.553 に答える
0 
inter <- 0 + t( apply(foo, 1, rank) >= 5)  # don't actually need it but shows how.

vals <- foo * ( 0+t(apply(foo, 1, rank)>=5))
#----------
   A  B  C  D  E  F  G
A  0 70  0 73 96  0  0
B  0  0  0 89 75 50  0
C 42 98  0  0  0  0 76
D 66  0  0  0 63  0 60
E 73  0 63 51  0  0  0
F 79  0 61  0  0  0 99
G  0  0 41 86 51  0  0

vals [ , rev( order( colSums(vals >0) )) ]  
 # Arrange in order of increasing numbers of top three
#-------
   E  D  A  G  C  B  F
A 96 73  0  0  0 70  0
B 75 89  0  0  0  0 50
C  0  0 42 76  0 98  0
D 63  0 66 60  0  0  0
E  0 51 73  0 63  0  0
F  0  0 79 99 61  0  0
G 51 86  0  0 41  0  0

A 列を最初にしたい場合は、次を使用できます vals [ , order( -colSums(vals >0) ) ]。そして、水平表示で解決できる場合：

> apply( vals2 , 2, function(col) paste0(rownames(foo)[col>0], " (", col[col>0], ")" ) )
$A
[1] "C (42)" "D (66)" "E (73)" "F (79)"

$D
[1] "A (73)" "B (89)" "E (51)" "G (86)"

$E
[1] "A (96)" "B (75)" "D (63)" "G (51)"

$C
[1] "E (63)" "F (61)" "G (41)"

$G
[1] "C (76)" "D (60)" "F (99)"

$B
[1] "A (70)" "C (98)"

$F
[1] "B (50)"
于 2013-01-02T23:58:26.317 に答える
0 
dat <- structure(c(0L, 47L, 42L, 66L, 73L, 79L, 25L, 70L, 0L, 98L, 20L, 
0L, 34L, 26L, 5L, 20L, 0L, 18L, 63L, 61L, 41L, 73L, 89L, 30L, 
0L, 51L, 56L, 86L, 96L, 75L, 30L, 63L, 0L, 12L, 51L, 46L, 50L, 
22L, 18L, 23L, 0L, 30L, 58L, 19L, 76L, 60L, 7L, 99L, 0L), .Dim = c(7L, 
7L), .Dimnames = list(c("A", "B", "C", "D", "E", "F", "G"), c("A", 
"B", "C", "D", "E", "F", "G")))

中間行列は、次のようにして計算できます。

in.top3 <- t(apply(-foo, 1, rank)) <= 3

次に、次のようにして、各列の上位 3 行のリストを作成できます。

apply(in.top3, 2, function(x)names(which(x)))
# $A
# [1] "C" "D" "E" "F"
# 
# $B
# [1] "A" "C"
# 
# $C
# [1] "E" "F" "G"
# 
# $D
# [1] "A" "B" "E" "G"
# 
# $E
# [1] "A" "B" "D" "G"
# 
# $F
# [1] "B"
# 
# $G
# [1] "C" "D" "F"
于 2013-01-02T23:24:07.497 に答える