劑量反應概率分析<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 值。 分析會將它們聯繫在一起。

你沒有提到的一件事是經常為這種類型的工作做的是記錄濃度。 我在上面顯示的模型對話框中這樣做了。 突出濃度模型效果,單擊轉換按鈕並選擇日誌。 結果張貼在這裡。

擔