米斯特黄-CSDN博客

原创聚类分析之K均值聚类

1.K均值简介麦克奎因1967年提出了K-均值聚类，该方法的思想是把每个样品聚集到其最近形心（均值）类中。在它的最简单说明中，这个过程由下列三步所组成：（1）把样品粗略分为K个初始类。（2）进行修改，逐个分派样品到其最近均值的类中（通常用标准化数据或者非标准化数据计算欧式距离）。重新计算接受新样品的类和失去样品的类的形心（均值）（3）重复第（2)步，直到各类无元素进出。2.K-均值法所需...

2019-07-14 11:08:11 2393

原创朴素贝叶斯的python实现

之前都在写关于R的一些内容，最近研究朴素贝叶斯比较多，就来更一个python的方法实现朴素贝叶斯吧，更的内容是关于离散型变量的分类，利用UCI数据集上的car数据来做，至于理论部分，请观看我的另一篇文章，>https://blog.csdn.net/heirenmin/article/details/94611260直接上代码，有什么问题，请在下方留言：#coding:utf-8#朴素...

2019-07-10 19:32:27 312

1.系统聚类简介系统聚类法是一种常用的聚类方法，假设有n个样品，每个样品有p个变量，系统聚类的基本步骤为：先将每个个体（样品或变量）各自看成一类，总共有r类（如果时Q型聚类，则r=n;如果时R型聚类，则r=p),根据个体间的相似程度（距离、相关系数等）将r类个体中最相似的一类合并称一类，得到r-1类，再在这r-1类中找出最相似的两类合并，得到r-2类，依次进行下去。对于Q性聚类情形，设Gs与G...

2019-07-07 13:49:27 7477 3

原创朴素贝叶斯的R实现

朴素贝叶斯分类器采用了“属性条件独立性假设”：对已知类别，假设属性独立。也就是说每个属性独立的对分类结果发生影响。根据贝叶斯公式有：又由于各个属性相互独立，所以有：这样朴素贝叶斯的后验概率就可以写成：由于对于所有的类别来说P(x)是相同的因此贝叶斯的判定准则是：这是朴素贝叶斯分类器的表达式令Dc表示训练集D中c类样本组成的集合，若有充足的独立同分布样本，则可容易地估计出类先验概...

2019-07-04 14:00:32 1575

转载贝叶斯统计R语言的包

http://cran.r-project.org/web/views/Bayesian.html一般模型========arm包: 包括使用lm,glm,mer,polr等对象进行贝叶斯推断的R函数BACCO: 随机函数的贝叶斯分析. 包含3个子包: emulator, calibrator, and approximator, 进行贝叶斯估计和评价计算机程序.bayesm: 市场与微经...

2019-07-03 18:44:55 11015

原创多元线性回归模型

多元线性回归模型通常用来研究一个应变量依赖多个自变量的变化关系，如果二者的以来关系可以用线性形式来刻画，则可以建立多元线性模型来进行分析。1.模型简介1.1模型的结构多元线性回归模型通常用来描述变脸y和x之间的随机线性关系，即：如果对y和x进行了x次观测，得到n组观察值yi,x1i,…,xki(i=1,2,…,n),他们满足一下关系式：1.2模型参数的检验在正态假定下，如果X是列...

2019-07-02 19:22:12 114436 3

原创误差修正模型

误差修正模型通常作为协整模型的补充模型出现。协整模型度量解释的是序列间的长期关系，而误差修正模型解释的是序列之间的短期关系。1.模型简介误差修正模型，简称ECM模型，构造原理如下：假设非平稳响应序列{yt}与非平稳输入序列{xt}之间具有协整关系，即：则回归残差序列为平稳序列：在第一个式子前同时减去yt-1，则有：将上面的式子整理一下得到：上式就是ECM 模型的表达式，说明误...

2019-07-02 10:14:18 21291 1

原创协整模型的R实现

在现实生活中，有些序列自身的变化虽然是非平稳的，但序列和序列之间却有着非常密切的长期均衡关系。比如农村家庭人均纯收入对数序列和人均生活消费支出对数序列，单整分析两个序列都是非平稳的，但是他们之间都有很强的相关性。这就可以建立协整模型对序列进行分析处理。1.模型简介假定自变量序列{X1},{x2},…{xk},响应变量序列为{yt},构造回归模型：假定回归残差序列平稳，我们就称响应序列与自...

2019-07-01 18:59:43 11182 6

原创平稳多元序列建模的R实现

