开始使用免费开始使用

评估正态性的图形方法

在视频中,您学习了如何用 20 个柱的直方图表示 FTSE 数据的概率密度,并将正态分布作为红色曲线添加到已有图形上:

> hist(ftse, nclass = 20, probability = TRUE)
> lines(ftse, dnorm(ftse, mean = mu, sd = sigma), col = "red")

可以看到,dnorm(x, mean, sd) 使用样本均值和标准差计算数据 x 的概率密度函数(PDF);这称为矩估计法(method-of-moments)

最后,若要计算数据 x 的密度估计,请使用 density(x)。它采用非参数方法构造所谓的核密度估计(KDE),无需对底层分布作出任何假设。

多种图形表明,该数据相较正态分布具有更肥尾。后续练习中,您将学习更好的图形与数值检验方法。

本练习中,您将为 2008–2009 年道琼斯指数的对数收益率拟合一个正态分布,并通过直方图与密度图将数据与拟合分布进行比较。包含道琼斯数据的对象 djx 已加载到您的工作空间。

本练习是课程的一部分

R 中的定量风险管理

查看课程

练习说明

  • 计算 djx 的均值与标准差(sd()),分别赋值给 musigma
  • 绘制 20 个柱的 djx 直方图来表示数据的概率密度。
  • lines()dnorm() 中补全参数,将 djx 的正态密度曲线以红色添加到直方图。
  • 使用 density()djx 绘制核密度估计。
  • 使用与上面相同的 lines() 命令,将 djx 的正态密度曲线以红色添加到 KDE。

交互式实操练习

通过完成这段示例代码来试试这个练习。

# Calculate average and standard deviation of djx
mu <- ___
sigma <- ___

# Plot histogram of djx
___(___)

# Add the normal density as a red line to histogram
lines(___, dnorm(___), col = ___)

# Plot non-parametric KDE of djx
___(___)

# Add the normal density as red line to KDE
lines(___, dnorm(___), col = ___)
编辑并运行代码