세 그룹의 평균 비교

# R

install.packages('moonBook')
require(moonBook)
data(acs)
acs
summary(acs)
str(acs)

moonBook::densityplot(age~sex, data=acs)
shapiro.test(acs$age[acs$sex=="Male"])
shapiro.test(acs$age[acs$sex=="Female"])
whicox.test(age-sex, data=acs)

var.test(age~sex, data=acs) # 등분산을 확인하기 위한 방법, P<0.05일 경우 등분산X
t.test(age~sex, data=acs, var.equal=TRUE) # Student's t test in R-studio
t.test(age~sex, data=acs, var.equal=FALSE) # Welch's test in R-studio

'''

세 그룹간의 평균 비교
1) 결과값이 연속변수
2) 각각의 결과값이 정규분포
3) 각각의 정규분포가 서로 등분산

1),2),3) 만족 => ANOVA

'''
data(acs)
View(acs)
moonBook::densityplot(LDLC~Dx, data=acs) # 그래프 확인

# 정규분포 : p-value가 0.05보다 크면 정규분포 이다.
# NSTEMI, Unstable Angina : 정규분포 아님
# STEMI : 정규 분포 (p-value=0.6066)
shapiro.test(acs$LDLC[acs$Dx=="NSTEMI"]) # 정규분포 확인/P>0.05 정규분포!!!
shapiro.test(acs$LDLC[acs$Dx=="STEMI"]) # p-value = 0.6066, P>0.05 정규분포
shapiro.test(acs$LDLC[acs$Dx=="Unstable Angina"]) # p-value = 2.136e-07

out = aov(LDLC~Dx, data=acs)
shapiro.test(resid(out)) # P < 0.05 : 세 그룹 중에서 하나이상이 정규분포를 하지 않음!

shapiro.test(resid(aov(LDLC~Dx, data=acs)))

 

 

<출처>

https://youtu.be/jcxhcQkM5rg

 

#통계 #ANOVA #t-test #ttest #정규분포 #shapiro #test #aov #densityplot #r #평균 #var #등분산 #연속변수