response<-rnorm(42, mean=10)
factor1<-factor(rep(c("A","B","C"), each=14))
factor2<-factor(rep(c(1,2,3,1,2,3,1,2,3,1,2,3,1,2), times=3))
factor3<-factor(rep(c("X","X","X","Y","Y","Y","X","X","X","Y","Y","Y","Y","Y"),times=3))
factor4<-factor(rep(c(rep(c("27 degrees","30 degrees"), each=6),"33 degrees","33 degrees"), times=3))

data1<-data.frame(response,factor1,factor2,factor3,factor4)


p1<- ggplot(data=data1, aes(x=factor2, y=response, fill=factor1))
p1<- p1+geom_bar(stat="identity")
p1<- p1+facet_grid(~factor4+factor3,scales = "free_x")+theme_bw()
p1<- p1+scale_y_continuous(expand = c(0,0))+theme(strip.background=element_rect( fill="white"))
p1<- p1 + scale_fill_brewer(guide = guide_legend(reverse=TRUE))
p1