El problema del árbol mágico del dinero


19

Pensé en este problema en la ducha, estaba inspirado en las estrategias de inversión.

Digamos que había un árbol de dinero mágico. Todos los días, puede ofrecer una cantidad de dinero al árbol de dinero y lo triplicará o lo destruirá con una probabilidad de 50/50. Inmediatamente nota que, en promedio, ganará dinero al hacer esto y está ansioso por aprovechar el árbol del dinero. Sin embargo, si ofreciera todo su dinero a la vez, tendría un 50% de perder todo su dinero. ¡Inaceptable! Eres una persona bastante reacia al riesgo, así que decides idear una estrategia. Desea minimizar las probabilidades de perderlo todo, ¡pero también desea ganar tanto dinero como pueda! Se le ocurre lo siguiente: todos los días, ofrece el 20% de su capital actual al árbol del dinero. Suponiendo que lo más bajo que puede ofrecer es 1 centavo, se necesitaría una racha de 31 pérdidas para perder todo su dinero si comenzara con 10 dólares. Y lo que es más, cuanto más efectivo ganes, más tiempo tendrá que ser la racha perdedora para que puedas perder todo, ¡increíble! Rápidamente comienza a ganar mucho dinero en efectivo. Pero entonces se te ocurre una idea: ¡puedes ofrecer el 30% cada día y ganar mucho más dinero! Pero espera, ¿por qué no ofrecer el 35%? 50%? Un día, con grandes signos de dólar en tus ojos, corres hacia el árbol del dinero con todos tus millones y ofreces hasta el 100% de tu efectivo, que el árbol del dinero quema rápidamente. Al día siguiente, obtienes un trabajo en McDonalds. que el árbol del dinero quema rápidamente. Al día siguiente, obtienes un trabajo en McDonalds. que el árbol del dinero quema rápidamente. Al día siguiente, obtienes un trabajo en McDonalds.

¿Existe un porcentaje óptimo de su efectivo que pueda ofrecer sin perderlo todo?

