I am trying to estimate the mean value of a lab test reading igg1_norm across different categories of a variable type (5 categories).
db <