贝叶斯估计、最大似然估计、最大后验概率估计

Bayes' Theorem

文章作者：Tyan
博客：noahsnail.com | CSDN | 简书

1. 引言

贝叶斯估计、最大似然估计(MLE)、最大后验概率估计(MAP)这几个概念在机器学习和深度学习中经常碰到，读文章的时候还感觉挺明白，但独立思考时经常会傻傻分不清楚(😭)，因此希望通过本文对其进行总结。

2. 背景知识

注：由于概率与数理统计需要了解的背景知识很多，因此这里只列出了部分内容，且写的较简略，许多概念的学习需要根据标题自己查找答案。

2.1 概率与统计

概率统计是很多人都学过的内容，但概率论与统计学的关系是什么？先看一下概率论与统计学在维基百科中的定义：

概率论是集中研究概率及随机现象的数学分支，是研究随机性或不确定性等现象的数学。
统计学是在数据分析的基础上，研究如何测定、收集、整理、归纳和分析反映数据数据，以便给出正确消息的科学。

下面的一段话引自LarrB Wasserman的《All of Statistics》，对概率和统计推断的研究内容进行了描述：

The basic problem that we studB in probabilitB is: 
Given a data generating process, what are the properities of the outcomes?

The basic problem of statistical inference is the inverse of probabilitB: 
Given the outcomes, what can we saB about the process that generated the data?

概率论是在给定条件（已知模型和参数）下，对要发生的事件（新输入数据）的预测。统计推断是在给定数据（训练数据）下，对数据生成方式（模型和参数）的归纳总结。概率论是统计学的数学基础，统计学是对概率论的应用。

2.2 描述统计和推断统计

统计学分为描述统计学和推断统计学。描述统计，是统计学中描绘或总结观察量基本情况的统计总称。推断统计指统计学中研究如何根据样本数据去推断总体数量特征的方法。

描述统计是对数据的一种概括。描述统计是罗列所有数据，然后选择一些特征量（例如均值、方差、中位数、四分中位数等）对总体数据进行描述。推断统计是一种对数据的推测。推断统计无法获取所有数据，只能得到部分数据，然后根据得到的数据推测总体数据的情况。

2.3 联合概率和边缘概率

假设有随机变量 $A$ 和 $B$ ，此时 $P(A=a,B=b)$ 用于表示 $A=a$ 且 $B=b$ 同时发生的概率。这类包含多个条件且所有条件同时成立的概率称为联合概率。请注意，联合概率并不是其中某个条件成立的概率，而是所有条件同时成立的概率。与之对应地， $P(A=a)$ 或 $P(B=b)$ 这类仅与单个随机变量有关的概率称为边缘概率。

联合概率与边缘概率的关系如下：

$P(A=a)=\sum_{b}P(A=a,B=b)$
$P(B=b)=\sum_{a}P(A=a,B=b)$

2.4 条件概率

条件概率表示在条件 $B=b$ 成立的情况下， $A=a$ 的概率，记作 $P(A=a|B=b)$ ，或者说条件概率是指事件 $A=a$ 在另外一个事件 $B=b$ 已经发生条件下的发生概率。为了简洁表示，后面省略a，b。

联合概率、边缘概率、条件概率的关系如下：

$P(A|B)=\frac {P(A,B)} {P(B)}$

转换为乘法形式：

$P(A,B)=P(B)*P(A|B)=P(A)*P(B|A)$

2.5 全概率公式

如果事件 $A_1，A_2，A_3，\ldots，A_n$ 构成一个完备事件组，即它们两两互不相容（互斥），其和为全集；并且 $P(A_i)$ 大于0，则对任意事件 $B$ 有 $P(B)=P(B|A_1)P(A_1)+P(B|A_2)P(A_2)+\ldots+ P(B|A_n)P(A_n)=\sum^n_{i=1}P(B|A_i)P(A_i)$ 上面的公式称为全概率公式。全概率公式是对复杂事件 $A$ 的概率求解问题转化为了在不同情况下发生的简单事件的概率的求和问题。

2.6 贝叶斯公式

