CN111984922A - 基于dna分子标记的拟烟品检验抽样方法及其应用设备 - Google Patents

基于dna分子标记的拟烟品检验抽样方法及其应用设备 Download PDF

Info

Publication number
CN111984922A
CN111984922A CN202010620585.2A CN202010620585A CN111984922A CN 111984922 A CN111984922 A CN 111984922A CN 202010620585 A CN202010620585 A CN 202010620585A CN 111984922 A CN111984922 A CN 111984922A
Authority
CN
China
Prior art keywords
cigarette
sampling
cigarettes
dna molecular
total
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN202010620585.2A
Other languages
English (en)
Inventor
张轲
龙杰
孙浩巍
陈丹
童治军
张冀武
杨青
潘威
池源
杨艺敏
韩小渊
李郸
李海燕
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Yunnan tobacco quality supervision and inspection station
Original Assignee
Yunnan tobacco quality supervision and inspection station
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Yunnan tobacco quality supervision and inspection station filed Critical Yunnan tobacco quality supervision and inspection station
Priority to CN202010620585.2A priority Critical patent/CN111984922A/zh
Publication of CN111984922A publication Critical patent/CN111984922A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
    • C—CHEMISTRY; METALLURGY
    • C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12M—APPARATUS FOR ENZYMOLOGY OR MICROBIOLOGY; APPARATUS FOR CULTURING MICROORGANISMS FOR PRODUCING BIOMASS, FOR GROWING CELLS OR FOR OBTAINING FERMENTATION OR METABOLIC PRODUCTS, i.e. BIOREACTORS OR FERMENTERS
    • C12M33/00—Means for introduction, transport, positioning, extraction, harvesting, peeling or sampling of biological material in or from the apparatus
    • C—CHEMISTRY; METALLURGY
    • C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12Q—MEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
    • C12Q1/00—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
    • C12Q1/68—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving nucleic acids
    • C12Q1/6876—Nucleic acid products used in the analysis of nucleic acids, e.g. primers or probes
    • C12Q1/6888—Nucleic acid products used in the analysis of nucleic acids, e.g. primers or probes for detection or identification of organisms
    • C12Q1/6895—Nucleic acid products used in the analysis of nucleic acids, e.g. primers or probes for detection or identification of organisms for plants, fungi or algae
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00—Administration; Management
    • G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063—Operations research, analysis or management
    • G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/04—Manufacturing
    • C—CHEMISTRY; METALLURGY
    • C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12Q—MEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
    • C12Q2600/00—Oligonucleotides characterized by their use
    • C12Q2600/156—Polymorphic or mutational markers
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
    • Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
    • Y02P90/30—Computing systems specially adapted for manufacturing

