Nilai yang diharapkan dari variabel acak pada pelemparan koin

8

Datang ke masalah yang menarik hari ini. Anda diberi koin dan uang x, Anda gandakan uang jika Anda mendapatkan kepala dan kehilangan setengah jika ekor pada setiap lemparan.

  1. Berapa nilai yang diharapkan dari uang Anda dalam percobaan
  2. Berapakah probabilitas mendapatkan lebih dari nilai yang diharapkan dalam (1)

Ini adalah bagaimana saya mendekatinya. Probabilitas kepala dan ekor sama (1/2). Nilai yang diharapkan setelah lemparan pertama = Jadi nilai yang diharapkan adalah setelah lemparan pertama. Demikian pula mengulangi harapan lemparan kedua pada 5x / 4, Nilai yang diharapkan setelah lemparan kedua =1/2(2x)+1/2(1/2x)=5x/45x/41/2(25x/4)+1/2(1/25x/4)=25x/16

Jadi Anda mendapatkan urutan nilai yang diharapkan: , , , ...5x/425x/16125x/64

Setelah mencoba, nilai yang Anda harapkan seharusnya .n(5n/4n)x

Jika cukup besar, nilai yang Anda harapkan harus mendekati rata-rata distribusi. Jadi probabilitas bahwa nilai lebih besar dari nilai yang diharapkan harus . Saya tidak yakin tentang yang ini.n0.5

pengguna862
sumber
1
Probabilitas untuk mengamati nilai di atas median adalah satu setengah, tetapi kuantitas Anda terdistribusi secara log-normal untuk n besar, yang memiliki rata-rata dan median yang berbeda, jadi Anda seharusnya tidak mengharapkan probabilitas untuk mendekati 1/2.
Yaroslav Bulatov

Jawaban:

9

Jika n cukup besar, nilai yang Anda harapkan harus mendekati rata-rata distribusi.

Ya itu benar.

Jadi probabilitas bahwa nilai lebih besar dari nilai yang diharapkan harus 0,5.

Ini hanya akan benar jika distribusinya simetris - yang dalam gim Anda tidak demikian. Anda dapat melihat ini dengan mudah jika Anda berpikir tentang apa nilai median dari kemenangan Anda setelahnyan melempar.


Anda dapat menganggap masalah Anda sebagai jalan acak . Jalan acak satu dimensi dasar adalah jalan pada garis nyata bilangan bulat, di mana pada setiap titik kita bergerak±1 dengan probabilitas p. Ini persis apa yang Anda miliki jika kita mengabaikan penggandaan / pembagian uang menjadi duap=0.5. Yang harus kita lakukan adalah memetakan ulang sistem koordinat Anda ke contoh ini. Membiarkanxmenjadi pot awal awal Anda. Lalu kami memetakan dengan cara berikut:

x*2^{-2} = -2
x*2^{-1} = -1 
  x = 0
 x*2 = 1  

yaitu 2kx=k. MembiarkanSn menunjukkan berapa banyak uang yang kami hasilkan dari game setelahnya n berbalik, lalu

Pr(Sn=2kx)=2n(n(n+k)/2)
untuk .n(n+k)/20

Ketika bukan kelipatan 2, maka . Untuk memahami ini, anggaplah bahwa kita mulai dengan £ 10. Setelah belokan, satu-satunya nilai yang mungkin adalah £ 5 atau £ 20, yaitu atau .(n+k)Pr(Sn)=0n=1k=1k=1

Hasil di atas adalah hasil standar dari jalan acak. Google random walks untuk info lebih lanjut. Juga dari teori random walk, kita dapat menghitung median return menjadi , yang tidak sama dengan nilai yang diharapkan.x

Catatan: Saya berasumsi bahwa Anda selalu dapat separuh uang Anda. Sebagai contoh, 1pence, 0.5pence, 0.25pence semuanya diizinkan. Jika Anda menghapus asumsi ini, maka Anda memiliki jalan acak dengan dinding penyerap.


Untuk kelengkapan

Berikut ini adalah simulasi cepat dalam R dari proses Anda:

#Simulate 10 throws with a starting amount of x=money=10
#n=10
simulate = function(){
  #money won/lost in a single game
  money = 10
  for(i in 1:10){
    if(runif(1) < 0.5)
      money = money/2
    else
      money = 2*money
  }
  return(money)
}

#The Money vector keeps track of all the games
#N is the number of games we play
N = 1000
Money = numeric(N)
for(i in 1:N)
  Money[i]= simulate()