由条件概率的乘法形式可得：

$P(A|B)=\frac {P(B|A)} {P(B)}*P(A)$

上面的式子称为贝叶斯公式，也叫做贝叶斯定理或贝叶斯法则。在贝叶斯定理中，每个名词都有约定俗成的名称：

$P(A|B)$ 是已知 $B$ 发生后 $A$ 的条件概率，也由于得自 $B$ 的取值而被称作 $A$ 的后验概率，表示事件 $B$ 发生后，事件 $A$ 发生的置信度。
$P(A)$ 是 $A$ 的先验概率或边缘概率，表示事件 $A$ 发生的置信度。
$P(B|A)$ 是已知 $A$ 发生后 $B$ 的条件概率，也由于得自 $A$ 的取值而被称作 $B$ 的后验概率，也被称作似然函数。
$P(B)$ 是 $B$ 的先验概率或边缘概率，称为标准化常量。
$\frac {P(B|A)} {P(B)}$ 称为标准似然比(这个叫法很多，没找到标准统一的叫法)，表示事件 $B$ 为事件 $A$ 发生提供的支持程度。

因此贝叶斯公式可表示为：后验概率=似然函数先验概率/标准化常量=标准似然比先验概率。根据标准似然比的大小，可分为下面三种情况：

如果标准似然比 $>1$ ，则先验概率 $P(A)$ 得到增强，事件 $B$ 的发生会增大事件 $A$ 发生的可能性；
如果标准似然比 $=1$ ，则先验概率 $P(A)$ 保持不变，事件 $B$ 的发生不影响事件 $A$ 发生的可能性；
如果标准似然比 $<1$ ，则先验概率 $P(A)$ 得到削弱，事件 $B$ 的发生会降低事件 $A$ 发生的可能性。

由全概率公式、贝叶斯法则可得：
$P(A_i|B)=\frac {P(B|A_i)P(A_i)} {P(B)}=\frac {P(B|A_i)P(A_i)} {\sum^n_{i=1}P(B|A_i)P(A_i)}$

2.7 似然与概率

在英文中，似然（likelihood）和概率（probability）是同义词，都指事件发生的可能性。但在统计中，似然与概率是不同的东西。概率是已知参数，对结果可能性的预测。似然是已知结果，对参数是某个值的可能性预测。

2.8 似然函数与概率函数

对于函数 $P(x|\theta)$ ，从不同的观测角度来看可以分为以下两种情况：

如果 $\theta$ 已知且保持不变， $x$ 是变量，则 $P(x|\theta)$ 称为概率函数，表示不同 $x$ 出现的概率。
如果 $x$ 已知且保持不变， $\theta$ 是变量，则 $P(x|\theta)$ 称为似然函数，表示不同 $\theta$ 下， $x$ 出现的概率，也记作 $L(\theta|x)$ 或 $L(x;\theta)$ 或 $f(x;\theta)$ 。

注：注意似然函数的不同写法。

2.9 推断统计中需要了解的一些概念

假设实际观测值与真实分布相关，试图根据观测值来推测真实分布
由于观测值取值随机，因此由它们计算得到的估计值也是随机值
估计方式多种多样，且不同估计方式得到的估计值也有所不同

样本、样本容量、参数统计、非参数统计、估计量、真实分布、经验分布。

2.10 频率学派与贝叶斯学派

注：频率学派与贝叶斯学派只是解决问题的角度不同。

频率学派与贝叶斯学派探讨「不确定性」这件事时的出发点与立足点不同。频率学派从「自然」角度出发，试图直接为「事件」本身建模，即事件 $A$ 在独立重复试验中发生的频率趋于极限 $p$ ，那么这个极限就是该事件的概率。

贝叶斯学派并不从试图刻画「事件」本身，而从「观察者」角度出发。贝叶斯学派并不试图说「事件本身是随机的」，或者「世界的本体带有某种随机性」，这套理论根本不言说关于「世界本体」的东西，而只是从「观察者知识不完备」这一出发点开始，构造一套在贝叶斯概率论的框架下可以对不确定知识做出推断的方法。

