西西河

主题:热烈祝贺中国农业部再次荣获“闭眼掩耳遗忘装傻盗铃”奖 -- 思想的行者

共:💬29 🌺37 🌵3 新:
全看分页树展 · 主题 跟帖
家园 您对统计(生统)不太熟悉吧

文中引用那篇文章《转cry1Ab/Ac基因抗虫水稻华恢1号在湖北省生产应用的安全证书》我看了,乍一看34页的内容很吓人。一共分了5组,control一组使用基础饲料,亲本稻谷饲料一组,高中低剂量的转基因稻谷饲料各一组。在长达106周试验期间,动物有死亡,各剂量组共死亡70只...总计达39%,所以剂量组共有70/0.39=180只sd大鼠,每组60只。

这篇论文没有给出control group和reference group的死亡比例确实不太严谨,这么写在国外统计,生统,生物期刊上根本发表不了。但是您认为没给出就是比例悬殊,猜测 70:9,甚至说“既然他一字不露只能认为对照组一只都没有死”,这个太想当然了。

那么我们来随便搜索一下sd大鼠.百度百科:1925年,美国斯泼累格·多雷(Sprague Dawley)农场用Wistar大鼠培育而成。其毛色白化。广泛用于药理、毒理、药效及GLP实验。寿命有人说是寿命为2.5-3年。我搜索了一篇在nih上被索引了38次的论文 The Laboratory Rat: Relating Its Age With Human's,其中索引目录里提到了Harlan, Inc. Sprague-Dawley growth chart.一个专门养实验动物公司的介绍 Life-span and Spontaneous Disease:

The median survival time was 24.5 months for virgins and 25 months for breeders (Cameron et al, 1982).

说明sd大鼠从可以喂食开始的寿命median是24.5-25个月,所以原文中说“按毒理程序规定,要求慢性毒性试验中24个月前动物死亡率不能大于50%”这句话没什么问题。单只sd大鼠活过106周的概率50%。

您自己猜测的70:0的概率是多少呢,control组,亲本稻谷饲料组共120只sd大鼠,0死亡的概率是0.5的120次方,大概是7.5x10的-37次方。如果只死了9只或以下的概率呢,用正态近似,总死亡数Y服从N(60,30)的正态分布,30是方差。所以P(Y<=9) =inverse Normal((9-60)/sqrt(30))=5.94 x 10的-20次方。

这篇论文有不严谨的地方,但是以我对统计粗浅的了解问题不大。

反转可以,但最好能用科学的方法去反,小崔前期出了不少笑话不就是这个原因吗?而且切忌一知半解盲目地得出结论。

我对转基因的看法也是认为还需要更多研究,目前的研究还看不出有毒的证据,不宜过早投入市场。

全看分页树展 · 主题 跟帖


有趣有益,互惠互利;开阔视野,博采众长。
虚拟的网络,真实的人。天南地北客,相逢皆朋友

Copyright © cchere 西西河