# install.packages ("RLAB") # Run First Time
требуют (RLAB) #used для функции распределения Бернулли rbern ()

#### Настройки
#Генеральная
live.plot = F #Watch Моделирование Live (работает медленнее)
max.bets<-1000 #NUMBER Ставка Конечного Simulation
итерации<-1000 #NUMBER имитаций для запуска
плата<-.0005 #Комиссия на перевод

Варианты #Bet Стратегия и Satoshi Dice
Win.Odds<-.5 #Odds выигрыша
Price.Multiplier<-1,957 #multiple Бет оплачена на Win
max.bet.size<-500 #Maximum допустимого размера Bet
start.wallet<-10 #Starting фонды
lose.multiplier<-2 #multiple потерять размер ставки использовать для следующей ставки
Настройки #End

## Calulate размеры ставок и выплаты
bet.num<-1: 20
bet.size<-.01 * lose.multiplier ^ (bet.num-1)
bet.size [который (bet.size>max.bet.size)]<-max.bet.size

WIN.SIZE<-Price.Multiplier * .01 * lose.multiplier ^ (bet.num-1) -fee
WIN.SIZE [который (WIN.SIZE>Price.Multiplier * max.bet.size)]<-Price.Multiplier * max.bet.size-плата



### Самая высокая плотность Interval калькулятор Функция (для графики)
get.HDI<-функции (sampleVec, credMass) {
  sortedPts = сортировки (sampleVec)
  ciIdxInc = этаж (credMass * длина (sortedPts))
  НСКР = длина (sortedPts) - ciIdxInc
  ciWidth = Rep (0, NCIS)
  для (я в 1: NCIS) {
    ciWidth [I] = sortedPts [г + ciIdxInc] - sortedPts [я]
  }
  HDImin = sortedPts [which.min (ciWidth)]
  HDImax = sortedPts [which.min (ciWidth) + ciIdxInc]
  HDIlim = с (HDImin, HDImax, credMass)
  вернуться (HDIlim)
}
####


###### Run Simulation
если (live.plot == Т) {
  dev.new ()
}

OUT2 = NULL,
пб<-txtProgressBar (мин = 0, макс = итерации, начальный = 0, стиль = 3)
для (J в 1: итераций) {
  
  бумажник<-start.wallet
  я<-1
  T<-1
  run.result = NULL,
  в то время как (кошелек>0 & T<= max.bets) {
    результат<-rbern (1, Win.Odds)
    если (результат == 0) {
      бумажник<-wallet-bet.size [я]
      цвет ="красный"
      я<-i + 1
    }
    
    если (результат == 1) {
      бумажник<-wallet-bet.size [я] + WIN.SIZE [я]
      цвет ="зеленый"
      я<-1
    }
    
    если (кошелек>0) {
      T<-t + 1
      run.result<-rbind (run.result, cbind (J, T, бумажник, я))
    }
    
    если (live.plot == Т) {
      участок (run.result [2], run.result [3],
           xlab ="Bet Номер", Ylab ="Монеты в Кошельке",
           цв = цвет,
           Основная = пасты ("Моделирование #", К)
      )
    }
  }
  
  out2<-rbind (OUT2, run.result)
  setTxtProgressBar (Pb, J)
}
закрыть (Рb)
#End Моделирование


#### Получить результаты
winning.iterations<-out2 [который (OUT2 [2] == max.bets), 1]
perc.wins<-100 * длина (которых (OUT2 [2] == max.bets)) / итераций
win.wallets<-out2 [который (OUT2 [2] == max.bets), 3]
HDI<-get.HDI (win.wallets, .95)

win.mode<-плотность (win.wallets) $ х [
    который (плотность (win.wallets) $ у ==
      макс (плотность (win.wallets) $ у)
    )]

fail.timepoints = матрица (nrow = (итерации длина (win.wallets)), Ncol = 1)
р<-1
для (я в 1: итераций) {
  температура<-out2 [который (OUT2 [1] == я),]
  если (макс (температура [2])    fail.timepoints [г]<-max (температура [2])
р<-r + 1
  }
}

#### Результаты Plot
dev.new ()
Схема (матрица (с (1,1,2,3), nrow = 2, Ncol = 2, byrow = 2))

Результаты моделирования #Plot
участок (0,0, тип ="N",
     xlab ="Bet Номер", Ylab ="Монеты в Кошельке",
     xlim = с (0, макс (OUT2 [2])), ylim = с (0, макс (OUT2 [3])),
     Основные = с (паста ("Начиная фонды =", Start.wallet, "  Win Odds =", Win.Odds),
            вставить("Bet Multiplier После потери = ", Lose.multiplier,"Икс", Сентябрь =""),
            вставить("Максимальные # Ставок =", Max.bets, "  # Имитаций =", итерация))
)
для (я в 1: итераций) {
  если (я% в% winning.iterations) {
    температура<-out2 [который (OUT2 [1] == я),]
    линии (температура [2], темп [3], Col = RGB (0,1,0, 0,1))
  } Еще {
    температура<-out2 [который (OUT2 [1] == я),]
    линии (температура [2], темп [3], Col = RGB (1,0,0, 0,1))
  }
}

#Plot Распределение ненулевых выигрыш
hist.counts<-hist (win.wallets, Col ="зеленый",
                  xlab ="Окончательный Кошелек Сумма",
                  перерывы = сл (0, (макс (win.wallets) + макс (win.wallets) / 10), с помощью = тах (win.wallets) / 10),
                  Основные = с (паста ("Процент Победа Симуляторы =", Perc.wins, "%"),
вставить("Mode =", Круглые (win.mode, 2), "  Среднее =", Круглые (среднее (win.wallets), 2))
),
                  суб = пасты (круглые (100 * HDI [3], 3), "%", "HDI (нижний, верхний):", Круглые (HDI [1], 1),",", Круглые (HDI [2], 1))
) $ Кол
Прямоугольник (HDI [1], 0, HDI [2] ,. 025 * макс (hist.counts), Col ="черный")

#Plot Распределение при банкротстве Происходило
HIST (fail.timepoints, Col ="красный",
     перерывы = сл (0, (макс (fail.timepoints) + макс (fail.timepoints) / 10), с помощью = тах (fail.timepoints) / 10),
     xlab ="Ставка Количество в банкротстве",
     Основные ="Ставка Количество в банкротстве")