频率学派的代表是最大似然估计；贝叶斯学派的代表是最大后验概率估计。

2.11 共轭先验

在贝叶斯统计中，如果后验分布与先验分布属于同类，则先验分布与后验分布被称为共轭分布，而先验分布被称为似然函数的共轭先验。

2.12 Beta分布

在概率论中，Beta分布也称Β分布，是指一组定义在 $(0,1)$ 区间的连续概率分布，有两个参数 $\alpha,\beta>0$ 。Beta分布的概率密度为：

$\begin{align}f(x;\alpha,\beta)&=\frac {x^{\alpha-1}(1-x)^{\beta-1}} {\int_{0}^1 \mu^{\alpha-1}(1-\mu)^{\beta-1}d\mu} \\\\&= \frac{\Gamma(\alpha+\beta)} {\Gamma(\alpha)\Gamma(\beta)}x^{\alpha-1}(1-x)^{\beta-1} \\\\&=\frac {1} {B(\alpha,\beta)}x^{\alpha-1}(1-x)^{\beta-1}\end{align}$ 其中， $\Gamma(z)$ 是 $\Gamma$ 函数。随机变量 $X$ 服从Beta分布写作 $X\sim Beta(\alpha,\beta)$ 。

3. 问题定义

以抛硬币为例，假设我们有一枚硬币，现在要估计其正面朝上的概率 $\theta$ 。为了对 $\theta$ 进行估计，我们进行了10次实验（独立同分布，i.i.d.），这组实验记为 $X=x_1，x_2，\ldots，x_{10}$ ，其中正面朝上的次数为6次，反面朝上的次数为4次，结果为 $(1,0,1,1,0,0,0,1,1,1)$ 。

4. 最大似然估计(MLE)

最大似然估计，英文为Maximum Likelihood Estimation，简写为MLE，也叫极大似然估计，是用来估计概率模型参数的一种方法。最大似然估计的思想是使得观测数据（样本）发生概率最大的参数就是最好的参数。

对一个独立同分布的样本集来说，总体的似然就是每个样本似然的乘积。针对抛硬币的问题，似然函数可写作： $L(X;\theta)=\prod_{i=0}^nP(x_i|\theta)=\theta^6(1-\theta)^4$ 根据最大似然估计，使 $L(X;\theta)$ 取得最大值的 $\theta$ 即为估计结果，令 $L(X;\theta)\prime =0$ 可得 $\hat{\theta}=0.6$ 。似然函数图如下：

MLE

由于总体的似然就是每个样本似然的乘积，为了求解方便，我们通常会将似然函数转成对数似然函数，然后再求解。可以转成对数似然函数的主要原因是对数函数并不影响函数的凹凸性。因此上式可变为： $lnL(X;\theta)=ln\prod_{i=0}^nP(x_i|\theta)=\sum_{i=0}^nln(P(x_i|\theta))=6ln(\theta)+4ln(1-\theta)$ 令 $ln(L(X;\theta)\prime) =0$ 可得 $\hat{\theta}=0.6$ 。

正态分布的最大似然估计

假设样本服从正态分布 $N\sim(\mu,\sigma^2)$ ，则其似然函数为 $L(\mu,\sigma^2)=\prod_{i=0}^n \frac {1} {\sqrt{2\pi} \sigma}e^{-\frac {(x_i-\mu)^2} {2\sigma^2}}$ 对其取对数得： $lnL(\mu,\sigma^2)=-\frac {n} {2}ln(2\pi) - \frac {n} {2} ln(\sigma^2) - \frac {1} {2\sigma^2} \sum_{i=0}^n(x_i-\mu)^2$
分别对 $\mu，\sigma^2$ 求偏导，并令偏导数为0，得： $\begin{cases} \frac {\partial lnL(\mu,\sigma^2)} {\partial \mu}= \frac {1} {\sigma^2} \sum_{i=0}^n(x_i-\mu) =0\\\\ \frac {\partial lnL(\mu,\sigma^2)} {\partial \sigma^2}= -\frac {n} {2\sigma^2} + \frac {1} {2\sigma^4}\sum_{i=0}^n(x_i-\mu)^2 =0 \end{cases}$

