等值檢驗中的P值計算公式一些客戶詢問了有關等效性測試的兩件事。
他們將安裝與VBA鏈接JMP的Antomotive數據關聯繫統(Visual Basic for application)
並希望在兩台設備內進行等效性測試。
問題1
如何在JMP的等效性測試中使用P值計算公式?
問題2
在等值測試的p值計算中,“中位數”是否可用而不是“均值”?
Re:等效性檢驗中的P值計算公式對於步驟4,突出顯示2列中位數,然後在列標題中單擊鼠標右鍵。進行如下所示的選擇以創建一系列差異。
Create a column for difference in medians為中位數創建一個差異列
對於步驟5,在步驟4中創建的列上運行分佈平台,並檢查分位數。例如,我們可以看到95%的bootstrap樣本的中位數差異為-5.18和-4.89。
因此,如果我們測試中位數差異超過5的假設,例如,我們將不能拒絕0.05水平的假設,因為我們的經驗置信區間並未被區間完全覆蓋(-5,5)。也就是說,我們沒有足夠的證據在0.05水平上得出結論,分佈的中位數差異小於5。
(事實上,我使用兩個均勻分佈來模擬這些數據:U(0,1)和U(5,6)。所以,這個結果就是我們所期望的。)
Examine distribution of difference in medians檢查中位數差異的分佈
Re:等效性檢驗中的P值計算公式嗨布拉迪,
感謝您的詳細解釋,如果我能理解更多,那就更好了。
你能解釋一下4)和5)
我不知道如何操縱JMP for 4)&5)
4)當引導表出現時,使用列變換創建一個新的列,這是中位數的差異。
5)使用此列上的分佈平台確定(alpha / 2)和(1-alpha / 2)分位數,它們形成中位數差異的經驗100(1-alpha)%置信區間
Re:等效性檢驗中的P值計算公式我應該補充一點,如果用戶有JMP Pro,我很容易執行我所描述的引導:
1)使用Fit Y by X將連續變量轉換為Y角色,將分類變量轉換為X角色,然後運行分析。
2)使用紅色三角形菜單打開分位數報告。
3)在分位數報告中,右鍵單擊中間列並選擇Bootstrap。
4)當引導表出現時,使用列變換創建一個新的列,這是中位數的差異。
5)使用此列上的分佈平台確定(alpha / 2)和(1-alpha / 2)分位數,它們形成中位數差異的經驗100(1-alpha)%置信區間。
Re:等效性檢驗中的P值計算公式嗨Ikju,
1)TOST程序的p值是每個單側測試的p值中的較高者。例如,如果實際差異為4,則2個測試具有以下假設:
H0(測試1):d> = 4
H0(測試2):d <= - 4
假設測試1的p值為0.00000003,測試2的p值為0.98983。然後報告的等效性測試的p值是0.98983。
只有當BOTH無效假設被拒絕時,才宣布實際上等同的手段。
更多信息在這裡: https://www.jmp.com/support/help/13-2/Example_of_an_Equivalence_Test.shtml
2)據我所知,JMP沒有提供基於中值的等價測試。理論上這應該是可能的,我已經看過一篇關於此的文章,但是文章不能免費閱讀: http://onlinelibrary.wiley.com/doi/10.1002/pst.384/abstract
如果樣本足夠大,我想你可以通過從數據中獲取bootstrap樣本(重複採樣替換)來構建非參數等價測試,計算每種情況下的中位數差異,並使用模擬結果形成置信區間中位數的差異。如果置信區間完全包含在區間(-k,k)中,其中k是被認為具有實際意義的差異,則可以聲明樣本實際上等效。
乾杯,
布雷迪