\begin{Schunk}
\begin{Sinput}
 library(hexbin)
 hbin <- hexbin(dat$educ, dat$sei, xbins = 40)
 plot(hbin, xlab = "Education", ylab = "Socio Economic Status", 
      main =  "Hexagon-binned Data Plot",
      lcex = 0.6)
\end{Sinput}
\end{Schunk}
