library(pastecs) library(BayesFactor) placebo<-c(3, 3, 4, 4, 4, 5, 6, 6, 6, 6, 6, 7, 7, 8, 8, 9, 9, 11, 11, 11, 12, 12, 12, 12, 12, 12) oxy<-c(3, 4, 4, 6, 6, 7, 8, 8, 8, 8, 9, 9, 10, 10, 10, 11, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12) round(stat.desc(placebo, basic = F), 1) round(stat.desc(oxy, basic = F), 1) placebo<-c(3, 3, 4, 4, 4, 5, 6, 6, 6, 6, 6, 7, 7, 8, 8, 9, 9, 10, 10, 10, 11, 11, 11, 12, 12, 12, 12, 12, 12) round(stat.desc(placebo, basic = F), 1) zak<-data.frame(gl(2, 29), c(placebo, oxy)) names(zak)<-c("Group", "MU") zak$Group<-factor(zak$Group, labels = c("Placebo", "Oxytocin")) wilcox.test(MU~Group, zak) yuenbt(placebo, oxy, tr=.2, alpha=.05, nboot = 2000, side = T) pb2gen(placebo, oxy, alpha=.05, nboot=2000, est=median) d<-(mean(oxy)-mean(placebo))/sd(placebo) ttestBF(oxy, y = placebo)