Apuesta

n=10
gana= 50
pierde= 25
prob=0.5
resultados <- rbinom(n=n, size=1, prob=prob); resultados
##  [1] 1 0 1 1 0 0 1 1 0 1
gan <- sum(resultados) * gana
saldo <-  gan-((n-sum(resultados))*pierde )

Pr =6/10 = 0.6

gan
## [1] 300
saldo
## [1] 200

Si 100 encuestados aceptan lanzar la moneda 10 veces ##########

ns=100
ind <- rbinom(n=ns, size=n, prob=prob)
gan <- gana*ind
saldo <- gan - (pierde*( n-ind ))
head(cbind(ind/n,saldo))
##          saldo
## [1,] 0.7   275
## [2,] 0.6   200
## [3,] 0.5   125
## [4,] 0.7   275
## [5,] 0.5   125
## [6,] 0.3   -25
hist(saldo, 10); abline(v=0, col=2, lty="dashed")

pierden <- sum(as.numeric(saldo<=0))
cbind(pierden,"Pr. que pierden"=(pierden/ns)) 
##      pierden Pr. que pierden
## [1,]      13            0.13