蒲公英 - 制药技术的传播者 GMP理论的实践者

搜索
查看: 4783|回复: 23
收起左侧

[其他] 悬浮粒子数分布的矛盾

  [复制链接]
药士
发表于 2013-5-2 09:55:50 | 显示全部楼层 |阅读模式

欢迎您注册蒲公英

您需要 登录 才可以下载或查看,没有帐号?立即注册

x
    由于有部分图片,且帖子排版太麻烦,所以直接上附件了,欢迎大家看了附件后,发表自己的观点。

悬浮粒子数分布的问题.rar

2.33 MB, 下载次数: 144

回复

使用道具 举报

大师
发表于 2013-5-2 10:01:48 | 显示全部楼层
关于悬浮粒子数分布的一个矛盾
无名
前几天看到论坛有人讨论,悬浮粒子在数据回顾的时候用什么图,进而讨论到,悬浮粒子数是属于什么分布的问题,以前我发过一篇帖子认为悬浮粒子属于正态分布,并推导得出UCL值不能用来进行回顾分析。但经过我的一些思考后,却发现了一个矛盾,特地发上来,大家一起讨论一下。
先罗列第一种看法:按控制图选择图来看,悬浮粒子属于单位体积内的数量属于计点值,应该属于泊松分布。
我们再来看另一种看法:在GB/T16292-2010或者ISO14644-1中提到了95%置信限的计算,其计算公式如下图:
1.jpg

注意上图中的t分布系数,下面我们会来做对比。
我们再去看看统计学中的单边置信区间的计算,在算置信区间的时候需要先知道样本所属总体的分布,其具体的选择依据及计算公式见下图:

2.jpg

我们可以发现该图中的“一个正态总体”的单侧置信区间其计算公式和95%UCL计算公式非常相似,那么他们前面的t分布系数是否相同呢?上图中的                              是t分布系数,详细见下图:(选自《概率论与数理统计》盛聚,谢式千,第三版)

3.jpg
我们可以发现该图中的“一个正态总体”的单侧置信区间其计算公式和95%UCL计算公式非常相似,那么他们前面的t分布系数是否相同呢?上图中的                              是t分布系数,详细见下图:(选自《概率论与数理统计》盛聚,谢式千,第三版)
回复

使用道具 举报

药徒
发表于 2013-5-2 10:11:10 | 显示全部楼层
回复

使用道具 举报

药生
发表于 2013-5-2 10:14:31 | 显示全部楼层
回复

使用道具 举报

药徒
发表于 2013-5-2 10:16:26 | 显示全部楼层
先下载学习了,然后再参与讨论
回复

使用道具 举报

药徒
发表于 2013-5-2 10:52:56 | 显示全部楼层
表示看不懂  概率没学过啊  
回复

使用道具 举报

药徒
发表于 2013-5-2 11:24:37 | 显示全部楼层
回复

使用道具 举报

药徒
发表于 2013-5-2 12:24:23 | 显示全部楼层
本帖最后由 blue8eyes 于 2013-5-2 12:27 编辑

下载了GB/T16292-2010查看相关的计算公式。发现如下其中有几个关键点需要注意。
第一:计算最后UCL数据的两个参数,M以及SE。
其中M以及SE的公式可以在上面大呆子的图片中看到。但是需要注意的是M以及SE公式中的“Ai”。
第二,Ai的来源依据
Ai的计算出自GB/T16292-2010中6.1小节,定义为采样点的平均悬浮粒子浓度,同时公式是Ai=ΣCi/n
其中Ci 为某一点采样点的粒子浓度。
也就是本身UCL计算的来源数据为平均值。
第三,UCL的意义
通过回顾UCL值的计算公式,发现其实GB/T16292-2010中所计算的UCL即控制上限实际上是各点平均值的控制上限。
我们平时所计算的控制上限=Xbar+3σ中是针对原始数据的,所以UCL计算的原始数据是“平均值”。
第四,中心极限定理
中心极限定理:设从均值为μ、方差为σ^2;(有限)的任意一个总体中抽取样本量为n的样本,当n充分大时,样本均值的抽样分布近似服从均值为μ、方差为σ^2/n的正态分布。
也就是说原本服从泊松分布的原始数据Ci,被计算出平均值Ai之后,Ai应服从正太分布。

我想可能是基于以上原因,UCL值的计算公式使用的是正态分布的公式。
以上为各位观点,欢迎各位指正及讨论,谢谢。
回复

使用道具 举报

药徒
发表于 2013-5-2 12:39:22 | 显示全部楼层
看了半天还糊涂
回复

使用道具 举报

药士
 楼主| 发表于 2013-5-2 13:14:33 | 显示全部楼层
