dff=readRDS("./dff.rds")

# 改一下網格點順序
SNdff=data.frame(S1.N86=dff$S1.N86 %>% unique %>% sort)
SNdff$reciprocal=1/(SNdff$S1.N86)
SNdff$rank=SNdff$reciprocal %>% rank

dff=merge(x=dff, y=SNdff, by=c('S1.N86'), all.x = T)
dff=data.frame(dff)

dff$S1.N86.new = factor(dff$S1.N86, levels=c(rev(c(dff$S1.N86 %>% unique %>% sort %>% as.character))))
dff$S1.N86.new %>% class

#確認一下病例分布
gggg2=aggregate(dff[, 15], list(dff$W1.E41,dff$S1.N86), sum) 


dff=dff[!is.na(dff$rainMA),]
dff=dff[!is.na(dff$t2mMA),]
dff$N %>% class

ggplot(dff)+
  facet_grid(S1.N86.new+S1.N86~W1.E41)+
  geom_line(data=dff, aes(x=date, y=rainMA/5), color='steelblue', size=0.1)+
  geom_point(data=dff, aes(x=date, y=t2mMA), color='red', size=0.05, alpha=0.5)+ 
  geom_bar(data=dff, aes(x=date, y=(dff$N)), stat="identity", color='#d9ef8b', alpha=0.1)+
  
  #scale_y_continuous(sec.axis = sec_axis(~.*5, name = "Relative humidity [%]"))+
  #geom_text(aes(label=N), vjust=1, color="white", size=3.5)+
  theme(axis.text.x = element_blank(),
        axis.ticks.x = element_blank(),
        axis.title.x=element_blank())+
  
  ggtitle(paste0("grid SickNum ",yy,' ',city))