跳到主要内容
访问密钥 NCBI主页 MyNCBI主页 主要内容 主导航
J应用统计。2022; 49(4): 782–802.
2020年10月31日在线发布。 数字对象标识:10.1080/02664763.2020.1840535
预防性维修识别码:PMC9041659
PMID:35707810

基于记录值的极值检验指数性

熊培汉, 威威庄,邱国欣a、,b条

关联数据

补充资料

摘要

在本文中,我们首先基于记录值的极值给出了指数分布的特征,然后介绍了指数性的良好性检验。蒙特卡罗模拟用于计算针对不同样本大小和显著性水平的拟议测试的临界值。为了证明所提出的测试的优势,我们采用了58个竞争测试,并计算了针对不同类型危险函数的不同备选方案的调整功率。功率结果表明,如果备选方案的故障率增加或浴缸的故障率减少-增加,特别是当样本量较小时,我们提出的测试具有更好的调整功率。最后,通过监控第页-测试值。

关键词:调整功率、极值、指数检验、蒙特卡罗模拟、记录值
数学学科分类(2010)初级62G30、次级63E10、62B10

1.简介

指数分布是统计工作中最常用的连续分布之一。指数分布由于其无记忆性,在可靠性理论、寿命试验和随机过程中得到了广泛的应用。在这种情况下,区分随机样本是否来自指数分布非常重要。那是因为n个独立同分布(iid)无序观测X(X)1,X(X)2,,X(X)n个从绝对连续的累积分布函数(cdf)F类,我们想测试

H(H)0:F类费用(λ)v.s.版本。H(H)1:F类费用(λ),

哪里费用(λ)表示具有cdf的指数分布类

F类(x个|λ)=1负极e(电子)负极λx个,x个0,0,x个<0,
(1)

λ>0是一个未知参数。为此,文献中有许多指数性的良好性测试,例如Zamanzade和Arghami[57]努格哈比和阿尔哈米[41]基于顺序统计矩,Noughabi和Balakrishnan[42]基于Phi-dixture、Grané和Fortiana[20]基于相关性和协方差,托拉比等。[52]基于经验分布函数等。在开发这些指数性良好性测试的技术中,指数分布的表征起着重要作用,可以作为一个基本成分。在这些测试中,我们有一些基于归一化间距的测试(例如爱泼斯坦测试[15],科查尔[27]),阶矩统计(例如Henze[22]),均匀分布的随机变量(例如Marange和Qin[34]),和相同分布的随机变量(例如努格哈比和阿尔哈米[39])以及其他许多人。但据我们所知,还没有人使用基于记录值的极值的特征来构建指数性测试。我们将在本文中填补这一空白。

对于随机变量X(X)具有概率密度函数(pdf)(f),其极值由Lad定义等。[30]作为

J型(X(X))=负极12负极(f)2(x个)d日x个.
(2)

学者们发现了许多有趣的性质,例如邱[45]、邱和贾[46,47]还有Jose和Sathar[25]. 此外,邱和贾[47]还有努格哈比和贾拉希费雷斯[43]提出了一些极值估计量,并将其应用于均匀性检验。我们的动机来自于我们之前的研究,其中给出了一个证明,表明指数分布可以由上记录值的极值唯一地表征(未发表的观察结果)。然后利用邱和贾给出的极值估计[46],我们提出了一种新的指数性检验方法。在这个数据丰富的时代,用最少的数据确定分布仍然很重要。例如,在许多寿命测试中,数据采集可能非常耗时且成本高昂。我们将进一步证明,当样本量较小时,我们的测试尤其有效。

本文的其余部分组织如下。在节中2,提出了一种新的检验统计量来检验指数性。此外,还显示了我们针对不同样本大小提出的测试的临界值。章节首先列出了现有的竞争对手测试和广泛使用的指数测试替代方法。然后,我们解释了为什么我们使用调整后的权力而不是惩罚权力或内在权力来比较我们提出的测试与竞争对手测试的表现。我们通过模拟研究完成了这一部分,这表明我们的测试对于增加或浴缸故障率的替代品具有很好的调整能力,并且在样本量较小时显著优于其他替代品。即使对于降低或单峰故障率的备选方案,我们建议的测试也不总是最差的测试。章节4给出了三个实例来说明所提出的测试和部分的优越性和鲁棒性5本文得出结论。读者可以在附录中看到Python测试统计的实现。完整的工作代码可以在补充文件中找到。

2.测试统计

X(X)1,X(X)2,,X(X)n个,是绝对连续cdf中的iid随机变量序列F类和pdf(f).表示j个第阶统计依据X(X)j个:n个,1j个n个,这是j个序列中第个最小值X(X)1,X(X)2,,X(X)n个.定义上限记录时间的序列U型(k个)作为U型(1)=1;U型(k个+1)=最小值{j个:j个>U型(k个);X(X)j个>U型k个X(X)},k个=1,2,k个第个上限记录值为U型k个X(X)=X(X)U型(k个),k个=1,2,。的pdfk个第个上记录值U型k个X(X)由提供

(f)U型k个X(X)(x个)=负极日志F类¯(x个)k个负极1Γ(k个)(f)(x个),负极<x个<,
(3)

哪里F类¯(x个)=1负极F类(x个)、和Γ(k个)=(k个负极1)!是完整的伽玛函数。读者可以参考阿诺德等。[2]用于定义和应用上限记录值。

根据(2)和(),的极值U型k个X(X)可以写成

J型(U型k个X(X))=负极12负极(f)U型k个X(X)2(x个)d日x个=负极12Γ2(k个)负极负极日志F类¯(x个)2k个负极2(f)2(x个)d日x个.
(4)

如果X(X)具有指数分布,cdf由(1),那么我们还有

J型(U型k个X(X))=负极12Γ2(k个)0(λx个)2k个负极2λ2e(电子)负极2λx个d日x个=负极λΓ(2k个负极1)22k个Γ2(k个).

使用(2)我们再次看到J型(X(X))=负极λ/4因此,如果X(X)具有指数分布,cdf由(1)那么,对所有人来说k个=1,2,,它认为

J型(U型k个X(X))=Γ(2k个负极1)22k个负极2Γ2(k个)J型(X(X)).

相反,如果上述等式适用于所有情况k个=1,2,,我们可以证明X(X)必须具有带速率参数的指数分布负极4J型(X(X))这在以下定理2.1中进行了说明。充分性的证明类似于邱中定理3.7的证明[45],此处省略。

定理2.1

非负随机变量X(X)具有速率参数的指数分布λ>0当且仅当

J型(U型k个X(X))=Γ(2k个负极1)22k个负极2Γ2(k个)J型(X(X)),k个=1,2,,

接下来,我们利用定理2.1介绍我们的指数性检验统计量。为此,我们让

E类k个=J型(U型k个X(X))负极Γ(2k个负极1)22k个负极2Γ2(k个)J型(X(X)).

从定理2.1可以看出E类k个=0为所有人k个=1,2,,当且仅当X(X)是指数型的。因此,E类k个可以用作指数的基本概念,可以用作指数性的测试统计量。假设E类ˆk个是的估计量E类k个通过随机抽样X(X)1,,X(X)n个基于X(X),则为的小值或大值E类ˆk个可以被视为非赞助人的症状,因此我们拒绝零假设H(H)0。为了简单起见,我们采用k个 = 2并考虑E类ˆ2作为我们在即将到来的讨论中对指数性的测试统计。

