[R 통계분석] Box's M test (공변량 행렬 동질성 검정) | MVTests 패키지

[R 통계분석] Box's M test (공변량 행렬 동질성 검정) | MVTests 패키지

#1. 방법

BoxM(data, group)

data에는 반드시 변수가 2개 이상인 data.frame이 들어가야 합니다. 공분산행렬을 비교하기 때문입니다.

#2. 예제

> 코드

#Box's M test를 하기 위해 MVTests 패키지를 먼저 설치합니다.
install.packages("MVTests")

#패키지 불러오기
library(MVTests)

#표본생성
set.seed(1)
height_male=rnorm(50,0,1)
height_female=rnorm(50,0,1)
weight_male=rnorm(50,0,1)
weight_female=rnorm(50,0,1)

#형식에 맞게 수정
height_all=c(height_male,height_female)
weight_all=c(weight_male,weight_female)
group=c(rep(1,50),rep(2,50))
DF=data.frame(height_all,weight_all)

#Box's M test 수행.
BoxM(data=DF[,1:2],group=group)

> 실행결과

> BoxM(data=DF,group=group)
$`Chisq`
[1] 1.801662

$df
[1] 3

$p.value
[1] 0.6145734

$Test
[1] "BoxM"

attr(,"class")
[1] "MVTests" "list"