You are describing a complex histogram, not a histogram. With ggplot
you can do it like this:
library(ggplot2) library(reshape2) x <- data.frame( Period = c(1,1,2,2,3,3,4,4), Sample = c("A","B","A","B","A","B","A","B"), Value1 = c(3,2,6,7,3,2,1,2), Value2 = c(1,0,5,2,2,0,2,5) ) mx <- melt(x, id.vars=1:2) ggplot(mx, aes(x=Period, y=value, fill=variable)) + geom_bar(stat="identity") + facet_grid(~Sample)

source share