推导E类ˆ2,我们首先回忆起邱和贾[46]引入了样本估计J型(X(X))通过J型2n个,即。

J型2n个=负极12n个=1n个c(c)/n个X(X)(+):n个负极X(X)(负极):n个,

哪里X(X)1:n个X(X)2:n个X(X)n个:n个订单统计是基于X(X)1,,X(X)n个.窗口大小是小于的正整数n个/2,X(X)(负极):n个=X(X)1:n个如果负极1,X(X)(+):n个=X(X)n个:n个如果+n个、和

c(c)=1+负极1,1,2,+1n个负极,1+n个负极,n个负极+1n个.
(5)

另一方面k个 = 每侧2个(4),我们有

J型(U型2X(X))=负极12负极日志2F类¯(x个)(f)2(x个)d日x个=负极1201日志2(1负极v(v))(f)(F类负极1(v(v)))d日v(v)=负极1201日志2(1负极v(v))d日F类负极1(v(v))d日v(v)负极1d日v(v).
(6)

因此,如果我们遵循瓦西塞克的想法[53],即替换分配函数F类通过经验分布函数F类n个,并使用差分算子代替微分算子,然后求F类负极1(v(v))在方程式中(6)可以通过订单统计的函数进行估计。形式上,对于随机样本X(X)1,,X(X)n个X(X),J型(U型2X(X))可以通过以下方式进行估算

J型(U型2X(X))ˆ=负极12n个=1n个日志21负极n个+12/n个X(X)(+):n个负极X(X)(负极):n个.

因此,对于E类2可以写成

E类ˆ2=J型(U型2X(X))图6负极J型2n个=负极12n个=1n个日志21负极n个+12/n个X(X)(+):n个负极X(X)(负极):n个+12n个=1n个c(c)/n个X(X)(+):n个负极X(X)(负极):n个=负极12n个=1n个2日志21负极n个+1负极c(c)/n个X(X)(+):n个负极X(X)(负极):n个,

哪里c(c)定义如下(5).是小于的正整数n个/2和X(X)(负极):n个=X(X)1:n个如果负极1,X(X)(+):n个=X(X)n个:n个如果+n个.

由于E类ˆ2取决于样本和窗口大小其渐近分布过于复杂,无法进行理论分析。通过Monte-Carlo方法,我们生成了10000个大小的样本n个 = 从零分布和(1负极α)将第个分位数作为显著性水平的临界值α.表格12给出不同样本大小在显著性水平上的临界值α=0.01和0.05。

表1。

的临界值|E类图62|在显著性水平α=0.01.
n个51020304050100
116.365410.31145.24112.94772.12462.04201.0722
22.18291.29120.68040.56530.46380.41170.3128
 0.77060.45640.38850.33700.31710.2642
4 0.63780.36770.33470.28770.28300.2420
5 0.57690.32280.28720.27590.25950.2284
6  0.28960.25750.26060.25270.2216
7  0.29620.24360.23590.23360.2150
8  0.30870.21940.23300.22720.2116
9  0.31600.21000.20940.21140.2068
10  0.34430.20680.19490.20480.2019
15   0.27060.17100.15540.1835
20    0.22640.14950.1649
30      0.1321
40      0.1011
50      0.1462

表2。

的临界值|E类ˆ2|在显著性水平α=0.05.
n个51020304050100
13.49622.03781.18650.88430.71610.65290.5083
20.83030.55880.38640.34690.31440.29580.2531
 0.38470.29060.26930.25220.24660.2207
4 0.33650.24680.23650.22440.22020.2044
5 0.33910.20480.21320.21260.20840.1953
6  0.18840.19420.19920.19850.1881
7  0.17610.17730.18650.18830.1835
8  0.19350.16130.17270.17890.1813
9  0.21210.14700.16130.17100.1772
10  0.23350.14120.14880.16300.1736
15   0.19010.11650.11850.1593
20    0.16860.10490.1445
30      0.1099
40      0.0759
50      0.1119

3.功率比较

许多学者开发了一系列现有的指数性测试。这里有几个,崔等。[9],易卜拉希米等。[13]努格哈比和阿尔哈米[38,39]Baratpour和Habibi Rad[4]提出了测试统计电视E,总费用,T型V(V)1,T型A类1C类K(K)L(左)n个通过使用不同的熵估计器来估计非对称Kullback-Leibler距离。扎曼扎德和阿尔哈米[57]努格哈比和阿尔哈米[41]推导了其估计量的非参数分布函数,并进一步提出了检验统计量T型V(V)2,T型E类2,T型A类2TZ公司Baringhaus和Henze研究了其他指数性测试[5]努格哈比和阿尔哈米[38,39]Ascher进行了功率比较[]、亨泽和美坦尼斯[23]和托拉比等。[52]. 我们总共采用了58项测试作为我们建议测试的合适竞争对手。据托拉比介绍等。[52]根据指数分布的特征,现有测试可分为以下12类。我们在表中总结了这58项竞争对手测试基于马兰奇和秦的类似工作[34]. 这些测试的详细信息可以在Ascher中找到[]、亨泽和美坦尼斯[23]和托拉比等。[52].

表3。

