可以先自定义函数
也可以用的时候再定义。
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
|
> mat <- matrix(c(1:3,7:9,4:6), byrow = T, nc = 3) > mat [,1] [,2] [,3] [1,] 1 2 3 [2,] 7 8 9 [3,] 4 5 6 > apply(mat, 2, function (x){order(x, decreasing=T)[1]}) # 查找每一列 [1] 2 2 2 > apply(mat, 1, function (x){order(x, decreasing=T)[1]}) # 查找每一行 [1] 3 3 3 > apply(mat, 1, function (x){ which .max(x)}) # 查找每一行 [1] 3 3 3 > n <- letters[1:5] > n [1] "a" "b" "c" "d" "e" > t <- apply(mat, 1, function (x){ which .max(x)}) > n[t] [1] "c" "c" "c" |
另一个例子:
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
|
MaxVar <- function (x, na. rm = FALSE) { ## compute `max` maxx <- max(x, na. rm = na. rm ) ## which equal the max wmax <- which (x == max(x)) ## how many equal the max nmax <- length(wmax) ## return out <- if (nmax > 1L) { c(999, NA) } else { c(maxx, wmax) } out } And use it like this: > new <- apply(Mydata[, -1], 1, MaxVar) > new [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [1,] 3 4 999 999 999 4 4 2 4 999 [2,] 1 4 NA NA NA 4 2 3 4 NA > Mydata <- cbind(Mydata, Max = new[1, ], Var = new[2, ]) > Mydata ID X1 X2 X3 X4 Max Var 1 1 3 1 1 1 3 1 2 2 1 2 1 4 4 4 3 3 1 1 1 1 999 NA 4 4 1 3 3 1 999 NA 5 5 2 2 2 1 999 NA 6 6 1 2 3 4 4 4 7 7 2 4 3 3 4 2 8 8 1 1 2 1 2 3 9 9 3 2 1 4 4 4 10 10 4 4 4 4 999 NA |
补充:使用R语言得到向量中所有的最大值或者最小值的下标
比如向量a:
1
|
a=c(1,2,4,3,4) |
代码如下:
1
|
which (a==max(a),arr.ind=TRUE) |
输出为:
[1] 3 5
以上为个人经验,希望能给大家一个参考,也希望大家多多支持服务器之家。如有错误或未考虑完全的地方,望不吝赐教。
原文链接:https://blog.51cto.com/matrix6ro/1899465