数据探索

  • 数据挖掘中特别的列
    • 只有一种值的列:
    • 列的值各不相同:
    • 派生变量:
  • 异常值分析
    • 3 σ 3{\sigma} 3σ原则
    • 箱线图分析
    • 一致性分析
  • 数据特征分析
    • 分布分析
      • 定量数据的分布分析
      • 定性数据分析
    • 统计量分析
      • 集中趋势度量
      • 离中趋势度量
      • 贡献度分析
    • 相关性分析
      • 计算相关系数
  • matlab主要数据探索函数
    • 统计特征函数
    • 统计作图函数
      • pie
      • boxplot
      • errorbar
        • 参考文献

数据挖掘中特别的列

列或者字段代表了每条记录里的数据。

只有一种值的列:

单值的列对于区分不同行不包含任何信息,它们缺乏任何信息内容,在数据挖掘中应当被忽略。

列的值各不相同:

另一个极端是列在每一行上几乎都不相同,由于每一个记录上这些数值几乎不同,所以它们妨碍了数据挖掘算法从不同行间找出规律。

派生变量:

如:电话号码和地址包含地理信息;
发动机的识别号码包含了生产年份、制造商、型号和原产地等;
客户号码也可能知道哪些是最新客户;
重要的特性应该作为派生变量提取出来,而忽略原始列。

异常值分析

异常值是指样本中的个别值,其数值明显偏离其余的观测值。异常值也称为离群点。异常值分析也称为离群点分析。

3 σ 3{\sigma} 3σ原则

如果数据服从正态分布,在 3 σ 3{\sigma} 3σ原则下,异常值被定义为一组测定值中与平均值的偏差超过三倍标准差的值。在正态分布的假设下,距离平均值 3 σ 3{\sigma} 3σ之外的值出现的概率为 P ( ∣ x − μ ∣ > 3 σ ) < = 0.003 P(|x-{\mu}|>3{\sigma})<=0.003 P(∣x−μ∣>3σ)<=0.003

箱线图分析

异常值通常被定义为小于 Q L − 1.5 I Q R Q_L-1.5IQR QL​−1.5IQR或大于 Q U + 1.5 I Q R Q_U+1.5IQR QU​+1.5IQR的值。 Q L Q_L QL​被称为下四分位数, Q U Q_U QU​被称为上四分位数。 I Q R IQR IQR表示四分位数间距。是上四分位数与下四分位数的插值。

function []=loss_detecting(sales)
%% 缺失值检测,并打印结果
rows = size(sales,1);
%% 缺失值检测 并打印结果
nanvalue = find(isnan( sales));
if isempty(nanvalue) %  没有缺失值disp('没有缺失值!');
elserows_ = size(nanvalue,1);disp(['缺失值个数为:' num2str(rows_) ',缺失率为:' num2str(rows_/rows) ]);
end

异常值检测函数:

function [] = abnormal_detecting(sales)
%% 异常值检测
% 箱形图上下界
q_= prctile(sales,[25,75]);
p25=q_(1,1);
p75=q_(1,2);
upper = p75+ 1.5*(p75-p25);
lower = p25-1.5*(p75-p25);
upper_indexes = sales(sales>upper);
lower_indexes = sales(sales<lower);
indexes =[upper_indexes;lower_indexes];
indexes = sort(indexes);
% 箱形图
figure
hold on;
boxplot(sales,'whisker',1.5,'outliersize',6);
rows = size(indexes,1);
flag =0;
for i =1:rowsif flag ==0text(1+0.01,indexes(i,1),num2str(indexes(i,1)));flag=1;elsetext(1-0.017*length(num2str(indexes(i,1))),indexes(i,1),num2str(indexes(i,1)));flag=0;end
end
hold off;
disp('数据缺失值及异常值检测完成!');

一致性分析

该问题主要集中在重复存放的数据未能进行一致性更新造成的。

数据特征分析

分布分析

能揭示数据的分布特征和分布类型。

定量数据的分布分析

对于定量、变量而言,选择组数和组宽是做频率分布分析时最主要的问题。
一般按照以下的步骤进行:
(1)求极差;
(2)决定组距和组数;
(3)决定分点;
(4)列出频率分布表;
(5)绘制频率分布直方图;
遵循的原则有:
(1)各组之间必须是相互排斥的;
(2)各组必须将所有的数据包含在内;
(3)各组的组宽最好相等。

定性数据分析