现有竞争对手测试。
测试统计的参考
熵估计量  
Choi公司等。电视EChoi公司等。[9]
Choi公司等。总费用Choi公司等。[9]
易卜拉欣等。 T型V(V)1易卜拉欣等。[13]
努格哈比和阿尔哈米 T型A类1努格哈比和阿尔哈米[38]
Baratpour和Habibi Rad C类K(K)L(左)n个Baratpour和Habibi Rad[4]
熵估计量的修正矩  
公园和公园 T型V(V)2公园和公园[44]
Park和Parki T型E类2公园和公园[44]
努格哈比和阿尔哈米 T型A类2Noughabi和Arghami[41]
扎曼扎德和阿尔哈米TZ公司扎曼扎德和阿尔哈米[57]
Phi发散  
努格哈比和巴拉克里什南TKL公司努格哈比和巴拉克里什南[42]
努格哈比和巴拉克里什南真实航向努格哈比和巴拉克里什南[42]
努加比和巴拉克里希南TJ公司努格哈比和巴拉克里什南[42]
努格哈比和巴拉克里什南TT公司努格哈比和巴拉克里什南[42]
努格哈比和巴拉克里什南 T型χ努格哈比和巴拉克里什南[42]
exp分布的特征  
Noughabi和Arghami(统一(0,1)) T型1努格哈比和阿尔哈米[39]
Noughabi和Arghami(费舍尔F类(2,2)) T型2努格哈比和阿尔哈米[39]
努格哈比和阿尔哈米(统一(-1,1)) T型努格哈比和阿尔哈米[39]
Noughabi和Arghami(籽粒密度) T型4努格哈比和阿尔哈米[40]
Volkova(罗斯伯格表征) 1沃尔科娃[54]
Volkova和Nikitin(Ahsanullah角色化) 2沃尔科娃和尼基丁[55]
基于间距  
科查尔 K(K)n个科查尔[27]
爱泼斯坦 E类P(P)S公司n个爱泼斯坦[15]
格内登科F检验 n个(第页)格涅坚科等。[19]、阿斯切尔[]
改进的格内登科F检验 n个(第页)哈里斯[21]
基于洛伦兹曲线和基尼指数  
洛伦兹 L(左)n个(第页)Gail和Gastwirth[17]
吉尼 G公司n个盖尔和加斯特沃思[18]
经验分布函数  
科尔莫戈罗夫·斯米尔诺夫 K(K)S公司n个利利福斯[32]
柯伊伯V(V)达戈斯蒂诺和斯蒂芬斯[11]
芬克尔斯坦和谢弗斯 S公司达戈斯蒂诺和斯蒂芬斯[11]
克拉梅尔·冯·米塞斯 W公司n个2达戈斯蒂诺和斯蒂芬斯[11]
安德森·达林 A类2达戈斯蒂诺和斯蒂芬斯[11]
托拉比等。 H(H)n个(1)托拉比等。[52]
托拉比等。 H(H)n个(2)托拉比等。[52]
基于残差平均函数  
克拉梅·冯·米塞斯 C类M(M)¯n个巴林豪斯和亨泽[6]
科尔莫戈罗夫·斯米尔诺夫 K(K)S公司¯n个巴林豪斯和亨泽[6]
Aboukhamseen和Aly T型1,n个Aboukhamseen和Aly[1]
Aboukhamseen和Aly T型2,n个Aboukhamseen和Aly[1]
Aboukhamseen和Aly T型1,n个(γ)Aboukhamseen和Aly[1]
Aboukhamseen和Aly T型2,n个(γ)Aboukhamseen和Aly[1]
经验拉普拉斯变换  
巴林豪斯和亨泽 B类H(H)n个巴林豪斯和亨泽[5]
亨泽 H(H)E类n个亨泽[22]
经验特征函数  
Epps和皮带轮 E类P(P)n个Epps和皮带轮[14]
Henze和Meintanis T型(n个,)(1)亨泽和美坦尼斯[23]
亨泽和美坦尼斯 T型(n个,)(2)亨泽和美坦尼斯[23]
基于相关性和协方差  
夏皮罗一威尔克W公司夏皮罗和威尔克[49]
斯蒂芬斯对SW的修改 W公司斯蒂芬斯[50]
Fortiana和Grané n个福蒂亚纳和格拉内[16]
蒙塔泽里和托拉比COV公司蒙塔泽里和托拉比[36]
其他测试  
考克斯和奥克斯 C类O(运行)n个考克斯和橡树[10]
莫兰 M(M)n个莫兰[37],奇里纳[51]
阿特金森 M(M)n个(第页)Mimoto和Zitikis[35]
Hollander和ProschanJ型Hollander和Proschan[24]
德什潘德 J型b条德什潘德[12]
金伯 D类第页金伯[26]
彼得拉 P(P)n个盖尔和加斯特沃思[18]
马朗和秦 E类E类L(左)R(右)1马朗和秦[34]
马朗和秦 E类E类L(左)R(右)2马朗和秦[34]

  1. 基于熵估计的测试。
  2. 基于熵估值器修正矩的测试。
  3. 基于Phi散度的测试。
  4. 基于指数分布特征的测试。
  5. 基于间距的测试。
  6. 基于洛伦兹曲线和基尼指数的测试。
  7. 基于经验分布函数的测试。
  8. 基于残差均值函数的测试。
  9. 基于经验拉普拉斯变换的测试。
  10. 基于经验特征函数的测试。
  11. 基于相关性和协方差的测试。
  12. 其他测试。

对于替代分布,我们使用Henze和Meintanis提出的相同分布[23]以及他们的参数选择。这些分布的定义如下。感兴趣的读者可以参考Torabi中的图2等。[52]对于几个参数值的故障率函数的形状θ正如Henze和Meintanis所指出的[23]这些概率分布族通常用作指数模型的替代品,包括故障率增加(IFR)、故障率减少(DFR)、单峰故障率增加-减少(UFR)以及浴缸故障率减少-增加(BFR)的密度。

  • 威布尔分布,表示为W公司(θ),pdf由提供
    (f)(x个)=θx个θ负极1经验(负极x个θ),θ>0, x个>0
  • 伽马分布,表示为Γ(θ),pdf格式由
    (f)(x个)=1Γ(θ)x个θ负极1经验(负极x个),θ>0, x个>0
  • 对数正态分布,表示为L(左)N个(θ),pdf由提供
    (f)(x个)=1θx个2π经验负极日志2x个2θ2,θ>0, x个>0
  • 半正态分布,表示为海南,pdf由提供
    (f)(x个)=2π经验负极x个22,x个>0
  • 均匀分布,表示为U型,pdf由提供(f)(x个)=1,0x个1
  • 修正极值分布,表示为E类V(V)(θ),cdf由给出
    F类(x个)=1负极经验{θ负极1(1负极e(电子)x个)},θ>0, x个>0
  • 线性递增失效率分布,表示为L(左)F类(θ),pdf由提供
    (f)(x个)=(1+θx个)经验负极x个负极θ2x个2,θ>0, x个>0
  • Dhillon定律,表示为D类L(左)(θ),cdf由给出
    F类(x个)=1负极经验负极日志(x个+1)θ+1,θ>0, x个>0
  • Chen分布,表示为C类H(H)(θ),cdf由给出
    F类(x个)=1负极经验{2(1负极e(电子)x个θ)},θ>0, x个>0

对于测试的功率比较,Cavus等。[8]注意到,为了比较测试的性能,必须同时考虑功率和I类错误概率。因此,他们引入了以下惩罚权力,其定义为:

γ=1负极β1+1负极α1α0,
(7)

哪里α1β分别为I类和II类错误概率,以及α0是标称水平。Cavus公司等。[8]一个例子表明,惩罚权力与调整权力相比具有优势(张和布斯[58])和内在动力(劳埃德[33]). 然而,在指数性测试的情况下,大多数测试统计数据,包括电视E,总费用,T型V(V)1,T型A类1,E类E类L(左)R(右)1,E类E类L(左)R(右)2,H(H)n个(1),H(H)n个(2)而我们提出的测试,其理论分布过于复杂,无法通过分析获得。由于无法计算这些测试的I型错误概率,因此惩罚功率和劳埃德固有功率在我们的研究中均无效。因此,我们转而使用张和布斯的调整后权力来比较我们的提案和竞争对手的测试。对于测试统计T型具有标称水平α0,其调整功率定义为:

第页ˆ=1k个=1k个T型1c(c)α,

哪里c(c)α是测试统计的临界值T型来自(1负极α0)真零假设下测试统计值的分位数,T型11,T型12,,T型1k个是测试统计值T型对于k个替代假设下的蒙特卡罗模拟,以及(A类)是集合的指示函数A类显然,在零假设下调整后的权力只是α因为这种方式C类α被选中。然后,可以公平地比较替代分布下不同测试的调整功率。根据这个定义,我们估计了E类ˆ2重复10000次大小样本n个 = 10%、20%、50%,标称水平为5%。如上所述E类ˆ2取决于窗口大小m。因此必须在计算其值之前指定,以使其获得良好(而非最佳)的调整功率。我们对不同样本大小进行了模拟,并建议=[n个/2]IFR和UFR分布和=[n个]用于DFR和BFR分布。结果如表所示49用粗体表示每列的最佳三次调整幂。竞争对手测试的一些调整后的力量来自Torabi等。[52]马兰奇和秦[34].

表5。

不同指数测试的调整功率,α0=0.05,n个 = 10
  UFR(超滤)DFR公司BFR公司
  LN(0.8)LN(1.5)DL(1)DL(1.5)W(0.8)G(0.4)瑞士(0.5)
