rm(list=ls(all=TRUE))  # clear all variables
graphics.off()  # clear all graphics
# Sternberg Search
# Greg Francis
# PSY 626
# 01 November 2020

# load full data file
SSdata<-read.csv(file="OneSubject.csv",header=TRUE,stringsAsFactors=FALSE)
# Treat MemorySetSize as a factor
SSdata$MemorySetSize = factor(SSdata$MemorySetSize)
levels(SSdata$MemorySetSize) = c("Small", "Medium", "Large")
# Treat Condition as a factor
SSdata$Condition = factor(SSdata$Condition)

# NHST ANOVA
summary(aov(RT ~ Condition*MemorySetSize, data=SSdata))


# load the BayesFactor library
library(BayesFactor)

# run the default BayesFactor ANOVA
bf = anovaBF(RT ~ Condition*MemorySetSize, data=SSdata)
# With more sampling iterations
newbf = recompute(bf, iterations = 500000)

# run all possible models
bf2 = anovaBF(RT ~ Condition*MemorySetSize, data=SSdata, whichModels="all", iterations=500000)

# Compare simpler models to "full" model (subtract one term)
bf3 = anovaBF(RT ~ Condition*MemorySetSize, data=SSdata, whichModels="top", iterations=500000)

# Compare more compliated models to "null" model (add one term)
bf4 = anovaBF(RT ~ Condition*MemorySetSize, data=SSdata, whichModels="bottom", iterations=500000)