R：在不同的參數上使用lapply的概率

Question

我有以下問題：

我有一個長度爲n的數組X，其中包含一些正數。 我現在想要構造第二個數組Y，它應該包含來自伯努利分佈的隨機數，但每次成功的概率都是m（X [i]，theta）。

我可以做一個for循環，但我想知道是否有可能做到一行。

我想是這樣的：

X = sample(c(seq(20)), replace = TRUE) 
m<-function(x,theta) return(x*theta) 
Y = lapply(lapply(X, m, theta=2), rbinom, n=1, size=1) 

，但它不能在一個形式我想工作。

我是一個REAL R初學者，所以也許答案很簡單，但我想不出該怎麼尋找找到的東西，可以幫助我......

Answer 1

我翻譯你的代碼

X <- sample(c(seq(20)), replace = TRUE) 
m <- function(x,theta) return(x*theta) 
prob <- lapply(X, m, theta=2) 
Y = unlist(lapply(prob, rbinom, n=1, size=1)) 

目前我們所有的NA爲Y。但是，如果prob是合適的，那麼代碼沒有任何問題。例如，

prob <- 1:10/10 
Y = unlist(lapply(prob, rbinom, n=1, size=1)) 

給出有意義的結果。

請注意，沒有必要使用lapply來獲得prob。只需做

prob <- m(X, theta) 

很好，因爲"*"運算符是向量化的。