Rock-Paper-Scissors
# This example illustrate how to find the binomial distribution of the Rock-Paper-Scissors example
numberTrial <- 12
parameter <- 1/3
probability <- function(n, x) {
factorial(n) / factorial(n-x) / factorial(x)*(parameter)^x*(1-parameter)^(n-x)
}
# create a dataframe for saving the probability of different number of scissors thrown out of 12 times under the null hypothesis
distvector <- vector('numeric',length = 13)
for (i in 0:12){
distvector[i+1] <- probability(12,i)
}
dis <- as.data.frame(cbind(seq(0,12),distvector))
dis
## V1 distvector
## 1 0 7.707347e-03
## 2 1 4.624408e-02
## 3 2 1.271712e-01
## 4 3 2.119520e-01
## 5 4 2.384460e-01
## 6 5 1.907568e-01
## 7 6 1.112748e-01
## 8 7 4.768921e-02
## 9 8 1.490288e-02
## 10 9 3.311751e-03
## 11 10 4.967626e-04
## 12 11 4.516023e-05
## 13 12 1.881676e-06
# Plot the distribution of throwing scissors
barplot(dis$distvector,ylim=c(0,0.3),names.arg = dis$V1)
