get_lunatic_data = function(){
    fn = '../../Data/lunatic.dat'
    DF = read.table(fn, header=FALSE, skip=14, col.names=c('COUNTY', 'NBR', 'DIST', 'POP', 'PDEN', 'PHOME'), fill=TRUE)
    DF[is.na(DF)] = 0
    DF
}