#required packages
library(psych)
library(Rcmdr)
library(granova)

#Two groups
attach(Guyer)
  #Descriptives
    describe.by(cooperation, condition)
    par(mfrow=c(2,1))
    hist(subset(cooperation,condition=="P"), xlim=c(20,85), ylim=c(0,7), xlab="", main="P")
    hist(subset(cooperation,condition=="A"), xlim=c(20,85), ylim=c(0,7), col="lightblue", xlab="", main="A")
    par(mfrow=c(1,1))
  #ANOVA
    AnovaModel.1 <- aov(cooperation ~ condition, data=Guyer)
    summary(AnovaModel.1)
    #Visualization
    granova.1w(cooperation, condition)
  #t-test
    t.test(cooperation~condition, var.equal=T, data=Guyer)
    cor(as.numeric(Guyer$condition),Guyer$cooperation)^2
  #Regression
    LinearModel.1 <- lm(cooperation ~ condition, data=Guyer)
    summary(LinearModel.1)
detach(Guyer)

#Three groups
attach(Friendly)
  #Descriptives
    describe.by(correct, condition)
    par(mfrow=c(3,1))
    hist(subset(correct,condition=="Meshed"), xlim=c(18,42), ylim=c(0,7), breaks=7, xlab="", main="Meshed")
    hist(subset(correct,condition=="Before"), xlim=c(18,42), ylim=c(0,7), breaks=7, col="lightblue", xlab="", main="Before")
    hist(subset(correct,condition=="SFR"), xlim=c(18,42), ylim=c(0,7), breaks=7, col="red", xlab="", main="SFR")
    par(mfrow=c(1,1))
  #ANOVA
    AnovaModel.2 <- aov(correct ~ condition, data=Friendly)
    summary(AnovaModel.2)
    plotMeans(Friendly$correct, Friendly$condition, error.bars="conf.int", level=0.95)
    stripchart(correct ~ condition, data=Friendly, vertical=T)
    granova.1w(correct, condition, v.rng=1)  #note there is overlap due to two groups having equal means