(sub-preguntas:

Si hay un porcentaje óptimo que debe ofrecer, ¿es estático (es decir, 20% todos los días) o el porcentaje debería crecer a medida que aumenta su capital?

Al ofrecer un 20% todos los días, ¿las probabilidades de perder todo su dinero disminuyen o aumentan con el tiempo? ¿Existe un porcentaje de dinero de donde las probabilidades de perder todo su dinero aumentan con el tiempo?


77
Esto parece una variación de la ruina del jugador
Robert Long

2
Gran parte de esta pregunta depende de si son posibles centavos fraccionales. Además, hay muchos objetivos posibles que alguien podría tener en esta situación. Diferentes objetivos tendrían diferentes estrategias óptimas.
Buge

Respuestas:


19

Este es un problema bien conocido. Se llama una apuesta de Kelly. La respuesta, por cierto, es 1/3. Es equivalente a maximizar la utilidad de registro de la riqueza.

Kelly comenzó tomando tiempo al infinito y luego resolviendo hacia atrás. Como siempre puede expresar retornos en términos de capitalización continua, también puede revertir el proceso y expresarlo en registros. Voy a usar la explicación de la utilidad de registro, pero la utilidad de registro es una conveniencia. Si está maximizando la riqueza como , terminará con una función que resulta ser la misma que la utilidad de registro. Si es la probabilidad de pago, y es la probabilidad de ganar, y es el porcentaje de la riqueza invertida, entonces la siguiente derivación funcionará.nbpX

Para una apuesta binaria, , para un período único y riqueza unitaria.E(log(X))=plog(1+bX)+(1p)log(1X)

ddXE[log(x)]=ddX[plog(1+bX)+(1p)log(1X)]
=pb1+bX1p1X

Poniendo la derivada a cero para encontrar los extremos,

psi1+siX-1-pag1-X=0 0

Multiplicación cruzada, terminas con

pb(1X)(1p)(1+bX)=0
pbpbX1bX+p+pbX=0
bX=pb1+p
X=bp(1p)b

En su caso,

X=3×12(112)3=13.

Puede expandir esto fácilmente a resultados múltiples o continuos resolviendo la utilidad esperada de la riqueza sobre una distribución de probabilidad conjunta, eligiendo las asignaciones y sujeto a cualquier restricción. Curiosamente, si lo realiza de esta manera, incluyendo restricciones, como la capacidad de cumplir con los pagos de la hipoteca, etc., ha contabilizado su conjunto total de riesgos y, por lo tanto, tiene un riesgo ajustado o al menos controlado. solución.

Desiderata El propósito real de la investigación original tenía que ver con cuánto apostar en función de una señal ruidosa. En el caso específico, cuánto apostar en una señal electrónica ruidosa donde indicaba el lanzamiento de armas nucleares por parte de la Unión Soviética. Ha habido varios lanzamientos cercanos tanto de Estados Unidos como de Rusia, obviamente por error. ¿Cuánto apuestas por una señal?


Esta estrategia daría un mayor riesgo de quiebra, creo que en comparación con las fracciones más bajas
probabilidad es

@probabilityislogic Solo en el caso de que existan centavos. En el caso discreto, se volvería cierto porque podría apostar su último centavo. No podías apostar un tercio de un centavo. En un mundo discreto, es intrínsecamente cierto que la probabilidad de quiebra debe estar aumentando en el tamaño de la asignación, independientemente del caso de pago. Una asignación del 2% tiene una mayor probabilidad de quiebra que un 1% en un mundo discreto.
Dave Harris

@probabilityislogic si comienza con 3 centavos, entonces es arriesgado. Si comienza con $ 550, entonces hay menos de una posibilidad en 1024 de quebrar. Para tamaños de maceta razonables, el riesgo de colapso discreto se vuelve pequeño a menos que realmente vaya al infinito, luego se convierte en una certeza a menos que se permitan préstamos.
Dave Harris

Esperaba que este fuera un problema conocido, pero no tenía idea de cómo buscarlo. Gracias por la mención de Kelly. Sin embargo, una pregunta: wikipedia sobre el criterio de Kelly menciona la siguiente fórmula para calcular el porcentaje óptimo: (bp-q) / b. Donde b es el #dólar que obtienes al apostar 1 $, p la probabilidad de ganar y q la posibilidad de perder. Si completo esto para mi escenario, obtengo: (2 * 0.5-0.5) /2=0.25. Lo que significa que el porcentaje óptimo para apostar sería del 25%. ¿Qué causa esta discrepancia con su respuesta de 1/3?
ElectronicToothpick

3
@ElectronicToothpick si completas b = 3 obtienes 1/3. La diferencia está en cómo considera el pago tres veces. Digamos que comienza con 1 dólar y apuesta 50 centavos, luego considera que el pago triple termina con cincuenta y cincuenta 50 centavos o 2 dólares (b = 2, es decir, menos 50 centavos o más 2 veces 50 centavos) frente a cincuenta y cincuenta 50 centavos o 2.50 dólares (b = 3, es decir, menos 50 centavos o más 3 veces 50 centavos).
Sextus Empiricus

5

Me gustó la respuesta dada por Dave harris. aunque trataría el problema desde una perspectiva de "bajo riesgo", en lugar de maximizar el beneficio

La caminata aleatoria que está haciendo, suponiendo que su apuesta de fracción es y la probabilidad de ganar tiene como donde . en promedio tiene Puede aplicar esto iterativamente para obtener con el valor esperado también puede expresar la cantidad en el tiempo en función de una sola variable aleatoria , pero observando que no es independiente deqp=0.5

Yt|Yt1=(1q+3qXt)Yt1
XtBernoulli(p)
E(Yt|Yt1)=(1q+3pq)Yt1
Yt|Y0=Y0j=1t(1q+3qXt)
E(Yt|Y0)=(1q+3pq)tY0
tZt=j=1tXtBinomial(t,p)ZtZt1
Yt|Y0=Y0(1+2q)Zt(1q)tZt

posible estrategia

podría usar esta fórmula para determinar un valor de "bajo riesgo" para . Supongamos que quiere asegurarse de que después de pérdidas consecutivas todavía tenga la mitad de su riqueza original. Luego configuraqkq=12k1

Tomando el ejemplo significa que establecemos , o con establecemos .k=5q=0.129k=15q=0.045

Además, debido a la naturaleza recursiva de la estrategia, este riesgo es lo que está tomando en cada apuesta. Es decir, en el momento , al continuar jugando, se asegura de que en el tiempo su riqueza sea al menossk+s0.5Ys

discusión

La estrategia anterior no depende de la recompensa de ganar, sino más bien de establecer un límite en la pérdida. Podemos obtener las ganancias esperadas al sustituir el valor de calculamos, y en el momento que se utilizó con el riesgo en mente.qk

sin embargo, es interesante observar la mediana en lugar de la rentabilidad esperada en el tiempo , que se puede encontrar asumiendo la . cuando tenemos la razón igual a . Esto se maximiza cuando y mayor que cuandotmedian(Zt)tp

Yk|Y0=Y0(1+2q)tp(1q)t(1p)
p=0.5(1+q2q2)0.5tq=0.251q<0.5

También es interesante calcular la probabilidad de adelantarse en el momento . para hacer esto necesitamos determinar el valor tal que reorganizando algo, encontramos que la proporción de victorias debería satisfacer Esto se puede conectar a una aproximación normal (nota: media de y error estándar de ) como tz

(1+2q)z(1q)tz>1
zt>log(1q)log(1q)log(1+2q)
0.50.5t
Pr(ahead at time t)Φ(tlog(1+2q)+log(1q)[log(1+2q)log(1q)])

que muestra claramente que el juego tiene muy buenas probabilidades. el factor multiplicador se minimiza cuando (valor maximizado de ) y disminuye monotónicamente en función de . entonces la estrategia de "bajo riesgo" es apostar una fracción muy pequeña de su riqueza y jugar una gran cantidad de veces.tq=013q

supongamos que comparamos esto con y . El factor para cada caso es y . Esto significa que después de juegos tendrías alrededor de un 95% de posibilidades de adelantarte con la apuesta pequeña, en comparación con un 75% de posibilidades con la apuesta más grande. Además, también tiene la posibilidad de ir a la quiebra con la apuesta más grande, suponiendo que tenga que redondear su apuesta a los 5 centavos o dólares más cercanos. Comenzando con esto podría ir . Esta es una secuencia de pérdidas de , y dado que el juego esperaríaq=13q=11000.110.32382013.35,8.90,5.95,3.95,2.65,1.75,1.15,0.75,0.50,0.35,0.25,0.15,0.1,0.05,0143819pérdidas, si tiene mala suerte con las primeras apuestas, incluso ganar puede no compensar una mala racha (por ejemplo, si la mayoría de sus ganancias se producen una vez que la mayor parte de la riqueza se ha ido) ir a la quiebra con la menor participación del 1% no es posible en juegos. La otra cara es que la apuesta más pequeña dará como resultado una ganancia mucho menor en promedio, algo así como un aumento de veces con la apuesta grande en comparación con un aumento de con la apuesta pequeña (es decir, espera tener 24 dólares después de 38 rondas con la apuesta pequeña). apuesta y 7000 dolares con la apuesta grande).383501.2


