R markdown file - iron deficiency for model 1