1电视E241917334117
2总费用254204621
T型V(V)1266204525
4 T型A类12612150100
5 C类K(K)L(左)n个1715143475
6 T型V(V)2231417361711
7 T型E类22320173542114
8 T型A类22324153352518
9TZ公司2329163172921
10TKL公司2411946100
11真实航向14710236655
12TJ公司14710236655
13TT公司14810236654
14 T型χ4720215947
15 T型1285204775136
16 T型230225043725
17 T型285204775036
18 T型42912148100
19 12161331104635
20 21991333115038
21 K(K)n个2011641100
22 E类P(P)S公司n个72059105136
23 n个(第页)631614124537
24 n个(第页)102061071412
25 L(左)n个(第页)14321231145543
26 G公司n个14391333144937
27 K(K)S公司n个16321330104131
28V(V)1725122993525
29 S公司19351436124634
30 W公司n个218351435114534
31 A类213371027176754
32 H(H)n个(1)2412046100
33 H(H)n个(2)24712216554
34 K(K)S公司¯n个1826153562820
35 C类M(M)¯n个16381332124534
36 T型1,n个1238924144635
37 T型2,n个14401029155040
38 T型1,n个(γ)1424112793525
39 T型2,n个(γ)1310112781816
40 B类H(H)n个17401334144737
41 H(H)E类n个17401335144737
42 E类P(P)n个15401435144737
43 T型n个,2.5(1)1951741220
44 T型n个,2.5(2)1919164052114
45 C类O(运行)n个16341537176853
46 M(M)n个1031722217462
47 M(M)n个(0)26191946116246
48 M(M)n个(0.99)16401232155038
49J型2311844100
50 J型0.52311740100
51 D类第页26271739105136
52 W公司n个836710112620
53 W公司15381130133527
54 n个12361025123427
55COV公司3122249110
56 P(P)n个14391128144536
57 E类E类L(左)R(右)12924215052115
58 E类E类L(左)R(右)21411331100
59 E类图6252526054

表6。

不同指数测试的调整功率,α0=0.05,n个 = 20
  国际财务报告
  G(1)W(1.4)G(2)海南U(0,1)信道(1)CH(1.5)左前(2)左前(4)企业价值(0.5)企业价值(1.5)
1电视E52742145210621827926
2总费用537502488188330421848
T型V(V)1537512387188230421848
4 T型A类1547633080208736522052
5 C类K(K)L(左)n个533382892218634482159
6 T型V(V)2531451983157625351540
7 T型E类2529421779137322311335
8 T型A类25243513699641725928
9TZ公司517277546501152145
10TKL公司549633184239141562259
11真实航向511215150214
12TJ公司51121515214
13TT公司50016111111
14 T型χ511122191214
15 T型1535551852137223341334
16 T型2536551855147425361435
17 T型535551751137224351334
18 T型4546622869198335501949
19 15274612289601726922
20 2529471439106720341128
21 K(K)n个544533388259142572564
22 E类P(P)S公司n个5810731621711712
23 n个(第页)522251862126522321339
24 n个(第页)569586116666
25 L(左)n个(第页)534462157147827401440
26 G公司n个535462170158429421546
27 K(K)S公司n个528401852136724351335
28V(V)526371666126822321235
29 S公司535492270168229431546
30 W公司n个2534472166147928421443
31 A类2530451762127624361237
32 H(H)n个(1)549603178249140552358
33 H(H)n个(2)66102182294728
34 K(K)S公司¯n个535462472187932441848
35 C类M(M)¯n个535472270168330431647
36 T型1,n个527361863137024361338
37 T型2,n个532442171148028421445
38 T型1,n个(γ)528362174167026381542
39 T型2,n个(γ)529382179157827411546
40 B类H(H)n个536482165148428421444
41 H(H)E类n个537482162148428421543
42 E类P(P)n个536482166158428421545
43 T型n个,2.5(1)545553386259142562563
44 T型n个,2.5(2)542503185239139542361
45 C类O(运行)n个537541950138125371337
46 M(M)n个5294613379701929926
47 M(M)n个(0)544642250168129411639
48 M(M)n个(0.99)534471963138327401343
49J型547622978228939532255
50 J型0.5543602665198435481949
51 D类第页533511751126623331331
52 W公司n个521201972136325341343
53 W公司532402275148428421449
54 n个532382386178530421754
55COV公司547662556188333441844
56 P(P)n个535462262148030431545
57 E类E类L(左)R(右)1547632767198236511951
58 E类E类L(左)R(右)2541453492269342562666
  E类ˆ2529237100839973898374

表7。

不同指数测试的调整功率,α0=0.05,n个 = 20
  UFR(超滤)DFR公司BFR公司
  LN(0.8)LN(1.5)DL(1)DL(1.5)W(0.8)G(0.4)瑞士(0.5)
1电视E5245306163118
2总费用4019266522414
T型V(V)1432627673219
4 T型A类14903581100
5 C类K(K)L(左)n个1843164882516
6 T型V(V)24342256265133
7 T型E类24550256095639
8 T型A类245592455136346
9TZ公司42642245176650
10TKL公司4103379100
11真实航向17011348979
12TJ公司17011338979
13TT公司17110348980
14 T型χ4702287764
15 T型151103174168673
16 T型25073175138267
17 T型5093174168570
18 T型45303680100
19 148112967188068
20 245162868208270
21 K(K)n个2502366100
22 E类P(P)S公司n个1030714116748
23 n个(第页)959933206856
24 n个(第页)1730101671311
25 L(左)n个(第页)27552264238270
26 G公司n个24671962247663
27 K(K)S公司n个30582056177156
28V(V)32441954136247
29 S公司33612066207561
30 W公司n个233622365207661
31 A类234622163268978
32 H(H)n个(1)3303074500
33 H(H)n个(2)871420348979
34 K(K)S公司¯n个28552262146247
35 C类M(M)¯n个27662163227561
36 T型1,n个23632052217359
37 T型2,n个25652160247764
38 T型1,n个(γ)22472350156147
39 T型2,n个(γ)19242252144236
40 B类H(H)n个26672164247763
41 H(H)E类n个29662064247963
42 E类P(P)n个25672064247663
43 T型n个,2.5(1)2718256743323
44 T型n个,2.5(2)22382161105037
45 C类O(运行)n个33602572289180
46 M(M)n个31522167319485
47 M(M)n个(0)47423481209077
48 M(M)n个(0.99)26672063257764
49J型4603579100
50 J型0.55003577100
51 D类第页56463471147962
52 W公司n个12601023154534
53 W公司20611653185342
54 n个18611550195643
55COV公司614384000
56 P(P)n个24642061227361
57 E类E类L(左)R(右)140563178155962
58 E类E类L(左)R(右)21918542781
  E类ˆ2402619595

表8。

不同指数测试的调整功率,α0=0.05,n个 = 50
  国际财务报告
  G(1)W(1.4)G(2)海南U(0,1)信道(1)CH(1.5)左前(2)左前(4)企业价值(0.5)企业价值(1.5)
