複雜的概率估算變得簡單 JMP 由於其精心設計的用戶界面,提供了許多易於使用的複雜統計分析。非線性估計並不那麼容易使用,因為該平台的難以置信的靈活性可能會挑戰一個統計上精明的人。但是,使用JMP腳本語言(JSL),可以輕鬆創建適合應用程序區域特定需求的簡單應用程序,並且適合快速重複使用。
估算非線性函數的參數首先要在數據表中創建包含公式的列,以指定作為回歸函數的一部分進行估計的參數,以及應當最小化以找到最終回歸估計的損失函數(如果不應該由最小二乘法完成)。
Probit回歸是二元結果(響應/無響應)的劑量反應實驗的常用分析。它的應用在JMP幫助和書籍中有詳盡的描述,但其用法至少需要基本的數學理解;此外,使用起來也很繁瑣,尤其是經常需要時。這裡的解決方案是編寫JSL腳本。您需要做的就是使用典型的命名法將數據集中的正確變量分配給回歸所需的角色,讓JMP完成剩下的工作。
數據集需要包含劑量列,待測物品總數和響應物品數量,如您所見圖1:數據。

圖1:數據
JSL腳本也很簡單。它由三部分組成:用戶對話框( 圖2:對話框 ),公式列的插入和非線性平台的應用。

圖2:對話框
您無需從頭開始編寫評估語句,因為您可以交互式地進行分析並將腳本保存到腳本窗口中。然後所有語句都存在,您只需要替換程序中的變量名稱。腳本解決方案的另一個優點是您可以覆蓋特定請求,例如各種逆預測的計算。它們是自動生成的,用戶無需單獨從上下文菜單中調用它們。

