R mauchly.test 莫奇利球形度檢驗

說明

測試 Wishart-distributed 協方差矩陣(或其變換)是否與給定矩陣成比例。

用法

mauchly.test(object, ...)
## S3 method for class 'mlm'
mauchly.test(object, ...)
## S3 method for class 'SSD'
mauchly.test(object, Sigma = diag(nrow = p),
             T = Thin.row(Proj(M) - Proj(X)), M = diag(nrow = p), X = ~0,
             idata = data.frame(index = seq_len(p)), ...)

參數

object	類 SSD 或 mlm 的對象。
Sigma	成正比的矩陣。
T	變換矩陣。默認情況下根據 M 和 X 計算。
M	說明外部投影的公式或矩陣(見下文)。
X	說明內部投影的公式或矩陣(見下文)。
idata	說明intra-block設計的數據幀。
...	要傳遞給其他方法或從其他方法傳遞的參數。

細節

這是一個通用函數，具有類 "mlm" 和 "SSD" 的方法。

基本方法是針對 SSD 類的對象，mlm 對象的方法隻是提取 SSD 矩陣並使用相同的選項和參數調用相應的方法。

T 參數用於在測試之前轉換觀察結果。這通常涉及到intra-block差異的轉換，但可能會遇到更複雜的within-block設計，從而需要更複雜的轉換。矩陣 T 可以直接給出，也可以指定為 M 和 X 所跨越的空間上的兩個投影之間的差，而這又可以作為矩陣或關於 idata 的模型公式給出(測試對於商空間的參數化是不變的 M/X )。

此測試的常見用途是在重複測量設計中，使用X = ~1。這與在未變換的協方差矩陣中測試複合對稱性幾乎但不完全相同。

請注意，默認值涉及 p ，它在內部計算為 SSD 矩陣的維度，以及 stats 命名空間中的幾個隱藏函數，即 proj ，它根據設計矩陣或模型公式計算投影矩陣和 Thin.row 從矩陣中刪除線性相關行，直到它具有完整的行秩。

值

"htest" 類的對象

注意

p 值與 SAS 的 p 值略有不同，因為二階項包含在以下公式中的漸近近似中：R.

例子

utils::example(SSD) # Brings in the mlmfit and reacttime objects

### traditional test of intrasubj. contrasts
mauchly.test(mlmfit, X = ~1)

### tests using intra-subject 3x2 design
idata <- data.frame(deg = gl(3, 1, 6, labels = c(0,4,8)),
                    noise = gl(2, 3, 6, labels = c("A","P")))
mauchly.test(mlmfit, X = ~ deg + noise, idata = idata)
mauchly.test(mlmfit, M = ~ deg + noise, X = ~ noise, idata = idata)

參考

T. W. Anderson (1958). An Introduction to Multivariate Statistical Analysis. Wiley.

也可以看看

SSD , anova.mlm , rWishart