주사위 1개 던지기
- runif(n) #0부터 1까지의 확률을 랜덤으로 n개 뽑고 싶을 때 (r:랜덤/unif)
runif(30)
[1] 0.763735464 0.437181863 0.877027896 0.809027348 0.921417037
[6] 0.849998626 0.715711981 0.222345358 0.212121592 0.250103505
[11] 0.363687967 0.589080022 0.553905812 0.922251850 0.572914371
[16] 0.669978790 0.467075639 0.782836641 0.373551849 0.144249398
[21] 0.317554528 0.383775095 0.852490050 0.946503954 0.163439841
[26] 0.822465883 0.477048066 0.001041733 0.526789160 0.242372675
- ceiling() #올림
ceiling(runif(10000)*6)
[1] 2 3 6 5 3 3 6 6 5 1 3 2 6 1 2 2 3 4 6 4 4 6 6 4 4 1 5 3 4 5 5 6 4
[34] 6 6 4 4 6 6 4 1 4 3 5 1 1 1 2 4 3 6 4 4 2 1 5 3 4 2 1 3 2 4 5 2 3
[67] 2 1 2 3 2 1 1 5 4 2 2 6 2 4 5 2 2 3 5 5 4 4 4 3 1 5 3 1 1 4 1 1 1
[100] 2 5 4 6 1 5 3 6 2 1 2 4 2 2 5 4 5 1 6 1 3 4 6 3 5 3 5 6 2 3 4 5 3
[133] 2 2 5 5 4 4 6 1 2 2 5 4 3 3 4 1 6 1 3 6 4 1 3 6 2 2 2 4 1 6 5 3 1
[166] 2 5 2 2 5 2 1 6 4 6 4 6 4 3 1 2 2 4 5 4 1 1 3 4 5 6 6 3 1 5 1 1 6
...
- table(variable)
table(die)
1 2 3 4 5 6
1709 1660 1637 1689 1640 1665
#1이 1709번 나옴 2가 1660번 나옴 ...~
- 변수에 넣어서 사용하는 것을 생활화하자
roll <- 100
die <-ceiling(runif(roll)*6)
b <-table(die)
barplot(b) #막대그래프를 그릴 수 있다.
주사위 2개 던지기
roll <- 1000000 #백만번 던져보자.
dice <- ceiling(runif(roll)6) + ceiling(runif(roll)6)
b <- table(dice)
barplot(b)
2 3 4 5 6 7 8 9 10 11
27486 55412 83023 110626 138214 166875 138932 112028 83709 56057
12
27638
#정규분포 모양 비슷하게 나온다.
