Watch this Author's files
Thank you for the interesting code! However, I suspect that there is an inconsistency in the interpretation of zone A and B. The totals and the plots are not compatible e.g. zoom in on clarke(200,243) and clarke(200, 160).