解得：
$\begin{cases} \hat{\mu}= \frac {1} {n} \sum_{i=0}^nx_i=\bar{x}\\\\ \hat{\sigma^2} = \frac {1} {n} \sum_{i=0}^n(x_i-\bar{x})^2 \end{cases}$

$\hat{\mu}，\hat{\sigma^2}$ 就是正态分布中 $\mu，\sigma^2$ 的最大似然估计。

最大似然估计的求解步骤：

确定似然函数
将似然函数转换为对数似然函数
求对数似然函数的最大值（求导，解似然方程）

5. 最大后验概率估计(MAP)

最大后验概率估计，英文为Maximum A Posteriori Estimation，简写为MAP。回到抛硬币的问题，最大似然估计认为使似然函数 $P(X|\theta)$ 最大的参数 $\theta$ 即为最好的 $\theta$ ，此时最大似然估计是将 $\theta$ 看作固定的值，只是其值未知；最大后验概率分布认为 $\theta$ 是一个随机变量，即 $\theta$ 具有某种概率分布，称为先验分布，求解时除了要考虑似然函数 $P(X|\theta)$ 之外，还要考虑 $\theta$ 的先验分布 $P(\theta)$ ，因此其认为使 $P(X|\theta)P(\theta)$ 取最大值的 $\theta$ 就是最好的 $\theta$ 。此时要最大化的函数变为 $P(X|\theta)P(\theta)$ ，由于 $X$ 的先验分布 $P(X)$ 是固定的（可通过分析数据获得，其实我们也不关心 $X$ 的分布，我们关心的是 $\theta$ ），因此最大化函数可变为 $\frac {P(X|\theta)P(\theta)} {P(X)}$ ，根据贝叶斯法则，要最大化的函数 $\frac {P(X|\theta)P(\theta)} {P(X)}=P(\theta|X)$ ，因此要最大化的函数是 $P(\theta|X)$ ，而 $P(\theta|X)$ 是 $\theta$ 的后验概率。最大后验概率估计可以看作是正则化的最大似然估计，当然机器学习或深度学习中的正则项通常是加法，而在最大后验概率估计中采用的是乘法， $P(\theta)$ 是正则项。在最大似然估计中，由于认为 $\theta$ 是固定的，因此 $P(\theta)=1$ 。

最大后验概率估计的公式表示： $\mathop{argmax}_{\theta}P(\theta|X)=\mathop{argmax}_{\theta}\frac {P(X|\theta)P(\theta)} {P(X)}\propto \mathop{argmax}_{\theta}P(X|\theta)P(\theta)$

在抛硬币的例子中，通常认为 $\theta=0.5$ 的可能性最大，因此我们用均值为 $0.5$ ，方差为 $0.1$ 的高斯分布来描述 $\theta$ 的先验分布，当然也可以使用其它的分布来描述 $\theta$ 的先验分布。 $\theta$ 的先验分布为： $\frac {1} {\sqrt{2\pi}\sigma}e^{-\frac {(\theta-\mu)^2} {2\sigma^2}} = \frac {1} {10\sqrt{2\pi}}e^{-50(\theta-0.5)^2}$ 先验分布的函数图如下：

Gaussian

在最大似然估计中，已知似然函数为 $P(X|\theta)=\theta^6(1-\theta)^4$ ，因此： $P(X|\theta)P(\theta)=\theta^6\times (1-\theta)^4\times \frac {1} {10\sqrt{2\pi}}\times e^{-50(\theta-0.5)^2}$ 转换为对数函数： $ln(P(X|\theta)P(\theta))=ln(\theta^6\times (1-\theta)^4 \times \frac {1} {10\sqrt{2\pi}}\times e^{-50(\theta-0.5)^2})=6ln(\theta)+4ln(1-\theta)+ln(\frac {1} {10\sqrt{2\pi}})-50(\theta-0.5)^2$

