主要内容

lillietest

Lilliefors测试

描述

例子

h= lillietest (x返回null假设的测试决定,即向量中的数据x来自于正态分布,而不是来自于这种分布,使用Lilliefors检验。结果h1如果检验在5%显著性水平上拒绝了原假设,并且0否则。

例子

h= lillietest (x名称,值返回带有由一个或多个名称-值对参数指定的附加选项的测试决定。例如,您可以针对不同的分布族测试数据、更改显著性级别或计算p-使用蒙特卡罗近似。

例子

hp) = lillietest (___还返回p价值p,使用前面语法中的任何输入参数。

例子

hpkstatcritval) = lillietest (___也返回测试统计信息kstat临界值critval的测试。

例子

全部折叠

加载示例数据。检验零假设:汽车里程,单位为英里每加仑(英里/加仑),在不同的车型上遵循正态分布。

负载carbig[h p、k、c] = lillietest (MPG)
警告:P小于表中最小值,返回0.001。
h = 1
p = 1.0000 e 03
k = 0.0789
c = 0.0451

测试数据k是否大于临界值c,所以lillietest返回h = 1表示在默认的5%显著性水平上拒绝原假设。警告指出返回的 p -value小于预计算值表中的最小值。找一个更准确的 p 值,使用MCTol进行蒙特卡罗近似。看到用蒙特卡罗近似法确定p值

加载示例数据。创建一个包含学生考试成绩数据第一列的向量。

负载examgradesx =成绩(:1);

检验样本数据来自1%显著性水平的正态分布的原假设。

(h p) = lillietest (x,“α”, 0.01)
h = 0
p = 0.0348

的返回值h = 0表明lillietest在1%显著性水平上不拒绝零假设。

加载示例数据。检验零假设:汽车里程,单位为英里每加仑(英里/加仑),在不同型号的汽车中呈指数分布。

负载carbigh = lillietest(英里/加仑,“分布”“指数”
h = 1

的返回值h = 1表明lillietest在默认的5%显著性水平上拒绝零假设。

生成两个样本数据集,一个来自威布尔分布,另一个来自对数正态分布。执行Lilliefors测试以评估每个数据集是否来自威布尔分布。通过使用威布尔概率图(wblplot).

从Weibull发行版生成示例。

rng (“默认”) data1 = wblrnd(0.5,2,[500,1]);

方法执行Lilliefors测试lillietest.若要测试数据的威布尔分布,请测试数据的对数是否具有极值分布。

h1 = lillietest(日志(data1),“分布”“极端值”
h1 = 0

的返回值h1 = 0表明lillietest未能在默认的5%显著性水平上拒绝原假设。使用威布尔概率图确定测试决策。

wblplot (data1)

图中包含一个坐标轴。标题为Weibull Probability Plot的坐标轴包含3个类型为line的对象。

图表明数据遵循威布尔分布。

从对数正态分布生成样本。

data2 = lognrnd(5 2[1] 500年);

执行Lilliefors测试。

h2 = lillietest(日志(data2),“分布”“极端值”
h2 = 1

的返回值h2 = 1表明lillietest在默认的5%显著性水平上拒绝零假设。使用威布尔概率图确定测试决策。

wblplot (data2)

图中包含一个坐标轴。标题为Weibull Probability Plot的坐标轴包含3个类型为line的对象。

图中显示数据不遵循威布尔分布。

加载示例数据。检验零假设:汽车里程,单位为英里每加仑(英里/加仑),在不同的车型上遵循正态分布。确定 p 用蒙特卡罗近似,最大蒙特卡罗标准误差为1的军医

负载carbig(h p) = lillietest(英里/加仑,“MCTol”1)的军医
h = 1
p = 8.3333 e-06

的返回值h = 1表明lillietest拒绝数据来自5%显著性水平的正态分布的零假设。

输入参数

全部折叠

样本数据,指定为向量。

数据类型:|

名称-值对的观点

指定可选的逗号分隔的对名称,值参数。的名字参数名和价值为对应值。的名字必须出现在引号内。可以以任意顺序指定多个名称和值对参数Name1, Value1,…,的家

例子:“分布”、“指数”,“阿尔法”,0.01在1%显著性水平下检验原假设,即总体分布属于指数分布族。

假设检验的显著性水平,指定为逗号分隔对组成“α”和范围(0,1)的标量值。

  • 如果MCTol没有使用,α必须在范围[0.001,0.50]内。

  • 如果MCTol使用,α必须在(0,1)范围内。

例子:“阿尔法”,0.01

数据类型:|

分布族为假设检验,指定为逗号分隔对组成“分配”下面是其中之一。

“正常” 正态分布
“指数” 指数分布
“极端值” 极端值分布

  • 测试x对于对数正态分布,测试是否日志(x)服从正态分布。

  • 测试x对于威布尔分布,测试是否日志(x)有一个极值分布。

例子:“分布”、“指数”

最大蒙特卡罗标准误差p,p-测试的值,指定为逗号分隔的对“MCTol”和范围(0,1)的标量值。

例子:“MCTol”,0.001

数据类型:|

输出参数

全部折叠

假设检验结果,返回为10

  • 如果h= 1,表示拒绝零假设α显著性水平。

  • 如果h= 0,这表示在α显著性水平。

p-测试的值,作为范围(0,1)内的标量值返回。p是观察到一个检验统计量与零假设下的观测值相同或更极端的概率。小的值p对无效假设的有效性产生了怀疑。

  • 如果MCTol没有使用,p对临界值表使用反插值计算,并返回为范围[0.001,0.50]的标量值。lillietest警告说,当p在列表范围内未找到,并返回最小或最大的列表值。

  • 如果MCTol使用,lillietest进行蒙特卡罗模拟,以计算更准确p值,p作为范围(0,1)中的标量值返回。

测试统计量,作为非负标量值返回。

假设检验的临界值,作为非负标量值返回。

更多关于

全部折叠

Lilliefors测试

Lilliefors检验是一种双边拟合优度检验,适用于零分布的参数是未知的,必须估计。这与单样本Kolmogorov-Smirnov检验相反,后者要求完全指定零分布。

Lilliefors检验统计量为:

D 马克斯 x | F x G x |

在哪里 F x 样本数据的实证CDF和 G x 为估计参数等于样本参数的假设分布的CDF。

lillietest是否可以用来测试数据向量x通过对数据向量应用转换并运行适当的Lilliefors测试,可以得到对数正态分布或Weibull分布:

  • 测试x对于对数正态分布,测试是否日志(x)服从正态分布。

  • 测试x对于威布尔分布,测试是否日志(x)有一个极值分布。

当原假设不是一个位置-规模的分布族时,Lilliefors检验不能使用。

蒙特卡罗标准误差

蒙特卡罗标准误差是由于模拟的误差p价值。

蒙特卡罗标准误差计算为:

年代 E p 1 p mcreps

在哪里 p 是估计的p-假设检验的值,并且mcreps为进行蒙特卡罗复制的次数。

蒙特卡罗复制的次数,mcreps的蒙特卡罗标准误差 p 小于指定的值MCTol

算法

为了计算假设检验的临界值,lillietest对小于1000的样本量和显著性水平在0.001和0.50之间的临界值,用蒙特卡罗模拟预先计算到一个表中。用的桌子lillietest比Lilliefors最初介绍的表更大,更准确。如果更准确的话p-值,或者如果所需显著性水平小于0.001或大于0.50,则MCTol输入参数可以用来运行蒙特卡罗模拟来计算p值更准确。

当检验统计量的计算值大于临界值时,lillietest在显著性水平上拒绝零假设α

lillietest对待x忽略丢失的值。

参考文献

w·J·科诺菲尔实际的非参数统计.霍博肯:约翰·威利父子公司,1980。

[2] Lilliefors, H. W.“关于均值未知指数分布的Kolmogorov-Smirnov检验”。美国统计协会杂志.1969年第64卷,387-389页。

[3] Lilliefors, H. W.“关于均值和方差未知的正态性的Kolmogorov-Smirnov检验”。美国统计协会杂志.1967年第62卷,399-402页。

之前介绍过的R2006a