R：調用函數

Question

我有一個矩陣：

x <- cbind(x1 = 3, x2 = c(4:1, 2:5)) 
dimnames(x)[[1]] <- letters[1:8] 

如何以下代碼工作？

cave <- function(x, c1, c2) c(mean(x[c1]), mean(x[c2])) 
apply(x,1, cave, c1="x1", c2=c("x1","x2")) 

特別是我不理解的說法，c(mean(x[c1]), mean(x[c2]))功能cave內。

另外我們用這種方式調用函數cave(x,a,b)。但是在apply函數中，當我們調用cave函數時，輸入是哪裏？

Answer 1

當你定義一個函數時，列出的最後一件事情是隱式返回的。因此，當定義爲cave

cave <- function(x, c1, c2) { 
    c(mean(x[c1]), mean(x[c2])) 
} 

這是相同的

cave <- function(x, c1, c2) { 
    return(c(mean(x[c1]), mean(x[c2]))) 
} 

在R，矢量是通過使用串聯c()元件限定。因此c(a, b)使用元素a和b作出長度爲2的向量。在這種情況下，a是在由c1給出的指數處的輸入x的平均值，並且b是在由c2給出的指數處的x的平均值。

當你apply()此功能x在尺寸1，你對x每一行計算cave()。

的apply()的語法如下：

apply(x, # object input to the function 
     1, # dimension of x over which to apply the function, 1=rows, 2=columns 
     cave, # function to apply 
     c1="x1", c2=c("x1", "x2") # any further arguments to the function 
    ) 

應用cave()第一行給了我們由x1平均的（這只是x1因爲它是一個單數的平均值）的矢量和平均值爲x1和x2。 3的平均值是3，3和4的平均值是3.5，所以輸出是c(3, 3.5)。

我們重複，每行x，並最終我們得到這樣的：

 a b c d e f g h 
[1,] 3.0 3 3.0 3 3.0 3 3.0 3 
[2,] 3.5 3 2.5 2 2.5 3 3.5 4 

你可以看到的x等a一行列a對應cave()。