added fit-hier.R