用JMP轻松计算半数致死量LD50

shenweis

案例：某研究者要测定太湖水中微囊藻毒素的毒力。分为不同的课题组，观察给药1h后各组小鼠的死亡情况，实验结果见表。求微囊藻毒素的LD50.
数据：
恢复数据
专家给出的结果：LD50为204.6mg/kg，其95%可信区间为183.1—228.8mg/kg.
注：此为专家手工计算。

shenweis

JMP的数据格式，还有其它，这里不介绍。

这里0代表死亡，1代表存活。

shenweis

变换剂量的数据，使其有原始数据，自然对数数据，常用对数数据。

shenweis

利用logistic回归计算，三种数据的分析结果汇总：
1、原始数据整体模型检验：


2、失拟检验：


3、逆预测设置死亡概率0.5：


4、逆预测结果：


5、自然对数数据结果：


6、常用对数数据结果：

shenweis

利用广义线性模型计算：


结果：


剂量数据进行自然对数或常用对数变换后，同理可以计算，这里不再赘述。

shenweis

本案例用logistic回归、广义线性模型两种方法，对剂量、剂量的自然对数、剂量的常用对数三种格式计算LD50，结果数值上有些稍有差异，但都是正确的计算方法。

用回归模型进行半数致死量的估计方法，由于对课题的设定没有特别要求，模型参数估计不受死亡率为0%或100%极端点的影响，而且可以通过考察模型的拟合度而选择更适合的模型，因而具有更广泛的适用性。

shenweis

JMP可以在逆预测中同时给出多种死亡概率的结果，可谓十分方便。

niksong

按照楼主的这种学习速度,用不了多久就要成为JMP专家了.

shenweis

兄弟才是专业，发现我的分析有问题就一定要指出来哦！

jappe2009

学习中，支持你，一直关注你，学习！