<< Chapter < Page | Chapter >> Page > |
If any one or more of these product expressions fail, the class is not independent. A similar situation holds for a class of four events: the product rule must hold for every pair,for every triple, and for the whole class. Note that we say “not independent” or “nonindependent” rather than dependent. The reason for this becomes clearer in dealing withindependent random variables.
We consider some classical exmples of nonindependent classes
As noted above, the replacement rule holds for any pair of events. It is easy to show, although somewhat cumbersome to write out, that if the rule holds for any finite number k of events in an independent class, it holds for any of them. By the principle of mathematical induction, the rule must hold for any finite subclass. We mayextend the replacement rule as follows.
Suppose the class is independent with respective probabilities , , and . Then
is independent and
is independent and
Similarly, the probabilities of the other six minterms, in order, are 0.21, 0.21, 0.06, 0.06, 0.09, and 0.09. With these minterm probabilities, the probability of any Booleancombination of , and C may be calculated
In general, eight appropriate probabilities must be specified to determine the minterm probabilities for a class of three events. In the independent case, three appropriate probabilities are sufficient.
Suppose is independent with , , and . Then and
With each of the basic probabilities determined, we may calculate the minterm probabilities, hence the probability of any Boolean combination of the events.
Frequently we have a large enough independent class that it is desirable to use MATLAB (or some other computational aid) to calculate the probabilitiesof various “and” combinations (intersections) of the events or their complements. Suppose the independent class has respective probabilities
It is desired to calculate (a) , and (b) .
We may use the MATLAB function prod and the scheme for indexing a matrix.
>>p = 0.01*[13 37 12 56 33 71 22 43 57 31];>>q = 1-p;>>% First case>>e = [1 2 4 7]; % Uncomplemented positions>>f = [3 5 6]; % Complemented positions>>P = prod(p(e))*prod(q(f)) % p(e) probs of uncomplemented factors
P = 0.0010 % q(f) probs of complemented factors>>% Case of uncomplemented in even positions; complemented in odd positions>>g = find(rem(1:10,2) == 0); % The even positions>>h = find(rem(1:10,2) ~= 0); % The odd positions>>P = prod(p(g))*prod(q(h))
P = 0.0034
In the unit on MATLAB and Independent Classes, we extend the use of MATLAB in the calculations for such classes.
Notification Switch
Would you like to follow the 'Applied probability' conversation and receive update notifications?