R Ping Pong的模擬循環

我需要找到概率Pr（X = i），i = 2，...。。。，6，通過使用R進行模擬，當兩個參與者A和B同意一個遊戲的獲勝者將獲得1分而失敗者0分;比賽結束，因爲其中一名球員超前2分或比賽數目達到6.假設A和B贏得比賽的概率分別爲2 3 y 1 3，並且每場比賽是獨立的。讓X表示結束遊戲所需的遊戲數量。R Ping Pong的模擬循環

我申請了下面的代碼：

juegos<-rbinom(6,1,2/3) 
juegos 
A<-cumsum(juegos) 
B<-cumsum(1-juegos) 
K<-abs(A-B)==2 

R<-rep(0,1000) 
for(i in 1:1000) 
{R[i]<-which.max(K)} 
R

但是我不知道什麼是發現的概率i = 2時，4和6

來源

2017-09-03 Willi Medina

這裏是下一步一種使用函數來模擬單個匹配的方式：

# Function to simulate one match 
one_match = function(p = 2/3){ 
    g = 0 
    score = 0 
    while (g < 6){ 
    g = g + 1 
    # Play one game & update score 
    if (runif(1) < p) 
     score = score + 1 
    else 
     score = score - 1 
    if (abs(score) == 2) break 
    } 
    return(g) 
} 

# Simulate matches 
n_sims = 100000 
outcomes = replicate(n_sims, one_match()) 
# Or, with a different winning probability, say p = 1/2 
# outcomes = replicate(n_sims, one_match(p = 1/2)) 

# Estimate probabilities 
probs = table(outcomes)/n_sims 
print(probs)

乾杯！

來源

2017-09-03 05:17:58 DodoD

R Ping Pong的模擬循環

回答

相關問題