1.模型简介对于平稳多元时间序列建模的思想是：假定响应序列{yt}和输入序列{x1t},{x2t},…,{xkt}均平稳，首先构造响应序列和输入序列的回归模型：因为响应序列和输入序列皆是平稳序列，那么其残差序列也是平稳序列：上述模型称为动态回归模型，简记ARIMAX.2.建模步骤（1）读入数据，并绘制其时序图，并判断其平稳性（2）绘制其输出序列自相关图，偏自相关图，并利用其自相关图...

2019-06-28 20:56:29 2151 2

原创季节加法模型

1.模型简介季节加法模型是指序列中季节效应和其他效应之间是加法关系，可以表示为：这时，各种效应信息的提取都非常容易。通常简单的周期步长差分即可将序列中的季节信息提取完毕，提取季节信息和趋势信息之后的残差序列就是一个平稳序列，可以用ARMA模型进行拟合。所以简单季节模型实际上就是通过趋势差分、季节差分将序列将序列转化为平稳序列，再对其进行拟合。它的模型结构通常如下：式中2.建模过程（...

2019-06-26 13:49:09 6735 1

原创 Holt两参数指数平滑法的R实现

指数平滑法是一种简单的高效的时间序列预测方法之一，简单指数平滑法只适合预测没有长期趋势和季节效应的时间序列，并且只能预测未来一期的值。对于具有长期趋势的，但是没有季节效应的时间序列，就比较适合利用Holt两参数指数平滑法来建模预测。1.模型简介Holton两参数指数平滑法适用于对含有线性趋势的序列进行修匀。他的基本结构为：2.建模步骤（1)读入时间序列（2)进行两参数指数平滑（3...

2019-06-25 13:44:48 10995 2

原创 ARCH模型的R语言实现

之前更新的模型都是假定方差齐性的条件，现在假定时间序列具有异方差性，那么之前的模型就不在适用。1.ARCH模型简介ARCH模型，全称为自回归条件异方差模型，有时简称条件异方差模型。该模型的构造原理如下：假设历史数据已知的情况下，零均值、纯随机序列残差序列具有异方差性：在正态分布的假定下，有：异方差等价与残差平方的均值：使用残差平方序列的自相关系数可以考虑异方差函数的自相关性。...

2019-06-24 16:59:11 26941 8

转载 R语言时间序列函数整理

【包】library(zoo) #时间格式预处理library(xts) #同上library(timeSeires) #同上library(urca) #进行单位根检验library(tseries) #arma模型library(fUnitRoots) #进行单位根检验library...

2019-06-24 08:57:30 6148 1

原创疏系数模型的R实现

疏系数模型是ARIMA模型的一种特殊情况，原来ARIMA模型中有一部分系数缺失了，那么这个模型就是疏系数模型。1.模型简介如果是自回归部分有缺失的系数，那么该疏系数模型就可以记为：ARIMA((p1,…pm),d,q)如果只是是移动平均部分有缺失的系数，那么该疏系数模型就可以简记为：ARIMA（p,d,(q1,…1n))如果自回归和移动平均部分都有系数缺失的话，该模型就可以简记为：A...

2019-06-23 14:48:24 6992 1

原创残差自回归模型的R实现

时间序列经常将非平稳时间序列分解为趋势部分，季节因素部分和随机部分，可以表示为：这里的Tt为长期趋势拟合，St为季节效应拟合。残差自回归是回归模型与ARMA模型的组合模型，由于回归模型对时间序列进行拟合是，序列中包含的信息可能不太充分，在拟合回归模型之后对其残差序列进行自相关性检验，若残差序列具有明显的自相关性，那么就需要对残差序列进行拟合。１．利用时间序列拟合回归模型对时间序列的趋势效...

2019-06-23 11:07:47 9566 4

原创 ARIMA模型的R实现

1.模型简介ARIMA模型，全称差分自回归移动平均模型。可以看作是差分运算和ARMA模型的结合。主要是针对具有线性趋势的非平稳的时间序列，来进行建模预测。ARIMA模型的形式如下：Φ(B)∇

2019-06-22 17:33:03 14713 8

原创关于ARMA模型的R语言实现

新手一枚，和大家一起学习R，以后基本每周都会更新1到2篇关于数据预测处理的模型和方法，希望和大家一起学习，一起成长。本周首先更新的是用R来实现ARMA模型。时间序列的模型，基本上都要建立在平稳的序列上，这里我们将来了解下ARIMA模型，以及其实现的R代码。　　ARMA（p,q)模型，全称移动平均自回归模型，它是由自回归（AR）部分和移动平均（MA）部分组成的，所以称之为ARMA模型。进行AR...

2019-06-22 15:03:07 36907 12

heirenmin的博客