我也再补充2点:
1.如果原始数据Ci符合泊松分布,则可以证明均值Ai也符合泊松分布;同样可证明若原始数据Ci符合正态分布,则均值Ai也符合正态分布。
2.中心极限定理应用的条件是n充分大,但貌似我们悬浮粒子的n不大。另:去年去听一个ISO里面专门搞悬浮粒子的专家曾今提到,以后再ISO的标准中可能会删除UCL的计算,理由1观n不大;理由2:由于悬浮粒子测试是由信号转化来计数,实际发现在转化过程中信号的的影响远大于外部环境的影响。

点评

我也可能当时建立这个指南就是基于中心极限定理吧,虽然n是不大。具体的理由和依据还需要询问制定该指南的相关组织了。  发表于 2013-5-2 13:22
回复

使用道具 举报

药士
发表于 2013-5-2 13:32:08 | 显示全部楼层
回复

使用道具 举报

药士
 楼主| 发表于 2013-5-2 13:35:40 | 显示全部楼层
人气要旺一点啊,欢迎大家踊跃发言啊!!
回复

使用道具 举报

药士
 楼主| 发表于 2013-5-3 10:42:32 | 显示全部楼层
blue8eyes 发表于 2013-5-2 12:24
下载了GB/T16292-2010查看相关的计算公式。发现如下其中有几个关键点需要注意。
第一:计算最后UCL数据的两 ...

现在看来,应该是悬浮粒子符合泊松分布,UCL值却是是用了中心极限定理。那在回顾分析时,有2个选择:
1.可以对每个点的平均值做C图,这个没啥争议。
2.对一个房间的平均值作图,这里就有争议了,按理泊松分布的均值还是泊松分布,该用C图,但依据UCL值计算   公式考虑了中心极限定理,又可以用单值I-MR图。这个究竟用那个呢?

点评

就像你之前说的,指南的编制人员也在考虑去除UCL值计算的要求,可能他们也意识到了相关问题吧。不过还是需要具体问问编写人员,不知道怎么联系他们。  发表于 2013-5-3 10:59
回复

使用道具 举报

药徒
发表于 2013-5-3 15:27:05 | 显示全部楼层
谢谢你辛苦了
回复

使用道具 举报

药士
发表于 2013-5-3 15:38:48 | 显示全部楼层
这个我一直在研究,目前没有水平喷
回复

使用道具 举报

发表于 2013-5-3 15:47:50 | 显示全部楼层
努力理解中!
回复

使用道具 举报

药士
 楼主| 发表于 2013-5-3 17:22:55 | 显示全部楼层
yuansoul 发表于 2013-5-3 15:38
这个我一直在研究,目前没有水平喷

欢迎开喷!
回复

使用道具 举报

药徒
发表于 2013-5-3 20:28:53 | 显示全部楼层
先学习再讨论
回复

使用道具 举报

药士
发表于 2013-5-4 06:06:03 | 显示全部楼层
glm1024 发表于 2013-5-3 17:22
欢迎开喷!

ok。我继续努力
回复

使用道具 举报

药士
 楼主| 发表于 2013-5-31 15:25:31 | 显示全部楼层
blue8eyes 发表于 2013-5-2 12:24
下载了GB/T16292-2010查看相关的计算公式。发现如下其中有几个关键点需要注意。
第一:计算最后UCL数据的两 ...

新的ISO14644-1已经出来了,已经取消了UCL值的计算,详细见https://www.ouryao.com/thread-136115-1-1.html
回复

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

×发帖声明
1、本站为技术交流论坛,发帖的内容具有互动属性。您在本站发布的内容:
①在无人回复的情况下,可以通过自助删帖功能随时删除(自助删帖功能关闭期间,可以联系管理员微信:8542508 处理。)
②在有人回复和讨论的情况下,主题帖和回复内容已构成一个不可分割的整体,您将不能直接删除该帖。
2、禁止发布任何涉政、涉黄赌毒及其他违反国家相关法律、法规、及本站版规的内容,详情请参阅《蒲公英论坛总版规》。
3、您在本站发表、转载的任何作品仅代表您个人观点,不代表本站观点。不要盗用有版权要求的作品,转贴请注明来源,否则文责自负。
4、请认真阅读上述条款,您发帖即代表接受上述条款。

QQ|手机版|蒲公英|ouryao|蒲公英 ( 京ICP备14042168号-1 )  增值电信业务经营许可证编号:京B2-20243455  互联网药品信息服务资格证书编号:(京)-非经营性-2024-0033

GMT+8, 2025-9-26 16:53

Powered by Discuz! X3.4

Copyright © 2001-2020, Tencent Cloud.

声明:蒲公英网站所涉及的原创文章、文字内容、视频图片及首发资料,版权归作者及蒲公英网站所有,转载要在显著位置标明来源“蒲公英”;禁止任何形式的商业用途。违反上述声明的,本站及作者将追究法律责任。
快速回复 返回顶部 返回列表