画出条形图和饼图。

统计量分析

集中趋势度量

(1)均值
一般采用加权均值;
截断均值是去掉高、低极端值后的平均数。
(2)中位数
将某一数据集 x : { x 1 , x 2 , . . . , x n } x:{\{x_1,x_2,...,x_n}\} x:{x1​,x2​,...,xn​}从小到大排序: { x ( 1 ) , x ( 2 ) , . . . , x ( n ) } \{x_{(1)},x_{(2)},...,x_{(n)}\} {x(1)​,x(2)​,...,x(n)​}。
当n为奇数时, M = x ( n + 1 2 ) M=x_{(\frac{n+1}{2})} M=x(2n+1​)​
当n为偶数时, M = 1 2 ( x ( n 2 ) + x ( n + 1 2 ) ) M=\frac{1}{2}(x_{(\frac{n}{2})}+x_{(\frac{n+1}{2})}) M=21​(x(2n​)​+x(2n+1​)​)
(3)众数
众数是指数据集中出现最频繁的值。

离中趋势度量

(1)极差
(2)标准差
(3)变异系数
变异系数度量标准差相对于均值的离中程度。
C V = s x ‾ ∗ 100 CV=\frac{s}{\overline{x}}*100% CV=xs​∗100
(4)四分位数间距:
其值越大,说明数据的变异程度越大,反之,说明数据的变异程度越小。

% 众数
mode_ = mode(sales);
% 极差
range_ = range(sales);
% 四分位数间距
q1=prctile(sales,25);
q2=prctile(sales,75);
distance=q3-q1;

贡献度分析

贡献度分析又称帕累托分析,它的原理是帕累托法则,又称20/80定律。比如:对于一个公司来讲,80%的利润常常来自于20%最畅销的产品,而其它80%的产品只产生了20%的利润。
帕累托分析也可以称为累计占比分析。

function [] = dish_pareto(num)
%% 帕累托图作图
rows = size(num,1);
hold on;
% 计算累计系数
yy_ = cumsum(num(:,end));
yy=yy_/yy_(end)*100;
[hAx,hLine1,hLine2]=plotyy(1:rows,num(:,end),1:rows,yy,'bar','plot');
set(hAx(1),'XTick',[])%去掉x轴的刻度
set(hLine1,'BarWidth',0.5);
set(hAx(2), 'XTick', 1:rows);
set(hAx(2),'XTickLabel',raw(2:end,2));
ylabel(hAx(1),'盈利:元') % left y-axis
ylabel(hAx(2),'累计百分比:%') % right y-axis
set(hLine2,'LineStyle','-')
set(hLine2,'Marker','d')
% 标记 80% 点
index = find(yy>=80);
plot(index(1),yy(index(1))*100,'d', 'markerfacecolor', [ 1, 0, 0 ] );
text(index(1),yy(index(1))*93,[num2str(yy(index(1))) '%'] );
hold off;
disp('餐饮菜品盈利数据帕累托图作图完成!');

相关性分析

计算相关系数

(1)Pearson相关系数
一般用于分析两个连续型变量之间的关系。

相关系数的取值范围:-1<=r<=1

(2)Spearman秩相关系数
Pearson线性相关系数要求连续变量的取值服从正态分布。不服从正态分布的变量,分类或等级变量至今的关联性可采用Spearman秩相关系数。
计算公式:

对两个变量成对的取值分别按照从小到大(或者从大到小的顺序)编秩。 R i R_i Ri​代表 x i x_i xi​的秩次, Q i Q_i Qi​代表 y i y_i yi​的秩次。 R i − Q i R_i-Q_i Ri​−Qi​为 x i , y i x_i,y_i xi​,yi​的秩次之差。

(3)判定系数
它是相关系数的平方,用 r 2 r^2 r2表示;用来衡量回归方程对y的解释程度。 r 2 r^2 r2越接近于1,表明x与y之间的相关性越强; r 2 r^2 r2越接近于0,表明两个变量之间几乎没有直线相关关系。

matlab主要数据探索函数

统计特征函数

MATLAB主要的统计特征函数

函数名 函数功能
mean() 计算数据样本的算术平均数
geomean() 计算数据样本的几何平均数
var() 计算数据样本的方差
std() 计算数据样本对的标准差
corr() 计算数据样本的spearman(Pearson)相关系数矩阵
cov() 计算数据样本的协方差矩阵
moment() 计算数据样本的指定阶中心距

