운동하는 사람들 중에서 체중의 평균이 그룹간 차이가 있나(anova)

 

문제:

귀무가설: 나이를 **세 그룹(30세 미만, 30~50세, 50세 이상)**으로 나누었을 때, 운동하는 사람들 중에서 체중의 평균이 그룹 간에 유의미한 차이가 없다

대립가설: 나이를 세 그룹으로 나누었을 때, 운동하는 사람들 중에서 체중의 평균이 그룹 간에 유의미한 차이가 있다.

 

가설 설정:

H0 (귀무가설): 나이별 그룹 간 체중 평균에 유의미한 차이가 없다

H1 (대립가설): 나이별 그룹 간 체중 평균에 유의미한 차이가 있다

근력운동 BE5_1   -> 1,8 비운동 2,3,4,5,6 운동

나이 age

몸무게 HE_wt

 

 

 

 

 

> s<-data[,c("age","BE5_1","HE_wt")]

 

> s<-s[s$BE5_1<9,]

 

> s<-s[complete.cases(s),]

 

> s$age_group<-cut(s$age,breaks=c(0,30,50,Inf),labels=c("under30","30to50","over50"))

 

> s$undong<-ifelse(s$BE5_1==1 | s$BE5_1==8,"notundong","undong")

 

> u<-s[s$undong=="undong",]

 

 

 

 

> bartlett.test(HE_wt~age_group, data=u)

 

        Bartlett test of homogeneity of variances

 

data:  HE_wt by age_group

Bartlett's K-squared = 72.463, df = 2, p-value < 2.2e-16

 

 

 

 

 

> one<-oneway.test(HE_wt~age_group, data=u, var.equal=FALSE)

> one

 

        One-way analysis of means (not assuming equal variances)

 

data:  HE_wt and age_group

F = 5.9957, num df = 2.00, denom df = 734.17, p-value = 0.002613

 

 

 

 

> pair<-pairwise.t.test(u$HE_wt,u$age_group,p.adjust.method="bonferroni")

> pair

 

        Pairwise comparisons using t tests with pooled SD 

 

data:  u$HE_wt and u$age_group 

 

       under30 30to50

30to50 0.6332  -     

over50 0.1694  0.0038

 

P value adjustment method: bonferroni