1电视E555822391139734541453
2总费用56481421003310054733385
T型V(V)156584431003210055743386
4 T型A类15809548993210064833383
5 C类K(K)L(左)n个56266541004410068854496
6 T型V(V)256485381002810051712982
7 T型E类25588231100219943642275
8 T型A类25326111986931834642
9TZ公司5517162046145
10TKL公司58797611004510076914693
11真实航向521458225917
12TJ公司521458225917
13TT公司510452125212
14 T型χ5000530161205
15 T型1567922985189842621959
16 T型2569923188209845642162
17 T型567922985189842621959
18 T型45819552953510064843584
19 1565902562159639601551
20 2567912872179843751757
21 K(K)n个58590711005710083955798
22 E类P(P)S公司n个513199678451216821
23 n个(第页)556654396309856733081
24 n个(第页)51124596256967
25 L(左)n个(第页)5769146933010062813081
26 G公司n个5799054993810069873890
27 K(K)S公司n个564833993269854732675
28V(V)559793599239950712476
29 S公司5769049993410066843488
30 W公司n个25759048983210064833286
31 A类25749245992910061813085
32 H(H)n个(1)5889465995010079924994
33 H(H)n个(2)5376213787942444646
34 K(K)S公司¯n个5718650993610065823688
35 C类M(M)¯n个5779053993710069873790
36 T型1,n个565814399309958783083
37 T型2,n个5758850993610068863690
38 T型1,n个(γ)5607443100319956753183
39 T型2,n个(γ)56880481003510063823588
40 B类H(H)n个5799153983710069873789
41 H(H)E类n个5799153983710069873789
42 E类P(P)n个5809154983810069873890
43 T型n个,2.5(1)58190641004810077914896
44 T型n个,2.5(2)58085641004910078924995
45 C类O(运行)n个5829645913010060803078
46 M(M)n个575953277219948672163
47 M(M)n个(0)5839741822910056742871
48 M(M)n个(099)5799052983610068873689
49J型5869657994110074894191
50 J型0.55849653973810070863887
51 D类第页5699135100239948692473
52 W公司n个56157551004110069864293
53 W公司57480571004210072894394
54 n个57379591004710074894796
55COV公司582964387309960782974
56 P(P)n个6758749963310064833385
57 E类E类L(左)R(右)15859457984110073894191
58 E类E类L(左)R(右)257783621004810076914796
  E类ˆ2556731009810095999899

表4。

不同指数测试的调整功率,α0=0.05,n个 = 10
  国际财务报告
  G(1)W(1.4)G(2)海南U(0,1)信道(1)瑞士(1.5)左前(2)左前(4)企业价值(0.5)企业价值(1.5)
1电视E5172310288361318818
2总费用527341853145922291431
T型V(V)1526331852145621281430
4 T型A类1529362051146124311433
5 C类K(K)L(左)n个523271860145622291433
6 T型V(V)2521271444114617221123
7 T型E类2519251341114316201021
8 T型A类25172311379391418919
9TZ公司5162210338361316817
10TKL公司527341850146023311433
11真实航向51015121112
12TJ公司51015121112
13TT公司50011100011
14 T型χ52221428224
15 T型1524331532124919251124
16 T型2627361635135220261326
17 T型524331532124919251124
18 T型4527351840135422291329
19 1513208156281014612
20 25162210247361217716
21 K(K)n个526321956156125321534
22 E类P(P)S公司n个5676175127768
23 n个(第页)510119318291013716
24 n个(第页)55659585656
25 L(左)n个(第页)5172211318431420921
26 G公司n个5162111379471420923
27 K(K)S公司n个5162111288361419819
28V(V)5152011358381419819
29 S公司519261340104817231024
30 W公司n个25182412369441622922
31 A类2513188287361217716
32 H(H)n个(1)527361847166023311432
33 H(H)n个(2)5111211111
34 K(K)S公司¯n个521251542114619251227
35 C类M(M)¯n个5172211368441722822
36 T型1,n个5131710307351217717
37 T型2,n个5162111358421420821
38 T型1,n个(γ)516191240103915201021
39 T型2,n个(γ)5162012439431521922
40 B类H(H)n个5172411338451421821
41 H(H)E类n个5172411338461421821
42 E类P(P)n个5162211359471319922
43 T型n个,2.5(1)526311852156023311535
44 T型n个,2.5(2)525301751145923301434
45 C类O(运行)n个5162510309451419819
46 M(M)n个5914617529811512
47 M(M)n个(0)524341432115118251125
48 M(M)n个(0.99)5162310318441319820
49J型525331747135722291331
50 J型0.5521281434104718241123
51 D类第页5192712249391520918
52 W公司n个510109327241114717
53 W公司5162210378461420822
54 n个5151911428441419822
55COV公司526361635125321281226
56 P(P)n个5162010318411319819
57 E类E类L(左)R(右)1527371737135521291329
58 E类E类L(左)R(右)2522241760135320271332
  E类ˆ2515141893568243655833

表9。

不同指数测试的调整功率,α0=0.05,n个 = 50
  UFR(超滤)DFR公司BFR公司
  LN(0.8)LN(1.5)DL(1)DL(1.5)W(0.8)G(0.4)瑞士(0.5)
1电视E93856596117855
2总费用7369449398566
T型V(V)179744895118972
4 T型A类18407099000
5 C类K(K)L(左)n个24852176196849
6 T型V(V)291895796209584
7 T型E类292925895279687
8 T型A类285944386349788
9TZ公司55932142339583
10TKL公司76065100000
11真实航向194015610098
12TJ公司19311559998
13TT公司1940057100100
14 T型χ89251399484
15 T型1942768993410098
16 T型2942268993110097
17 T型942768993410098
18 T型4925731008667
19 193256898379997
20 292296899389997
21 K(K)n个3803995000
22 E类P(P)S公司n个19541129139069
23 n个(第页)19912179389587
24 n个(第页)4457244371312
25 L(左)n个(第页)60904998489997
26 G公司n个47953997489894
27 K(K)S公司n个71914695389892
28V(V)75844694289687
29 S公司73935298439995
30 W公司n个276945298449995
31 A类2869358995210099
32 H(H)n个(1)5105098000
33 H(H)n个(2)479523875910098
34 K(K)S公司¯n个62924396359790
35 C类M(M)¯n个60954497469994
36 T型1,n个55933894409892
37 T型2,n个56954297469995
38 T型1,n个(γ)51853588299383
39 T型2,n个(γ)39583291328677
40 B类H(H)n个47954097489995
41 H(H)E类n个46954097489994
42 E类P(P)n个45953997489994
43 T型n个,2.5(1)49744396259283
44 T型n个,2.5(2)32833293319384
45 C类O(运行)n个669255995610099
46 M(M)n个778660995610099
47 M(M)n个(0)8482701004810099
48 M(M)n个(0.99)44953996499995
49J型8606999000
50 J型0.59107299000
51 D类第页97837199319995
52 W公司n个16931761318470
53 W公司29932888348875
54 n个26932586358977
55COV公司96580100000
56 P(P)n个43943896469893
57 E类E类L(左)R(右)163935299399798
58 E类E类L(左)R(右)239303592145599
  E类ˆ201120100100