几点说明:
var
使用格式:v=var(X),计算样本X的方差v。若X为向量,则计算向量的样本方差;若X为矩阵,则v为X的各列向量的样本方差构成的行向量。
corr
使用格式:R=corr(x,y,‘name’,‘value’),计算列向量x,y的相关系数矩阵R。

name value 说明
Type pearson 皮尔森相关系数
Type Kendall 卡德尔系数
Type Spearman 斯皮尔曼系数
Rows all 全部数据,默认选项
Rows complete 只使用没有缺失值的行
Rows pairwise 计算R(i,j)只使用第i和j列中没有缺失值的数据

(6)cov
使用格式:R=cov(X),计算样本X的协方差矩阵R。样本X可为向量或矩阵。当X为向量时,R表示X的方差;当X为矩阵时,cov(X)计算方差矩阵。
R=cov(x,y)。函数等价于cov([x,y])。参数x,y为长度相等的列向量。
(7)moment
使用格式:m=moment(X,order)计算样本X的order阶次的中心矩m,参数order为正整数。样本X可为向量、矩阵或多维数组。
说明:一阶中心矩为0,二阶中心矩为用除数n得到的方差。其中n为向量X的长度或矩阵X的行数。

统计作图函数

作图函数名 作图函数功能
pie() 绘制饼形图
hist() 绘制二维条形直方图
boxplot() 绘制样本数据的箱线图
errorbar() 绘制误差条形图

pie

使用格式:pie(X).绘制矩阵X中非负数据的饼形图。若X中非负元素和小于1,则函数仅画出部分饼形图,且非负元素X(i,j)的值直接限定饼形图中扇形大小;若X中非负元素和大于等于1,则非负元素X(i,j)代表饼形图中扇形大小通过X(i,j)/Y的大小来决定,其中Y为矩阵X非负元素和。
例题:

x=[1,3,1.5,4,1.5];
explode=[1,0,0,0,0];% 第一个元素为1表示饼图中第一部分分离出来
pie(x,explode); % 画出饼形图

boxplot

使用格式:boxplot(X,notch),绘制矩阵样本X的箱型图。notch=1时,绘制矩阵样本X的带刻槽的凹盒图。参量notch=0时,绘制矩阵样本X的无刻槽的矩形箱型图。

x1=normrnd(4,5,100,1);
x2=normrnd(8,6,100,1);
boxplot([x1,x2],1);

errorbar

使用格式:errorbar(X,Y,L,U),绘制误差条形图。参量X,Y,L,U必须为同型向量或矩阵。若同为向量则在点(X(i),Y(i))画出向下长为L(i),向上长为U(i)的误差棒;若同为矩阵,则在点(X(i,j),Y(i,j))处画出向下长为L(i,j),向上长为U(i,j)的误差棒。

x=0:pi/10:pi;
y=2*x.*sin(x);
e=std(y).*ones(size(x)); % 产生误差棒长度
errorbar(x,y,e,e);

参考文献

张良均等. MATLAB数据分析与挖掘实战. 机械工业出版社,2015.

