773_zh-TW
尋求有關使用 JMP 計算 LD50/ED50 的建議JMP 用戶您好,
我試圖弄清楚這一點,但統計指南非常模糊,我無法完全找到完成這項工作所需的答案。
我正在嘗試讓 JMP 根據觀察到的死亡率比例計算 LD50 值。
JMP 中的文件和例程似乎依賴邏輯回歸,並暗示 0.5 的逆預測將產生 LD50。 問題是,這種方法依賴於每個[log] 劑量水平下群體中每個個體的二項式(活/死)數據......我沒有個體動物數據,只有每個動物挑戰數量的死亡數量劑量。
因此,例如,我對藥物進行對數稀釋系列,並對十隻動物中的每隻施用每個劑量等級的藥物。 第一級,沒有動物存活,2->無,3->75%死亡,4->35%,5,6,7,所有動物都存活。
要使用邏輯回歸平台,可以為每隻動物的一行輸入 1 或 0(或“活”/“死”),對每個劑量水平重複此操作(表中的 70 行),這樣就可以正常工作,但我不知道哪個動物是哪個。
我正在尋找一種輸入 1, 1, 0.75, 0.35, 0, 0, 0 作為對數稀釋水平 [負] 1,2,3,4,5,6,7 的響應的方法。
有人知道如何做到這一點嗎?回覆:尋求有關使用 JMP 計算 LD50/ED50 的建議您需要知道每個級別的實際死亡或生存人數,而不僅僅是比例。 然後有一個二進制結果變量,包含 1 和 0 或“死”和“活”,並添加一個帶有計數的額外列。 將計數列放置在邏輯對話方塊的 FREQ 欄位中。 然後逆向預測將為您提供帶有置信區間的 LD50(或您需要的任何 LD)。
許多生物測定使用機率而不是 Logit。 JMP可以使用擬合最小平方法下的廣義線性模型進行Probit,只需選擇二項式分佈和Probit Link。 平台沒有反向預測選項(也許是下一個版本?)但您可以使用帶有置信帶的預測刻畫器來讀取反向預測和近似置信限。 只需滑動垂直紅色虛線,直到垂直線達到 0.5(或您要尋找的任何 LD)。 然後您可以從圖表中讀出對數稀釋度和信賴區間。 此 CI 可能與另一個套件產生的基準限制不同,而且我不確定它是否與 SAS Proc Probit 為您提供的 CI 相符。 如果您有 SAS,您可以向 SAS 提交報表,看看會得到什麼。
玩得開心
甘特
Comments
There are no comments yet