2010-08-03 191 views
1

我想在R中複製SPSS的線性判別分析輸出,而且我很難找到執行m-box測試的方法。R Box'M測試homocedasticity

我發現的唯一的事情是在論壇上發佈了一些代碼來手動實現這個過程,但我想知道是否沒有爲此目的已經包含在語言本身中的任何東西。

+0

我不熟悉「m-box」測試。 SPSS專用術語?我聽說過Box-Pierce,Box-Jenkins和Box-Q測試,但這些都與自相關有關。 – 2010-08-03 23:22:30

+0

我拼錯了名字,它被稱爲Box'm測試。無論如何,它不像我那樣常用(我剛開始學習這個主題),可能是因爲它的高敏感性,正如約翰在下面評論的那樣。 它用於測試MANOVA假設之一,看起來大多數人要求它是SPSS用戶,因爲該軟件將其作爲標準輸出之一。 只是爲了防止別人感興趣,這裏是SPSS算法的規範: ftp://ftp.spss.com/pub/spss/statistics/spss/algorithms/ap14boxm.pdf (維基百科有一個定義,但只有在德語版本中) – Sam 2010-08-04 14:15:54

回答

1

code,可以發現一個簡單的rseek搜索。這通常不會完成,因爲它非常高的靈敏度會導致顯着的p值,這可能沒有多大意義。

編輯: 那舊的鏈接不工作了,但事實證明,試驗是在Biotools公司包的功能boxM實現。這仍然是一個相對容易的搜索。而且,你可能不應該像所有這些測試一樣使用它。你應該仔細檢查你的協方差矩陣和你的假設。

+0

你說得對,那是我所指的代碼,但它是用矩陣手動操作的,我想知道是否有一個例程已經在語言中實現。 – Sam 2010-08-04 13:22:41

+0

不是正確答案...代碼是代碼。你在網上有什麼功能。將它添加到R,現在它是你的R環境的一部分。 – John 2010-08-04 14:46:04

+1

我知道這已經有一段時間了,但是鏈接被破壞了嗎?你可以在這裏發佈代碼嗎? – 2015-11-04 22:09:58

1

在包Biotools公司你可以找到函數boxM(數據分組)。它根據一個分類因子對從多元正態數據獲得的協方差矩陣的同質性進行Box的M檢驗。該測試基於卡方近似。