Как выработать условные вероятности в системе Mathematica. Возможно ли это?

Может ли система Mathematica выполнять вычисления условной вероятности по правилу Байеса, не выполняя вычисления вручную? Если да, то как?

Я искал подсказку и в документах Mathemtaica, и в Интернете, но ничего не нашел. Меня не интересует, как сделать правило Байеса вручную через систему Mathematica, я хочу знать, есть ли способ определить условные вероятности и вычислить другие автоматически.

Итак, чтобы использовать игрушечный пример, предполагающий распределение Бернулли

P(Cancer+) = 0.01
P(Cancer-) = 0.99

P(Test+|Cancer+) = 0.9
P(Test-|Cancer+) = 0.1
P(Test+|Cancer-) = 0.2
P(Test-|Cancer-) = 0.8

Можно ли вычислить

P(Cancer+|Test+) = 0.0434

Итак, используя следующее.

Print["P(C+) = ", PCancerT=BernoulliDistribution[0.01]];
Print["P(C-) = ", PCancerF=BernoulliDistribution[0.99]];
Print[]
Print["P(T+|C+) = ", PTestTGivenCancerT=BernoulliDistribution[0.9]];
Print["P(T-|C+) = ", PTestFGivenCancerT=BernoulliDistribution[0.1]];
Print["P(T+|C-) = ", PTestTGivenCancerF=BernoulliDistribution[0.2]];
Print["P(T-|C-) = ", PTestFGivenCancerF=BernoulliDistribution[0.8]];
Print[]
Print["P(T+,C+) = ", PTestTAndCancerT = Probability[vCT&&vTTCT,{vCT\[Distributed]PCancerT,vTTCT\[Distributed]PTestTGivenCancerT}]];
Print["P(T-,C+) = ", PTestFAndCancerT = Probability[vCT&&vTFCF,{vCT\[Distributed]PCancerT,vTFCF\[Distributed]PTestFGivenCancerT}]];
Print["P(T+,C-) = ", PTestTAndCancerF = Probability[vCF&&vTTCF,{vCF\[Distributed]PCancerF,vTTCF\[Distributed]PTestTGivenCancerF}]];
Print["P(T-,C-) = ", PTestFAndCancerF = Probability[vCF&&vTTCF,{vCF\[Distributed]PCancerF,vTTCF\[Distributed]PTestFGivenCancerF}]];
Print[]
Print["P(C+|T+) = ?"];
Print["P(C+|T-) = ?"];
Print["P(C-|T+) = ?"];
Print["P(C-|T-) = ?"];

Я могу вычислить совместные вероятности, определяя все таблицы вероятностей вручную, но есть ли способ заставить систему Mathematica сделать эту тяжелую работу? Есть ли способ определить и вычислить такие условные вероятности?

Большое спасибо за любую помощь, даже если это "Вы не можете... прекратите попытки" :)

PS: это была попытка сделать что-то в этом роде? Символическое условное ожидание в системе Mathematica

11
задан Community 23 May 2017 в 11:52
поделиться