I have a longitudinal dataset of radiation exposures of an occupational cohort. A percentage of the exposure values are missing and I would like to multiply impute the missing values (it is one option of several we are comparing). The data are recorded in long format (one row for each exposure entry) and there are multiple exposure measurements per worker. However, the data are time-unstructured (different data collection schedules for each worker) and unbalanced.
I want to account for the correlation between repeated measurements on the same worker. However, because of the time-unstructured nature of the dataset, I am unable to convert my dataset into wide format and impute that way. I have begun reading about about using multilevel imputation for such a scenario, but I rather unfamiliar with this approach, including within R. Is this an appropriate method to investigate?