• Home
  • About R Trix
  • Teaching Resources
  • RTricks Library
  • Interactive Apps
  • Training Network
RTFR landingdist

R Output for Discrete Probability Example 6

> N <- 1000

> result <- rep(0, N)

> for (i in 1:N) result[i] <- sample(1:6, 1) + sample(1:6, 1)

> table(result)

23456789101112
28 58 74 120 132 181 143 107 85 49 23

> prob1 <- table(result)/N

> prob1

23456789101112
0.028 0.058 0.074 0.120 0.132 0.181 0.143 0.107 0.085 0.049 0.023

> dice2 <- outer(1:6, 1:6, "+")

> rownames(dice2) <- 1:6

> colnames(dice2) <- 1:6

> dice2

123456
1 2 3 4 5 6 7
2 3 4 5 6 7 8
3 4 5 6 7 8 9
4 5 6 7 8 9 10
5 6 7 8 9 10 11
6 7 8 9 10 11 12

> prob2 <- table(dice2)/36

> prob2

23456789101112
0.028 0.056 0.083 0.111 0.139 0.167 0.139 0.111 0.083 0.056 0.028