在R和Python中求峰度和偏度
文章目录
- 一、 偏度
- 1.1 计算公式
- 1.2 自定义函数以及调用
- 二、峰度
- 2.1 计算公式
- 2.2 自定义函数以及调用
- 三、参考资料
一、 偏度
1.1 计算公式
k-阶中心矩为:
uk=1n∑i=1n(xi−x‾)ku_k=\frac1n\sum^n_{i=1}{(x_i-\overline{x})^k} uk=n1i=1∑n(xi−x)k
偏度公式:
g1=n2u3(n−1)(n−2)s3g_1=\frac{n^2u_3}{(n-1)(n-2)s^3} g1=(n−1)(n−2)s3n2u3
1.2 自定义函数以及调用
(一)R语言中实现
编写skew函数
skew <- function(x){n <- length(x) #长度meanx <- mean(x) #均值s <-sd(x) #均方差u <- sum((x-meanx)^3 )/n #3阶中心矩g1 <- (n^2*u)/((n-1)*(n-2)*s^3) #偏度return(g1) #输出结果
}
调用函数以及测试结果
setwd("E:/R")
#测试偏度函数
source('skew.R')
datastudent <- c(11,2,4,5,8,9,10)
skew(datastudent)
得到的结果为: -0.366937
(二)Python中实现
在Python有现成的函数,导入相应的库之后使用skew函数即可求出偏度。
import pandas as pd
import numpy as np
da = [11,2,4,5,8,9,10]
das = pd.Series(da)
das.skew()
得到的结果为:-0.36693703063253225
二、峰度
2.1 计算公式
峰度公式:
g2=n2(n+1)u4(n−1)(n−2)(n−3)s4−3(n−1)2(n−2)(n−3)g_2=\frac{n^2(n+1)u_4}{(n-1)(n-2)(n-3)s^4}-3\frac{(n-1)^2}{(n-2)(n-3)} g2=(n−1)(n−2)(n−3)s4n2(n+1)u4−3(n−2)(n−3)(n−1)2
2.2 自定义函数以及调用
(一)R语言中实现
kurt <- function(y)
{s<- sd(y) #均方差n<- length(y) #长度uk4 <- sum((y-mean(y))^4)/n #4阶中心矩g2 <- (n^2*(n+1)*uk4)/((n-1)*(n-2)*(n-3)*s^4)-((3*(n-1)^2)/((n-2)*(n-3))) #峰度return(g2) #输出结果
}
调用函数以及测试结果
#测试峰度函数
source('kurt.R')
datastudent <- c(11,2,4,5,8,9,10)
kurt(datastudent)
得到的结果为:-1.490657
(二)Python中实现
import pandas as pd
import numpy as np
da = [11,2,4,5,8,9,10]
das = pd.Series(da)
das.kurt()
得到的结果为:-1.490657439446367
三、参考资料
[1].「量学堂-3」统计矩
在R和Python中求峰度和偏度相关推荐
- python随机森林库_随机森林库:R和Python中的不同结果
下面的代码用R和python训练一个随机森林模型.正如您所注意到的,R(1-0.27=0.73)中的精度要比Python中的(0.69)好.此外,特性在R和Python中的重要性是不同的.在 [EDI ...
- Python中求最大值和最小值max()函数、min()函数
[小白从小学Python.C.Java] [Python全国计算机等级考试] [Python数据分析考试必会题] ● 标题与摘要 Python中求最大值和最小值 max()函数.min()函数 ● 选 ...
- 【数模】数据统计中的峰度与偏度
数据统计中的峰度与偏度 峰度: 峰度(peakedness:kurtosis)又称峰态系数.表征概率密度分布曲线在平均值处峰值高低的特征数. 峰度包括正态分布(峰度值=3),厚尾(峰度值>3), ...
- 先马后看!详解线性回归、朴素贝叶斯、随机森林在R和Python中的实现应用!(附代码)...
来源| analyticsvidhya 编译| 火火酱,责编| Carol 出品 | AI科技大本营(ID:rgznai100) 谷歌的自动驾驶汽车和机器人得到了媒体的广泛关注,但是公司真正的未来是在 ...
- python中求差的函数_Python编程基础11:函数和模块
一.函数和模块概述 (一)函数概述 函数可以看成是语句的集合,通过函数调用来执行其包含的语句.函数可以返回一个计算结果,根据每次函数调用的参数,可以返回不同的计算结果.Python利用函数提高代码的重 ...
- python中求最小公约数,python求最大公约数和最小公倍数的简单方法
python求最大公约数和最小公倍数的简单方法 python怎么求最大公约数和最小公倍数 一.求最大公约数 用辗转相除法求最大公约数的算法如下: 两个正整数a和b(a>b),它们的最大公约数等于 ...
- python中求斐波那契数偶数和_[译] 斐波那契数列中的偶数 (Python vs. JavaScript)
斐波那契数列中的偶数 (Python vs. JavaScript) 对于雇主来说,用某种方式来生成斐波那契数列是一道热门的面试题.而求斐波那契数列中的所有偶数便是其热门的变体之一.这里,我将用 Py ...
- Python中求列表list的平均数
Python中如何求列表list的平均数 当列表list中只包含数字时,如何求取它的平均数: from numpy import * a = [52,69,35,65,89,15,34] b = me ...
- python中求差的函数_python3函数取差是什么?如何写代码?
我们在数学中进行加减运算,两个数相减就可以得到差.今天小编就带来一个新问题,函数取差.有的小伙伴会很困惑,函数是一群英文使用,又不是数字怎么还能进行取差呢?这里也算是学习python有趣的地方,函数取 ...
- python中求余_python中取余
Python的基础语法 Python的基础语法 在对Python有了基础的认识之后,下面来了解一下Python的基础语法,看看它和C语言.java之间的基础语法差异. 一.变量.表达式和语句 Pyth ...
最新文章
- Apache Spark源码走读之8 -- Spark on Yarn
- 终端服务器有多种运行模式,云终端的三种工作模式你都知道的吗
- kettle根据参数动态派生列
- 洛谷 P1377 [TJOI2011]树的序 解题报告
- 牛掰本机限速软件appband
- python输出最大值与最小值求成绩平均值_Python第3次作业
- 【SpringCloud】Spring Cloud Config 配置中心
- 代理ip如何使用_为什么在使用代理IP爬虫时会出现超时?
- 来自 Dropbox 的可扩展性设计经验
- 用递归的方式处理数组 把递归方法方法定义到数组的原型上 (这是一次脑洞大开的神奇尝试)...
- react-testing-library
- 【大数据面试宝典】 第一篇 Hadoop 面试题
- 荣耀路由器外界硬盘与win10
- kdj买卖指标公式源码_KDJ胜率极高的买入机会,死记一个指标(附KDJ背离指标源码)...
- python算日期_python中的日期计算
- java自行车怎么调档_变速的自行车如何调档
- 迭代最小二乘拟合椭圆
- Guitar Pro8苹果mac最新版本下载安装教程
- 神威超级计算机配置,神威太湖轻型超级计算机配置和性能列表
- MIDle生命周期详解,以及工作原理
热门文章
- 计算机与信息科学相关教材,AdobeIllustrator实例教程/计算机与信息科学系列规划教材...
- 与报文交换比较,分组交换有哪些优点?
- ios打包发布testflight
- assimp android build,Android assimp编译及引用
- SpringMVC快速上手教程及SSM整合案例
- Python兼职:300-800元/天,各行各业都能运用!
- h5php大转盘抽奖,微信公众号H5抽奖幸运大转盘,PHP带后端,微信抽奖大转盘源码可以独立使用...
- wpsppt设置页码和总页数_wps的ppt页码怎么设置
- 数据结构——图的基本操作
- 谷歌应用内购神器Freedom原理解析