`
jobar
  • 浏览: 341386 次
  • 性别: Icon_minigender_1
  • 来自: 北京
社区版块
存档分类
最新评论

R与t分布(1) 概率密度函数

 
阅读更多
在概率论和统计学中,学生t-分布(Student's t-distribution),可简称为t分布。应用在估计呈正态分布的母群体之平均数。它是对两个样本均值差异进行显著性测试的学生t检定的基础。学生t检定改进了Z检定(Z-test),因为Z检定以母体标准差已知为前提。虽然在样本数量大(超过30个)时,可以应用Z检定来求得近似值,但Z检定用在小样本会产生很大的误差,因此必须改用学生t检定以求准确。

在母体标准差未知的情况下,不论样本数量大或小皆可应用学生t检定。在待比较的数据有三组以上时,因为误差无法压低,此时可以用变异数分析(ANOVA)代替学生t检定。

t分布的推导由英国人威廉·戈塞特(Willam S. Gosset)于1908年首先发表,当时他还在爱尔兰都柏林的吉尼斯(Guinness)啤酒酿酒厂工作。酒厂虽然禁止员工发表一切与酿酒研究有关的成果,但允许他在不提到酿酒的前提下,以笔名发表t分布的发现,所以论文使用了“学生”(Student)这一笔名。之后t检定以及相关理论经由罗纳德·费雪(Sir Ronald Aylmer Fisher)的发扬光大,为了感谢戈塞特的功劳,费雪将此分布命名为学生t分布(Student's t)。
参见维基百科
T的概率密度函数是
v  等于n − 1。 T的分布称为t-分布。参数v 一般被称为自由度。

t分布 - 特征
1.以0为中心,左右对称的单峰分布;
2.t分布是一簇曲线,其形态变化与n(确切地说与自由度ν)大小有关。自由度ν越小,t分布曲线越低平;自由度ν越大,t分布曲线越接近标准正态分布(u分布)曲线
t分布,在概率统计中,在置信区间估计、显著性检验等问题的计算中发挥重要作用。
set.seed(1000)
x<-seq(-5,5,length.out=1000)
y<-dt(x,1,0)

plot(x,y,col="red",xlim=c(-5,5),ylim=c(0,0.5),type='l',
     xaxs="i", yaxs="i",ylab='density',xlab='',
     main="The T Density Distribution")

lines(x,dt(x,5,0),col="green")
lines(x,dt(x,5,2),col="blue")
lines(x,dt(x,50,4),col="orange")

legend("topleft",legend=paste("df=",c(1,5,5,50)," ncp=", c(0,0,2,4)), lwd=1, col=c("red", "green","blue","orange"))


结果如下:


  • 大小: 86.2 KB
分享到:
评论

相关推荐

    mvnfast:多元正态分布的快速方法

    mvnfast :快速多元正态分布和学生t分布mvnfast R软件包提供了与多元正态分布和Student t分布有关的高效计算工具。 由于通过Rcpp \ RcppArmadillo软件包使用了C ++代码并通过OpenMP API进行了并行化,因此这些工具...

    beta分布的概率密度的matlab代码-Matlab2R_Pitfalls:将Matlab代码转换为R代码的陷阱

    beta分布的概率密度的matlab代码将Matlab代码转换为R的注意事项和陷阱 这是我将Matlab代码转换为R代码时遇到的陷阱。 这是一个很粗糙的版本。 稍后我会添加详细信息 var1 &lt;-var2 Matlab:var1 &lt;(-var2),...

    可靠性复习精简版+老师给的计算题总结

    1. 可靠性的定义产品(系统)在规定的条件下和规定的时间内,完成规定功能的能力。 2. 可靠度:产品在规定条件下和规定时间内完成规定功能的概率。...4. 失效概率密度函数 :不可靠度函数F(t)随时间的变化率。

    R语言与网站分析

    sample 542.10.4 重复序列 rep 552.10.5 概率分布 552.11 控制流 582.11.1 分支语句 582.11.2 循环语句 592.12 运算符、函数和过程 612.12.1 常用运算符 612.12.2 函数 622.12.3 过程 632.13 数据的读写等...

    Python实现直方图均衡基本原理解析

    通过一个变换,将输入图像的灰度级转换为`均匀分布`,变换后的灰度级的概率密度函数为 $$P_s(s) = \frac{1}{L-1}$$ 直方图均衡的变换为 $$s = T(r) = (L-1)\int_0^r {P_r(c)} \,{\rm d}c $$ $s$为变换后的灰度级...

    狭义数据分析1.pdf

    array(0.), array(0.)) 利⽤公式 ss.norm.pdf(0.0)# 带⼊具体的x值 # Out[14]: 0.3989422804014327 ss.norm.pdf(0.0) Out[14]: 0.3989422804014327 ss.norm.ppf(0.9)# 概率密度函数从-⽆穷到+⽆穷的积分是1,...

    Excel公式与函数大辞典.宋翔(带书签高清文字版).pdf

    1.2.1 A1引用样式与R1C1引用样式 6 1.2.2 输入公式 8 1.2.3 修改公式 9 1.2.4 移动和复制公式 9 1.2.5 删除公式 12 1.2.6 改变公式的计算方式 13 1.3 在公式中使用函数 14 1.3.1 为什么使用函数 14 1.3.2 ...

    核密度非参数估计的matlab代码-Non_Parametric_Kernel_Density_Estimation:我们建议使用基于核密度估

    核密度非参数估计的matlab代码Non_Parametric_Kernel_Density_Estimation 我们建议使用基于核密度估计 (KDE) 的方法进行分类。 这种非参数方法本质上以一种有原则的方式为每个类提供了成员资格的可能性。 该实现用于...

    R语言经典实例(中+英)

    R语言经典实例(中+英) 第1章 R入门和获得帮助 7  1.1 下载和安装R软件 8  1.2 开始运行R软件 10  1.3 输入R命令 13  1.4 退出R 15  1.5 中断R正在运行的程序 16  1.6 查看帮助文档 17  1.7 获取函数的帮助...

    现代通信原理作业答案_曹志刚版

    5、 窄带高斯噪声的一维随机包络服从_________ 分布,其概率密度函数p(x)= ;OOK信号加窄带高斯噪声的 一维随机包络服从 分布。 6、 在0—-T时刻内高度为A的矩形信号,与之相匹配的滤波器的单位冲激响应h(t)...

    DSP算法大全C语言版本-完整版

    均匀分布的概率密度函数为 ≤x≤b fCx) 其它 通常用U(b表示均匀分布的均值为“士方差为2 产生均匀分布随机数的方法如下 首先,由给定的初值x0,用混合同余法 T-1+c)(mod M) y; /M 产生(0,1)区间上的随机数y。...

Global site tag (gtag.js) - Google Analytics