#df = data.frame(subdata,mu.tilde$muhat,mu.hat$muhat)

library(maps,mapproj)
library(ggplot2)
FL <- map_data('county','florida')
names(FL)[6] <- 'County'

mydat <- merge(FL,subdat)

fig <- qplot(long,lat,data=mydat,group=group,fill=mu.tilde$muhat,geom='polygon')+
#  facet_wrap(~Date,ncol=3)+
  geom_path(colour='grey10',lwd=0.5)+
  scale_fill_gradient2('',low='cyan',mid='white',high='navy',
                       guide='colorbar',limits=c(0,0.0003),na.value='orange', 
                       breaks=c(0,0.0001,0.0002,0.0003), 
                       labels=c('0','1e-4','2e-4','3e-4'))+
  guides(fill=guide_colorbar(barwidth=75,barheight=1,direction='horizontal'),
         cex=1.2)+theme_bw(base_size=30)+xlab('longitude')+ylab('latitude')+
  theme(legend.position='bottom',
        axis.ticks=element_blank(),
        line=element_blank(),
        axis.text=element_blank(),
        panel.border=element_rect(color='black',size=1.2),
        axis.line=element_line(colour='black'),
        legend.margin=margin(-10,0,0,-50),
        legend.box.margin=margin(0,0,0,0))
 ggsave('fig92.pdf',width=16,height=12)