生统笔记3-常用的离群之检测方法及R语言代码实现

定量方法

Kurtosis measure(峰度)

  • 统计学中,峰度(Kurtosis)衡量的是随机变量概率分布的峰态。峰度越大,说明该数据系列中的极端值越多。
峰度的计算公式
  • 缺点:这个指标会根据样本量的大小而变化,一般样本量越大峰度也就越大,不同样本量大小的峰度是不可比的。

  • 代码

kurtosis <- function(x){
  scale_x <- scale(x)
  sum(scale_x**4)/length(scale_x)
}

统计检验

  • Grubbs' Test: 这是测试单个异常值时推荐的检验。
  • Tietjen-Moore Test: 这是 Grubbs 检验对多个异常值情况的推广。它的限制是必须准确指定异常值的数量。
  • Generalized Extreme Studentized Deviate (ESD) Test : 此检验仅需要可疑异常值数量的上限,并且是在不知道异常值的确切数量时推荐的测试。

Grubbs' Test

  • 需要使得数据服从正态分布
  • H0 (null hypothesis): There is no outlier in the data.
    HA (alternative hypothesis): There is an outlier in the data.
  • 代码
library(outliers)
grubbs.test(x, type = 10, opposite = FALSE, two.sided = FALSE)

基于近邻度的模型

  • 基于近邻度的模型(proximity-based models),其核心观念是:基于相似度和距离函数,离群点是那些与其他节点孤立的数据点。近邻模式是离群值分析这一领域中最流行的一类方法。其中包括聚类方法(clustering)、基于密度的方法(density-baed)和最近邻方法(nearest-neighbor)。
  • 在最近邻方法中,计算每个数据点到其k个最近邻的距离,并将其视为离群得分。
  • 选定k (k>1),在该方法中,会将那些与大部队相去甚远的小节点群体视为离群点。这是情有可原的,因为在数据生成过程中,如若发生了异常,也可能会产生一系列比较少的相关节点,它们虽然相互之间离的很近,但也应视作异常点。

Local Outlier Factor (LOF, 局部离群因子)

  • LOF算法的基本思想是,根据数据点周围的数据密集情况,首先计算每个数据点的一个局部可达密度,然后通过局部可达密度进一步计算得到每个数据点的一个离群因子,该离群因子即标识了一个数据点的离群程度。LOF值越大,表示离群程度越高,LOF值越小,表示离群程度越低。

  • 代码

wget http://www2.uaem.mx/r-mirror/src/contrib/DMwR_0.3.1.tar.gz
devtools::install_local('./software/DMwR_0.3.1.tar.gz')
DMwR::lofactor(data, k)      # k为根据几个最近邻居计算
  • 缺点:只能根据LOF的大小筛选前n个离群点,不好定cutoff

定性方法

3 sigma原则

如果任何数据点超过标准偏差的3倍,那么这些点很可能是异常值。

3 sigma原则

四分位间距(IQR)

四分位间距(IQR)是我们常画的boxplot中用于定义异常值,低于(Q1 - 1.5 * IQR)或高于(Q3+1.5 * IQR)的观测值定义为异常值。


IQR

参考

https://blog.csdn.net/weixin_36274103/article/details/112717897
https://zhuanlan.zhihu.com/p/385238291
https://www.itl.nist.gov/div898/handbook/eda/section3/eda35h.htm

©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 203,937评论 6 478
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 85,503评论 2 381
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 150,712评论 0 337
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 54,668评论 1 276
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 63,677评论 5 366
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 48,601评论 1 281
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 37,975评论 3 396
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 36,637评论 0 258
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 40,881评论 1 298
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 35,621评论 2 321
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 37,710评论 1 329
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 33,387评论 4 319
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 38,971评论 3 307
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 29,947评论 0 19
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 31,189评论 1 260
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 44,805评论 2 349
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 42,449评论 2 342

推荐阅读更多精彩内容

  • 一、基本概念 异常对象被称作离群点。异常检测也称偏差检测和例外挖掘。 常见的异常成因:数据来源于不同的类(异常对象...
    王尔德的小人书阅读 1,947评论 0 0
  • 一、什么是异常检测 数据中的异常数据通常被认为是异常点、离群点或孤立点,这些数据的特征与大多数数据不一致,呈现出“...
    全村人民的希望阅读 848评论 0 0
  • 1、什么是异常检测 异常检测(Outlier Detection),顾名思义,是识别与正常数据不同的数据,与预期行...
    noob鸽阅读 471评论 0 0
  • 一、异常检测简介 异常检测是通过数据挖掘方法发现与数据集分布不一致的异常数据,也被称为离群点、异常值检测等等。 1...
    算法进阶阅读 2,144评论 0 9
  • 一、离群点是什么? 离群点,是一个数据对象,它显著不同于其他数据对象,与其他数据分布有较为显著的不同。有时也称非离...
    堂堂正正的大号阅读 2,428评论 0 2