2010-11-30 19 views
1

St Petersburg paradox est un jeu de hasard où vous payez un montant fixe pour entrer dans la partie. Vous lancez une pièce à plusieurs reprises jusqu'à ce qu'une pile soit lancée. Votre gain est la somme de 1 à n de 2^n où n est le nombre de têtes avant les premières queues. Si cela n'a pas de sens essayez l'article wikipediarésoudre le paradoxe de St Peterburg en R

Je faisais un article sur la théorie de l'utilité attendue et j'écrivais sur le paradoxe de Saint-Pétersbourg et je pensais qu'il serait intéressant (mais pas pertinent pour mon article) d'essayer et de faire un carlo monte en R pour combien vous attendez à gagner après 10000 essais

Je veux essentiellement faire http://www.mathematik.com/Petersburg/Petersburg.html en R avec 10.000 essais

+2

Cela ressemble à un super projet. Pourquoi avez-vous poster ici? As-tu une question? – abelenky

+0

Il est étroitement lié aux enchères progressives dans un jeu de casino 50/50, où vous commencez avec une somme fixe X, si vous avez perdu votre offre X * 2, puis X * 4, etc. Si vous avez un tampon d'argent infini, vous va sûrement gagner ce X initial à la fin. – ruslik

+0

Le site lié ne fonctionne pas correctement, car si vous payez 20 $ chaque fois que vous ne gagnez pas. cela ne fonctionne que lorsque le salaire est d'environ 2^3. – John

Répondre

3

C'est facile dans R. le jeu suit la distribution géométrique avec p = 1/2:

N <- 1e+4 
out <- replicate(N, mean(2^rgeom(1000, .5))) 

Parce que le gain espéré du jeu est infini, vous obtiendrez une distribution empirique extrêmement biaisée que vous même pas être en mesure de représenter correctement:

hist(out) 

échelle logarithmique pourrait être une meilleure idée.

hist(log(out)) 
+0

Merci beaucoup vitoshka! Je savais que quelqu'un saurait comment le faire facilement. – Max