改进的T-SPRT在电子商务产品质量检验中的应用
Application of Improved T-SPRT in E-Commerce Product Quality Inspection
DOI: 10.12677/ecl.2024.132168, PDF, HTML, XML, 下载: 55  浏览: 171 
作者: 龙荣进, 叶茂越:贵州大学数学与统计学院,贵州 贵阳;胡思贵:贵州医科大学生物与工程学院,贵州 贵阳
关键词: 序贯检验截尾序贯概率比检验截尾序贯检验Sequential Test Truncated Sequential Probability Ratio Test Truncated Sequential Test
摘要: 为保证消费者合法权益以及商家和企业的长久利益,本文对质量指标服从正态分布的电子商务产品在均值参数已知时,对标准差σ参数的截尾序贯检验方案进行了研究。给出了σ参数关于截尾序贯概率比检验(T-SPRT)的截尾序贯检验方案T的定义,方案犯两类错误的概率、平均试验样本量以及综合平均实验样本量等评价指标的计算表达式。并从截距项对T-SPRT方案T进行改进,给出了改进方案Tt的搜索步骤,并利用仿真实验验证方案Tt和评价指标的精确性。仿真结果表明,其仿真实验评价指标值与统计计算得到的评价指标值基本保持一致,表明本文建立的统计指标计算表达式和Tt方案是正确的。与T-SPRT方案T相比,改进方案Tt能更好的控制犯两类错误的概率,保证了实验双方的利益。同时,在样本量截尾值相等时,Tt方案拥有更少的综合平均实验样本量,其下降的比率在10%左右,使得实验成本得到了大幅度减少。
Abstract: In order to ensure the legitimate rights and interests of consumers and the long-term interests of merchants and enterprises, this paper studies the truncation sequential test scheme of the standard deviation σ parameter of e-commerce products whose quality indicators obey normal distribution when the mean parameter is known. The definition of the σ-parameter about the truncated sequential test plan of the truncated sequential probability ratio test (T-SPRT), the probability of the plan making two types of errors, the average test sample size, and the comprehensive average test sample size and other evaluation indicators are given. And the T-SPRT scheme T is improved from the intercept term, and the search steps for the improved scheme Tt are given. Simulation experiments are used to verify the accuracy of the scheme Tt and the evaluation index. The simulation results show that the evaluation index value of the simulation experiment is basically consistent with the evaluation index value obtained by statistical calculation, indicating that the statistical index calculation expression and Tt scheme established in this article are correct. Compared with T-SPRT scheme T, the scheme Tt can better control the probability of making two types of errors, ensuring the interests of both parties in the experiment. At the same time, when the sample size cut-off values are equal, the plan has a smaller comprehensive average experimental sample size, and the reduction rate is about 10%, which greatly reduces the experimental cost.
文章引用:龙荣进, 胡思贵, 叶茂越. 改进的T-SPRT在电子商务产品质量检验中的应用[J]. 电子商务评论, 2024, 13(2): 1368-1377. https://doi.org/10.12677/ecl.2024.132168

1. 引言

电子商务让消费者通过网络实现了足不出户可网上购物和网上支付。与传统线下购物相比,电子商务节省了客户与企业、消费者和商家的时间和空间,大大提高了交易效率 [1] 。现如今,随着网络渠道(如淘宝、京东、抖音电商、天猫、拼多多以及多多买菜等)层出不穷,中国电商交易规模呈线性增长,其中,2016年到2022年中国电商交易规模统计见图1。然而,在电子商务迅猛发展的同时也存在一系列的问题,如电子商务产品质量参差不齐 [2] 、假货率高 [3] 、直播退货率高 [4] 、质量满足不了客户和消费者要求以及产品信息虚假等问题,其实质为产品质量不过关。电子商务产品作为电子商务的核心内容,在如此“卷”的局面中,满足消费者需求才是核心竞争力 [5] 。为保证消费者合法权益不受侵害以及从商家和企业长久利益来看,在产品售卖前和产品验收时对其进行质量检验实验是十分必要的。因此,如何设计一个优良的抽样检验方案成为降低实验成本的关键。

在实际生活中,很多产品的质量指标服从或近似服从正态分布,非常普遍,如机械零部件的尺寸,药品的有效性指标以及设备的误差等。对于该类产品标准差 σ 参数的抽样检验,常用的检验方法有传统固定实验样本量的抽样检验(SN检验)和截尾序贯概率比检验(T-SPRT) [6] [7] [8] [9] 。序贯检验采用“试试看,看试试”的实验策略,由于在实验过程中就能做出统计决策,因此,与SN检验相比,序贯检验方案能大幅度地缩减产品抽样检验的平均试验样本量或平均试验时间 [10] [11] [12] [13] 。

Figure 1. Statistics on the scale of e-commerce transactions in China (100 million yuan)

图1. 中国电商交易规模统计(亿元)

众所周知,在第二次世界大战后期,为提高军工产品的抽样效率和降低抽样成本,Wald提出了一种一般性的序贯检验方法,即序贯概率比检验(SPRT) [14] 。SPRT检验示意图如图2所示,其检验方案由两条平行的直线构成。由于SPRT继续实验区为开放式区域,使其缺乏最大实验样本量的上限值 [10] 。这对试验预算试验成本和控制都带了困难。因此,在实际应用中采用的序贯检验方案一般都是截尾的,例如T-SPRT、2-SPRT、SMT等等。T-SPRT因其构造简单,使用范围广,例如国际标准IEC61124、俄罗斯国准GOSR27.402、我国国家标准GBT2828、国际标准IEC1123等均采用T-SPRT,T-SPRT检验示意图如图3所示。需要说明的是,T-SPRT仍有不足,该检验方案存在检验水平利用不充分或超用的现象。同时,当给定检验水平相同时,方案犯两类错误的概率相差值较大,平均实验样本量达不到最优。

为此,本文对质量指标服从正态分布的产品在均值参数已知时,对标准差 σ 参数的抽样截尾序贯检验进行了研究,给出了改进T-SPRT方案的方法步骤,使得在误差允许范围内对应的截尾序贯检验方案拥有更少的综合平均实验样本量,从而达到提高产品抽样效率,减少抽样检验试验成本的效果。

本文余下内容安排如下:第2节介绍正态分布下标准差 σ 参数的SN检验、SPRT检验以及T-SPRT检验;第3节给出正态分布下标准差 σ 参数的T-SPRT方案的定义,以及建立方案的评价指标体系,在T-SPRT的基础上,对T-SPRT进行改进,给出改进方案的搜索步骤,并进行仿真实验验证其方案的可靠性;第4节将改进的T-SPRT方案与T-SPRT方案进行对比分析;第5节为总结部分。

2. 基本检验方法

设X为产品的质量指标,服从均值为 μ ,标准差为 σ 的正态分布,即 X ~ N ( μ , σ ) 。其中 μ 已知, σ 未知,对于该类产品标准差 σ 参数的计量型抽样检验,常讨论的假设检验为:

H 0 : σ = σ 0 vs H 1 : σ = σ 1 = λ σ 0 ( λ > 1 ) . (1)

如(1)式的计量型检验问题,常采用的检验方法有固定样本量的检验法(SN检验)和截尾序贯概率比检验(T-SPRT)这两种。下面对这两种检验方法进行介绍。

2.1. SN检验

若有检验样本 X 1 , X 2 , , X k ,记检验统计量 S k = i = 1 k ( x i μ ) 2 ,则 S k ~ G a ( k 2 , 1 2 σ 2 ) ,密度函数如(3)式所示:

f ( s k ) = { 1 ( 2 σ 2 ) k / 2 Γ ( k / 2 ) e s k 2 σ 2 s k k 2 1 , s k >0 , 0 , . (2)

设拒绝原假设的临界点为 U k 0 ,则拒绝域为 { S k 0 U k 0 } ,给定检验水平 ( α 0 , β 0 ) ,可建立方程组:

{ α 0 = U k 0 + f ( s k 0 ; σ 0 ) d s k 0 β 0 = 0 U k 0 f ( s k 0 ; σ 1 ) d s k 0 (3)

取定检验水平 ( α 0 , β 0 ) 后,求解(3)式即可得到 ( k 0 , U k 0 ) 。此时,必须对 k 0 个样本进行检验,才能做出实验决策,故该方法称为固定样本量的检验,即SN检验。

2.2. SPRT检验

f ( x i ; σ ) 为序贯检验第i次实验时样本 X i 的密度函数, L R k 为似然比,计算公式如下:

L R k = i = 1 k f ( x i ; σ 1 ) i = 1 k f ( x i ; σ 0 ) = λ k E x p { S k ( 1 2 σ 0 2 1 2 σ 1 2 ) } (4)

取满足条件 0 < A < 1 < B 的常数A和B,对不等式 A < L R k < B 进行对数变换,即可得到序贯概率比检验(SPRT)方案接受线l和拒绝线u的表达式,其表达如下:

{ l = a k + b 1 u = a k + b 2 (5)

其中,

A = β 0 1 α 0 , B = 1 β 0 α 0 , a = 2 ln λ 1 / σ 0 2 1 / σ 1 2 b 2 = 2 ln B 1 / σ 0 2 1 / σ 1 2 , b 1 = 2 ln A 1 / σ 0 2 1 / σ 1 2 (6)

此时检验方法的判断规则为:进行到第k次实验,若 S k a k + b 1 ,则接受原假设,实验结束;若 a k + b 1 < S k < a k + b 2 ,则继续实验;若 S k a k + b 2 ,则拒绝原假设,实验结束。该方法称为序贯概率比检验(SPRT),其检验示意图见图2

Figure 2. Schematic diagram of SPRT inspection plan

图2. SPRT检验方案示意图

2.3. T-SPRT

T-SPRT是在SPRT的基础上进行了强制截尾。当给定样本量截尾值 r 0 和成功判别值 u r 0 * 后,作直线 k = r 0 与接受线l和拒绝线u相交,作直线 S k = u r 0 * 与拒绝线u和直线 k = r 0 相交。此时,SPRT继续实验区域得到了限制,得到的检验方案即为T-SPRT检验方案,其检验示意图见图3

Figure 3. Schematic diagram of T-SPRT inspection plan

图3. T-SPRT检验方案示意图

当给定 r 0 后, u r 0 * 的取值可取为 [15] :

u r 0 * = l r 0 + u r 0 2 (7)

T-SPRT检验方法的序贯检验准则为,当 k = 1 , 2 , , r 0 1 时:若 S k a k + b 1 ,则接受原假设,实验结束;若 S k min { a k + b 2 , u r 0 * } ,则拒绝原假设,实验结束。若 a k + b 1 < S k < min { a k + b 2 , u r 0 * } ,则继续实验;当 k = r 0 时:若 S r 0 < u r 0 * ,接受原假设,实验结束;若 S r 0 u r 0 * ,拒绝原假设,实验结束。

3. 改进的T-SPRT

3.1. 方案评价指标

为了更好地评价T-SPRT算法,根据T-SPRT检验示意图和检验准则,定义关于T-SPRT的截尾序贯检验方案如下:

定义1设 l 1 , l 2 , , l r 0 u 1 , u 2 , , u r 0 为(5)式计算出的两列单调不减的实数列。当 l i < 0 ,令 l i = 0 。设存在非零常数 u r 0 * 满足 l r 0 < u r 0 * < u r 0 。若时,令 u i = u r 0 * ,若 u i u r 0 * ,令 u i = u i ,则称

T ( r 0 ) = ( u 1 u 2 u r 0 * u r 0 * 0 0 l r 0 1 u r 0 * ) (8)

为样本量截尾值为 r 0 ,实验成功判别值为 u r 0 * 时,关于T-SPRT的截尾序贯检验方案,简记为T。

为第k次实验时,接受原假设的情况。同理可得为第k次实验时,拒绝原假设的情况。此时,方案T犯两类错误的实际概率计算公式如下:

(9)

其中, k f ( s k ; σ 0 ) 表示截尾序贯检验方案上的概率转移密度函数,其计算公式如下:

(10)

M = min { k | S k l k or S k u k , k = 1 , , r 0 } 。此时,方案T在 σ 0 σ 1 成立时的平均实验样本量(平均实验次数) E σ 0 ( M | T ) E σ 1 ( M | T ) 的计算公式如下所示:

(11)

E σ 0 ( M | T ) E σ 1 ( M | T ) 给予相同的权重得到平均实验样本量的综合指标 E ( M ¯ | T ) ,计算公式为:

E ( M ¯ | T ) = E σ 0 ( M | T ) + E σ 1 ( M | T ) 2 (12)

3.2. 改进T-SPRT的步骤

图3和公式(9)可知,在不改变斜率a和成功判别值 u r 0 * 的情况下,当截距项b2减小(增大)时,拒绝线u往下(上)移,此时拒绝域增大(减小),犯第I类错误的概率会增大(减少);当截距项b1减小(增大)时,接受线l往下(上)移,此时接受域减小(增大),犯第II类错误的概率会减少(增大)。据这一规律,本文考虑对截距项b1和b2分别进行上调和下调处理,具体步骤如下:

第1步:给定检验水平 ( α 0 , β 0 ) ,根据(3)式和(6)式确定初始值k0、a、b1、b2

第2步:令 r 0 = k 0

第3步:令 δ 1 = 1 为截距项b1初始的变化率;

第4步:令 b 1 * = δ 1 b 1 ,此时接受线l的表达式为 l * = a + b 1 * k ,根据该表达式计算出 k = 1 , 2 , , r 0 时对应的接受边界点;

第5步:令 δ 2 = 1 为截距项b2初始的变化率;

第6步:令 b 2 * = δ 2 b 2 ,此时拒绝线的表达式为 u * = a + b 2 * k ,根据该表达式计算出 k = 1 , 2 , , r 0 时对应的拒绝边界点;

第7步:按式(13)式计算成功判别值 u r 0 *

u r 0 * = l r 0 * + u r 0 * 2 (13)

第8步:根据定义1计算截尾序贯概率比检验方案T;

第9步:根据(9)式计算方案T犯第I类错误的概率和第II类错误的概率;

第10步:若 α ( T ) α 0 , β ( T ) β 0 (当给定的检验水平 α 0 = β 0 时,为保证生产方和使用方双方承担的风险相同,应再增加条件 | α ( T ) β ( T ) | < 0.005 )时,计算方案T所需的综合平均实验样本量,并储存于集合C,之后令 δ 2 = δ 2 0 .005 ,重复第6步到第10步,直到 δ 2 = 0 。若 α ( T ) > α 0 β ( T ) > β 0 ,则令 δ 2 = δ 2 0 .005 ,重复第6步到第10步,直到 δ 2 = 0

第11步:若 δ 2 = 0 ,则令 δ 1 = δ 1 0 .005 ,重复第4步到第11步,直到 δ 1 = 0

第12步:若 δ 1 = 0 ,计算集合C的长度L,若 L = 0 ,则令 k 0 = k 0 + 1 ,重复第2步到第12步,若 L > 0 ,则跳到第13步;

第13步:找到满足 α ( T ) α 0 , β ( T ) β 0 条件下综合平均实验样本量最小的方案,记为 T t

3.3. 应用及仿真验证

例1. 自动机床的加工精度,设计误差指标的均方差 σ 0 = 1 λ = 2.1 μ = 99.978 ,取 β 0 = 0.200 α 0 = 0.233 。给出下式检验问题的检验方案。

H 0 : σ 0 = 1 vs H 1 : σ 1 = 2.1 (14)

解:由(3)式计算SN检验所需的实验样本量 k 0 = 3 ,临界值 U k 0 = 4 .64 。再根据(6)式有 A = 0 .2913 B = 3 .8350 ,斜率 a = 1 .9190 ,截距项 b 1 = 3 .1906 b 2 = 3 .4767 。当 r 0 = k 0 = 3 时,根据改进的T-SPRT计算可得 r 0 = 3 b 1 * = 3 .1109 b 2 * = 1 .1125 ,对应的 T t ( 3 ) 方案如下:

T t ( 3 ) = ( 3 .03 4 .76 4 .76 0 .00 0 .73 4 .76 ) (15)

此时, T t ( 3 ) 方犯两类错误的实际概率以及综合平均实验样本量(综合平均实验次数)依次为 α ( T t ( 3 ) ) = 0 .199689 β ( T t ( 3 ) ) = 0 .232469 E ( M ¯ | T t ( 3 ) ) = 2 .188346

r 0 = 3 ,根据SPRT检验理论,计算可得SPRT检验方案 T ( 3 ) 为:

T ( 3 ) = ( 5 .40 5 .90 5 .90 0 .00 0 .65 5 .90 ) (16)

例2. 自动机床的加工精度,设计误差指标的均方差 σ 0 = 1 λ = 1.6536 μ = 0 ,取 β 0 = 0.1 α 0 = 0.1 。给出下式检验问题的检验方案。

H 0 : σ 0 = 1 vs H 1 : σ 1 = 1.6536 (17)

根据(3)式解出的 ( k 0 , U k 0 ) = ( 14 , 21.152 ) 。当 α 0 = β 0 = 0.1 ,取 r 0 = 14 时,由(6)式有 a = 1.5859 b 1 = 6 .9282 b 2 = 6 .9282 。公式(7)得此时的成功判别值 u 14 * = 22 .20 。根据改进的T-SPRT检验法求解出的 r 0 = 18 b 1 * = 6.8835 b 2 * = 4.8151 ,此时对应的成判别值 u 18 * = 25 .71 ,方案 T t ( 18 ) 表1所示。

Table 1. Scheme T t ( 18 )

表1. T t ( 18 ) 方案

表2展示的 T t ( 18 ) T t ( 3 ) 方案的评价指标值及其仿真实验结果,其中SD表示50组仿真实验值的标准差。由该表可知,改进的T-SPRT方案 T t ( 18 ) T t ( 3 ) 在犯两类错误的概率和综合平均实验样本量的统计计算值和仿真实验值基本一致,表明本文建立的评价指标计算表达式和改进的T-SPRT方案是正确的。此外, T t ( 18 ) T t ( 3 ) 方案犯两类错误的概率小于给定检验水平,保护了商家和消费者权益。

Table 2. The evaluation index values and simulation results of T '

表2. T ' 方案的评价指标值及其仿真结果

4. 对比分析

为了更全面直观的比较 T ( 3 ) T t ( 3 ) 方案以及 T ( 6 ) T t ( 6 ) 方案,绘制了这两组方案综合平均实验样本量 E ( M ¯ ) 随着待检验参数值 σ 变化的曲线图,其结果如图4所示。

Figure 4. Comparison between T and T t on E ( M ¯ )

图4. T、 T t 方案关于 E ( M ¯ ) 的比较

图4可知,当待检验参数值 σ 的真实取值在原假设和备择假设的之间时,与方案 T ( 3 ) T ( 18 ) 相比,方案 T t ( 18 ) 所需的综合平均实验样本量(综合平均实验次数)更少,节省的成本更多。

Table 3. Comparison between T ' and T

表3. T、 T ' 方案的比较

表3展示的是在不同检验问题和检验水平下,当样本量(实验次数)截尾值 r 0 相等时,求解出的 T ( r 0 ) T t ( r 0 ) 方案。根据该表可知,与 T ( r 0 ) 方案相比, T t ( r 0 ) 方案能更好的控制犯两类错误的概率,同时所需的综合平均实验样本量也较少,其下降的幅度在10%左右。

5. 总结

为保证消费者合法权益不受侵害以及从商家和企业长久利益来看,本文对质量指标服从正态分布的电子商务产品在均值参数已知时,对标准差参数的截尾序贯检验方案进行了研究。给出了关于标准差 σ 参数T-SPRT的截尾序贯检验方案的定义、方案犯两类错误的概率以及综合平均实验样本量等统计指标的计算表达式。并给出了改进的T-SPRT方案截距项的搜索步骤,并对得到满足要求的截尾序贯检验方案 T t 进行仿真实验。研究表明,其仿真实验评价指标值与统计计算得到的评价指标值基本保持一致,表明本文建立的评价指标计算表达式和改进的T-SPRT方案是正确的。同时,与原有的T-SPRT方案相比,改进的方案 T t 不仅满足了犯两类错误的概率小于给定检验水平的要求,在样本量截尾值相等时,其拥有的综合平均实验样本量(综合平均实验次数)也得到了下降,下降的幅度在10%左右。

参考文献

[1] 陈宇希, 朱豪. 电子商务产品质量研究文献可视化分析[J]. 品牌与标准化, 2024(1): 102-104.
[2] 信用中国. 协同治理理念下电子商务信用失范现象及对策[EB/OL].
https://credit.yinchuan.gov.cn/zxzlhlwly/20210805/1395.html, 2024-01-21.
[3] 吴靖文. 浅谈我国电商平台产品的质量问题及对策[EB/OL].
https://www.docin.com/p-2875358258.html, 2024-01-21.
[4] 信用中国. 重磅: 2021年电子商务领域“双十一”信用环境评价和提示来了[EB/OL].
http://credit.mhk.gov.cn/web/news/non/detail/dbf0edde-a36c-409d-ba88-7e8d62098db9, 2024-01-21.
[5] 昝妍. “产品不是核心竞争力, 满足用户需求才是”——访山西金蝉电子商务有限公司总经理韩智勇[J]. 中国招标, 2023(11): 69-72.
[6] 闫章更, 濮晓龙. 现代军事抽样检验方法及应用[M]. 北京: 国防工业出版社, 2008.
[7] 高明坤. 试验设备误差的统计假设检验法[J]. 数学的实践与认识, 1987(4): 55-62.
[8] 王燕飞. 正态分布产品的Bayes可靠性抽样检验方案[J]. 吉林化工学院学报, 2020, 37(5): 88-92.
https://doi.org/10.16039/j.cnki.cn22-1249.2020.05.020
[9] 吕建华, 吴启光. 正态分布下产品可靠性抽样检验方案——方差未知综合双侧情形[J]. 系统工程理论与实践, 2002(9): 63-69 96.
[10] 方茂达, 胡思贵, 李秋德, 陈慧娟. 指数分布下贝叶斯截尾序贯检验研究[J]. 运筹与模糊学, 2023, 13(1): 157-167.
[11] 陈慧娟, 胡思贵, 方茂达. 关于序贯检验相关研究的文献综述[J]. 运筹与模糊学, 2022, 12(4): 1251-1264.
https://doi.org/10.12677/ORF.2022.124132
[12] 胡思贵. 计数型截尾序贯检验的样本空间排序法[J]. 应用概率统计, 2013, 29(4): 200-213.
[13] 杨欣昱. 威布分布下高质量产品的序贯验收方法与成本优化[D]: [硕士学位论文]. 成都: 电子科技大学, 2021.
https://doi.org/10.27005/d.cnki.gdzku.2021.000846
[14] Wald, A. (1947) Sequential Analysis. Wiley, New York.
[15] 金振中. 二项分布的Bayes截尾序贯检验方法及其应用[J]. 战术导弹技术, 1996(1): 21-27.
https://doi.org/10.16358/j.issn.1009-1300.1996.01.005