############ csv ????????? ???????????? ?????? ########################## KY <- read.csv(file = "/Users/jinmedi/Desktop/KY.csv",header = TRUE, sep = ",") ############ pwr ????????? ?????? ?????? ???????????????. ################## library(pwr) ############ ??? ????????? ????????? ????????? ?????? ##################### mean.1 <- mean(KY$score[KY$group == 1]) mean.2 <- mean(KY$score[KY$group == 2]) ############ ??? ????????? ??????????????? ????????? ?????? ################# sd.1 <- sd(KY$score[KY$group == 1]) sd.2 <- sd(KY$score[KY$group == 2]) ############ ?????? ???????????? ???????????? ?????? ####################### summary(KY) ############ ?????? ???????????? ???(row) ????????? ???????????? ?????? ######## nrow(KY) ############ ?????? ???????????? ???(columns) ????????? ???????????? ?????? #### length(KY) ############?????? ?????? ????????? ?????? ############################ sqrt(2) ############ ???????????? ????????? ?????? ############################## abs(1-3) ############ effect size??? ????????? ?????? ######################### effectSize <- abs(mean.1-mean.2)/(sqrt((sd.1^2+sd.2^2)/2)) effectSize <- abs(mean.1-mean.2)/(sqrt((sd.1^2*(nrow(group.1)-1)+sd.2^2*(nrow(group.2)-1))/((nrow(group.1)-1)+(nrow(group.2)-1)))) ############ power analysis??? ?????? ?????? n ?????? ????????? ?????? ############################### pwr.t.test(d= effectSize,power=.8,sig.level=.05,type="two.sample",alternative="two.sided") ############ ????????? ????????? ???????????? ???????????? ????????? ?????? ?????? ############################ group.1 <- subset(KY, KY$group==1) group.2 <- subset(KY, KY$group==2) ############ ??? ????????? ???????????? ???????????? ?????? ?????? ####################################### shapiro.test(group.1$score) shapiro.test(group.2$score)