JMP数据分析论坛

标题: 用JMP轻松计算半数致死量LD50 [打印本页]

作者: shenweis 时间: 2014-4-23 14:56
标题: 用JMP轻松计算半数致死量LD50
本帖最后由 shenweis 于 2014-4-23 16:00 编辑

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

作者: shenweis 时间: 2014-4-23 14:58
本帖最后由 shenweis 于 2014-4-23 15:55 编辑

JMP的数据格式，还有其它，这里不介绍。
[attach]1854[/attach]
这里0代表死亡，1代表存活。

作者: shenweis 时间: 2014-4-23 14:59
本帖最后由 shenweis 于 2014-4-23 15:07 编辑

变换剂量的数据，使其有原始数据，自然对数数据，常用对数数据。
[attach]1856[/attach]

作者: shenweis 时间: 2014-4-23 15:01
本帖最后由 shenweis 于 2014-4-23 15:39 编辑

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

2、失拟检验：
[attach]1859[/attach]

3、逆预测设置死亡概率0.5：
[attach]1861[/attach]

4、逆预测结果：
[attach]1857[/attach]

5、自然对数数据结果：
[attach]1862[/attach]

6、常用对数数据结果：
[attach]1865[/attach]

作者: shenweis 时间: 2014-4-23 15:08
本帖最后由 shenweis 于 2014-4-23 15:43 编辑

利用广义线性模型计算：
[attach]1866[/attach]

结果：
[attach]1867[/attach]

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

作者: shenweis 时间: 2014-4-23 15:11
本帖最后由 shenweis 于 2014-4-23 15:50 编辑

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

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

作者: shenweis 时间: 2014-4-23 15:53
JMP可以在逆预测中同时给出多种死亡概率的结果，可谓十分方便。

作者: niksong 时间: 2014-4-25 21:26
按照楼主的这种学习速度,用不了多久就要成为JMP专家了.

作者: shenweis 时间: 2014-4-25 21:32
兄弟才是专业，发现我的分析有问题就一定要指出来哦！

作者: jappe2009 时间: 2014-6-18 15:33
学习中，支持你，一直关注你，学习！

作者: juna 时间: 2014-6-23 13:26
喜欢这种带具体例子的讲解，多谢楼主分享

欢迎光临 JMP数据分析论坛 (http://www.jmpforum.net/)