令 $ln(P(X|\theta)P(\theta))\prime=0$ ，可得： $100\theta^3-150\theta^2+40\theta+6=0$ 由于 $0\le\theta\le1$ ，解得： $\hat{\theta}\approx0.529$ 。 $P(X|\theta)P(\theta)$ 的函数图像如下，基本符合 $\theta$ 的估计值 $\hat{\theta}$ ：

MAP

如果我们用均值为 $0.6$ ，方差为 $0.1$ 的高斯分布来描述 $\theta$ 的先验分布，则 $\hat{\theta}=0.6$ 。由此可见，在最大后验概率估计中， $\theta$ 的估计值与 $\theta$ 的先验分布有很大的关系。这也说明一个合理的先验概率假设是非常重要的。如果先验分布假设错误，则会导致估计的参数值偏离实际的参数值。

先验分布为Beta分布

如果用 $\alpha=3,\beta=3$ 的Beta分布来描述 $\theta$ 的先验分布，则 $P(X|\theta)P(\theta)=\theta^6\times (1-\theta)^4\times \frac {1} {B(\alpha,\beta)}\times \theta^{\alpha-1}(1-\theta)^{\beta-1}$ 令 $P(X|\theta)P(\theta)\prime=0$ 求解可得： $\hat{\theta}=\frac {\alpha+5} {\alpha + \beta +8}=\frac {8} {3 + 3 +8}\approx 0.57$

$Beta(3,3)$ 的概率密度图像如下图：

Beta(3,3)

最大后验概率估计的求解步骤：

确定参数的先验分布以及似然函数
确定参数的后验分布函数
将后验分布函数转换为对数函数
求对数函数的最大值（求导，解方程）

6. 贝叶斯估计

贝叶斯估计是最大后验估计的进一步扩展，贝叶斯估计同样假定 $\theta$ 是一个随机变量，但贝叶斯估计并不是直接估计出 $\theta$ 的某个特定值，而是估计 $\theta$ 的分布，这是贝叶斯估计与最大后验概率估计不同的地方。在贝叶斯估计中，先验分布 $P(X)$ 是不可忽略的。回到抛硬币的例子中，在已知 $X$ 的情况下，描述 $\theta$ 的分布即描述 $P(\theta|X)$ ， $P(\theta|X)$ 是一种后验分布。如果后验分布的范围较窄，则估计值的准确度相对较高，反之，如果后验分布的范围较广，则估计值的准确度就较低。

贝叶斯公式： $P(\theta|X)=\frac {P(X|\theta)P(\theta)} {P(X)}$

在连续型随机变量中，由于 $P(X)=\int_{\Theta}P(X|\theta)P(\theta)d\theta$ ，因此贝叶斯公式变为： $P(\theta|X)=\frac {P(X|\theta)P(\theta)} {\int_{\Theta}P(X|\theta)P(\theta)d\theta}$

从上面的公式中可以看出，贝叶斯估计的求解非常复杂，因此选择合适的先验分布就非常重要。一般来说，计算积分 $\int_{\theta}P(X|\theta)P(\theta)d\theta$ 是不可能的。对于这个抛硬币的例子来说，如果使用共轭先验分布，就可以更好的解决这个问题。二项分布参数的共轭先验是Beta分布，由于 $\theta$ 的似然函数服从二项分布，因此在贝叶斯估计中，假设 $\theta$ 的先验分布服从 $P(\theta)\sim Beta(\alpha, \beta)$ ，Beta分布的概率密度公式为： $f(x;\alpha,\beta)=\frac {1} {B(\alpha,\beta)}x^{\alpha-1}(1-x)^{\beta-1}$ 因此，贝叶斯公式可写作： $\begin{aligned} P(\theta|X)&=\frac {P(X|\theta)P(\theta)} {\int_{\Theta}P(X|\theta)P(\theta)d\theta} \\\\ &=\frac {\theta^6(1-\theta)^4 \frac {\theta^{\alpha-1}(1-\theta)^{\beta-1}} {B(\alpha,\beta)} } {\int_{\Theta}\theta^6(1-\theta)^4 \frac {\theta^{\alpha-1}(1-\theta)^{\beta-1}} {B(\alpha,\beta)}d\theta} \\\\&=\frac {\theta^{\alpha+6-1}(1-\theta)^{\beta+4-1}} {\int_{\Theta}\theta^{\alpha+6-1}(1-\theta)^{\beta+4-1}d\theta} \\\\ &=\frac {\theta^{\alpha+6-1}(1-\theta)^{\beta+4-1}} {B(\alpha+6-1,\beta+4-1)} \\\\ &=Beta(\theta|\alpha+6-1,\beta+4-1) \\\\&=Beta(\theta|\alpha+6,\beta+4)\end{aligned}$ 从上面的公式可以看出， $P(\theta|X) \sim Beta(\theta|\alpha+6,\beta+4)$ 。其中 $B$ 函数，也称 $Beta$ 函数，是一个标准化常量，用来使整个概率的积分为1。 $Beta(\theta|\alpha+6,\beta+4)$ 就是贝叶斯估计的结果。

