[K-MOOC 통계학의 이해 1] 5-4 확률계산 실습 및 퀴즈

#과제
#생일문제: 상황 Ⓑ에서 k=5, 10, 25, 50, 100일 때 통계적 확률을 추정하여라. 단 반복수는 100,000번
n <- 365
daily <- scan()
2000 2000 1500 1500 1500 1000 1000 1000 1500 1500 1500 2000

monthday <- c(31,28,31,30,31,30,31,31,30,31,30,31)
births <- rep(daily,monthday)
birthprob <- births/sum(births)

# k=50일 때 표본
k <- 50
x <- sample.int(n,k,replace=T,prob=birthprob) #1~365(생년월일을 일로 표현한 것) 숫자 중 무작위로 k만큼 복원추출해서 뽑기, 확률 반영
result <- length(unique(x)) == k
result

# k=50일 때 100000개 표본
noSim <- 100000
result <- 0
for (i in 1:noSim)
{
  x <- sample.int(n,k,replace=T,prob=birthprob)
  result <- result + (length(unique(x)) == k) #True이면 1
}
cat("생일이 모두 다른 경우:",result,"\n") #cat = print
cat("생일이 모두 다를 확률:",result/noSim,"\n")

 

 

 

 

*강의의 모든 내용을 적은게 아닌, 정리하며 눈에 보이는 것만 후딱 요약해 적은겁니다.

강의도 무료이니 여인권 교수님의 강의를 직접 들어보길 추천해요!

 

모든 내용 출처 : http://www.kmooc.kr/courses/course-v1:SookmyungK+SM_sta_004k+2019_03SM_02/course/#block-v1:SookmyungK+SM_sta_004k+2019_03SM_02+type@chapter+block@eb0a1c109e664458999982d10e934a63