圖3:結果
將包裝為加載項時,可以將整個腳本作為一鍵式操作進行分發和安裝。
您可以下載這個JSL應用程序來自JMP文件交換。
Arati Mejdal寫道:嗨,那裡。是的,它可用。貝科... Arati Mejdal 中寫道:
嗨,您好。是的,它可用。因為我們移動了文件交換,所以文件的位置發生了變化。這裡是: https://community.jmp.com/docs/DOC-6974
抱歉給你帶來不便。
請幫助寫道:它仍然可以下載嗎?一世...請幫忙中寫道:
它還可以下載嗎?我似乎無法找到它。謝謝。
Bernd Heinen寫道:大衛你好,我為因為...而道歉...貝恩德海寧中寫道:
大衛你好,
我為有點不專心而道歉。現在它再次可用,關鍵字搜索“probit”將帶您到那裡。
最好的祝福
貝恩德
大衛寫道:嗨貝恩德,這個JSL腳本仍然可用......大衛中寫道:
嗨貝恩德,
此JSL腳本是否仍可從JMP文件交換中獲得?我找不到它搜索“probit”或“回歸”
乾杯
Bernd Heinen寫道:你好Joe,我無法複製......貝恩德海寧中寫道:
你好喬,
我無法重現該錯誤,請將您的JMP數據表發送給我bernd.heinen@jmp.com
喬寫道:我和史蒂夫有同樣的問題。我有...喬中寫道:
我和史蒂夫有同樣的問題。我有3列;一個用於劑量,一個用於總生物,一個用於總事件。當我使用probit附加組件時,我收到的相同錯誤消息是“無法為列'Dose'設置值,因為行號(-1)無效。”
Bernd寫道:你好Steve,probit計算本身......貝恩德中寫道:
史蒂夫你好
概率單位計算本身有點繁瑣,但各個實驗的結果可以用非常簡單的結構報告。它可能看起來像這個表:
劑量|被測單位數量|觀察到的事件數量
10,000 | 10 | 10
1,000 | 10 | 10
100 | 10 | 9
10 | 10 | 6
1 | 10 | 3
0.1 | 10 | 0
這些列對應於腳本對話框中的字段(日誌)劑量,#總計和#響應。無論您的數據最初採用何種形式,JMP都擁有在其表格中重新整形的所有工具。
史蒂夫寫道:您好,我需要幫助估算...史蒂夫中寫道:
嗨,我需要幫助與probit估計JSL,所以有些原因我無法想像每次我嘗試運行我得到一個錯誤消息“無法為列'Dose'設置值,因為行號(-1)無效“
有誰知道我做錯了什麼?請幫忙急需幫助我的論文做一些劑量反應曲線
Bernd Heinen寫道:你好Tisha,首先你的觀察...貝恩德海寧中寫道:
你好Tisha,
首先,您的觀察是正確的,與LDParm相關的值是回歸線的斜率。
另外一個提示:在解決方案大綱之上有一個置信限制的關鍵。您可以指定置信度,然後獲得兩個參數的置信度限制。每個置信度都需要單獨的迭代。結果顯示在參數表的其他列中。它們也可能對您的比較有所幫助。
遺憾的是,協方差矩陣未在報告中打印出來。得到它的唯一方法是去熱點(非線性擬合和保存公式前面的紅色三角形 - >預測的節省標準誤差。這會使用公式將新列插入到數據表中。您可以打開此公式,其中的第一部分是協方差矩陣。您可以復制並粘貼它,因此可以將其用於計算。
我希望,我能幫助解釋一下。
最好的祝福
貝恩德
Tisha King-Heiden寫道:我已成功使用此sc ... Tisha King-Heiden 中寫道:
我已成功使用此腳本計算JMP11中的LC50,並希望使用LC50比率來比較不同處理的不同LC50。要做到這一點,我需要協方差。我確定LDParm是斜率,我想知道協方差是否也有不同的標記。如果沒有,我如何確定協方差?
Gail Massari寫道:Siamak,您可以提交來自...的SAS代碼蓋爾馬薩里中寫道:
Siamak,
您可以從JMP提交SAS代碼。如果您正在尋找JMP中的等效功能,請考慮參與JMP軟件社區論壇以獲取特定問題的問題和答案。JMP用戶發現這是從其他用戶那裡獲得幫助的最快捷方式。JMP開發人員和技術支持人員也參與其中。看到 www.jmp.com/forum 。
凱倫克萊德寫道:我想用這個超過200分開......凱倫克萊德中寫道:
我想將它用於> 200個獨立的數據組。沒有“by”選項,似乎沒有一種簡單的方法可以做到這一點(與使用擬合模型/廣義線性模型啟動分析不同)。有任何想法嗎?
Sajeewa Amaradasa寫道:嗨Bernd;你能解釋一下... Sajeewa Amaradasa 中寫道:
嗨貝恩德;
你能解釋一下如何使用JMP來計算殺菌劑的EC50(50%生長抑制的有效濃度)嗎?我想為不同的真菌分離物找到EC50。對於每種分離物的不同殺真菌劑劑量,我對4次重複具有不同的%抑制。可以使用上面公佈的相同應用進行分析。
謝謝
Sajeewa
Bernd Heinen寫道:嗨Siamak,你回應我的...貝恩德海寧中寫道:
嗨Siamak,
你回答我的博客是關於JMP中的概率分析,並且沒有Proc Probit,而是一個非線性平台。如果運行此選項,結果報告將始終具有紅色三角形(熱點),您可以從中選擇“估計” - >“反向預測”。它會打開一個表格,您可以在其中輸入要計算的多個ED值。
另一方面,如果您在那裡使用SAS Stat和Proc Probit,則需要使用以下語法輸入INVERSECL選項:INVERSECL <(PROB = rates)>。
我希望我能提供幫助
貝恩德
Siamak Roshandel寫道:嗨,我想知道我可以使用proc ... Siamak Roshandel 中寫道:
你好
我想知道我可以使用proc probit來計算LT50嗎?如果有可能在哪個部分我可以找到它