如果使用贝叶斯估计得到的 $\theta$ 分布存在一个有限均值，则可以用后验分布的期望作为 $\theta$ 的估计值。假设 $\alpha=3,\beta=3$ ，在这种情况下，先验分布会在 $0.5$ 处取得最大值，则 $P(\theta|X) \sim Beta(\theta|9,7)$ ， $Beta(\theta|9,7)$ 的曲线如下图：

Beta(9,7)

从上图可以看出，在 $\alpha=3,\beta=3$ 的情况下， $\theta$ 的估计值 $\hat{\theta}$ 应该在 $0.6$ 附近。根据Beta分布的数学期望公式 $E(\theta)=\frac {\alpha} {\alpha+\beta}$ 可得： $\hat{\theta}=\int_{\Theta} \theta P(\theta|X)d\theta=E(\theta)=\frac {\alpha} {\alpha+\beta}=\frac {9} {9+7}=0.5625$

注：二项分布参数的共轭先验是Beta分布，多项式分布参数的共轭先验是Dirichlet分布，指数分布参数的共轭先验是Gamma分布，⾼斯分布均值的共轭先验是另⼀个⾼斯分布，泊松分布的共轭先验是Gamma分布。

贝叶斯估计要解决的不是如何估计参数，而是用来估计新测量数据出现的概率，对于新出现的数据 $\tilde{x}$ ：

$P(\tilde{x}|X) = \int_{\Theta}P(\tilde{x}|\theta)P(\theta|X)d\theta=\int_{\Theta}P(\tilde{x}|\theta)\frac {P(X|\theta)P(\theta)} {P(X)}d\theta$

贝叶斯估计的求解步骤：

确定参数的似然函数
确定参数的先验分布，应是后验分布的共轭先验
确定参数的后验分布函数
根据贝叶斯公式求解参数的后验分布

7. 总结

从最大似然估计、最大后验概率估计到贝叶斯估计，从下表可以看出 $\theta$ 的估计值 $\hat{\theta}$ 是逐渐接近 $0.5$ 的。从公式的变化可以看出，使用的信息是逐渐增多的。最大似然估计、最大后验概率估计中都是假设 $\theta$ 未知，但是确定的值，都将使函数取得最大值的 $\theta$ 作为估计值，区别在于最大化的函数不同，最大后验概率估计使用了 $\theta$ 的先验概率。而在贝叶斯估计中，假设参数 $\theta$ 是未知的随机变量，不是确定值，求解的是参数 $\theta$ 在样本 $X$ 上的后验分布。

注：最大后验概率估计和贝叶斯估计都采用Beta分布作为先验分布。

Type	MLE	MAP	BE
$\hat{\theta}$	0.6	0.57	0.5625
$f$	$P(X \| \theta)$	$P(X\|\theta)P(\theta)$	$\frac {P(X\|\theta)P(\theta)} {P(X)}$

参考资料

最后编辑于：2021.08.03 16:01:17