que es si tenemos en cuenta que se elige de una manera de bajo riesgo y que no está calculando que para t > > k , no se trata tan mal una aproximación. Por lo tanto, probablemente exagere el beneficio de la gran estrategia de apuestas. qt>>k
probabilidadislogic

Su enfoque para maximizar la mediana de es en realidad el mismo enfoque de Dave Harris que maximiza la media de Z t (que es la misma que la mediana de Z t ). Sería diferente cuando uno está maximizando la media de Y t que está distribuida de forma lognormal y para la cual la media y la mediana no son iguales. ZtZtZtYt
Sextus Empiricus

5

No creo que esto sea muy diferente de la Martingala. En su caso, no hay apuestas duplicadas, pero el pago ganador es 3x.

Codifiqué una "réplica viva" de tu árbol. Corro 10 simulaciones. En cada simulación (traza), comienza con 200 monedas e intenta con el árbol, 1 moneda cada 20,000 veces.

Las únicas condiciones que detienen la simulación son la bancarrota o haber "sobrevivido" 20k intentos

enter image description here

Creo que cualesquiera que sean las probabilidades, tarde o temprano la bancarrota te espera.


El código es JavaScript improvisado pero sin dependencia: https://repl.it/@cilofrapez/MagicTree-Roulette

Le muestra los resultados de inmediato. El código es fácil de modificar: ejecutar cuantas simulaciones, cantidad de apuesta, cuantos intentos ... ¡No dude en jugar!

