警戒限和行动限评估报告——模板.docxVIP

PAGE6 / NUMPAGES6 警戒限和行动限评估报告 纯化水系统 职责 姓名/部门 签字 起草人 审核人 批准人 目录 TOC \o 1-1 \h \z \u 1 背景 3 2 范围 3 3 评估原理 3 4 统计原理 3 5 命令调用 4 6 数据分析 4 7 结论 5 8 原始数据 5  背景 为了制定一个合适的警戒限和行动限，本报告结合风险评估及统计学来制定一个科学合理的行动限和警戒限。 数据是通过日常环境监测获得，涉及偏差的数据未包括在内，以避免已查明原因的特殊原因造成的变差影响计算结果，降低警戒灵敏度。 原始数据包括总需氧菌（TAMC），电导（Conductivity），总有机碳（TOC）。 范围 数据评价周期为yyyy.mm.dd~yyyy.mm.dd 纯化水系统 评估的项目包括总需氧菌，电导，总有机碳。 参数选择 显著性水平0.05 95% quantile would be used as alert limit and 99% quantile would be used as action limit. 控制策略 Baseline（可选） 微生物，警戒限的基线，行动限的基线 电导 TOC Update rule 微生物，计算值和当前限度相差9倍以上，需要更新限度。 电导/TOC，计算值和当前限度相关10%以上，需要更新限度。 控制策略 警戒限和行动限可以由两部分组成。基础限度为基于风险制定的固有限度。 如果计算的限度低于基础限度，我们可以使用因有限度作为警戒限或行动限。 如果计算的限度高于基础限度，但是低于法规限度，我们可以使用计算的限度作为行动限和警戒限。 如果计算的限度高于法规限，须采取措施进行整改。并使用基础限度作为警戒限和行动限。 新厂房应使用基础限度作为警戒限和行动限。 根据数据分布特点，把95%分位数处的值作为警戒限，把99.73%分位数处的值作为行动限。 本报告中使用0.05作为显著性水平。 基础限度请参见表1 表1: 基础限度 项目 限度 基础限度 总需氧菌 警戒限 1000cfu/100ml 电导 警戒限 2.1μs/cm 总有机碳 警戒限 50ppb 统计原理 使用JMP 版本 10.0.0 ，Minitab 19， R 控制台 REF _Re\r \h 5 的函数进行计算。 在限度计算之前需要对数据进行正态性检验。如果数据服从正态分布，我们可以使用 REF _Ref429422637 \r \h 4.3所示的正态方法计算警戒限或行动限。如果不服从正态分布，我们可以使用负二项分布或者伽马分布，但是如果数据的均值大于方差，即使不服从正态分布，我们仍可以使用切比雪夫定律来进行计算，可以使用均值标准差的方法计算警戒限和行动限，或使用其它方法。 正态分布又称高斯分布，均值加减1个标准差范围内的概率密度为68.27%，均值加减2个标准差范围内的概率密度为95.45%，均值加减3个标准差范围内的概率密度为99.73%。分布表达式为。 X~N(u, δ) 负二项分布可以用来计算限度，来避免低估限度。本函数有3个参数。X代表一系列伯努利实验后的第x次成功实验出现。N代表伯努利实验的总次数。P代表每次实验的成功概率。公式如下： pnbinom(x, n, p) 因为负二项分布为非闭合的函数，计算过程需要进行迭代计算。泊松-伽马分布为一个闭合的函数以计算限度。P为概率，西格玛是过度分散参数，lambda为均值参数。函数如下： Gamma Poisson Quantile (lambda, sigma, p) eq \o\ac(○,1) 伽马分布同样可用于计算限度，式中，p为事件出现的概率。alpha为形状参数，beta为尺度参数。函数公式如下： Gamma Quantile (p, alpha, beta) 命令调用 R控制台 方法一 输入 source(“limit.R”) LC(mean, variance, t, way, AL, AC) 这里mean为均值，variance为方差，t为正态性检验的P值，way为方法:1为伽玛分布，其它值为负二项分布，默认为1，AL为警戒限的累积概率。AC为行动限的累积概率。 注意：如果t大于等于0.05，将会使用正态方法进行计算，无论方法参数m为任何值。 方法二 输入 source(“limit1.R”) 复制数据，第一列为日期，第二列为数据 调用LC(way, AL, AC, print)，way为方法:1为伽玛分布，其它值为负二项分布，默认为1，AL为警戒限的累积概率。AC为行动限的累积概率。Print为是否输出结果到excel，默认为FLASE。 JMP 使用XXX模型拟合数据