Landscapes

  • Engineering & Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Business, Economics & Management (AREA)
  • Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Human Resources & Organizations (AREA)
  • Organic Chemistry (AREA)
  • Wood Science & Technology (AREA)
  • Theoretical Computer Science (AREA)
  • Zoology (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Analytical Chemistry (AREA)
  • Strategic Management (AREA)
  • Biotechnology (AREA)
  • Economics (AREA)
  • Proteomics, Peptides & Aminoacids (AREA)
  • General Engineering & Computer Science (AREA)
  • Data Mining & Analysis (AREA)
  • General Health & Medical Sciences (AREA)
  • Tourism & Hospitality (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Genetics & Genomics (AREA)
  • Computational Mathematics (AREA)
  • Biochemistry (AREA)
  • Development Economics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Educational Administration (AREA)
  • Pure & Applied Mathematics (AREA)
  • Molecular Biology (AREA)
  • Mathematical Physics (AREA)
  • Microbiology (AREA)
  • Marketing (AREA)
  • General Business, Economics & Management (AREA)
  • Operations Research (AREA)
  • Algebra (AREA)
  • Game Theory and Decision Science (AREA)
  • Botany (AREA)

Abstract

本发明公开了基于DNA分子标记的拟烟品检验抽样方法及其应用设备，属于烟草卷烟产品鉴别检验技术领域。包括以下步骤：S1:确定“弃真”错误α、“存伪”错误β和总体判定阈值(总体卷烟所占比例)；S2:定义总体卷烟支数比例为p，则从总体中取得卷烟的支数X符合0～1分布；S3:从总体中抽取容量为n的样本，Yn近似服从正态分布，即Yn～N(0，1)；S4:定义容量为n的样本中检出卷烟支数为m，计算出样本的平均值；S5:计算出两个实根
Figure DDA0002565031440000011
总体中卷烟支数比例为p的置信概率为1‑α，置信区间为
Figure DDA0002565031440000012
S6：确定最小取样量nmin；S7：确定最大取样数nt。该方法可有效保证DNA分子标记鉴别检验的时效性，节省检验成本，仅检验较少的样品的即可判定总体样本情况。

Description

基于DNA分子标记的拟烟品检验抽样方法及其应用设备
技术领域
本发明涉及烟草卷烟产品鉴别检验技术领域，具体涉及基于DNA分子标记的拟烟品鉴别检验抽样方法及其应用设备。
背景技术
近年来，除了电子雾化烟和加热不燃烧烟草制品的流行，所谓的花烟、茶烟、药烟等外形与卷烟相差无己的拟烟产品也出现在了消费者视线中。此类产品采用卷烟包装和产品样式，原料为茶丝(茶烟)、花丝(花烟)、中药叶丝(药烟)，部分产品甚至直接采用烟用辅料进行卷制和包装。以戒烟、卷烟替代品、健康保健等为宣传口号，更甚者将产品宣传为可以戒烟的“烟”，可像卷烟产品一样吸食，但不具有卷烟含有的烟碱，也可冲泡作为花饮、茶饮、药饮。大量充斥于旅游产品、伴手礼，甚至卷烟零售户的商铺进行销售，种类繁杂，品牌多样，有一定的回购率和消费群体。
但是，此类拟烟产品尚缺乏明确的行业标准，产品宣传的健康安全也存在争议。并且借卷烟打擦边球营销，扩大产品影响力，对卷烟市场正常的经营秩序已经产生了一定影响。近期又有未成年人吸食和茶烟、药烟中检出烟碱(尼古丁)成分的报道，缺乏有效的市场监管。
目前，鉴别检验拟烟产品主要原理是测定样品中的烟碱和N-亚硝胺含量，来检验判定此类产品是否为烟草专卖品。但烟碱并非烟草所独有，很多茄科植物也含有，并且烟草特有的N-亚硝胺的含量极低，基本都在纳克/克的水平，常出现卷烟样品中检不出烟草特有N-亚硝胺的情况，因此，利用该方法用于鉴别检验“拟烟产品”存在较大误判风险。
而对于抽样方法，现行一般采用《卷烟产品鉴别检验规程》进行抽样，先根据条装或盒装的商标品牌、品牌规格、包装形式等进行分类。再根据表1中规则进行抽样。
若采用该抽样方法进行基于DNA分子标记的“拟烟产品”鉴别检验，在 95％置信水平下，对于现行的抽样方案给出的结论为抽样样本的卷烟比例，而实际上该样本推断批次总体卷烟比例及其置信区间如表1所示。由表1可知现行方案最小一次抽样数量为80支，随着产品基数增加，最小抽样量增加至400 支，不利于控制抽样检测费用。并且，当抽样量较小时，置信区间过宽，说明该抽样方式对总体的估计能力较差，无法有效保障抽样精度。再者，由于此类产品供货方式多为送货制，销售点基本没有成件的存货，也不具备上述的抽样条件，因此，迫切需要一种科学合理的改进抽样方法。
表1《卷烟产品鉴别检验规程》抽样估计置信区间
Figure RE-GDA0002718289260000021
*注：设定当检出超过5％的抽检烟支为卷烟时，不能拒绝总体批次为卷烟。
针对上述已有技术状况，本发明申请人做了大量反复而有益的探索，最终产品取得了有效的成果，并且形成了下面将要介绍的技术方案。
发明内容
为此，本发明提供了基于DNA分子标记的拟烟品检验抽样方法及其应用设备，以解决现有技术中在对拟烟产品的检验抽样中存在所需样品量大、时效性低、抽样精度较差的问题。
为了实现上述目的，本发明提供如下技术方案：
基于DNA分子标记的拟烟品检验抽样方法，包括以下步骤：
S1：确定“弃真”错误α、“存伪”错误β和总体判定阈值(总体卷烟所占比例)；
S2：定义总体卷烟的支数比率为p，则从总体中取得卷烟的支数X符合 0～1分布；可得公式如下：
P(X＝x)＝px(1-p)1-x，x＝0或1；
S3:从总体中抽取容量为n的样本；Yn近似服从正态分布，即Yn～N(0，1)；
S4:定义容量为n的样本中检出卷烟支数为m，则样本的平均值为
Figure RE-GDA0002718289260000031
设定置信概率为1-α，则有
Figure RE-GDA0002718289260000032
其中，满足P(Yn>Zα)＝α的点Zα为服从正态分布的Yn的上α分位点；
S5:由
Figure RE-GDA0002718289260000033
得
Figure RE-GDA0002718289260000034
展开，可得到 ap2+bp+c＜0，其中，
Figure RE-GDA0002718289260000035
xi＝0或1(i＝1，2，…，n)，由于
Figure RE-GDA0002718289260000036
于是可得，
Figure RE-GDA0002718289260000037
设二次三项式ap2+bp+c的两个实根为
Figure RE-GDA0002718289260000038
则总体中卷烟支数比率为p的置信概率为1-α的置信区间为
Figure RE-GDA0002718289260000039
S6：确定最小取样量nmin；
S7：确定最大取样数nt。
进一步地，步骤S6的具体过程包括如下：
设定α＝0.05，β＝0.10，在置信水平95％(1-α)时，抽检方认为产品总体含低于5％的卷烟支数时判断该批次为拟烟产品；由概率统计理论知，在同样的可靠性和精确度下，样本越小，合格的允许样本中含有卷烟的概率越小；当使用序贯抽样方案时，最小累积卷烟允许值为0，即逐个检验没有发现卷烟，所预测的产品总体卷烟比率的单尾置信上限
Figure RE-GDA0002718289260000041
为5％，此时Za＝1.645，根据数学模型反算得到此时的样本量为51.4053，向上取整数，则最小抽样量nmin＝52；
此时根据
Figure RE-GDA0002718289260000042
求得两个实根分别为0和4.95％，即当抽取n＝52 只拟烟产品，经DNA分子标记检测，若未检出卷烟烟支，可以有95％的把握确定该批产品总体卷烟烟支比率在0～4.95％之间。
进一步地，步骤S7的具体过程包括如下：
在实际鉴别检验过程中很可能会出现多次检验均未达到判定阈值的情况，即累计的卷烟支数总是介于序贯抽样表所列的上下限之间而难以确定其总体专卖品比例，使检验不能终止；因此，需要计算最大抽样数nt来终止抽样过程；
根据设定的α＝0.05，β＝0.10，在GB/T8051-2008《计数序惯抽样检验方案》附表1中查得截尾拒收数为25，即最大累积卷烟支数拒收值为25，此时所预测的产品总体卷烟比率的双尾置信下限
Figure RE-GDA0002718289260000043
为5％，由于Zβ＝1.654，根据S2-S5 数学模型反算得到此时的样本量为363.3349，向下取整得到最大抽样数 nt＝363.3349，向下取整数得到最大抽样数nt＝363。
进一步地，还包括S8，步骤S8的具体过程包括如下：
将抽取的若干烟支同时切断，并使其内烟丝同步落入对应取样瓶内，将取样瓶内烟丝与烟草基因组数据做烟草特异性DNA分子标记的鉴别检验。
应用基于DNA分子标记的拟烟品检验抽样方法的设备，包括底座体、操作台、取样室及取样管；其中，
所述底座体与所述操作台之间留设有放置取样管的空间，以形成取样室；
所述操作台的顶端面开设有至少一个烟支容纳槽；
所述烟支容纳槽的外沿均开设有一个切断容置槽；
所述烟支容纳槽内在对应所述切断容置槽的位置开设有落烟口，所述操作台的底端对应所述落烟口的外沿设有连接凸起；
所述连接凸起与所述取样管的管口相固接；
所述操作台的一侧端转动设有与所述切断容置槽相对应的切刀。
进一步地，所述底座体固接设有向上延伸的支撑柱，所述支撑柱在远离所述底座体的一端与操作台相固接。
进一步地，所述烟支容纳槽呈长条弧状槽且具有一个延伸方向；
所述烟支容纳槽沿其延伸方向的一端与操作台的边沿对应设有槽口，所述烟支容纳槽沿其延伸方向的另一端位于操作台且对应设有限位端。
进一步地，所述限位端与所述切断容置槽之间的长度不小于3cm，切刀能够完全切除烟支滤嘴。
进一步地，所述切刀的侧端均匀固接设有至少一个侧刃，所述侧刃与所述烟支容纳槽相对应。
进一步地，所述连接凸起与取样管之间的固接方式为螺合或过盈配合。
本发明具有如下优点：
1、该方法可有效保证时效性，节省检测成本。
2、鉴于拟烟产品实际销售的方式和DNA分子标记鉴别检验技术的特点，可检验较少的样品的就可以判定总体样本情况，提升了抽样精度。
附图说明
为了更清楚地说明本发明的实施方式或现有技术中的技术方案，下面将对实施方式或现有技术描述中所需要使用的附图作简单地介绍，本说明书所绘示的结构、比例、大小等，均仅用以配合说明书所揭示的内容，以供熟悉此技术的人士了解与阅读，任何结构的修饰、比例关系的改变或大小的调整，在不影响本发明所能产生的功效及所能达成的目的下，均应仍落在本发明所揭示的技术内容得能涵盖的范围内。
图1为本发明实施例提供的方法流程示意图。
图2为本发明实施例提供的应用设备的整体结构图。
图3为本发明实施例提供的以ncum为52，110，142以及对应的《卷烟产品鉴别检验规程》一次抽样的n为40，200，400的OC曲线图。
图4为本发明实施例提供的初次抽样(nmin＝52)烟草特异性SSR引物扩增结果示意图(M为DNA Marker分子量为100-1500bp，S为卷烟标准样品扩增结果，n1～n52为初次抽样“拟烟产品”样品1～52扩增结果)。
附图中，各标号所代表的部件列表如下：
底座体1、操作台2、取样室3、烟支容纳槽4、限位端41、落烟口42、切断容置槽5、切刀6、侧刃7、取样管8。
具体实施方式
以下由特定的具体实施例说明本发明的实施方式，熟悉此技术的人士可由本说明书所揭露的内容轻易地了解本发明的其他优点及功效，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
如图1所示，本发明实施例提供了一种基于DNA分子标记的拟烟品检验抽样方法，基于DNA分子标记的鉴别检验成本较高的技术特性，本实施例选择可以节省平均样本量的序贯抽样方案。但既然是抽样检验，必然会产生两类错误，一类错误是把“拟烟产品”判定为卷烟，称之为第一类错误或“弃真”错误，其概率用α表示。另一类错误是把含烟草特征性成分的“拟烟产品”判定为非卷烟，称之为第二类错误或“存伪”错误，其概率用β表示。
具体方法步骤包括如下：
S1:为兼顾抽样精度和成本，本发明实施例中采用GB/T8051-2008《计数序惯抽样检验方案》6.1常规组合α＝0.05和β＝0.10。从统计学角度出发只有在一定置信水平，推断出总体一定数量(比例)的产品为卷烟或“拟烟产品”才能判定该鉴别检验案例为卷烟或“拟烟产品”。本发明设定在1-α＝95％的置信水平下，推断出总体的卷烟比例为5％该批产品即判断为烟草专卖品。
S2:定义某待检的拟烟产品中总体卷烟的支数比率为p，则从总体中任取一支“拟烟产品”取得卷烟的支数X符合0-1分布；可得公式如下：
P(X＝x)＝px(1-p)1-x，x＝0或1。
S3:从总体中抽取容量为n的样本X1，X2…Xn；当n较大时，由中心极限定理可知，随机变量
Figure RE-GDA0002718289260000071
Yn近似服从正态分布，即Yn-N(0，1)。
S4:定义容量为n的样本中检出卷烟支数为m，则样本的平均值为
Figure RE-GDA0002718289260000072
用此来估计总体的卷烟支数m所占比率，由于设定置信概率为1-α，则有
Figure RE-GDA0002718289260000073
其中，满足P(Yn>Zα)＝α的点Zα为服从正态分布的随机变量Yn的上α分位点。
S5:由
Figure RE-GDA0002718289260000074
得
Figure RE-GDA0002718289260000075
将该式展开，可得到 ap2+bp+c＜0，其中，
Figure RE-GDA0002718289260000076
xi＝0或1(i＝1，2，…， n)，由于
Figure RE-GDA0002718289260000078
于是可得，
Figure RE-GDA0002718289260000077
设二次三项式ap2+bp+c的两个实根为
Figure RE-GDA0002718289260000081
因此，样品总体中卷烟支数比率为p的置信概率为1-α的置信区间为
Figure RE-GDA0002718289260000082
Figure RE-GDA0002718289260000083
S6：确定最小取样量
设定在置信水平95％时，抽检方认为产品总体含低于5％的卷烟支数时判断该批次为拟烟产品。由概率统计理论知，在同样的可靠性和精确度下，样本越小，合格的允许样本中含有卷烟的概率越小。当使用序贯抽样方案时，最小累积卷烟允许值为0，即逐个检验没有发现卷烟，所预测的产品总体卷烟比率的单尾置信上限
Figure RE-GDA0002718289260000084
应为5％，此时Za＝1.645，根据数学模型反算得到此时的样本量为51.4053，向上取整数，则最小抽样量nmin＝52。因此，最小取样量(初始抽样量)nmin＝52。
此时根据
Figure RE-GDA0002718289260000085
求得两个实根分别为0和4.95％，即当抽取n＝52 只拟烟产品，经DNA分子标记检测，若未检出卷烟烟支，可以有95％的把握确定该批产品总体卷烟烟支比率在0～4.95％之间。
以上述设定的置信概率1-α＝95％，总体5％的卷烟支数来推断总体的卷烟支数。当第一样本数量为52时，由概率统计理论知，在同样的可靠性和精确度下，样本越小，合格的允许样本中含有卷烟制品的概率越小，此样本中有很大可能全为“拟烟产品”，没有卷烟，此时可以求出样品中卷烟制品的比率的置信上限
Figure RE-GDA0002718289260000086
可以推断出其总体的卷烟制品的比率在(0，
Figure RE-GDA0002718289260000087
)区间内。
由正态分布表查得，当α为0.5时，Zα为1.645，由S5中下列公式得出：
Figure RE-GDA0002718289260000088
c＝0，
根据
Figure RE-GDA0002718289260000089
求得两个实根分别为0和4.95％，即，当抽取n＝52 只拟烟产品，经DNA分子标记检测，若未检出卷烟烟支，可以有95％的把握确定该批产品总体卷烟烟支比率在0-4.95％之间。
其余情况该批“拟烟产品”总体的卷烟产品比率的推断详见附表，计算方法与举例置信区间的推断方法相同。表中不合格批次的判定是指如果达到相应的累计抽样数，其中卷烟支数达到表格中的阈值即判定该产品属于卷烟。
S7：确定最大取样数
在实际鉴别检验过程中可能会出现多次检验均未达到判定阈值的情况，即累计的专卖品支数总是介于序贯抽样表所列的上下限之间而难以确定其总体专卖品比例，使检验不能终止。因此，需要计算最大抽样数nt来终止抽样过程。根据设定的α＝0.05，β＝0.10，在GB/T8051-2008《计数序惯抽样检验方案》附表1中查得截尾拒收数为25，即最大累积卷烟支数拒收值为25，此时所预测的产品总体卷烟比率的双尾置信下限
Figure RE-GDA0002718289260000091
为5％，由于Zβ＝1.654，根据《计数序惯抽样检验方案》及步骤S2-S5数学模型得到此时的样本量为363.3349，向下取整得到最大抽样数nt＝363.3349，向下取整数得到最大抽样数nt＝363。
表2序贯抽样参考表(双尾)
Figure RE-GDA0002718289260000092
S8：将抽取的若干烟支同时切断，并使其内烟丝同步落入对应取样瓶内，将取样瓶内烟丝与烟草品种基因组数据做烟草特异性DNA分子标记鉴别检验。
抽样特性曲线(OC曲线)是基于抽样方案(ncum，c)绘制的，即从整体为N的产品中抽取样本数量为ncum的样本进行检验，ncum处于nmin与nt之间，。若符合判定标准c，则判定该批产品总体合格，否则判定不合格，通过OC曲线与抽样方案变化的关系很可以进一步验证抽样方案的合理性。以ncum为52， 110，142以及对应的《卷烟产品鉴别检验规程》一次抽样的n为40，200，400 的OC曲线为例，如图3。
由图3可知当ncum逐渐增大时，OC曲线有向左移动的趋势，曲线急剧下降，斜度变大，说明本研究中OC曲线具有相对较高的确定性和灵敏度，即当“拟烟产品”相对一致性较高且总体卷烟比例未超过规定界限时，大概率判定为“拟烟产品”，不会将其判定为卷烟。当其总体卷烟比例超过规定界限时也能准确地判定为卷烟，具有较少中间概率判定阶段。反观，一次抽样曲线，随着抽样量的增加先有向右移动的趋势，OC曲线交叉后才有向左移的趋势，随着抽样量的增大，OC曲线灵敏度和确定度也明显不如本研究的抽样曲线。
为了验证抽样方法在检测中的实际效果，将待检材料中的两条“拟烟产品”破除包装后，将烟支充分混匀，随机抽取52支烟支为初始样品，对照为随机抽取的2019年某品牌成品卷烟标准样品，如图4。
卷烟标准样品(S)在250bp分子量处扩增出了特异性条带。若初次抽样的 52个“拟烟产品”(n1～n52)样品均未在250bp分子量处扩增出特异性条带，即累计允许卷烟支数Ac＝0，初次抽样判定为“拟烟产品”(合格)的条件，即判定该产品为“拟烟产品”。
从现有的一次抽样方案来看，若将该方案应用于基于DNA标记的“拟烟产品”鉴别检验，较大的抽样数量，不适合个体检验准确，检测成本较高的 DNA分子标记鉴别检验。并且，本研究的目的是要通过样本推算总体的卷烟比例，一次抽样的方案当样本较小时，置信区间过宽，对总体的估计能力较差，无法有效保障抽样精度。同时，OC曲线也表明，本研究的抽样方案在较小抽样量的基础上，随着样品数量的增加比一次抽样具有更佳的灵敏度和确定度。
序贯抽样的较其他抽样方法的明显优势就是可降低平均样本量，从而降低检验成本和时间。但该抽样方法需要每次抽样结束后进行检验，判断是否进行下进一步的抽检。与随机抽样、分层抽样等方法相比，存在一定误判风险，对检验人员的专业技术熟练度要求相对较高。当检验接近判定阈值的批次时，会导致整个鉴别检验周期延长，需平衡时间和经济成本。但本研究中无论是卷烟还是“拟烟产品”均属于经工艺和机械加工后的工业产品，其质量有一定的稳定性，也就是说产品的相对一致性是较好的。所以结合应用实例，在实际检鉴别检验过程中，仅部分烟支含烟草特征性成分的“拟烟产品”产品很少，即混有卷烟的“拟烟产品”在生产过程中被严格控制，大多数的被检样品相对一致性较高，只需检验很小的样本就能通过。
如所示，本发明还提供一种应用基于DNA分子标记的拟烟品检验抽样方法的设备，包括底座体1、操作台2、取样室3、烟支容纳槽4、切断容置槽5、切刀6、侧刃7及取样管8，用以实现将待检测的烟支滤嘴进行切断，同时能使不同烟支取样后的烟丝在滤嘴切断时可以自动落入不同的取样瓶内，功能实用性强。具体设置如下：
所述底座体1的边沿固接设有向上延伸的支撑柱，所述支撑柱在远离所述底座体1的一端固接设有操作台2，用以借助支撑柱使得底座体1与操作台2 之间能够形成取样室3，取样室3可用于放置取样管8。
所述操作台2的顶端面均匀开设有若干个烟支容纳槽4，所述烟支容纳槽 4呈长条弧状槽且具有一个延伸方向；烟支容纳槽4沿其延伸方向的一端与操作台2的边沿对应设有槽口，烟支容纳槽4沿其延伸方向的另一端位于操作台2且对应设有限位端41，用以限制烟支滤嘴端的位置。
每个所述烟支容纳槽4的外沿均开设有一个切断容置槽5，若干个切断容置槽5之间呈线性对应，用以可通过切刀6同时作用于若干个切断容置槽5。
具体的是，所述限位端41与所述切断容置槽5之间的长度不小于3cm，使用时可将烟支滤嘴紧贴限位端41，,以保证切刀6能够完全切除烟支滤嘴。
所述烟支容纳槽4内在对应所述切断容置槽5的位置开设有落烟口42，所述操作台2的底端对应所述落烟口42设有连接凸起，用以可借助连接凸起与取样管8的管口之间相固接。
所述连接凸起与取样管8之间的固接方式包括但不限于螺合、过盈配合。
作为本实施例的一种优选方案，所述操作台2的一侧端转动设有切刀6，所述切刀6与呈线性对应的若干所述切断容置槽5分别对应；所述切刀6的侧端均匀固接设有若干个侧刃7，相邻两个所述侧刃7之间的距离与相邻两个烟支容纳槽4之间距离相等，用以通过侧刃7可在切断烟支滤嘴时，在烟支上同时切出一个纵开口，并可带动烟丝从纵开口下落，大大提升了功能实用性。
以上所述仅为本发明的较佳实施例，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (10)

1.基于DNA分子标记的拟烟品检验抽样方法，其特征在于，包括以下步骤：
S1：确定“弃真”错误α、“存伪”错误β和总体判定阈值(总体卷烟所占比例)；
S2：定义总体卷烟的支数比率为p，则从总体中取得卷烟的支数X符合0～1分布；可得公式如下：
P(X＝x)＝px(1-p)1-x，x＝0或1；
S3:从总体中抽取容量为n的样本；Yn近似服从正态分布，即Yn～N(0，1)；
S4:定义容量为n的样本中检出卷烟支数为m，则样本的平均值为
Figure FDA0002565031410000011
设定置信概率为1-α，则有
Figure FDA0002565031410000012
其中，满足P(Yn>Zα)＝α的点Zα为服从正态分布的Yn的上α分位点；
S5:由
Figure FDA0002565031410000013
得
Figure FDA0002565031410000014
展开，可得到ap2+bp+c＜0，其中，
Figure FDA0002565031410000015
xi＝0或1(i＝1，2，…，n)，由于
Figure FDA0002565031410000016
于是可得，
Figure FDA0002565031410000017
设二次三项式ap2+bp+c的两个实根为
Figure FDA0002565031410000018
则总体中卷烟支数比率为p的置信概率为1-α的置信区间为
Figure FDA0002565031410000019
S6：确定最小取样量nmin；
S7：确定最大取样数nt。
2.如权利要求1所述的基于DNA分子标记的拟烟品检验抽样方法，其特征是，步骤S6的具体过程包括如下：
设定α＝0.05，β＝0.10，在置信水平95％(1-α)时，抽检方认为产品总体含低于5％的卷烟支数时判断该批次为拟烟产品；由概率统计理论知，在同样的可靠性和精确度下，样本越小，合格的允许样本中含有卷烟的概率越小；当使用序贯抽样方案时，最小累积卷烟允许值为0，即逐个检验没有发现卷烟，所预测的产品总体卷烟比率的单尾置信上限
Figure FDA0002565031410000021
为5％，此时Za＝1.645，根据数学模型反算得到此时的样本量为51.4053，向上取整数，则最小抽样量nmin＝52；
此时根据
Figure FDA0002565031410000022
求得两个实根分别为0和4.95％，即当抽取n＝52只拟烟产品，经DNA分子标记检测，若未检出卷烟烟支，可以有95％的把握确定该批产品总体卷烟烟支比率在0～4.95％之间。
3.如权利要求2所述的基于DNA分子标记的拟烟品检验抽样方法，其特征是，步骤S7的具体过程包括如下：
在实际鉴别检验过程中很可能会出现多次检验均未达到判定阈值的情况，即累计的卷烟支数总是介于序贯抽样表所列的上下限之间而难以确定其总体专卖品比例，使检验不能终止；因此，需要计算最大抽样数nt来终止抽样过程；
根据设定的α＝0.05，β＝0.10，在GB/T8051-2008《计数序惯抽样检验方案》附表1中查得截尾拒收数为25，即最大累积卷烟支数拒收值为25，此时所预测的产品总体卷烟比率的双尾置信下限
Figure FDA0002565031410000023
为5％，由于Zβ＝1.654，根据S2-S5数学模型反算得到此时的样本量为363.3349，向下取整得到最大抽样数nt＝363.3349，向下取整数得到最大抽样数nt＝363。
4.如权利要求1所述的基于DNA分子标记的拟烟品检验抽样方法，其特征是，还包括S8，步骤S8的具体过程包括如下：
将抽取的若干烟支同时切断，并使其内烟丝同步落入对应取样瓶内，将取样瓶内烟丝与烟草基因组数据做烟草特异性DNA分子标记的鉴别检验。
5.应用权利要求4所述的基于DNA分子标记的拟烟品检验抽样方法的设备，其特征是，包括底座体、操作台、取样室及取样管；其中，
所述底座体与所述操作台之间留设有放置取样管的空间，以形成取样室；
所述操作台的顶端面开设有至少一个烟支容纳槽；
所述烟支容纳槽的外沿均开设有一个切断容置槽；
所述烟支容纳槽内在对应所述切断容置槽的位置开设有落烟口，所述操作台的底端对应所述落烟口的外沿设有连接凸起；
所述连接凸起与所述取样管的管口相固接；
所述操作台的一侧端转动设有与所述切断容置槽相对应的切刀。
6.如权利要求5所述的设备，其特征是，所述底座体固接设有向上延伸的支撑柱，所述支撑柱在远离所述底座体的一端与操作台相固接。
7.如权利要求5所述的设备，其特征是，所述烟支容纳槽呈长条弧状槽且具有一个延伸方向；
所述烟支容纳槽沿其延伸方向的一端与操作台的边沿对应设有槽口，所述烟支容纳槽沿其延伸方向的另一端位于操作台且对应设有限位端。
8.如权利要求7所述的设备，其特征是，所述限位端与所述切断容置槽之间的长度不小于3cm，切刀能够完全切除烟支滤嘴。
9.如权利要求5所述的设备，其特征是，所述切刀的侧端均匀固接设有至少一个侧刃，所述侧刃与所述烟支容纳槽相对应。
10.如权利要求5所述的设备，其特征是，所述连接凸起与所述取样管之间的固接方式为螺合或过盈配合。
CN202010620585.2A 2020-07-01 2020-07-01 基于dna分子标记的拟烟品检验抽样方法及其应用设备 Pending CN111984922A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010620585.2A CN111984922A (zh) 2020-07-01 2020-07-01 基于dna分子标记的拟烟品检验抽样方法及其应用设备

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010620585.2A CN111984922A (zh) 2020-07-01 2020-07-01 基于dna分子标记的拟烟品检验抽样方法及其应用设备

Publications (1)

Publication Number Publication Date
CN111984922A true CN111984922A (zh) 2020-11-24

Family

ID=73437703

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010620585.2A Pending CN111984922A (zh) 2020-07-01 2020-07-01 基于dna分子标记的拟烟品检验抽样方法及其应用设备

Country Status (1)

Country Link
CN (1) CN111984922A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114897349A (zh) * 2022-05-09 2022-08-12 中国人民解放军海军工程大学 一种成败型序贯抽样试验方案确定系统及方法

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114897349A (zh) * 2022-05-09 2022-08-12 中国人民解放军海军工程大学 一种成败型序贯抽样试验方案确定系统及方法
CN114897349B (zh) * 2022-05-09 2023-09-05 中国人民解放军海军工程大学 一种成败型序贯抽样试验方案确定系统及方法

Similar Documents

Publication Publication Date Title
CN104122179B (zh) 一种评价电子烟烟雾量的方法
CN106018168B (zh) 一种在线测量卷烟吸阻的方法
CN103844344A (zh) 一种调控不同批次卷烟烟丝质量均匀性的方法及其应用
CN110196265B (zh) 基于分布概率评价批内卷烟烟支内烟丝分布稳定性的方法
Alpert et al. Free nicotine content and strategic marketing of moist snuff tobacco products in the United States: 2000–2006
CN108241003B (zh) 一种非恒温式加热电子烟具一致性检测方法及系统
CN111984922A (zh) 基于dna分子标记的拟烟品检验抽样方法及其应用设备
CN105445047B (zh) 除湿机除湿量的检测方法及装置
CN109892699A (zh) 电子烟抽吸次数统计方法、电子烟、存储介质及设备
Di Miceli et al. Instrumental and sensory evaluation of eating quality of peaches and nectarines
CN101992862A (zh) 卷烟包装卷制质量缺陷分类判定方法
CN106324181B (zh) 一种用于卷烟燃吸过程中瞬时燃烧速率的测定方法
CN212647465U (zh) 基于dna分子标记的拟烟品检验的取样设备
Wadkin et al. E‐cigarette puffing topography: The importance of assessing user behaviour to inform emissions testing
CN117172723B (zh) 一种基于人工智能的卷烟生产加工智能管理系统
CN208030250U (zh) 用于烟叶物理特性检测及切丝的装置
CN106918525A (zh) 一种烟梗体积变化率检测方法
CN109100264A (zh) 一种快速预测细支卷烟抽吸均匀性的方法
CN106645530B (zh) 一种基于烟叶致香成分的多模型评价烟叶原料相似度的方法
CN111667176A (zh) 一种基于卷烟品牌需求的烟叶风格特征彰显程度评价方法
US20230160845A1 (en) Detection of contamination of fluids
CN114062537B (zh) 一种准确检测加热卷烟中施加料液扩散速率的方法
CN110162732A (zh) 一种设计卷烟单支重量的方法
CN110333177B (zh) 一种基于卷烟纸透气度稳定的卷烟分段吸阻无损快速检测评价方法
CN114877844B (zh) 卷烟抽吸过程稳定性的检测方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination