剂量反应概率分析<meta http-equiv="Content-Type" content="text/html; charset=UTF-8" />嗨组,
我正在尝试重现我在 SPSS 中执行的简单剂量反应分析。 所以在 JMP 中,我使用“拟合模型”模块,选择广义线性模型、二项分布和概率链接。
我有三个变量,进行的化验次数(10 次化验),阳性结果的化验次数和浓度(从 0 到 36)
0 10 2
1 10 4
2 10 9
6 10 16
8 10 25
9 10 36
0 10 0
我选择执行的化验数量作为“Y”变量,结果为阳性的化验数量作为“Freq”,日志中的浓度作为变量来模拟影响。 但是当我点击运行时,我会在结果中收到这条消息:
“警告:Hessian 不是正定的。 奇异性可能是由于模型拟合不佳或模型协变量之间的线性相关性所致。 拟合和结果的价值值得怀疑”
并且不执行分析。
任何人都可以分享一些评论吗?
谢谢,
戴夫
回复:剂量反应概率分析这确实是输入变量的正确顺序。 谢谢你指出这一点,迈克尔。
戴夫
回复:剂量反应概率分析出于某种原因,原始数据的最后一行 0 10 0 使分析不需要正确的顺序。 如果您将其排除在外,那么顺序就很重要。
迈克尔
回复:剂量反应概率分析一直感兴趣地关注这个线程。
来自 JMP 文档:
如果您的数据汇总为事件和试验的频率,请指定两个连续的
列按此顺序:成功次数的计数,以及成功次数的计数
试验。或者,您可以指定失败次数而不是成功次数。
(我的重点)。
因此,如果您切换输入 2 个 y 变量的顺序,我认为您会看到更多您期望的内容,尤其是如果您将 x 轴更改为对数。
也就是说,我无法解释为什么订单似乎与您的原始数据无关。
迈克尔
回复:剂量反应概率分析亲爱的丹诺,
抱歉再次打扰您,但是将数字更改为演示,我似乎得到了错误的结果。 请查看此数据集,变量相同但数字不同。 对于三个最高浓度,我得到 100% 阳性
阳性结果数量 进行的检测数量 浓度
0 1140 0.01
241 1140 0.1
590 1140 0.5
963 1140 1
1134 1140 5
1140 1140 10
1140 1140 50
1140 1140 100
我按照你的建议运行 Probit

运行后我得到一个奇怪的结果,这对我来说没有意义。 它不是在 Y 上做正确的百分比,而是在做一些不同的事情。 我究竟做错了什么?

非常感谢你的帮助
戴夫
回复:剂量反应概率分析非常感谢 DanO,这是一个非常明确的答案。 有趣的是,JMP 知道它们是相关联的。
关于日志转换,是的,我会考虑到这一点(它确实在我原来的帖子中提到过,尽管我必须承认非常模糊!)
最好的,
戴夫
回复:剂量反应概率分析戴夫,
您的“阳性结果的检测数量”也需要作为“Y”变量输入。 Frequency 作用是表示数据表中一行实际出现了多少次。 在这种情况下,GLM 模型期望执行的检测次数和阳性结果的数量作为“单独的”Y 值。 分析会将它们联系在一起。

你没有提到的一件事是经常为这种类型的工作做的是记录浓度。 我在上面显示的模型对话框中这样做了。 突出浓度模型效果,单击转换按钮并选择日志。 结果张贴在这里。

担