R语言笔记 自定义峰度 和 偏度 函数
admin
2023-02-16 05:40:04
0次

> mtcars<- function(x,na.omit=FALSE){

+ if(na.omit)

+ x<-x[!is.na(x)]

+ m<-mean(x)

+ n<-length(x)

+ s<-sd(x)

+ skew<-sum((x-m)^3/s^3)/n

+ kurt<-sum((x-m)^4/s^4)/n-3

+ return(c(n=n,mean=m,stdev=s,skew=skew,kurtosis=kurt))

+ }


vars<- c("mpg","hp","wt")


> sapply(mtcars[vars],mtstat)

               mpg          hp          wt

n        32.000000  32.0000000 32.00000000

mean     20.090625 146.6875000  3.21725000

stdev     6.026948  68.5628685  0.97845744

skew      0.610655   0.7260237  0.42314646

kurtosis -0.372766  -0.1355511 -0.02271075

> 


上一篇：R语言笔记 计算描述性统计量的函数

下一篇：go语言的循环语句

相关内容

热门资讯

邓煜，诗人，数学家，老二次元 来源 视觉中国菲尔兹奖的通知邮件抵达时，普林斯顿大学数学博士、芝加哥大学数学系教授邓煜正在看一本恋爱...
安全生产许可证状态不是“有效”... 工程公司的安全生产许可证在官方系统显示为“其他”状态，有一栏则显示“整改”，却成为高标准农田项目的第...
深化产教融合 推进数智育人 哈... 7月23日，由阿里国际人工智能人才孵化中心（以下简称“阿里国际AITIC”）主办的“智启未来·数智赋...
中国工程院举办学习陈立泉、贲德... 北京7月24日电 (记者 孙自法)中国工程院7月24日在北京举办学习陈立泉、贲德院士科学家精神座谈会...
从成都看低空经济，专访杨军：万... 电子科技大学教授、深思实验室主任杨军接受封面新闻记者专访 封面新闻记者 王一理 摄影 陈光旭 近年来...
蚂蚁集团副总裁张俊杰：希望与成... 7月24日下午，APEC人工智能高级别论坛在成都举行。会议期间，蚂蚁集团副总裁、健康事业群总裁张俊杰...
北大博雅塔为王虹邓煜亮灯 极目新闻记者 陈洋洋7月24日晚，多位北京网友发帖称，在北京大学校园内的博雅塔上，看到了庆祝中国数学...
要和马六甲争高下的陆桥黄了，暴... （文/观察者网 郭光昊、张菁娟）绕开马六甲海峡，在泰国南部的狭窄地峡上“劈开一个口子”。这是每个中国...
AI做科研，走过概念期 2024年，诺贝尔化学奖颁给了谷歌DeepMind的两位人工智能科学家德米斯·哈萨比斯和约翰·杰珀，...
FLUX 3多模态AI模型登场... IT之家 7 月 24 日消息，Black Forest Labs 昨日（7 月 23 日）发布博文...