mean(Money);median(Money)
#Probabilities
#Simulated
table(Money)/1000
#Exact
2^{-10}*choose(10,10/2)

#Plot the simulations
plot(Money)
csgillespie
sumber
Karena (n + k) / 2 tidak harus integral, pertimbangkan menulis ulang probabilitas sebagai Pr (S_n = 2 ^ {2k-n}) = 2 ^ -n Comb (n, k), 0 <= k <= n. (Ada juga sesuatu yang mencurigakan tentang penyamaan 2 ^ kx = k.)
whuber
haruskah variabel Anda moneymenjadi x?
Jeromy Anglim
@Jeromy: Ya, saya sudah mengubahnya. @whuber: Anda benar, saya sudah mencoba membuat probabilitas sedikit lebih jelas. BTW, integral-> integer dalam komentar Anda.
csgillespie
1

Biarkan menjadi kekayaan setelah memainkan game ini, di mana kita menganggap Godaan di sini adalah untuk mengambil , dan mempelajari sebagai jalan acak simetris, dengan inovasi ukuran . Ternyata, ini akan baik-baik saja untuk pertanyaan kedua, tetapi bukan yang pertama. Sedikit pekerjaan akan menunjukkan bahwa, asimptotik kita memiliki . Dari sini Anda tidak dapat menyimpulkan bahwa adalah log asimtotik yang biasanya didistribusikan dengan Operasi log tidak bolak-balik dengan batas. Jika ya, Anda akan mendapatkan nilai yang diharapkan dari sebagaiSkkS0=1.Xk=logSkXk±log2XkN(0,k(log2)2)Skμ=0,σ=log2k.Skexp(klog2log2/2) , yang hampir benar, tetapi tidak cukup.

Namun, metode ini baik-baik saja untuk menemukan kuantil , dan pertanyaan probabilitas lainnya, seperti pertanyaan (2). Kami memilikiKuantitas di sisi kiri dari ketidaksetaraan terakhir adalah, asimptotik, standar normal, dan dengan demikian probabilitas bahwa melebihi rata-rata mendekati mana adalah CDF dari standar normal. Ini mendekati nol dengan cukup cepat.SkSk(54)kXkklog(5/4)Xk/klog2klog(5/4)/log2.Sk1Φ(klog(5/4)/log2),Φ

Kode Matlab untuk memeriksa ini:

top_k = 512;
nsamps = 8192;
innovs = log(2) * cumsum(sign(randn(top_k,nsamps)),1);
s_k = exp(innovs);
k_vals = (1:top_k)';
mean_v = (5/4) .^ k_vals;
exceed = bsxfun(@ge,s_k,mean_v);
prob_g = mean(double(exceed),2);

%theoretical value
%(can you believe matlab doesn't come with normal cdf function!?)
nrmcdf = @(x)((1 + erf(x / sqrt(2)))/2);
p_thry = 1 - nrmcdf(sqrt(k_vals) * log(5/4) / log(2));

loglog(k_vals,prob_g,'b-',k_vals,p_thry,'r-');
legend('empirical probability','theoretical probability');

grafik yang dihasilkan: teks alternatif

shabbychef
sumber
-1

Anda benar tentang harapan itu.

Anda sebenarnya juga memiliki jawaban yang tepat untuk kemungkinan mendapatkan lebih dari taruhan asli Anda, meskipun bukan bukti yang tepat. Pertimbangkan, alih-alih jumlah uang mentah yang Anda miliki, logaritma basis-2-nya. Ini ternyata menjadi berapa kali Anda menggandakan uang Anda, dikurangi berapa kali Anda membagi dua. Ini adalah jumlah dari variabel acak independen, masing-masing sama dengan atau dengan probabilitas . Probabilitas yang Anda inginkan adalah probabilitas bahwa ini positif. Jika aneh, maka dengan simetri itu persis ; jika adalah even (sebut saja ) maka ituSnn+111/2n1/2n2k1/2minus setengah probabilitas bahwa . Tapi , yang mendekati sebagai .Sn=0P(S2k=0)=(2kk)/22k0k

Michael Lugo
sumber
Demonstrasi ini tampaknya mengasumsikan jumlah uang yang diharapkan adalah nol, tetapi tidak. Selain itu, Anda tampaknya menganggap bahwa jaring penggandaan dikurangi separuh harus tidak negatif, yang juga tidak benar. Secara kolektif kesalahan ini menghasilkan batas yang benar, tapi itu kecelakaan.
whuber