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

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

#1. 방법

BoxM(data, group)

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

#2. 예제

> 코드

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

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

#표본생성
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)

    Box's M-test for Homogeneity of Covariance Matrices

data:  DF
Chi-Sq (approx.) = 1.8017, df = 3, p-value = 0.6146