En la parte inferior del código, los resultados de cada simulación (por defecto 10) se guardan en un archivo CSV con dos columnas: número de giro y dinero. Lo hice para poder enviarlo a un trazador en línea para los gráficos.

Sería fácil tenerlo todo automatizado localmente usando la biblioteca de Gráficos de Google, por ejemplo. Si solo desea ver los resultados en la pantalla, puede comentar esa última parte como mencioné en el archivo.

EDITAR

Código fuente:

/**
 * License: MIT
 * Author: Carles Alcolea, 2019
 * Usage: I recommend using an online solution like repl.it to run this code.
 * Nonetheless, having node installed, it's as easy as running `node magicTree.js`.
 *
 * The code will run `simulations` number of scenarios, each scenario is equal in settings
 * which are self-descriptive: `betAmount`,`timesWinPayout`, `spinsPerSimulation`, `startingBankRoll`
 * and `winningOdds`.
 *
 * At the end of the code there's a part that will generate a *.csv file for each simulation run.
 * This is useful for ploting the resulting data using any such service or graphing library. If you
 * wish the code to generate the files for you, just set `saveResultsCSV` to true. All files will
 * have two columns: number of spin and current bankroll.
 */

const fs = require('fs'); // Only necessary if `saveResultsCSV` is true

/**
 * ==================================
 * You can play with the numbers of the following variables all you want:
 */
const betAmount          = 0.4,   // Percentage of bankroll that is offered to the tree
      winningOdds        = 0.5,
      startingBankRoll   = 200,
      timesWinPayout     = 2,
      simulations        = 5,
      spinsPerSimulation = 20000,
      saveResultsCSV     = false;
/**
 * ==================================
 */

const simWins = [];
let currentSim = 1;

//* Each simulation:
while (currentSim <= simulations) {
  let currentBankRoll = startingBankRoll,
      spin            = 0;
  const resultsArr  = [],
        progressArr = [];

  //* Each spin/bet:
  while (currentBankRoll > 0 && spin < spinsPerSimulation) {
    if (currentBankRoll === Infinity) break; // Can't hold more cash!
    let currentBet = Math.ceil(betAmount * currentBankRoll);
    if (currentBet > currentBankRoll) break;  // Can't afford more bets... bankrupt!

    const treeDecision = Math.random() < winningOdds;
    resultsArr.push(treeDecision);
    if (treeDecision) currentBankRoll += currentBet * timesWinPayout; else currentBankRoll -= currentBet;
    progressArr.push(currentBankRoll);
    spin++;
  }

  const wins = resultsArr.filter(el => el === true).length;
  const losses = resultsArr.filter(el => el === false).length;
  const didTheBankRollHold = (resultsArr.length === spinsPerSimulation) || currentBankRoll === Infinity;

  const progressPercent = didTheBankRollHold ? `(100%)` : `(Bankrupt at aprox ${((resultsArr.length / parseFloat(spinsPerSimulation)) * 100).toPrecision(4)}% progress)`;

  // Current simulation summary
  console.log(`
  - Simulation ${currentSim}: ${progressPercent === '(100%)' ? '✔' : '✘︎'}
    Total:      ${spin} spins out of ${spinsPerSimulation} ${progressPercent}
    Wins:       ${wins} (aprox ${((wins / parseFloat(resultsArr.length)) * 100).toPrecision(4)}%)
    Losses:     ${losses} (aprox ${((losses / parseFloat(resultsArr.length)) * 100).toPrecision(4)}%)
    Bankroll:   ${currentBankRoll}
  `);

  if (didTheBankRollHold) simWins.push(1);

  /**
   * ==================================
   * Saving data?
   */
  if (saveResultsCSV) {
    let data = `spinNumber, bankRoll`;
    if (!fs.existsSync('CSVresults')) fs.mkdirSync('CSVresults');
    progressArr.forEach((el, i) => {
      data += `\n${i + 1}, ${el}`;
    });
    fs.writeFileSync(`./CSVresults/results${currentSim}.csv`, data);
  }
  /**
   * ==================================
   */

  currentSim++;
}