对于IFR替代方案,我们建议的测试是最强大的测试之一。除了W公司(1.4),G公司(2)G公司(1)(即指数分布),我们的建议是所有样本大小的所有IFR分布的最佳工作测试之一。特别是考虑小样本时(n个 = 10) ,我们提案的执行情况U型(0,1),C类H(H)(1),C类H(H)(1.5),L(左)F类(2),L(左)F类(4)E类V(V)(0.5)比竞争对手高出20%-30%。虽然随着样本量的增加,大多数测试都有更好的性能,但我们提出的测试相对于其他测试仍有相当大的优势C类H(H)(1),L(左)F类(2),E类V(V)(0.5)对于中等样本量(n个 = 20, 50). 虽然建议的测试在W公司(1.4)G公司(2),这不是最糟糕的。IFR替代方案的其他强大测试包括T型A类1,TKL公司,K(K)n个,H(H)n个(1)T型n个,2.5(1)。功能最差的测试是真实航向,TJ公司,TT公司,T型χH(H)n个(2).

对于UFR和DFR分布,建议的测试是功能最差的测试之一,因为它只适用于G公司(0.4)什么时候n个 = 20, 50. 我们观察到,对于所有四个UFR分布,没有一个单独的测试能够很好地稳定地工作,我们只能说,对于UFR来说,相对最强大的测试是T型4COV公司所有样本大小。对于DFR发行版,最强大的测试是真实航向,TJ公司,TT公司,A类2,H(H)n个2,C类O(运行)n个M(M)n个,适用于所有样本大小。DFR功能最差的测试是T型A类1,TKL公司,T型4,K(K)n个,H(H)n个(1),J型,J型0.5COV公司.

对于BFR分布,我们的测试具有最强大的功能真实航向,TJ公司,TT公司,A类2,H(H)n个(2)M(M)n个什么时候n个 = 10.何时n个 = 20分,50分,我们提出的测试优于所有竞争对手。我们还注意到E类E类L(左)R(右)1E类E类L(左)R(右)2适用于中等样本量(n个 = 20、50),当样本量较小时,他们没有表现得比其他人更好(n个 = 10). BFR最不强大的测试是T型A类1,TKL公司,K(K)n个,H(H)n个(1),J型,J型0.5COV公司.

一般来说,我们从表中看到49测试的性能取决于替代分布,没有任何一个测试能够击败其他所有替代测试。我们提出的测试对IFR和BFR分布具有很好的调整能力,特别是在样本量较小的情况下,效果很好。即使对于DFR和UFR替代分布,我们提出的测试也不总是最差的。因此,我们希望拟议的测试能够在实际应用中得到很好的应用。

4.真实数据应用

在本节中,我们使用三个实际数据集来展示我们提出的测试统计的适用性和稳健性。第一个数据集符合指数分布,而其余两个数据集则符合Chen分布和均匀分布。当样本呈指数分布时,所有竞争对手的测试都表现良好,第一个示例将简单演示我们的测试统计是如何工作的。为了进行比较,我们采用了一些经典统计数据,如Kolmogorov-Smirnov的K(K)S公司n个,柯伊伯氏V(V)克拉梅·冯·米塞斯W公司n个2和Anderson-Darling的A类2它被从业者广泛使用,并在许多统计软件中可用,并且是IFR替代方案的一些最强大的竞争对手,即。T型A类1,H(H)n个(1)、和K(K)n个。在剩下的两个示例中,我们将把我们的建议与这些测试统计数据进行比较。所有测试均在5%标称水平下进行,所有模拟均采用10000次重复。

  • 示例1:年清洁梯度地下水监测井的氯乙烯数据(μ/L(左))

在第一个例子中,我们借用了Bhaumik和Gibbons的数据[7]它代表了从清理梯度监测井获得的氯乙烯数据,如下所示:

5.1,1.2,1.3,0.6,0.5,2.4,0.5,1.1,8.0,0.8,0.4,0.6,0.9,0.4,2.0,0.5,5.3,.2,2.7,2.9,2.5,2.3,1.0,0.2,0.1,0.1,1.8,0.9,2.0,4.0,6.8,1.2,0.4,0.2.

该数据集已用于许多研究,例如Krishnamoorthy等。[28]还有库尔卡尼和鲍瓦尔[29]用于建模寿命分布,如Shanker所示等。[48]用指数分布拟合得很好。这一结果得到了马兰奇和秦的验证[34]通过使用测试统计数据E类E类L(左)R(右)1E类E类L(左)R(右)2。我们展示了图中数据的QQ图1(a) 。该图还表明指数分布很适合此数据集。蒙特卡罗模拟结果表明,我们的测试基于E类ˆ2用一个第页-值0.7105无法拒绝零假设,因此表明氯乙烯数据确实符合指数分布。这表明我们的测试证实了马兰奇和秦的发现[34]和Shanker等。[48].

保存图片、插图等的外部文件。对象名称为CJAS_A_1840535_F0001_c.jpg

(a)氯乙烯数据,(b)电器故障数据的QQ图。

  • 例2:在寿命试验中,电器的1000次循环失效次数

在这个例子中,我们考虑了来自Lawless的数据[31]它表示电器在寿命试验中发生故障的1000 s周期数。该数据包括50个观察结果,如下所示:

0.014,0.034,0.059,0.061,0.069,0.080,0.123,0.142,0.165,0.210,0.381,0.464,0.479,0.556,0.574,0.839,0.917,0.969,0.991,1.064,1.088,1.091,1.174,1.270,1.275,1.355,1.397,1.477,1.578,1.649,1.702,1.893,1.932,2.001,2.161,2.292,2.326,2.337,2.628,2.785,2.811,2.886,2.993,3.122,3.248,3.715,3.790,3.857,3.912,4.100.

此示例的主要重点是监视第页-测试值,以比较拟议测试和本节开头提到的测试的性能。1(b) 给出了QQ图,这表明指数分布不适合此数据集。事实上,Yousaf等。[56]表明陈分布比指数分布更适合于拟合该数据集。我们可以从表中看到10我们建议的测试是最低的第页-值,这明确拒绝了5%显著性水平上的指数性无效假设。我们还注意到,除了K(K)n个,其他所有竞争对手都未能在此数据集上检测到指数分布和Chen分布之间的差异。如果我们选择1%作为显著性水平,即使是基于K(K)n个不能拒绝零假设。这进一步显示了我们测试的优势。

  • 示例3:转换后的氯乙烯数据

表10。

电器故障数据的指数检验结果。
测试 |E类ˆ2| K(K)S公司n个V(V) W公司n个2 A类2 T型A类1 H(H)n个(1) K(K)n个
第页-价值0.00550.18420.12850.09790.08990.20280.09700.0356

由于示例1中的数据集已被证明是指数分布的良好拟合,因此我们使用概率积分变换将数据转换为均匀分布。转换后的数据如下:

0.0518,0.0518,0.1009,0.1009,0.1917,0.1917,0.1917,0.2336,0.2336,0.2336,0.2733,0.2733,0.3467,0.3805,0.3805,0.4126,0.4431,0.4719,0.4719,0.4993,0.6162,0.6550,0.6550,0.7059,0.7211,0.7356,0.7623,0.7863,0.8178,0.8810,0.9337,0.9404美元,0.9732,0.9858

在最后一个例子中,我们有兴趣进行一项自举研究,以评估所提出的测试统计量在使用均匀分布数据检测偏离指数性方面的能力和稳健性。我们首先使用基于密度的经验似然比检验(dbEmpLikeGOF)来检查转换数据的一致性。dbEmpLikeGOF给出了一个第页-值为0.7443,这表明数据与均匀分布一致。为了进行功率比较,我们选择了表中使用的所有测试统计数据10。通过随机删除24、14和4个观察结果,我们进行了引导研究,并将结果显示在表中11结果清楚地表明,我们提出的测试比其竞争对手具有更大的优势和相当大的鲁棒性,因为我们提出的试验统计在所有模拟场景下都保持了显著的高功率。当样本量较小时,这种优势更加明显(n个 = 10).

