library(psych)
setwd("/Users/Pria/Desktop/Data/")
##Load data
latency = read.csv("23138xTRP1.csv", header = TRUE)
##head shows first few lines of data
head(latency)
##Summary statd
summary(latency)
##More info on stats
describe(latency)
res.aov2 <- aov(Gal4.UAS ~ Gal4 + UAS, data = latency)
summary(res.aov2)
res.aov3 <- aov(Gal4.UAS.1 ~ Gal4.1 + UAS.1, data = latency)
res.aov3 <- aov(Gal4.UAS.1 ~ Gal4.1 + UAS.1 + Gal4:UAS, data = latency)
summary(res.aov3)