來自希望在運行效力測試之前進行 logit 函數轉換的客戶的生物測定數據大家好，
我有一位客戶試圖模仿下面所附的使用 Combistats 軟體製作的報告，對劑量反應數據（總數中的死亡人數）執行效力、並行性分析等。請參閱所附資料表。
這裡的問題是，在這份報告中，Y 是用 logit 函數轉換的。我們如何在 JMP 中執行此操作？ JMP 能否與此報告中的結果相符？
客戶熟悉用於劑量反應擬合的擬合曲線中的邏輯回歸。這裡的問題是需要模仿報表並使用 logit 函數執行轉換。
任何見解將不勝感激。
回覆：來自希望在運行效力測試之前進行 logit 函數轉換的客戶的生物測定數據沒有完整的答案，但它一定是某種調整後的邏輯，然後是標準變異數分析和效力計算。比例上的樸素 logit 不起作用，因為 logit(0) 是無限的。在https://www.unistat.com/guide/bioassay-analysis-quantal-response-method/的第 10.3.3 節中找到了相同的數據和其他線索
回覆：來自希望在運行效力測試之前進行 logit 函數轉換的客戶的生物測定數據右鍵單擊比例並選擇新列Formula>specialty>logit 將給出所需的轉換，但如Russ 所說，0 比例將給出-inf 值（1 將給出+inf 值），並將導致缺失值。在這種情況下，我看到人們使用廣義線性模型（JMP 可以做到這一點），其中響應和總數在 y 中，劑量和樣本在 x 中。然後用 logit 連結進行二項式分佈。這裡的問題是，沒有方便的方法來測試並行性、效力等，您必須手動計算所有這些。
回覆：來自希望在運行效力測試之前進行 logit 函數轉換的客戶的生物測定數據基於CombiStats 網頁，他們使用歐洲藥典作為分析基礎（分位數回應法）。我在 2019 年 Discovery EU 上就此做了一次演講，但主要關注相對效力的計算。嘗試在 JMP 中重現這一點存在一些挑戰。首先，您需要使用廣義線性模型或廣義迴歸來進行 Logit 擬合。您的輸出將包括卡方統計數據，但不包括平方和和 F 檢定。當我嘗試使用演講中的過程來執行此操作時，我無法匹配結果。我不確定他們所說的非線性是否與我使用的相同。他們使用 MLE 還是 LS？這有什麼不同嗎？用於收集資料的設計也會產生影響，我不確定我所說的話是否與範例中的設計一致。
這不是我第一次看到 EP 5.3 的出現。我有興趣收到@SamGardner的來信
回覆：來自希望在運行效力測試之前進行 logit 函數轉換的客戶的生物測定數據事實證明這是一個有趣的小問題@yasmine_hajar感謝@DonMcCormack和@Chris_Kirchberg 的評論。雖然我們處理這個問題的方式要好得多並且符合最佳現代統計實踐，但這不是OP。我記得多年來嘗試從其他軟體包重現結果的一些困難和令人沮喪的案例。
CombiStats 文件中的變異數分析表是老式的，帶有經典的平方和等，但根本不清楚他們是如何計算它的。該文件似乎沒有提供足夠的詳細資訊來說明他們的具體方法。我只是快速嘗試了幾件事，包括使用 0.5 和其他一些值的連續性校正的 logits、各種權重，甚至 GLMM 插件https://community.jmp.com/t5/JMP-Add-Ins/Generalized- Linear-Mixed-Model-Add-in/ta-p/284627 @MeichenDong在@ElizabethClaassen 在擬合模型（計算偽邏輯和權重）中實現它之前寫過，但沒有任何東西與他們的方差分析相匹配。
如果有人對數字/統計難題感興趣，我相信可以根據方差分析表的結果對他們用作 Y 的確切 8 個數字進行逆向工程，但我還沒有嘗試進行設定。它將是一個非線性方程組，具有平方和的經典變異數分析公式和 8 個 Y 作為未知數。他們也有可能使用權重－文件提到使用逆方差，對於二項式來說，它是 1/(n*p*(1-p))。
@yasmine_hajar想知道客戶是否願意聯繫 CombiStats 並詢問有關如何計算示例 5.3.2（Logits）的詳細信息，表明他們希望獨立重現結果。一旦我們看到，我確信我們可以以某種方式重現它。