表11。

引导电源结果。
删除的观察结果 |E类ˆ2| K(K)S公司n个V(V) W公司n个2 A类2 T型A类1 H(H)n个(1) K(K)n个
240.88970.17670.19370.23550.17710.39970.41700.4375
1410.35280.62820.55430.52100.86030.94470.9525
410.77350.94780.97830.98900.99640.98650.9943

5.结论

在本文中,我们提出了一种新的基于记录值外部性的指数性测试,并通过一组具有不同类型故障率函数的替代分布,将其性能与58种经典和最新的指数性测试进行了比较。使用蒙特卡罗方法进行了模拟研究,结果表明,我们提出的测试在IFR和BFR备选方案方面优于竞争对手的测试。即使在DFR和UFR替代方案的情况下,我们建议的测试也不总是最差的测试。我们还提供了三个实际例子,以表明所提出的测试可以自信地应用于实践。

补充材料

工作代码:

致谢

我们感谢两位匿名评审和一位副编辑的重要和建设性意见,他们改进了演示。本研究得到了国家自然科学基金(71871208,71971204)、安徽省自然科学基金会(1908085MG236,gxbjZD54)的资助。

附录。

导入数学

将numpy导入为np

定义calE2(样本,n,m):

sample.sort()

Junx=范围(1,n+1)内i的总和([np.log(1-i/(n+1))**(2*k-2)/(样本[min(i+m-1,n-1)]-样本[max(i-m-1,0)]))

Jlnx=范围(1,m+1)内i的总和([(1+(i-1)/m)/(样本[最小(i+m-1,n-1)]-样本[最大(i-m-1,0)])

Jlnx+=范围(m+1,n-m+1)内i的总和([2/(样本[最小(i+m-1,n-1)]-样本[最大(i-m-1,0)])

Jlnx+=i在范围(n-m+1,n+1)内的总和([(1+(n-i)/m)/(样本[最小(i+m-1,n-1)]-样本[最大(i-m-1,0)])

E2=-m/(n**2*math.gama(k)**2)*Junx+m/(2*n**2)*Jlnx

返回E2

资金筹措表

本研究得到了国家自然科学基金(71871208,71971204)、安徽省自然科学基金会(1908085MG236,gxbjZD54)的资助。

披露声明

提交人没有报告任何潜在的利益冲突。

工具书类

1Aboukhamseen S.M.和Aly E.E。,基于平均剩余寿命函数的指数检验,科威特科学杂志。 43(2016),第185-198页。[谷歌学者]
2Arnold B.C.、Balakrishnan N.和Nagaraja H.N。,记录约翰·威利父子公司,纽约,1988年。[谷歌学者]
三。Ascher S.公司。,指数性检验综述,通信统计。理论方法 19(1990年),第1811-1825页。doi:10.1080/03610929008830292[交叉参考][谷歌学者]
4Baratpour S.和Habibi Rad A。,基于累积剩余熵的指数分布优良性检验,通信统计。理论方法 41(2012),第1387–1396页。doi:10.1080/03610926.2010.542857[交叉参考][谷歌学者]
5Baringhaus L.和Henze N。,基于经验拉普拉斯变换的指数性一致性检验,Ann.Inst.Stat.数学。 43(1991年),第551–564页。doi:10.1007/BF00053372[交叉参考][谷歌学者]
6Baringhaus L.和Henze N。,基于平均剩余寿命函数表征的指数拟合检验,统计师。论文 41(2000),第225–236页。doi:10.1007/BF02926105[交叉参考][谷歌学者]
7Bhaumik D.K.和Gibbons R.D。,至少为第页属于在每个第页位置,技术计量学 48(2006),第112-119页。doi:10.1198/004017005000000355[交叉参考][谷歌学者]
8Cavus M.、Yazici B.和Sezer A。,当I型错误概率不同时,比较测试功率的惩罚功率法,通信统计。模拟。计算。(2019). 可于10.1080/03610918.2019.1588310获取。[交叉参考][谷歌学者]
9Choi B.、Kim K.和Song S.H。,基于Kullback-Leibler信息的指数性优度检验,通信统计。模拟。计算。 33(2004),第525-536页。doi:10.1081/SAC-120037250[交叉参考][谷歌学者]
10考克斯·D.R.和奥克斯·D。,生存数据分析查普曼和霍尔出版社,伦敦,1984年。[谷歌学者]
11D’Agostino R.B.和Stephens M.A。,fit Goodness技术马塞尔·德克尔,纽约,1986年。[谷歌学者]
12Deshpande有限公司。,一类针对递增失效率平均方案的指数检验,生物特征 70(1983年),第514-518页。doi:10.1093/biomet/70.2.514[交叉参考][谷歌学者]
13易卜拉欣·N、哈比布拉·M和苏菲·E·S。,基于Kullback-Leibler信息的指数性检验,J.R.统计社会服务。B统计方法。 54(1992年),第739–748页。[谷歌学者]
14Epps T.和滑轮L。,由经验特征函数导出的指数性与单调风险方案的检验,J.R.统计社会服务。B Methodol公司。 48(1986),第206–213页。[谷歌学者]
15爱泼斯坦B。,假设生命的基本分布是指数分布的有效性检验,技术计量学 2(1960),第83–101页。[谷歌学者]
16Fortiana J.和GranéA。,基于最大相关性的指数分布的无标度拟合优度统计量,J.统计计划。推断 108(2002),第85-97页。doi:10.1016/S0378-3758(02)00272-0[交叉参考][谷歌学者]
17Gail M.H.和Gastwirth J.L。,基于基尼统计量的指数分布无标度优良性检验,J.R.统计社会服务。B Methodol公司。 40(1978年),第350–357页。[谷歌学者]
18Gail M.H.和Gastwirth J.L。,基于Lorenz曲线的指数分布无标度优良性检验,美国统计协会。 73(1978年),第787-793页。[谷歌学者]
19Gnedenko B.V.,Belyayev Yu。K.和Solovyev A.D。,可靠性理论的数学模型,学术出版社,纽约,1969年。[谷歌学者]
20GranéA.和Fortiana J。,基于最大相关的指数分布的无位置和无标度优良度统计量,统计 43(2009),第1-12页。doi:10.1080/02331880802187667[交叉参考][谷歌学者]
21哈里斯·C。,关于指数性测试的注释,海军后勤部。 23(1976年),第169-175页。doi:10.1002/nav.3800230116[交叉参考][谷歌学者]
22亨泽N。,指数性的一类新的灵活综合检验,通信统计。理论方法 22(1992年),第115-133页。doi:10.1080/03610929308831009[交叉参考][谷歌学者]
23Henze N.和Meintanis S.G。,指数性的最新和经典测试:部分综述和比较,Metrika公司 61(2005),第29-45页。doi:10.1007/s001840400322[交叉参考][谷歌学者]
24Hollander M.和Proschan F。,测试新产品是否优于旧产品,安。数学。统计师。 43(1972年),第1136-1146页。doi:10.1214/网址/1177792466[交叉参考][谷歌学者]
25Jose J.和Sathar E.A。,的剩余极值k个-记录值,统计师。普罗巴伯。莱特。 146(2019年),第1-6页。doi:10.1016/j.spl.2018.10.109[交叉参考][谷歌学者]
26金伯A.C。,基于稳定概率图的指数分布、Weibull分布和Gumbel分布测试,生物特征 72(1985),第661-663页。doi:10.1093/biomet/72.3.661[交叉参考][谷歌学者]
27科查尔股份有限公司。,针对单调故障率平均值的指数性测试,通信统计。理论方法 14(1985),第381-392页。doi:10.1080/03610928508828919[交叉参考][谷歌学者]
28Krishnamoorthy K.、Mathew T.和Mukherjee S。,伽马分布的基于正态分布的方法:预测和公差区间以及应力强度可靠性,技术计量学 50(2008),第69-78页。doi:10.1198/00401700000353[交叉参考][谷歌学者]
29Kulkarni H.V.和Powar S.K。,威布尔分布的简单正态逼近及其在预测上限估计中的应用,J.概率。斯达。 2011(2011),第1-10页。doi:10.115/2011/863274[交叉参考][谷歌学者]
30Lad F.、Sanfilippo G.和Agro G。,外部性:熵的互补对偶,统计科学。 30(2015),第40-58页。doi:10.1214/14-STS430[交叉参考][谷歌学者]
31Lawless J.F。,寿命数据的统计模型和方法,卷。362霍博肯·威利(Wiley,Hoboken),2011年。[谷歌学者]
32Lilliefors H.W.公司。,均值未知指数分布的Kolmogorov-Smirnov检验,美国统计协会。 64(1969年),第387-389页。doi:10.1080/01621459.1969.10500983[交叉参考][谷歌学者]
33劳埃德·C.J。,根据尺寸调整的估计测试功率,J.统计计算。模拟。 75(2005),第921-933页。doi:10.1080/00949650412331321160[交叉参考][谷歌学者]
34Marange C.S.和Qin Y。,基于概率积分变换的指数性矩经验似然比检验,J.应用。斯达。 46(2019年),第2786–2803页。doi:10.1080/02664763.2019.1614542[交叉参考][谷歌学者]
35Mimoto N.和Zitikis R。,阿特金森指数、莫兰统计和指数检验,J.日本统计师。Soc公司。 38(2008),第187-205页。doi:10.14490/jjss.38.187[交叉参考][谷歌学者]
36Montazeri N.H.和Torabi H。,基于样本协方差的指数性测试,伊朗J.Stat.Res 10(2014),第147-161页。doi:10.18869/acadpub.jsri.10.2.147[交叉参考][谷歌学者]
37莫兰P.A.P。,间隔的随机划分——第二部分,J.R.统计社会服务。B Methodol公司。 13(1951年),第147-150页。[谷歌学者]
38Noughabi H.A.和Arghami N.R。,使用不同熵估计的五种指数检验的蒙特卡罗比较,J.统计计算。模拟。 81(2011),第1579-1592页。doi:10.1080/00949655.2010.496368[交叉参考][谷歌学者]
39Noughabi H.A.和Arghami N.R。,基于指数分布特征的指数性检验,J.统计计算。模拟。 81(2011),第1641-1651页。doi:10.1080/0949655.2010.498373[交叉参考][谷歌学者]
40Noughabi H.A.和Arghami N.R。,使用转换数据测试指数性,J.统计计算。模拟。 81(2011),第511-516页。doi:10.1080/00949650903348171[交叉参考][谷歌学者]
41Noughabi H.A.和Arghami N.R。,基于熵估计修正矩的优良性检验,通信统计。模拟。计算。 42(2013),第499–513页。doi:10.1080/03610918.2011.634535[交叉参考][谷歌学者]
42Noughabi H.A.和Balakrishnan N。,基于Phi散度的拟合优度检验,J.应用。斯达。 43(2016),第412-429页。doi:10.1080/02664763.2015.063116[交叉参考][谷歌学者]
43Noughabi H.A.和Jarrahiferiz J。,关于极值的估计,J.非参数。斯达。 31(2019年),第88–99页。doi:10.1080/10485252.2018.1533133[交叉参考][谷歌学者]
44Park S.和Park D。,基于两个熵估计的拟合优度检验矩修正,J.统计计算。模拟。 73(2003),第685-694页。doi:10.1080/0094965031000070367[交叉参考][谷歌学者]
45邱刚。,订单统计和记录值的极值,统计师。普罗巴伯。莱特。 120(2017),第52-60页。doi:10.1016/j.spl.2016.09.016[交叉参考][谷歌学者]
46邱刚、贾凯。,Extopy估计及其在一致性检验中的应用,J.非参数。斯达。 30(2018),第182-196页。doi:10.1080/10485252.2017.1404063[交叉参考][谷歌学者]
47邱刚、贾凯。,顺序统计量的剩余极值,统计师。普罗巴伯。莱特。 133(2018),第15-22页。doi:10.1016/j.spl.2017.09.014[交叉参考][谷歌学者]
48Shanker R.、Hagos F.和Sujatha S。,使用指数分布和Lindley分布对寿命数据建模,生物。生物学家。国际期刊。 2(2015),第140-147页。[谷歌学者]
49Shapiro S.S.和Wilk M.B。,指数分布的方差检验分析(完全样本),技术指标 14(1972年),第355-370页。网址:10.1080/00401706.1972.10488921[交叉参考][谷歌学者]
50斯蒂芬斯文学硕士。,关于原点已知的指数性的W检验,技术计量学 20(1978年),第33-35页。网址:10.1080/00401706.1978.10489614[交叉参考][谷歌学者]
51Tchirina股份有限公司。,基于Moran统计量的指数检验的Bahadur效率和局部最优性,数学杂志。科学。 127(2005),第1812-1819页。文件编号:10.1007/s10958-005-0143-x[交叉参考][谷歌学者]
52Torabi H.、Montazeri N.H.和GranéA。,指数检验的广泛回顾和两个竞争性建议及其在可靠性方面的应用,J.统计计算。模拟。 88(2018),第108–139页。doi:10.1080/00949655.2017.1379522[交叉参考][谷歌学者]
53瓦西塞克·欧。,基于样本熵的正态性检验,J.R.统计社会服务。B类 38(1976年),第54–59页。[谷歌学者]
54沃尔科娃·K·。,基于Rossbergs特征的指数检验的渐近效率,数学杂志。科学。 167(2010),第486–494页。doi:10.1007/s10958-010-9934-9[交叉参考][谷歌学者]
55Volkova K.Y.和Nikitin Y.Y。,基于Ahsanullah特征的指数检验及其效率,数学杂志。科学。 204(2015),第42-54页。doi:10.1007/s10958-014-2185-4[交叉参考][谷歌学者]
56Yousaf F.、Ali S.和Shah I。,基于上记录值的Chen分布的统计推断,Ann.数据科学。 6(2019年),第831-851页。doi:10.1007/s40745-019-00214-7[交叉参考][谷歌学者]
57Zamanzade E.和Arghami N.R。,基于修正熵估计修正矩的优良性检验,J.统计计算。模拟。 81(2011),第2077–2093页。doi:10.1080/00949655.2010.517533[交叉参考][谷歌学者]
58张杰和Boos D.D。,蒙特卡罗实验中的调整功率估计,通信统计。模拟。计算。 23(1994),第165-173页。doi:10.1080/03610919408813162[交叉参考][谷歌学者]

文章来自应用统计学杂志由以下人员提供泰勒和弗朗西斯