五、数据科学导论——数学基础之概率第1关概率基础之贝叶斯任务描述本关任务根据相关知识完成与数据概率相关的选择题。相关知识我们可以把概率论视为从事件空间中抽取的事件的不确定性进行量化的一种方式。空间是指所有可能的结果的集合这些结果的任何一部分就是一个事件。比如“骰子掷出的点数为1”“骰子掷出的点数是偶数”等都是事件。不独立和独立如果 E 发生意味着 F 发生或者 F 发生意味着 E 发生我们就称事件 E 与事件 F 为不相互独立dependent。反之E 与 F 就相互独立。举个例子如果两次掷起一枚均匀的硬币那么我们无法根据第一次掷硬币的结果是否是正面朝上来判断第二次的结果是否正面朝上。第一次掷硬币的结果和第二次掷硬币的结果这两个事件就是独立的。相反如果第一次掷硬币的结果是正面朝上那么我们能很明显能判断出两次掷硬币是否都是反面朝上。如果第一次掷硬币的结果是正面朝上那么很明显两次掷硬币的结果不可能都是反面朝上。因此第一次结果正面朝上和两次结果不可能都是反面朝上这两个事件就是不独立的。从数学角度讲事件 E 和事件 F 独立意味着两个事件同时发生的概率等于它们分别发生的概率的乘积P(E,F)P(E)P(F)在上例中“第一次掷硬币正面朝上”的概率为 1/2“两次掷硬币都是背面朝上”的概率为 1/4但“第一次掷硬币正面朝上并且两次掷硬币都是背面朝上”的概率为 0。条件概率如果事件 E 与事件 F 独立那么定义式如下P(E,F)P(E)P(F)如果两者不一定独立并且 F 的概率不为零那么 E 关于 F 的条件概率式如下P(E,F)P(E∣F)P(F)如果 E 和 F 独立则上式应该表示为P(E∣F)P(E)上式表示F 是否发生并不会影响 E 是否发生的概率。如果我们假设每个孩子是男孩和是女孩的概率相同第二个孩子的性别概率与第一个孩子的性别概率独立那么事件“没有女孩”的概率是 1/4事件“一个男孩一个女孩”的概率为 1/2事件“两个女孩”的概率为 1/4。现在我们的问题是事件 B“两个孩子都是女孩”关于事件 G“大孩子是女孩”的条件概率是多少用条件概率的定义式进行计算如下P(B∣G)P(B,G)/P(G)P(B)/P(G)1/2事件B与G的交集“两个孩子都是女孩并且大孩子是女孩”刚好是事件 B 本身。一旦你知道两个孩子都是女孩那大孩子必然是女孩。贝叶斯定理贝叶斯定理是条件概率的某种逆运算假设我们需要计算事件 E 基于已发生的事件 F 的条件概率但我们已知的条件仅仅是事件 F 基于已发生的事件 E的条件概率。两次利用条件概率的定义可以得到下式P(E∣F)P(E,F)/P(F)P(E∣F)P(E)/P(F)事件 F 可以分割为两个互不重合的事件“F 和 E 同时发生”与“F 发生 E 不发生”。我们用符号 !E 指代“非 E”即“E 没有发生”有下式P(F)P(F,E)P(F,!E)因此P(E∣F)P(F∣E)P(E)/[P(F∣E)P(E)P(F∣!E)P(!E)假设有这样一种病10000 个人中会有一个得这个病。还假设有种针对该病的测试具有 99% 的可能性能给出正确判断如果患病测试显示“有病”如果健康则显示“无病”。阳性的测试结果意味着什么呢我们用T表示“测试结果阳性”用 D 表示“你患有该病”。那么根据贝叶斯定理如果测试结果为阳性那么你患有该病的概率是P(D∣T)P(T∣D)P(D)/[P(T∣D)P(D)P(T∣!D)P(!D)我们知道P(T|D)即一个人测试结果为阳性并且本人实际患病的概率为 0.99。P(D)即一个人实际患病的概率是 1/100000.0001。P(T |J D)即一个不患病的人检测结果呈阳性的概率是 0.01。P(J D)即一个人实际上不患该病的概率为 0.9999。如果将以上数据代入贝叶斯定理可得P(D∣T)0.98结果表示测试结果为阳性的人实际患病的概率不到 1%。编程要求根据相关知识按照要求完成右侧选择题任务。测试说明平台会对你选择的答案进行判断全对则通过测试。开始你的任务吧祝你成功第2关概率基础之数据分布任务描述本关任务根据下列的相关知识完成与数据概率相关的选择题。相关知识模型的基线取决于数据的好坏数据的好坏取决与你对数据的理解。所以为了更加懂数据就先理解一下数据有哪些分布。随机变量随机变量指这样一种变量其可能的取值有一种联合概率分布。定义一个简单的随机变量掷起一枚硬币如果正面朝上随机变量等于1如果背面朝上随机变量等于0。可以再定义更复杂些的随机变量如掷起一枚硬币 10 次查看正面朝上的次数或者从range(10) 取出的一个值每个数值被取到的可能性都相等。联合分布对变量实现每种可能值都赋予了概率。通过掷硬币得到的随机变量等于 0 的概率为 0.5等于 1 的概率也为 0.5。从 range(10) 中生成随机变量的分布为从 0 到 9 之间的每个数值赋予 0.1 的概率。随机变量也可以基于某些条件事件产生就像其他事件一样。比如上关中提到的双生子例子如果X 是表示女孩个数的随机变量那么 X 等于 0 的概率为 1/4 等于1 的概率为 1/2 等于 2 的概率为 1/4。在已知至少一个孩子是女孩后接着再定义一个新的随机变量 Y表示基于这个条件的女孩的个数。Y 等于 1 的概率为 2/3等于 2 的概率为 1/3 。还定义另一个新的随机变量 Z表示基于大孩子是女孩的条件之上的女孩的个数Z 等于 1 的概率为1/2 等于 2 的概率为 1/2。大部分情况下我们会隐式使用随机变量的概念即在使用时并没有对此加以特别关注。连续分布掷硬币对应的是离散分布对离散的结果赋予正概率。我们常常希望对连续结果的分布进行建模。例如均匀分布函数对 0 到 1 之间的所有值都赋予相同的权重。因为 0 和 1 之间有无数个数字因而对每个点而言赋予的权重几乎是零。因此我们用带概率密度函数的连续分布来表示概率一个变量位于某个区间的概率等于概率密度函数在这个区间上的积分。均匀分布的密度函数如下def uniform_pdf(x):return 1 if x 0 and x 1 else 0一个服从均匀分布的随机变量落在 0.2 和 0.3 之间的概率为 1/10 。在 Python 中random.random() 是按均匀分布生成的伪随机数的函数。累积分布函数如下def uniform_cdf(x):if x 0: return 0 # 均匀分布的随机变量不会小于0elif x 1: return x # e.g. P(X 0.4) 0.4else: return 1 # 均匀分布的随机变量总是小于1这个函数给出了一个随机变量小于等于某一特定值的概率。正态分布正态分布是典型的钟型曲线形态分布函数可以完全由两个参数决定均值 μmu和标准差 σ sigma。均值描述钟型曲线的中心标准差描述曲线有多“宽”。正态分布的分布函数如下实现代码def normal_pdf(x, mu0, sigma1):sqrt_two_pi math.sqrt(2 * math.pi)return (math.exp(-(x-mu) ** 2 / 2 / sigma ** 2) / (sqrt_two_pi * sigma))让我们来绘制一些概率密度函数来看看它们的形状如何xs [x / 10.0 for x in range(-50, 50)]plt.plot(xs,[normal_pdf(x,sigma1) for x in xs],-,labelmu0,sigma1)plt.plot(xs,[normal_pdf(x,sigma2) for x in xs],--,labelmu0,sigma2)plt.plot(xs,[normal_pdf(x,sigma0.5) for x in xs],:,labelmu0,sigma0.5)plt.plot(xs,[normal_pdf(x,mu-1) for x in xs],-.,labelmu-1,sigma1)plt.legend()plt.title(多个正态分布的概率密度函数)plt.show()如果 μ0 并且σ 1这个分布称为标准正态分布。如果 Z 是服从标准正态分布的随机变量则有如下转换式XσZμ其中 X 也是正态分布但均值是 μ标准差是 σ 。相反如果 X 是均值为 μ 标准差为σ的正态分布那么Z(X−μ)/σ是标准正态分布的随机变量。标准正态分布的累积分布函数无法用“基本”的解析形式表示但在 Python 中可以用函数 math.erf 描述def normal_cdf(x, mu0,sigma1):return (1 math.erf((x - mu) / math.sqrt(2) / sigma)) / 2再绘出一系列概率累积分布函数:xs [x / 10.0 for x in range(-50, 50)]plt.plot(xs,[normal_cdf(x,sigma1) for x in xs],-,labelmu0,sigma1)plt.plot(xs,[normal_cdf(x,sigma2) for x in xs],--,labelmu0,sigma2)plt.plot(xs,[normal_cdf(x,sigma0.5) for x in xs],:,labelmu0,sigma0.5)plt.plot(xs,[normal_cdf(x,mu-1) for x in xs],-.,labelmu-1,sigma1)plt.legend(loc4) # 底部右边plt.title(多个正态分布的累积分布函数)plt.show()中心极限定理正态分布的运用如此广泛很大程度上归功于中心极限定理。这个定理说一个定义为大量独立同分布的随机变量的均值的随机变量本身就是接近于正态分布的。特别地如果 x1 ,…,xn 都是均值为 μ、标准差为 σ 的随机变量且 n 很大那么近似正态分布且均值为 μ。等价于:​上式近似正态分布均值为 0 标准差为 1。举一个例子带有n和p两个参数的二项式随机变量。一个二项式随机变量Binonimal(n,p) 是 n 个独立伯努利随机变量 Bernoulli(p) 之和每个伯努利随机变量等于 1 的概率是 p等于 0 的概率是 1-pdef bernoulli_trial(p):return 1 if random.random() p else 0def binomial(n, p):return sum(bernoulli_trial(p) for _ in range(n))每个伯努利随机变量 Bernoulli(p) 的均值为 p标准差为:​根据中心极限定理当 n 变得很大一个二项式随机变量 Binonimal(n,p)近似于一个正态分布的随机变量其中均值为 μnp标准差为:如果把两个分布都在图上绘出来很容易看出相似性def make_hist(p, n, num_points):data [binomial(n, p) for _ in range(num_points)]# 用条形图绘出实际的二项式样本histogram Counter(data)plt.bar([x - 0.4 for x in histogram.keys()],[v / num_points for v in histogram.values()],0.8,color0.75)mu p * nsigma math.sqrt(n * p * (1 - p))# 用线形图绘出正态近似xs range(min(data), max(data) 1)ys [normal_cdf(i 0.5, mu, sigma) - normal_cdf(i - 0.5, mu, sigma) for i in xs]plt.plot(xs,ys)plt.title(二项分布与正态近似)plt.show()make_hist(0.75, 100, 10000)近似表达的意义在于如果你想知道掷起一枚均匀的硬币 100 次中正面朝上超过 60 次的概率那么可以用一个正态分布 Normal(50, 5) 的随机变量大于 60 的概率来估计。这比计算二项式分布 Binonimal(100, 0.5) 的累积分布函数更容易。编程要求根据相关知识按照要求完成右侧选择题任务。测试说明平台会对你选择的答案进行判断全对则通过测试。开始你的任务吧祝你成功有任何问题都可以随时关注私信