帶有刪失數據的配對 t 檢驗您好,我有經過審查的微生物計數數據(日誌 10)。希望將測試期與基線進行比較,並確定微生物數量是否在統計上顯著減少。據我所知,只有 Gen Reg 可以進行審查數據並使用阻塞因子。阻塞因子是我能想到的將配對測試和基線結果聯繫起來的唯一方法。但是我似乎沒有足夠的自由度來進行測試?這些塊對我來說沒有意義,JMP 解釋它們的方式?有什麼建議?
分析回复:帶有刪失數據的配對 t 檢驗讓我們從您的數據開始。 它有問題。
根據您的願望陳述,我對兩個時期的非參數分佈估計進行了比較。

這是結果:

你看到問題了嗎? 基期沒有非參數估計。
這將您的數據按期間與審查員製表。 似乎你所有的基期觀察都被審查了。 測試期間的觀察不是。

這是您的數據的點圖,按週期著色。 所有基本觀測值都具有相同的值。

假設數據不存在上述問題。 現在的問題是,如果要比較兩個時期,為什麼要添加 Block? 它是一個解釋變量,可以描述兩個時期之間的差異嗎? 將其放入模型的決定不應取決於是否可以產生統計數據。 所以我建議從你的模型中刪除 Block,看看結果是否對你有意義,假設數據是好的。
除了 GenReg,另外兩個平台可以幫助您研究帶有分組變量的刪失觀察數據:帶有分組的 Life Distribution(上面的第一個屏幕截圖)和 Fit Life by X。 如果您只需要比較觀察組,則兩者可能更合適。
但最終,您的數據存在嚴重問題。 請不要通過運行軟件推土機數據,並對結果感到滿意。 可能很危險。
回复:帶有刪失數據的配對 t 檢驗作為 Peng,謝謝你 - 我的邏輯是錯誤的,因為在 Censored 列中將觀察標記為“Y”或“N”。 “Y”應該出現在 Period = Base 且 TSA.Count.L10 = 5.54 的地方。我使用的公式如下。但現在我意識到它不起作用,所以我只是手動編碼了審查數據。另外,當我堆疊數據時,我以某種方式使所有基值 = 5.54,這是不正確的。我可能在嘗試編寫 Censor 公式時這樣做了。
塊因子的原因是這實際上是一個配對 t 檢驗。比較的是基礎與測試期間的 TSA.Count.L10,在每件設備內,例如在製冰機 + 蓋子內。我會直接使用 Paired T-test 平台,但它不考慮審查數據。Gen Reg 可以。然後塊指定配對比較,至少這是我的希望。
問題:數據被右審查——Gen Reg 平台如何知道審查在哪一方?
If( :Period == "Base" & (:TSA.Count.L10 == 5.54),
"Y",
"N"
)
回复:帶有刪失數據的配對 t 檢驗感謝您更正數據。 我想我更好地理解了這個問題。
首先,讓我回答最後一個問題。 如果您提供審查列,則軟件假定它使用您在啟動對話框中選擇的審查代碼指示正確審查。 這是所有涉及審查的 JMP 平台所遵循的慣例。 如果您想表達不同類型的審查,例如左審查,您必須遵循約定以不同方式創建數據。 這個談話可能是最有幫助的:審查數據分析簡介.
我在文檔中找不到 GenReg 支持配對測試。 如果您能指出該位置,我將不勝感激。 另一方面,我會使用其他工具來解決這個問題。
首先,我拆分數據以形成一個新數據:

結果如下:

如果一個人會進行配對 t 檢驗,如果沒有數據被刪失,那麼首先會計算各個行的成對差異。 那麼問題就變成了差異是否與零沒有區別。 但我們審查了數據。 我假設正確的審查。 否則,相應地更改計算。
首先,我將創建一個列作為測試 - 基礎,就像對配對 t 檢驗所做的那樣。

現在看看這些值,即差異,並思考它們的含義。 查看第一行:4.91 - 5.54 = -0.635583025471626。 這裡 5.54 被刪失,這意味著實際值大於 5.54。 所以實際差異會小於-0.635583025471626。 因此,這使得 -0.635583025471626 成為左刪失觀察。
同樣的邏輯適用於所有基本觀察值被審查的行。 如果未刪失基本觀測值,則差異是精確值。
現在創建兩個新列,就像這裡所說的:審查數據分析簡介.

對於單個行,如果 Left 缺失,Right 不是,這是一個左刪失觀察。 如果 Left 和 Right 都沒有丟失,並且它們相等,那麼這是一個精確的觀察。 使用壽命分佈如下:

然後擬合一個分佈,Normal 對我來說看起來不錯。

現在問題變成了位置參數估計是否顯著不同於零。 如果答案是肯定的,那麼這兩組是顯著不同的。 不,否則。 在這個例子中,假設我猜對了數據中的審查意味著正確審查,那麼兩組看起來沒有顯著差異。 如果我猜錯了,差異的審查類型將是正確審查。 然後相應地更改步驟。
回复:帶有刪失數據的配對 t 檢驗謝謝 Peng,我了解您的解決方案,這正是我想要的。此外,我之前的措辭有點含糊:我試圖說 Gen Reg 適用於審查數據,但不適用於配對數據。問候,保羅