2011-10-19 6 views
9

Czy możliwe jest obliczenie prawdopodobieństwa 3 lub więcej głowic z 4 rzutów monetą przy użyciu funkcji prawdopodobieństwa lub funkcji NProbability.Korzystanie z opcji NProbability [] lub Probability [] w celu obliczenia prawdopodobieństwa 3 lub więcej głowic z 4 rzutów monetą

To nie jest pytanie o trywialną odpowiedź na ten problem, bardziej jest zrozumienie, jak rozwiązać ten problem z Mathematica za pomocą dystrybucji.

Więc za pomocą 4 zmiennych losowych z podziału P

Miałem nadzieję, coś takiego by rade, ale to nie działa. Otrzymuję 0.

P = BernoulliDistribution[0.5]; 
vars = List[Distributed[a,P],Distributed[b,P],Distributed[c,P],Distributed[c,P]]; 
NProbability[Count[ {a,b,c,d}, 1] >= 3, vars] 

Wszelkie pomysły będą mile widziane.

+0

dwa błędy tam: c zadeklarowane dwa razy i nigdy nie używając Count zamiast Total. –

+0

Być może @Sjoerd jest zbyt grzeczny. Po jego modyfikacjach kod OP działa –

+0

Ha! A więc tak. Podwójna deklaracja to błąd wklejania kopii na stronie. Używanie Total zamiast Counta dało podstęp. Wielkie dzięki chłopaki. – Bart

Odpowiedz

10

Nie ekspertem używając MMA dla statystyki tutaj, ale to wydaje się działać:

l = TransformedDistribution[ 
     x + y + w + z, {x \[Distributed] BernoulliDistribution[0.5], 
         y \[Distributed] BernoulliDistribution[0.5], 
         z \[Distributed] BernoulliDistribution[0.5], 
         w \[Distributed] BernoulliDistribution[0.5]}]; 

Table[NProbability[x > i, x \[Distributed] l], {i, -1, 4}] 
(* 
{1, 0.9375, 0.6875, 0.3125, 0.0625, 0.} 
*) 
+0

Wow. Dzięki wielkie! Wyjaśniłem to doskonale :) – Bart

+0

Na konferencji mma tech właśnie teraz z leonidem. Nie mamy zbyt wiele zębów do robienia czegokolwiek innego. –

+0

@Sjoerd Conference? Z piwem czy bez piwa? –

8
In[10]:= Probability[a + b + c + d >= 3, vars] 

Out[10]= 0.3125 

rzut monetą jest łatwiejsze opisane z BinomialDistribution:

In[12]:= Probability[m >= 3, m \[Distributed] BinomialDistribution[4, 0.5]] 

Out[12]= 0.3125 
+0

Tak, to dużo czystsze. Zauważyłem, że TransformedDistribution, zgodnie z sugestią belisariusa, jest oceniane jako BinomialDistribution. Thankyou :) – Bart

+0

@Bart Myślałem, że zrobiłeś to celowo, ponieważ chciałeś wiedzieć, jak dodać losowe vary! Ha! Głupi ja. –

+0

@belisarius Chciałem wiedzieć, jak dodać przypadkowe vary, a nie tylko ten problem z klapką monet. Zatem ocenianie na dwumianowy dysk jest interesującym punktem w tym konkretnym przypadku. To, co dostarczyłeś, jest tym, czego naprawdę chciałem, ogólną wiedzą na temat tego rodzaju obliczeń. Więc byłeś na miejscu! :) – Bart

Powiązane problemy