확률변수와 모수(parameter)의 관계

예제 코드

rbernoulli <- function(n, p){
    # 결과값을 받을 변수 result
    result <- rep(0, n)
    
    # 스탑와치 돌린 결과
    stop_watch <- runif(n)
    
    # 필터링을 이용한 1값을 부여
    result[stop_watch < p] <- 1
    
    return(result)
}

우리가 만든 함수를 실행을 시켜 봅시다! 1이 나오는 비율이 30%인 베르누이 확률변수에서 숫자 10개 뽑기.

# 베르누이 p = 0.3 숫자 10개 뽑기
rbernoulli(10, 0.3)

##  [1] 1 0 0 0 1 0 0 0 1 1