I have a data set that contains 40 firms over 15 years with firm and individual level data, with the individuals being sample randomly from each firm over 15 years. There are time invariant variables at the firm level and the panel design is unbalanced. How should I be looking to treat my mixed model to account for the random sampling? Any help or advice would be great.

Thank you.