####Dependent Variable Codes
#Frontal vs. Rostral and Caudal
Sero_long$FrontalVsRosCaud <- 2*(Sero_long$Region=="Rostral") -1*(Sero_long$Region=="Middle") -1*(Sero_long$Region=="Caudal")
#Rostral vs. Caudal
Sero_long$MiddlevsCaudal <- 0*(Sero_long$Group=="Rostral") + 1*(Sero_long$Group=="Middle") -1*(Sero_long$Group=="Caudal")

####Independent Variable Codes
#Linear contrast for Homecage vs. Shock
Sero_long$HCvsShock <- 2*(Sero_long$Group=="HC") -1*(Sero_long$Group=="IS") -1*(Sero_long$Group=="ES")
#Quadratic contrast code for IS vs. ES
Sero_long$ISvsES <- 0*(Sero_long$Group=="HC") + 1*(Sero_long$Group=="IS") -1*(Sero_long$Group=="ES")

modelFR1 <- lm(FrontalVsRosCaud~HCvsShock+ISvsES, data=Sero_long)