使用R拟合分布
几个常用的概率函数介绍
这里,参考R语言实战,以及[Fitting Distribution with R]的附录。
一.认识各种分布的形态
1.1 连续型随机变量的分布
分布,F−分布,Wishart
分布。
先上个图,一睹为快。
以上几个分布之间的关系如以下结构图所示。
[广义线性模型导论3rd edition,p10]
1.1.1 正态分布
正态分布N(μ,σ2)
的密度函数:
正态分布的形态如图。
library(ggplot2)ggtitle
("正态分布密度函数")
正态分布可以衍生出如下的分布。
若Zi∼ i.i.dN(0,1)
,则有
set.seed(123)
data_chisq<-data.frame(x1 = rchisq(200, 10, ncp = 0),
x2 = rchisq(200, 50, ncp =0),
x3 = rchisq(200, 100, ncp = 0))
data_chisq_long<-melt(data_chisq)
## No id variables; using all as measure variables
ggplot(data = data_chisq_long, aes(x = value, colour = variable))+
geom_density()+
ggtitle
("卡方分布密度函数")
1.1.3t−
分布
若,并且Z和X2独立,则有
set.seed(123)
data_t<-data.frame(x1 = rt(200, 10, ncp = 0),
x2 = rt(200, 50, ncp =0),
x3 = rt(200, 100, ncp = 0))
data_t_long<-melt(data_t)
## No id variables; using all as measure variables
ggplot(data = data_t_long, aes(x = value, colour = variable))+
geom_density()+
ggtitle
("t-分布密度函数")
1.1.4F−
分布
若,并且X1和X2相互独立,则有
set.seed(123)
data_f<-data.frame(x1 = rf(200,df1 = 10, df2 = 10, ncp = 0),
x2 = rf(200,df1 = 5, df2 = 3, ncp =0),
x3 = rf(200, df1 = 3, df2 = 5, ncp = 0))
data_f_long<-melt(data_f)
## No id variables; using all as measure variables
ggplot(data = data_f_long, aes(x = value, colour = variable))+
geom_density()+
ggtitle
1.1.5Wishart
分布
Wishart
分布是χ2分布在p
维正态情况下的推广。这里对多维情况暂不展开介绍。
1.2 离散型随机变量的分布
1.2.1 伯努利(Bernoulli
)分布
伯努利分布记为Bernoulli(p)
,只有0和1两种取值。概率测度函数如下:
1.2.2 二项(Binomial
二项分布B(n,p)
的可能取值范围为0,1,...,n。其概率测度函数如下:
ggtitle
("二项分布概率分布图")
1.2.3 负二项(NegativeBinomial
)分布
负二项分布。
1.2.4 几何(Geometric
)分布
1.2.5 泊松(Poission
)分布
1.3 指数分布族及其相互联系
1.3.1 指数分布族
[广义线性模型导论3rd edition,p58]
1.3.2 指数分布
1.3.3 Weibull分布
1.3.4 Beta分布
1.3.5 Gama分布
1.3.6 双指数(DoubleExponential
)分布
1.4 其他分布
1.4.1 均匀(Uniform
)分布
1.4.2 柯西(Cauchy
)分布
1.4.3 对数正态(Lognormal
)分布
1.5 可视化探索的步骤举例
首先,通过直方图,经验累积分布形态等来观察数据的分布形态。
#产生一组服从N(10,2)分布的随机数
二.模型选择
三.参数估计
模拟估计
矩估计
极大似然估计
四.拟合优度指标
五.拟合优度检验
数据分析咨询请扫描二维码
若不方便扫码,搜微信号:CDAshujufenxi
持证人简介:贺渲雯 ,CDA 数据分析师一级持证人,互联网行业数据分析师 今天我将为大家带来一个关于用户私域用户质量数据分析 ...
2025-04-18一、CDA持证人介绍 在数字化浪潮席卷商业领域的当下,数据分析已成为企业发展的关键驱动力。为助力大家深入了解数据分析在电商行 ...
2025-04-17CDA持证人简介:居瑜 ,CDA一级持证人,国企财务经理,13年财务管理运营经验,在数据分析实践方面积累了丰富的行业经验。 一、 ...
2025-04-16持证人简介: CDA持证人刘凌峰,CDA L1持证人,微软认证讲师(MCT)金山办公最有价值专家(KVP),工信部高级项目管理师,拥有 ...
2025-04-15持证人简介:CDA持证人黄葛英,ICF国际教练联盟认证教练,前字节跳动销售主管,拥有丰富的行业经验。在实际生活中,我们可能会 ...
2025-04-14在 Python 编程学习与实践中,Anaconda 是一款极为重要的工具。它作为一个开源的 Python 发行版本,集成了众多常用的科学计算库 ...
2025-04-14随着大数据时代的深入发展,数据运营成为企业不可或缺的岗位之一。这个职位的核心是通过收集、整理和分析数据,帮助企业做出科 ...
2025-04-11持证人简介:CDA持证人黄葛英,ICF国际教练联盟认证教练,前字节跳动销售主管,拥有丰富的行业经验。 本次分享我将以教培行业为 ...
2025-04-11近日《2025中国城市长租市场发展蓝皮书》(下称《蓝皮书》)正式发布。《蓝皮书》指出,当前我国城市住房正经历从“增量扩张”向 ...
2025-04-10在数字化时代的浪潮中,数据已经成为企业决策和运营的核心。每一位客户,每一次交易,都承载着丰富的信息和价值。 如何在海量客 ...
2025-04-09数据是数字化的基础。随着工业4.0的推进,企业生产运作过程中的在线数据变得更加丰富;而互联网、新零售等C端应用的丰富多彩,产 ...
2025-04-094月7日,美国关税政策对全球金融市场的冲击仍在肆虐,周一亚市早盘,美股股指、原油期货、加密货币、贵金属等资产齐齐重挫,市场 ...
2025-04-08背景 3月26日,科技圈迎来一则重磅消息,苹果公司宣布向浙江大学捐赠 3000 万元人民币,用于支持编程教育。 这一举措并非偶然, ...
2025-04-07在当今数据驱动的时代,数据分析能力备受青睐,数据分析能力频繁出现在岗位需求的描述中,不分岗位的任职要求中,会特意标出“熟 ...
2025-04-03在当今数字化时代,数据分析师的重要性与日俱增。但许多人在踏上这条职业道路时,往往充满疑惑: 如何成为一名数据分析师?成为 ...
2025-04-02最近我发现一个绝招,用DeepSeek AI处理Excel数据简直太爽了!处理速度嘎嘎快! 平常一整天的表格处理工作,现在只要三步就能搞 ...
2025-04-01你是否被统计学复杂的理论和晦涩的公式劝退过?别担心,“山有木兮:统计学极简入门(Python)” 将为你一一化解这些难题。课程 ...
2025-03-31在电商、零售、甚至内容付费业务中,你真的了解你的客户吗? 有些客户下了一两次单就消失了,有些人每个月都回购,有些人曾经是 ...
2025-03-31在数字化浪潮中,数据驱动决策已成为企业发展的核心竞争力,数据分析人才的需求持续飙升。世界经济论坛发布的《未来就业报告》, ...
2025-03-28你有没有遇到过这样的情况?流量进来了,转化率却不高,辛辛苦苦拉来的用户,最后大部分都悄无声息地离开了,这时候漏斗分析就非 ...
2025-03-27