R markdown file - iron deficiency (model 2)