sapply in R, способ использования?

Question

Я программист на C++, и я новичок в R. Кто-то сказал мне, что использование цикла for в R - плохая идея и что лучше использовать sapply. Я написал следующий код, чтобы вычислить вероятность birthday coincidence:

prob <- 1   # prob of no coincidence 
days <- 365 
k <- 50    # how many people 
probability <- numeric() #probability vector (empty right now) 
for(i in 1:k){ 
    prob <- (days - i + 1)/days * prob # Formula for no coincidence 
    probability[i] <- 1 - prob 
} 

Как я могу сделать то же самое с sapply? Я хочу сделать что-то вроде:

1 - sapply(1:length(m), function(x) prod(m[1:x])) 

Но как использовать формулу для несовместимости с днем ​​рождения?

Answer 1

Вы можете сделать:

m <- (days - seq_len(k) + 1)/days 
probability <- 1 - sapply(seq_along(m), function(x) prod(m[1:x])) 

но будет отсутствовать на полезный cumprod функции:

probability <- 1 - cumprod(m) 

, который будет намного быстрее.

(Также дал вам пик seq_along и seq_len, которые являются более надежными, чем : при работе с векторами нулевой длины.)

Answer 2

Для вас конкретный вопрос, он, вероятно, лучше всего использовать встроенный в день рождения вероятностный калькулятор

sapply(1:50, pbirthday)