人面猴
序言：七十年代末，一起剥皮案震惊了整个滨河市，随后出现的几起案子，更是在滨河造成了极大的恐慌，老刑警刘岩，带你破解...
沈念sama阅读 159,015评论 4赞 362
死咒
序言：滨河连续发生了三起死亡事件，死亡现场离奇诡异，居然都是意外死亡，警方通过查阅死者的电脑和手机，发现死者居然都...
沈念sama阅读 67,262评论 1赞 292
救了他两次的神仙让他今天三更去死
文/潘晓璐我一进店门，熙熙楼的掌柜王于贵愁眉苦脸地迎上来，“玉大人，你说我怎么就摊上这事。” “怎么了？”我有些...
开封第一讲书人阅读 108,727评论 0赞 243
道士缉凶录：失踪的卖姜人
文/不坏的土叔我叫张陵，是天一观的道长。经常有香客问我，道长，这世上最难降的妖魔是什么？我笑而不...
开封第一讲书人阅读 43,986评论 0赞 205
港岛之恋（遗憾婚礼）
正文为了忘掉前任，我火速办了婚礼，结果婚礼上，老公的妹妹穿的比我还像新娘。我一直安慰自己，他们只是感情好，可当我...
茶点故事阅读 52,363评论 3赞 287
恶毒庶女顶嫁案：这布局不是一般人想出来的
文/花漫我一把揭开白布。她就那样静静地躺着，像睡着了一般。火红的嫁衣衬着肌肤如雪。梳的纹丝不乱的头发上，一...
开封第一讲书人阅读 40,610评论 1赞 219
城市分裂传说
那天，我揣着相机与录音，去河边找鬼。笑死，一个胖子当着我的面吹牛，可吹牛的内容都是我干的。我是一名探鬼主播，决...
沈念sama阅读 31,871评论 2赞 312
双鸳鸯连环套：你想象不到人心有多黑
文/苍兰香墨我猛地睁开眼，长吁一口气：“原来是场噩梦啊……” “哼！你这毒妇竟也来了？” 一声冷哼从身侧响起，我...
开封第一讲书人阅读 30,582评论 0赞 198
万荣杀人案实录
序言：老挝万荣一对情侣失踪，失踪者是张志新（化名）和其女友刘颖，没想到半个月后，有当地人在树林里发现了一具尸体，经...
沈念sama阅读 34,297评论 1赞 242
护林员之死
正文独居荒郊野岭守林人离奇死亡，尸身上长有42处带血的脓包…… 初始之章·张勋以下内容为张勋视角年9月15日...
茶点故事阅读 30,551评论 2赞 246
白月光启示录
正文我和宋清朗相恋三年，在试婚纱的时候发现自己被绿了。大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
茶点故事阅读 32,053评论 1赞 260
活死人
序言：一个原本活蹦乱跳的男人离奇死亡，死状恐怖，灵堂内的尸体忽然破棺而出，到底是诈尸还是另有隐情，我是刑警宁泽，带...
沈念sama阅读 28,385评论 2赞 253
日本核电站爆炸内幕
正文年R本政府宣布，位于F岛的核电站，受9级特大地震影响，放射性物质发生泄漏。R本人自食恶果不足惜，却给世界环境...
茶点故事阅读 33,035评论 3赞 236
男人毒药：我在死后第九天来索命
文/蒙蒙一、第九天我趴在偏房一处隐蔽的房顶上张望。院中可真热闹，春花似锦、人声如沸。这庄子的主人今日做“春日...
开封第一讲书人阅读 26,079评论 0赞 8
一桩弑父案，背后竟有这般阴谋
文/苍兰香墨我抬头看了看天上的太阳。三九已至，却和暖如春，着一层夹袄步出监牢的瞬间，已是汗流浃背。一阵脚步声响...
开封第一讲书人阅读 26,841评论 0赞 195
情欲美人皮
我被黑心中介骗来泰国打工，没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留，地道东北人。一个月前我还...
沈念sama阅读 35,648评论 2赞 274
代替公主和亲
正文我出身青楼，却偏偏与公主长得像，于是被迫代替她去往敌国和亲。传闻我的和亲对象是个残疾皇子，可洞房花烛夜当晚...
茶点故事阅读 35,550评论 2赞 270