// Total summary
console.log(`We ran ${simulations} simulations, with the goal of ${spinsPerSimulation} spins in each one.
Our bankroll (${startingBankRoll}) has survived ${simWins.length} out of ${simulations} simulations, with ${(1 - winningOdds) * 100}% chance of winning.`);
```

1
¿Puede publicar el código que escribió para esto también, por favor?
baxx

1
Esto es apostar con una apuesta constante, pero apostando una proporción fija de su riqueza, como aquí, cada vez produciría un resultado diferente. Es posible que deba adaptar esto para evitar monedas fraccionarias (por ejemplo, redondear hacia abajo a menos que esto produzca un valor inferior a1moneda, en cuyo caso apostará1moneda)1411
Henry

@baxx Claro, acabo de actualizar la publicación. Henry, no estoy seguro de entenderte. Puedo adaptar el código para adaptarse a diferentes necesidades si lo desea.
Carles Alcolea

@CarlesAlcolea Solo decía que sería bueno que el código que usaste para la publicación estuviera contenido en la publicación misma. No estoy seguro de si el enlace a las respuestas que ha publicado morirá en algún momento o no
baxx

1
@baxx ¡Claro! Después de escribir este programa improvisado, pensé que debería hacer una pequeña aplicación en línea, para poder explorar casi cualquier situación de este tipo fácilmente. No encontré ninguno. Ahora estoy ahogado en el trabajo, así que por el momento dejo el código en la publicación y la aplicación en mi lista de tareas pendientes
Carles Alcolea

4

Planteamiento del problema

Yt=log10(Mt)Mtt

q

Y0=1YL=2YWYW

Caminata aleatoria

Yt

Yt=Y0+i=1tXi

dónde

P[Xi=aw=log(1+2q)]=P[Xi=al=log(1q)]=12

Probabilidad de quiebra

Martingala

La expresion

Zt=cYt

c

caw+cal=2
c<1q<0.5

E[Zt+1]=E[Zt]12caw+E[Zt]12cal=E[Zt]

Probabilidad de terminar en bancarrota

Yt<YLYt>YWYWYLaw

E[Zτ]τE[Z0]

Así

cY0=E[Z0]=E[Zτ]P[Yτ<L]cYL+(1P[Yτ<L])cYW

y

P[Yτ<YL]cY0cYWcYLcYW

YW

P[Yτ<YL]cY0YL

Conclusiones

¿Existe un porcentaje óptimo de su efectivo que pueda ofrecer sin perderlo todo?

Cualquiera que sea el porcentaje óptimo dependerá de cómo valore las diferentes ganancias. Sin embargo, podemos decir algo sobre la probabilidad de perderlo todo.

Solo cuando el jugador está apostando una fracción cero de su dinero, ciertamente no irá a la quiebra.

qqgambler's ruin

qgambler's ruin=11/b
cawal

b=2

¿Las probabilidades de perder todo su dinero disminuyen o aumentan con el tiempo?

La probabilidad de ir a la quiebra depende de la distancia desde la cantidad de dinero donde el jugador quiebra. Cuandoq<qruina del jugador el dinero del jugador aumentará, en promedio, y la probabilidad de quiebra disminuirá, en promedio.

Probabilidad de quiebra cuando se utiliza el criterio de Kelly.

Cuando utiliza el criterio de Kelly mencionado en la respuesta de Dave Harris, q=0,5(1-1/ /si), para si siendo la relación entre pérdida y ganancia en una sola apuesta, entonces independiente de si El valor de C será igual a 0.1 y la probabilidad de quebrar será 0.1S-L.

Es decir, independiente del parámetro de asimetría. sidel árbol mágico, la probabilidad de quiebra, cuando se utiliza el criterio de Kelly, es igual a la proporción de la cantidad de dinero donde el jugador quiebra y la cantidad de dinero con la que el jugador comienza. Por diez dólares y 1 centavo, esta es una probabilidad de 1: 1000 de ir a la quiebra, cuando se utiliza el criterio de Kelly.

Simulaciones

Las simulaciones a continuación muestran diferentes trayectorias simuladas para diferentes estrategias de juego. Las trayectorias rojas son las que terminaron en bancarrota (golpear la líneaYt=-2)

simulations

Distribución de beneficios después del tiempo. t

Para ilustrar mejor los posibles resultados de las apuestas con el árbol del dinero, puede modelar la distribución de Ytcomo un proceso de difusión unidimensional en un campo de fuerza homogéneo y con un límite absorbente (donde el jugador se arruina). La solución para esta situación ha sido dada por Smoluchowski

Smoluchowski, Marian V. "Über Brownsche Molekularbewegung unter Einwirkung äußerer Kräfte und deren Zusammenhang mit der verallgemeinerten Diffusionsgleichung". Annalen der Physik 353.24 (1916): 1103-1112. (disponible en línea a través de: https://www.physik.uni-augsburg.de/theo1/hanggi/History/BM-History.html )

Ecuación 8:

W(x0,x,t)=ec(xx0)2Dc2t4D2πDt[e(xx0)24Dte(x+x0)24Dt]

This diffusion equation relates to the tree problem when we set the speed c equal to the expected increase E[Yt], we set D equal to the variance of the change in a single steps Var(Xt), x0 is the initial amount of money, and t is the number of steps.

The image and code below demonstrate the equation:

  • The histogram shows the result from a simulation.

  • The dotted line shows a model when we use a naive normal distribution to approximate the distribution (this corresponds to the absence of the absorbing 'bankruptcy' barrier). This is wrong because some of the results above the bankruptcy level involve trajectories that have passed the bankruptcy level at an earlier time.

  • The continuous line is the approximation using the formula by Smoluchowski.

illustration as diffusion in force field

Codes

#
## Simulations of random walks and bankruptcy:
#

# functions to compute c
cx = function(c,x) {
  c^log(1-x,10)+c^log(1+2*x,10) - 2
}
findc = function(x) {
  r <- uniroot(cx, c(0,1-0.1^10),x=x,tol=10^-130)
  r$root
}


# settings
set.seed(1)
n <- 100000
n2 <- 1000
q <- 0.45

# repeating different betting strategies
for (q in c(0.35,0.4,0.45)) {
  # plot empty canvas
  plot(1,-1000,
       xlim=c(0,n2),ylim=c(-2,50),
       type="l",
       xlab = "time step", ylab = expression(log[10](M[t])) )

  # steps in the logarithm of the money
  steps <- c(log(1+2*q,10),log(1-q,10))

  # counter for number of bankrupts
  bank <- 0

  # computing 1000 times
  for (i in 1:1000) {
    # sampling wins or looses
    X_t <- sample(steps, n, replace = TRUE)
    # compute log of money
    Y_t <- 1+cumsum(X_t)
    # compute money
    M_t <- 10^Y_t
    # optional stopping (bankruptcy)
    tau <- min(c(n,which(-2 > Y_t)))
    if (tau<n) {
      bank <- bank+1
    }
    # plot only 100 to prevent clutter
    if (i<=100) {
      col=rgb(tau<n,0,0,0.5)
      lines(1:tau,Y_t[1:tau],col=col)
    }
  }
  text(0,45,paste0(bank, " bankruptcies out of 1000 \n", "theoretic bankruptcy rate is ", round(findc(q)^3,4)),cex=1,pos=4)
  title(paste0("betting a fraction ", round(q,2)))
}

#
## Simulation of histogram of profits/results
#

# settings
set.seed(1)
rep <- 10000  # repetitions for histogram
n   <- 5000   # time steps
q   <- 0.45    # betting fraction
b   <- 2      # betting ratio loss/profit
x0  <- 3      # starting money

# steps in the logarithm of the money
steps <- c(log(1+b*q,10),log(1-q,10))

# to prevent Moiré pattern in
# set binsize to discrete differences in results
binsize <- 2*(steps[1]-steps[2]) 

for (n in c(200,500,1000)) {

  # computing several trials
  pays <- rep(0,rep)
  for (i in 1:rep) {
    # sampling wins or looses
    X_t <- sample(steps, n, replace = TRUE)
      # you could also make steps according to a normal distribution
      # this will give a smoother histogram
      # to do this uncomment the line below
    # X_t <- rnorm(n,mean(steps),sqrt(0.25*(steps[1]-steps[2])^2))

    # compute log of money
    Y_t <- x0+cumsum(X_t)
    # compute money
    M_t <- 10^Y_t
    # optional stopping (bankruptcy)
    tau <- min(c(n,which(Y_t < 0)))
    if (tau<n) {
      Y_t[n] <- 0
      M_t[n] <- 0
    }
    pays[i] <- Y_t[n]
  }

  # histogram
  h <- hist(pays[pays>0],
            breaks = seq(0,round(2+max(pays)),binsize), 
            col=rgb(0,0,0,0.5),
            ylim=c(0,1200),
            xlab = "log(result)", ylab = "counts",
            main = "")
  title(paste0("after ", n ," steps"),line = 0)  

  # regular diffusion in a force field (shifted normal distribution)
  x <- h$mids
  mu <- x0+n*mean(steps)
  sig <- sqrt(n*0.25*(steps[1]-steps[2])^2)
  lines(x,rep*binsize*(dnorm(x,mu,sig)), lty=2)

  # diffusion using the solution by Smoluchowski
  #   which accounts for absorption
  lines(x,rep*binsize*Smoluchowski(x,x0,0.25*(steps[1]-steps[2])^2,mean(steps),n))

}

"That is, independent from the assymetry parameter b of the magic tree, the probability to go bankrupt, when using the Kelly criterion, is equal to the ratio of the amount of money where the gambler goes bankrupt and the amount of money that the gambler starts with. For ten dollars and 1 cent this is a 1:1000 probability to go bankrupt" Im a bit surprised about this. So this means the probability to go bankrupt will be 1:1000 even if the payout is 10 times the offered money per round? How is this possible when the odds of going bankrupt decrease as your money grows?
ElectronicToothpick

1
@ElectronicToothpick If the payout is larger, and if you do not change the fraction that you gamble, then the probability to go bankrupt will be smaller. However, when you increase the fraction that you gamble, then this may not be true anymore. With the Kelly criterion, you will increase the fraction to gamble when the payout is higher. This will increase the expected value of the logarithm of the money, but as a consequence, the probability to go bankrupt will remain the same.
Sextus Empiricus

1
Actually, when the gambler is not using the Kelly criterion, which optimizes E[logMt], but instead chooses to optimize E[Mt], then the consequence is that a higher fraction of the amount of money is being gambled. Possibly this might lead to an increase in the risk of bankruptcy when the payout is made larger. I could add an analysis of this, but I am afraid that my answer is already too long and/or complex.
Sextus Empiricus
Al usar nuestro sitio, usted reconoce que ha leído y comprende nuestra Política de Cookies y Política de Privacidad.
Licensed under cc by-sa 3.0 with attribution required.