DOE 0323_일원분산분석(R 실습)
#< 0325 R 코드> # R 4.3.3
group<-c(rep(1,4),rep(2,4),rep(3,4),rep(4,4),rep(5,4))
y<-c(2.4,2.7,3.1,3.1,
0.7,1.6,1.7,1.8,
2.4,3.1,5.4,6.1,
0.3,0.3,2.4,2.7,
0.5,0.9,1.4,2.0)
sol<-cbind(group,y)
group<-as.factor(group)
aov1<-aov(y~group)
plot(y~group)
summary(aov1)
summary.lm(aov1)
tapply(y,group,mean)
par(mfrow=c(2,2))
plot(aov1)
bartlett.test(y~group)
pairwise.t.test(y,group,p.adjust="none",pool.sd=TRUE)
pairwise.t.test(y,group,p.adjust="bonferroni",pool.sd=FALSE)
a.tukey<-TukeyHSD(aov1, ordered=TRUE)
a.tukey
tapply(y,group,mean)
tapply(y,group,var)
tapply(y,group,sd)
by(y, group, FUN=function(x){c(mean(x), sd(x))}) #-"그룹별 반응값의 평균과 표준편차를 보여줌"
model.tables(aov1)- "그룹별 수준효과 보여줌"
model.tables(aov1,type="mean") #- #그룹별 평균 보여줌#
summary(aov1) #-"일원분산분석결과"
summary.lm(aov1)#-"일원분산분석결과+그룹간 비교결과 제시"
----------------------------
L<-LSD.test(aov1, "group")
L
---------------------------------
install.packages("agricolae")
library(agricolae) #-agricolae패키지를 이용한 다중비교
L1<-LSD.test(aov1,"group",p.adj="none",main="Y of different groups")
L1
-------------------------------
library(agricolae)#@-agricolae패키지를 이용한 다중비교
LSD.test(aov1,"group",p.adj="none",main="Y of different groups")
L<-LSD.test(aov1, "group")
L
댓글
댓글 리스트-
작성자storyteller 작성자 본인 여부 작성자 작성시간 22.05.01 install.packages("agricolae")
library(agricolae)-)-agricolae패키지를 이용한 다중비교
L1<-LSD.test(aov1,"group",p.adj="none",main="Y of different groups")
L1 -
작성자storyteller 작성자 본인 여부 작성자 작성시간 25.09.26 library(agricolae)
#@-agricolae패키지를 이용한 다중비교
LSD.test(aov1,"group",p.adj="none",main="Y of different groups")
L<-LSD.test(aov1, "group")
L