= round(rand(1000000,3) – 0.2);
for i= 1:6
b = a(1:10^i,3;
c(i) =sum(b)/(10^i);
end
c
for i =1:6
b = (a(1:10^i,1))&(a(1:10^i,2))&a(1:10^i,3);
d(i) = sum(b)/(10^i);
end
d
for I = 1:6
b=a(1:10^i,1)&a(1:10^i,2)&a(1:10^i,3);
e(i) = sum(b)/(10^i);
end
e
这段代码是用来计算条件概率的。具体来说,它计算了给定三个随机变量的条件下,另外两个随机变量的条件概率。
首先,通过调用rand
函数生成一个 1000000 x 3 的随机矩阵,并使用round
函数四舍五入到最接近的整数。这样做是为了得到某个随机变量的一组观测值。
然后,通过循环计算,分别计算了在每个观测子集的条件下,第三个随机变量的概率。具体而言,首先根据观测子集大小10^i
获取相应的观测值子集b
,然后通过求和并除以观测子集大小,计算出第三个随机变量的条件概率,并将结果存储在向量c
中。
接下来,使用类似的方法分别计算第二个和第一个随机变量的条件概率,将结果分别存储在向量d
和e
中。
这段代码的目的是通过模拟计算条件概率来估计每个随机变量的条件概率。