数据挖掘与分析(基于MATLAB)——数据探索相关推荐

  1. 《对冲基金建模与分析基于MATLAB》简介及PDF下载

    转 <对冲基金建模与分析--基于MATLAB>简介及PDF下载 内容简介 本书是关于用MATLAB对对冲基金进行建模和分析的入门读物.在对对冲基金的基本概念.分类.相关工具和指标系统介绍的 ...

  2. 决策树模型回归可视化分析_商业分析过程与数据探索

    通常,商业模型以历史数据构建,一旦对模型精确度感到满意,就可以将模型部署到生产中进行实时数据分析. 第一阶段,理解商业问题 在项目初始阶段,重点是从商业角度理解商业问题,商业目标和商业要求,然后将这种 ...

  3. 基于matlab的ofdm系统仿真及性能分析,基于MATLAB的OFDM系统性能仿真

    基于MATLAB的OFDM系统性能仿真(论文15000字) 摘要:本文介绍OFDM系统的基本原理与实现方法,以及移动多径环境下使用OFDM的优势,详细阐述了OFDM系统中符号间干扰ISI.载波间干扰I ...

  4. 扫描电镜图片的Matlab分析,基于Matlab的图像增强与复原技术在SEM图像中的应

    0引  言 根据国内外的相关文献,研究和发展图像处理工具,改善图像质量是当今研究的热点.图像增强与复原是一种基本的图像处理技术.其按照特定的需要突出一幅图像中的某些信息或强化某些感兴趣的特征,将原来不 ...

  5. 基于matlab的光学薄膜特性分析,基于matlab的光学薄膜特性分析

    <基于matlab的光学薄膜特性分析>由会员分享,可在线阅读,更多相关<基于matlab的光学薄膜特性分析(17页珍藏版)>请在金锄头文库上搜索. 1.目录1 技术指标 .12 ...

  6. 基于matlab的光学薄膜特性分析,基于matlab的光学薄膜特性分析.doc

    目录 1 技术指标1 2 基本原理1 2.1 单层光学薄膜1 2.2 多层光学薄膜及矩阵理论2 3 基于matlab的光学薄膜特性分析4 4 结论9 5 参考文献9 6 附录:基于matlab的光学薄 ...

  7. 基于matlab的对称振子天线分析,基于Matlab的对称振子天线分析

    天线方面文献 第27卷第2期 2009年3月泉州师范学院学报(自然科学)JournalofQuanzhouNormalUniversity(NaturalScience)Vol.27No.2Mar.2 ...

  8. matlab随机信号分析,基于MATLAB的随机信号分析方法.ppt

    <基于MATLAB的随机信号分析方法.ppt>由会员分享,可在线阅读,更多相关<基于MATLAB的随机信号分析方法.ppt(31页珍藏版)>请在人人文库网上搜索. 1.基于MA ...

  9. matlab传热分析,基于matlab的超临界流体对流传热可视化分析方法

    基于matlab的超临界流体对流传热可视化分析方法 [技术领域] [0001]本发明涉及一种基于MATLAB的管内超临界流体对流传热可视化分析方法,属于MATLAB应用技术领域. [背景技术] [00 ...

  10. matlab方波假频现象分析,基于MATLAB 的信号时域采样及频率混叠现象分析

    龙源期刊网 http://www.doczj.com/doc/7d04774eae45b307e87101f69e3143323968f58a.html 基于MATLAB 的信号时域采样及频率混叠现象 ...

最新文章

  1. bootstrap-fileinput 使用
  2. getOwnPropertyNames() keys()
  3. rdf mysql持久化l_Redis进阶(数据持久化RDF和AOF)
  4. .NET中多线程调试的一个办法
  5. 13-Mybatis 注解开发
  6. Qt——P13 Q4版本信号槽连接
  7. mysql5.7.13编译安装_编译安装MySQL-5.7.13
  8. gcc 源码下载地址
  9. cactiez v11使用配置mysql_安装cactiez v11对windows和linux系统进行监控
  10. 正态性检验(Normality test)
  11. Win 10 + Ubuntu 18.04双系统 卸载Ubuntu
  12. 什么是组合(Composition),有什么作用
  13. UFS UIC Layer: MIPI UNIPRO
  14. outlook客户端不显示table表格边框线的问题
  15. 高效能人士的七个习惯——由内而外全面造就自己
  16. miui10android9,MIUI 10稳定版刚到 基于安卓9的MIUI10 8.9.7开发版也来
  17. 数据结构几类排序的总结和完整代码 待续。。
  18. javac.exe 、 java.exe、 javaw.exe 、 javaws.exe有什么区别?
  19. 年化超额收益计算matlab,超额年化收益率的计算 看了一下雪球上一些雪友贴出的基金或大V的历年收益率,发现一些人统计的不对,想说一下: 1.超额收益率 某基金当年收益率为... - 雪球...
  20. 利用html创建新闻页面

热门文章

  1. 已知空间中的三点 求三角形面积_三角形的面积公式八叙
  2. 电子警察原理大揭秘 注意!不被电子警察拍照的秘密(图解)
  3. 安装过千月(BlueSoleil)蓝牙软件后,网络图标上就一直显示红叉,但却能上网
  4. 黑莓9780(BB9780)UD Caller 显示坐标
  5. IST:Iterative Shrinkage/Thresholding和Iterative Soft Thresholding
  6. 明远智睿A40i嵌入式核心板助力智能公交系统构建
  7. 遥感影像正射矫正及高分二号遥感影像获取
  8. Metasequoia for Mac(水杉3D建模器)
  9. 三星电子股东大会今日召开 股东或借机要求把李在镕干掉?
  10. 算法竞赛入门(2)学习笔记——循环结构程序设计