R中利用apply、tapply、lapply、sapply、mapply

來源:互聯網
上載者:User

標籤:

  • apply函數(對一個數組按行或者按列進行計算):
使用格式為:
apply(X, MARGIN, FUN, ...)其中X為一個數組;MARGIN為一個向量(表示要將函數FUN應用到X的行還是列),若為1表示取行,為2表示取列,為c(1,2)表示行、列都計算。
範例程式碼:
> ma <- matrix(c(1:4, 1, 6:8), nrow = 2)
> ma
     [,1] [,2] [,3] [,4]
[1,]    1    3    1    7
[2,]    2    4    6    8
> apply(ma, c(1,2), sum)
     [,1] [,2] [,3] [,4]
[1,]    1    3    1    7
[2,]    2    4    6    8
> apply(ma, 1, sum)
[1] 12 20
> apply(ma, 2, sum)
[1]  3  7  7 15
  • 函數tapply(進行分組統計):
使用格式為:
tapply(X, INDEX, FUN = NULL, ..., simplify = TRUE)其中X通常是一向量;INDEX是一個list對象,且該list中的每一個元素都是與X有同樣長度的因子;FUN是需要計算的函數;simplify是邏輯變數,若取值為TRUE(預設值),且函數FUN的計算結果總是為一個標量值,那麼函數tapply返回一個數組;若取值為FALSE,則函數tapply的傳回值為一個list對象。需要注意的是,當第二個參數INDEX不是因子時,函數 tapply() 同樣有效,因為必要時 R 會用 as.factor()把參數強制轉換成因子。
範例程式碼:
> fac <- factor(rep(1:3, length = 17), levels = 1:5)
> fac
 [1] 1 2 3 1 2 3 1 2 3 1 2 3 1 2 3 1 2
Levels: 1 2 3 4 5
> tapply(1:17, fac, sum)
 1  2  3  4  5
51 57 45 NA NA
> tapply(1:17, fac, sum, simplify = FALSE)
$`1`
[1] 51

$`2`
[1] 57

$`3`
[1] 45

$`4`
NULL

$`5`
NULL
> tapply(1:17, fac, range)
$`1`
[1]  1 16

$`2`
[1]  2 17

$`3`
[1]  3 15

$`4`
NULL

$`5`
NULL
#利用tapply實作類別似於excel裡的樞紐分析表的功能:
> da
   year province sale
1  2007        A    1
2  2007        B    2
3  2007        C    3
4  2007        D    4
5  2008        A    5
6  2008        C    6
7  2008        D    7
8  2009        B    8
9  2009        C    9
10 2009        D   10
> attach(da)
> tapply(sale,list(year,province))
 [1]  1  4  7 10  2  8 11  6  9 12
> tapply(sale,list(year,province),mean)
      A  B C  D
2007  1  2 3  4
2008  5 NA 6  7
2009 NA  8 9 10
  • 函數table(求因子出現的頻數):
使用格式為:
table(..., exclude = if (useNA == "no") c(NA, NaN), useNA = c("no",    "ifany", "always"), dnn = list.names(...), deparse.level = 1)其中參數exclude表示哪些因子不計算。
範例程式碼:
> d <- factor(rep(c("A","B","C"), 10), levels=c("A","B","C","D","E"))
> d
 [1] A B C A B C A B C A B C A B C A B C A B C A B C A B C A B C
Levels: A B C D E
> table(d)
d
 A  B  C  D  E
10 10 10  0  0
> table(d, exclude="B")
d
 A  C  D  E
10 10  0  0

 

  • 函數lapply與函數sapply:
lapply的使用格式為:
lapply(X, FUN, ...)lapply的傳回值是和一個和X有相同的長度的list對象,這個list對象中的每個元素是將函數FUN應用到X的每一個元素。其中X為List對象(該list的每個元素都是一個向量),其他類型的對象會被R通過函數as.list()自動轉換為list類型。函數sapply是函數lapply的一個特殊情形,對一些參數的值進行了一些限定,其使用格式為:
sapply(X, FUN,..., simplify = TRUE, USE.NAMES = TRUE)sapply(*, simplify = FALSE, USE.NAMES = FALSE) 和lapply(*)的傳回值是相同的。如果參數simplify=TRUE,則函數sapply的傳回值不是一個list,而是一個矩陣;若simplify=FALSE,則函數sapply的傳回值仍然是一個list。
範例程式碼:
> x <- list(a = 1:10, beta = exp(-3:3), logic = c(TRUE,FALSE,FALSE,TRUE))
> lapply(x, quantile)
$a
   0%   25%   50%   75%  100%
 1.00  3.25  5.50  7.75 10.00

$beta
         0%         25%         50%         75%        100%
 0.04978707  0.25160736  1.00000000  5.05366896 20.08553692

$logic
  0%  25%  50%  75% 100%
 0.0  0.0  0.5  1.0  1.0

> sapply(x, quantile,simplify=FALSE,use.names=FALSE)
$a
   0%   25%   50%   75%  100%
 1.00  3.25  5.50  7.75 10.00

$beta
         0%         25%         50%         75%        100%
 0.04978707  0.25160736  1.00000000  5.05366896 20.08553692

$logic
  0%  25%  50%  75% 100%
 0.0  0.0  0.5  1.0  1.0
#參數simplify=TRUE的情況
> sapply(x, quantile)
         a        beta logic
0%    1.00  0.04978707   0.0
25%   3.25  0.25160736   0.0
50%   5.50  1.00000000   0.5
75%   7.75  5.05366896   1.0
100% 10.00 20.08553692   1.0
  • 函數mapply:
函數mapply是函數sapply的變形版,mapply 將函數 FUN 依次應用每一個參數的第一個元素、第二個元素、第三個元素上。函數mapply的使用格式如下:
mapply(FUN, ..., MoreArgs = NULL, SIMPLIFY = TRUE,USE.NAMES = TRUE)其中參數MoreArgs表示函數FUN的參數列表。範例程式碼:
> mapply(rep, times=1:4, x=4:1)
[[1]]
[1] 4

[[2]]
[1] 3 3

[[3]]
[1] 2 2 2

[[4]]
[1] 1 1 1 1

#直接使用函數rep的結果:
> rep(1:4,1:4)
 [1] 1 2 2 3 3 3 4 4 4 4

R中利用apply、tapply、lapply、sapply、mapply

聯繫我們

該頁面正文內容均來源於網絡整理,並不代表阿里雲官方的觀點,該頁面所提到的產品和服務也與阿里云無關,如果該頁面內容對您造成了困擾,歡迎寫郵件給我們,收到郵件我們將在5個工作日內處理。

如果您發現本社區中有涉嫌抄襲的內容,歡迎發送郵件至: info-contact@alibabacloud.com 進行舉報並提供相關證據,工作人員會在 5 個工作天內聯絡您,一經查實,本站將立刻刪除涉嫌侵權內容。

A Free Trial That Lets You Build Big!

Start building with 50+ products and up to 12 months usage for Elastic Compute Service

  • Sales Support

    1 on 1 presale consultation

  • After-Sales Support

    24/7 Technical Support 6 Free Tickets per Quarter Faster Response

  • Alibaba Cloud offers highly flexible support services tailored to meet your exact needs.