# bias correction in H2Omg

five <- with(all_empty ,all_empty[ hour( StartTime ) >= 5 & hour( StartTime ) < 6 , ] )
six <- with(all_empty ,all_empty[ hour( StartTime ) >= 6 & hour( StartTime ) < 7 , ] )
seven <- with(all_empty ,all_empty[ hour( StartTime ) >= 7 & hour( StartTime ) < 8 , ] )
eight <- with(all_empty ,all_empty[ hour( StartTime ) >= 8 & hour( StartTime ) < 9 , ] )
nine <- with(all_empty ,all_empty[ hour( StartTime ) >= 9 & hour( StartTime ) < 10 , ] )
ten <- with(all_empty ,all_empty[ hour( StartTime ) >= 10 & hour( StartTime ) < 11 , ] )
eleven <- with(all_empty ,all_empty[ hour( StartTime ) >= 11 & hour( StartTime ) < 12 , ] )
twelve <- with(all_empty ,all_empty[ hour( StartTime ) >= 12 & hour( StartTime ) < 13 , ] )

median(five$H2Omg)
median(six$H2Omg)
median(seven$H2Omg)
median(eight$H2Omg)
median(nine$H2Omg)
median(ten$H2Omg)
median(eleven$H2Omg)
median(twelve$H2Omg)



> median(five$H2Omg)
[1] 0.003249901
> median(six$H2Omg)
[1] 0.02742861
> median(seven$H2Omg)
[1] 0.06084972
> median(eight$H2Omg)
[1] 0.09536454
> median(nine$H2Omg)
[1] 0.02219424
> median(ten$H2Omg)
[1] 0.03273413
> median(eleven$H2Omg)
[1] 0.02152175
> median(twelve$H2Omg)
[1] 0.006491296

### here is plot using the bias correction indicated above