- 博客(38)
- 资源 (2)
- 收藏
- 关注
RSS订阅原创 R——factor()、gl（） 构造因子
因子——表示名义型变量或有序变量。名义变量一般表示类别，如性别，种族等等。有序变量是有一定排序顺序的变量，如职称，年级等等。使用factor()函数f <- factor(x=charactor(), levels, labels=levels, exclude = NA, ordered = is.ordered(x), namax = NA)其中：x 为创建因子的数据，是一个...
2019-12-30 21:09:39
2262
1
原创 数据标准化处理
Z-score 标准化适用范围1.数据的最大值和最小值未知2.有超出取值范围的离群数据经过处理的数据符合标准正态分布，即均值为0，标准差为1，其转化函数为：x∗=（x−μ）/σ标准化后的变量值围绕0上下波动，大于0说明高于平均水平，小于0说明低于平均水平。缺点估算Z-Score需要总体的平均值与方差，但是很难得到，大多数情况下是用样本的均值与标准差替代。对于数据的分布有一定的...
2019-12-28 21:51:15
3848
原创 均方误差——MSE 和标准差 的区别
均方误差——MSE参数估计值与参数真值之差平方的期望值;均方根误差——RMSE平均绝对误差——MAE绝对误差的平均值标准差 ——SDMSE 机器学习 预测时 用，是看看预测的准不准。SD 是看它与平均值之间的差距，看看你这个数据的离散程度。...
2019-12-28 12:43:55
3494
原创 Shapiro-Wilk检验（一般配合QQ图）
检验数据是否符合正态分布类似于线性回归的方法一样，是检验其于回归曲线的残差。使用范围：样本量N<20软件实现：R语言的stats包中的shapiro.test函数该检验原假设：数据服从正态分布（所以要想证明数据不服从正态分布必须拒绝原假设，一般是p值小于0.05）具体数学表达式：一般先用QQ图观察一下，再用这个检验一下。...
2019-12-27 14:26:48
4695
原创 R——输出数据和结果
sink("E:\\桌面\\data.txt",append=FALSE,split=TRUE)a<-c(1,2,3,4,5,6,7,8,9,10)print(a) 用print去输出append=TRUE代表的意思就是在源文件上面追加写入split=TRUE的意思就是不仅仅在文件中输出，也在控制台输出...
2019-12-27 14:00:17
2206
原创 excel————&符号
连接文本时， 文本需要用英文 " " 包含住：=“A”&“B” 返回的是AB=A1&A2 返回的是在表格A1和A2位置连接起来=A1&“-”&A2 返回的是在表格A1的位置的数—A2位置的数...
2019-12-26 12:54:57
623
原创 协方差、协方差矩阵
在机器学习中，理解协方差矩阵的关键在于牢记它计算的是同一个样本不同特征维度之间的协方差，而不是不同样本之间。拿到样本矩阵之后，我们首先要明确一行是样本还是特征维度。一般来说，样本矩阵中一行是一个样本，一列为一个特征维度。所以要按列计算均值（期望），再按行计算出协方差矩阵，把每一行的协方差矩阵相加再除以行数（即样本数），得到样本矩阵的协方差矩阵一、协方差从公式上看，协方差是两个变量与自身...
2019-12-24 15:47:21
3785
原创 共轭 、 共轭转置、共轭矩阵、酉矩阵、正定矩阵、半正定矩阵
共轭复数实数部分相同而虚数部分互为相反数的两个复数。矩阵的共轭转置把矩阵转置后，再把每一个数换成它的共轭复数。自共轭矩阵矩阵中每一个第i 行第j 列的元素都与第j 行第i 列的元素的共轭相等。酉矩阵AH 是A 的共轭转置A叫做酉矩阵...
2019-12-24 13:36:43
25313
原创 matlab————imshow（）
例如有一张图queen是double类型的：（1）直接用imshow(queen)结果——将queen中为0的显示为黑色，大于等于1的显示为白色（2）imshow(queen,[low,high])灰度低于low的像素被显示为黑色，高于high的被显示为为白色，介于low和high之间的像素按比例显示为各个等级的灰色（3） imshow(queen,[ ])函数会将图像矩阵中的min(...
2019-12-23 11:16:28
233
原创 奇异值分解——SVD
奇异值分解(Singular Value Decomposition)一、可以将一个矩阵A分解为三个矩阵的乘积一个正交矩阵U(orthogonal matrix)一个对角矩阵S(diagonal matrix)一个正交矩阵V的转置。（1）正交矩阵：方阵，列（行）向量是标准正交向量（正交的单位向量，内积为0）（2）对角矩阵：除主对角线外的值都是0。可以是方阵、长矩阵（3）实对称...
2019-12-22 23:59:16
1513
原创 R语言基本操作
运行q()——退出R程序tab——自动补全ctrl+L——清空consolerm(list=ls())——删除工作空间的所有对象ESC——中断当前计算文件目录设置setwd（）——设置工作文件目录getwd（）——获取当前工作文件目录list.files()——查看当前文件目录中的文件ls() 和 objects（）——查看当前工作空间中存在的对象（变量）methods(x)...
2019-12-22 16:28:27
2682
原创 R——读取excel
一、本身文件的属性是用的\ ， 但是在R里 要么用\ 要么用/C:\Users\Desktop\sth.txtC:/Users/Desktop/sth.txt二、提取方式（1）从剪贴板里读取选中EXCEL里所需要的内容，复制，在R中输入：data<- read.table(“clipboard”,header=TRUE,sep=’\t’)(2) 转换成文本txt读取data&...
2019-12-22 16:11:58
281
1
原创 R语言——par时一些字母
plot 二维坐标绘图dotchart 点图barplot 条形图hist 直方图pie 饼图points 添加点lines 添加线text 添加文字title 添加标题boxplot 箱线图##画图plot(x~y,xlab="",ylab="",main="",xlim=c(0,45),ylim=c(0,45),pch=18,col=2,cex=5)#xlab/ylab...
2019-12-22 16:02:23
1124
原创 the moran effect ————莫兰效应
环境随机性中的空间相关性可以使大范围内的种群同步，这种现象被称为莫兰效应。在种群生态学中，莫兰定理(或莫兰效应)指出，同一物种的两个不同种群的时间相关性等于它们所处环境的多样性之间的相关性。...
2019-12-21 23:59:43
841
原创 单纯形算法
LP问题： linear programming ——线性规划问题单纯形法的一般解题步骤①把线性规划问题的约束方程组表达成典范型方程组,找出基本可行解作为初始基本可行解.②若基本可行解不存在,即约束条件有矛盾,则问题无解.③若基本可行解存在,从初始基本可行解作为起点,根据最优性条件和可行性条件,引入非基变量取代某一基变量,找出目标函数值更优的另一基本可行解.④按步骤3进行迭代,直到对应...
2019-12-19 14:16:49
136
原创 经验分布函数
例题随机观察总体X，得到一个容量为10的样本值：3.2, 2.5, -2， 2.5, 0, 3, 2, 2.5, 2, 4 。 求X经验分布函数【解 】：把样本值按从小到大的顺序排列为 于是得经验分布函数为说明(1) 经验分布函数Fn(x)是一个阶梯形函数, 当样本容量增大时, 相邻两阶梯的跃度变低, 阶梯宽度变窄, 容易想像, 这样的阶...
2019-12-18 21:58:21
30579
3
原创 R————时间序列————指数预测模型
指数预测模型单指数模型: 拟合的是只有常数水平项和时间点ii处随机项的时间序列，这时认为时间序列不存在趋势项和季节效应。双指数模型 / Holt指数平滑:拟合的是有水平项和趋势项的时序。三指数模型/Holt-Winters指数平滑:拟合的是有水平项、趋势项以及季节效应的时序。可使用以下函数进行拟合指数模型：R中自带的HoltWinters()函数forecast包中的ets()函数 ...
2019-12-17 14:46:02
1613
原创 R语言————par函数
为什么很多代码一开始先写一个par( ) ？【这个函数会“永久性”改变作图设置，而我们有时并不想要这种功能，我们可能希望画下一幅图的时候，之前的参数可以被“还原”回来】【1】在一幅图开始之前先把作图参数保存到一个对象中，比如op = par()【2】在作这幅图的过程中用par()函数任意更改设置以适合需要【3】作完图后，再用par(op)语句把之前保存的参数设置“释放”出来，【4】这样，...
2019-12-16 23:00:00
13138
原创 分位数详解+QQ图
二分位数中位数四分位数把所有数值由小到大排列并分成四等份，处于三个分割点位置的数值就是四分位数。1）第一四分位数(Q1)，等于该样本中所有数值由小到大排列后第25%的数字；2）第二四分位数(Q2)，又称“中位数”，等于该样本中所有数值由小到大排列后第50%的数字；3）第三四分位数(Q3)，等于该样本中所有数值由小到大排列后第75%的数字。第三四分位数与第一四分位数的差距又称四分位距。...
2019-12-16 22:10:52
18668
2
原创 R语言————管道符————%>%
计算向量x的log值的指数函数值的绝对值的整数近似值平常：round(abs(exp((log(x)))))利用管道符号：library(dplyr) #使用dplyr包功能x %>% log( ) %>% exp( ) %>% abs( ) %>% round( )...
2019-12-16 18:41:09
2821
原创 R——Error in plot.new() : figure margins too large——解决方法
使用 windows（） 弹出新的窗口来画。或者拉大一下画图的区域的长和宽。
2019-12-16 15:59:15
1762
原创 R语言————%in%
a<-c(3,6,4)b<-c(1,3,4,5,7,8)a%in%b % a是否在b里面，出来的是逻辑值true or false。which(a%in%b) %a里面的哪些值，第几个值，在b里面。a[which(a%in%b)] %取出a中所有在b中出现过的元素match(a,b) %想要看看a出现在b中的哪个位置，指的是b的位...
2019-12-15 14:22:45
1451
原创 R语言————用箱线图检查异常值
[stats, n, conf, out]= boxplot.stats(x, coef=1.5, do.conf=TRUE, do.out=TRUE) x：所要分析的数值向量coef：盒子的长度，也就是说上图中盒须的长度需要是几倍的盒长（IQR）do.conf: 是否要输出confdo.out： 是否要输出outstats：返回5个元素的值：盒须最小值，盒最小值，中位数，盒最大值，盒...
2019-12-15 14:15:44
12113
原创 excel——如何选中奇数行/偶数行/某一列只选择其数值为偶数/奇数
可以在所要进行筛选的那一列后面添加一列然后用=mod（row（），2）/ =mod（ ，2）row（A2）：返回A2的行数。mod（B,2）：用B的值去除以2 的余数。这样就会出现0101010101…然后，再用筛选功能，筛选出这一列只是1或者只是0的就好了。...
2019-12-14 17:24:21
5332
原创 matlab————导出图片不清楚如何解决
在画出来的图中左上方找：文件————导出设置————左方属性里选渲染——分辨率设为300/600————右侧导出
2019-12-13 22:59:35
8525
1
原创 matlab————在图中标数值
将二维矩阵转化成随着数值大小而变换颜色的图——imagescimagesc(c) %将c这个矩阵转化成图colormap(parula);colorbar %使用的是matlab自带的颜色parula %colorbar表示显示颜色变化的图ax...
2019-12-13 22:53:06
1678
原创 excel技巧——时间复制粘贴后变成数字
我们经常需要用到复制粘贴 而且需要只粘贴文本这时候就会发现 日期变为了数字 如2001-1-1 变成了35555 …解决方案复制粘贴后， 选中粘贴的那一列，选中,并且去找这里：然后就会发现更改过来了。或者：直接选择，右键，设置单元格格式：日期 or 自定义...
2019-12-13 17:53:43
14285
原创 excel技巧——删除空白行
选择数据区域内容最多的一列（避免删除有用信息，因为有的列有部分值没有，但是这列数据依然有效）按F5调出定位功能对话框定位条件选择空值然后在开始菜单中，找到删除删除工作表行...
2019-12-13 17:34:49
105
原创 excel处理两列时间序列
时间序列剔除不一致时间的数据的做法：1.针对时间列用vlookup将前两列的时间序列与后两列匹配的 第二列的数据填入其后。2. 对未匹配值进行删除对于没有匹配上的，自然会出现 #NA，那么就要进行删除。用筛选，选出 #NA，然后选中，删除。 然后去掉筛选，剩下的就是都匹配上的了。这样两个时间序列就把互相不对应的时间删除，并且可以用一列相同的时间了。...
2019-12-13 17:28:34
1641
原创 R语言——STL时间序列分解
stl(x, s.window, s.degree=0, t.window=NULL, t.degree=1, robust=FALSE, na.action=na.fail)s.window————提取季节性时的loess算法 时间窗口的宽度，需要是奇数；s.degree————提取季节性时局部拟合多项式的阶数，需要是0或者1；t.window————提取趋势性时的loess算法 时间窗口...
2019-12-12 22:31:49
13460
5
原创 R语言笔记——ts函数（不定时更新）
1.ts（）函数将其转化为时间序列时间序列对象是一种专为时间序列分析而设计的对象类型，其中包括两个维度，一个是描述指标的数值，还有一维是时间。时间序列对象和一般数值型向量类似，只不过是加了一个时间的描述。在R语言中可以使用ts(数据向量，frequency=表示将时间分开的时间间隔，start=c(第一个数据所表示的年，月))。sales.data<-read.csv(“citysa...
2019-12-12 19:17:00
21087
1
原创 STL笔记
时间序列数据的变化是众多复杂因素共同作用的结果，影响因素主要包括长期趋势、季节变动、周期变动和不规则变动，其中不规则变动又包含随机变动和突发变动，这一影响因素往往难以测定，一般作为干扰项处理。时间序列分解能够帮助分析者去除其他因素的影响，单纯分析某一确定性因素影响下的序列分布规律。目前常用的分解模型有加法模型和乘法模型，如果季节变动的幅度以及趋势周期的波动不随时间变化或者变化幅度不大，适合采用...
2019-12-12 16:50:51
465
原创 switch-case结构
switch value %value为需要进行判决的标量或字符串 case test1 commands1 %如果value的值等于test1，执行commands1，并**结束此结构** case test2 commands2 %如果value的值等于test2，执行commands2，并**结束此结构** .....
2019-12-06 22:01:33
221
原创 交叉验证crossvalind
例1Data = rand(9,3);%创建维度为9×3的随机矩阵样本indices = crossvalind(**'Kfold'**, 9, 3);%将数据样本随机分割为3部分for i = 1:3 %循环3次，分别取出第i部分作为测试样本，其余两部分作为训练样本 test = (indices == i); %这里用==生成的是逻辑类型，（0,1表示的） train...
2019-12-06 14:25:18
424
原创 交叉验证——CrossValidation
CrossValidation【CV】基本思想将原始数据(dataset)分组：一部分为训练集(train set)一部分为验证集(validation set)用训练集对分类器进行训练,再利用验证集来测试训练得到的模型(model),以此来做为评价分类器的性能指标.常见方法Hold-Out Method将原始数据随机分为两组,一组做为训练集,一组做为验证集,利用训练集训练分类器,...
2019-12-06 13:31:09
189
原创 递归思想
递归的思想以此类推递归的两个条件可以通过递归调用来缩小问题规模，且新问题与原问题有着相同的形式。（自身调用）存在一种简单情境，可以使递归在简单情境下退出。（递归出口）举例factorial(n)=n*factorial(n-1)，其中n为非负整数,且0!=1,1!=1public static long factorial(int n) throws Exception { ...
2019-12-06 10:11:14
107
原创 插值函数interp1
插值用函数f (x)在已知的一些点的函数值，作出适当的特定函数，来估计其他点的值。注意事项所有的插值方法都要求x是单调的，并且xi不能够超过x的范围。用法yi= interp1(x,y,xi,‘method’)其中x，y为插值点，yi为在被插值点xi处的插值结果；x,y为向量，**‘method’**表示采用的插值方法：**‘nearest’**是最邻近插值**‘linear’*...
2019-12-05 19:27:02
1855
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人