概率分析 - 擬合檢驗(回歸)和重複我正在嘗試確定我在體外(在培養皿中)嘗試過的某些殺菌劑與不同分離株的 ec50(與對照相比,有 50% 徑向真菌生長的濃度)。 我用不同的殺菌劑測試了相同的 9 種不同的分離株作為放射狀菌絲體的生長。 對於每種殺真菌劑,有 5、6 或 7(取決於殺真菌劑)濃度測試,每個濃度重複 3 次。 對於每種殺真菌劑,所有培養皿均在同一天測量。
概率 - 二項式分佈似乎很好地擬合了數據,至少在大多數情況下是這樣。
我的問題是如何在分析中包含重複以及如何確定標準誤差。
附上數據截圖。
datadata數據
fit modelfit model適合模特
regression plotregression plot回歸圖
ec 50ec 50EC 50
回复:概率分析 - 適合測試(回歸)和復制所以我像上圖一樣設置了分析(下面的快速腳本來設置它)
Fit Model(
Y( :inhibition, :N ),
By( :isolate ),
Effects( :log ppm, :Replicate, :Replicate * :log ppm ),
Personality( "Generalized Linear Model" ),
GLM Distribution( "Binomial" ),
Link Function( "Logit" ),
Overdispersion Tests and Intervals( 0 ),
Name( "Firth Bias-Adjusted Estimates" )(0),
Run
)
按住控制鍵(OSx 上的命令鍵),同時單擊最高級別的紅色三角形,從菜單中選擇反向預測。 然後對於彈出的每個對話框,輸入 0.5 作為預測概率。 這需要一段時間,因為系統會提示您輸入每個隔離級別。
完成後,右鍵單擊其中一個分離株的逆向預測表之一,然後選擇製作組合數據表。 該表是我在上面的示例中進行操作的表,因此我可以繪製圖表並進行分析。
回复:概率分析 - 適合測試(回歸)和復制我在擬合模型後使用紅色三角形下的逆預測區間選項來獲得拐點/ec50s
回复:概率分析 - 適合測試(回歸)和復制謝謝。 你能告訴我你是如何計算平均值的嗎? 您是否將它們導出到表格並從那裡準備圖表? 我似乎無法弄清楚。 先感謝您。
回复:概率分析 - 適合測試(回歸)和復制在這種情況下,一組重複是在一起的。 這不是一個壞方法,因為它基本上符合均值複製,從而更好地適應響應模式。 請記住,重複不是獨立的觀察值(取決於實驗單位),並且 JMP 假定每一行都是獨立的觀察值。 如果數據是堆疊的(每個重複一列,每個分離物一行)
您可能還想嘗試為每個分離株分別擬合每個複制。像這樣:
請注意,您必須為每個分離物手動添加每個逆向預測。 (注意:超級乏味,所以請在 jmp dot com 上發送功能建議以支持)當您選擇逆向預測選項時按住 CTRL,然後將出現用於計算每個分離物的 0.5 概率的對話框。
在這種情況下,將計算隔離對每個複制的適合度。 錯誤不會為 replcates 合併。 這就像用獨立的 sigma 估計做一個混合模型。
在下圖中,“平均”水平與您在帖子中所做的一樣,所有重複項都在一起。 注意置信區間的差異。
注意:JMP 14有區間作用,所以這個很容易製作。
所以要處理的問題是,如果表格像發布的文件一樣採用高格式,則行中的重複項將被視為獨立觀察。 我將做出一個安全的假設,即它們不是獨立的。 重複取決於實驗單元。 所以這裡是如何使用它。 另一種方法是像這樣拆分錶:
current Datatable() << Split(
Split By( :Replicate ),
Split( :inhibition ),
Group( :isolate, :ppm ),
Sort by Column Property
);
不確定哪個是最好的處理重複的方法,平均值還是幾何平均值?
但無論哪種方式,很容易看出強制均值的置信區間(具有獨立重複的 tall 格式的表。)比拆分錶的均值區間窄很多。 這是因為通過將重複視為獨立來擴大自由度。

此圖的數據來自以不同方式擬合的數據的 EC50,它不是默認報告圖。