I\'m analysing data from this kaggle dataset. I\'m trying to analyse the effect of bmi from smokers on medical cost, but also I want to fit a hierarchical model with cluster