I'm using a linear mixed model (R,lme4) to evaluate the effect of drug dose on auditory thresholds across a range of 14 frequencies: lmer(Threshold ~ Frequency + Condition + (1|Subject)
We are missing observations at the higher frequencies for a variety of reasons. A few cases are because we tested the first few subjects with a narrower range of frequencies. The majority of missing observations are due to the drug treatment increasing thresholds outside the range at which we can test them (our speaker can only get so load). Missing observations increase with frequency. At the highest 3 frequencies, we have 0 observations for the drug conditions. The attached plot compares our controls to one dose of the drug. The numbers in parentheses indicate proportion of subjects for which we have a response at the given frequency (we have 100% of responses below 500 Hz for all subjects).
Is there a way for the LMM to handle this data and allow us to make an inference about the difference in thresholds between conditions at the higher frequencies or would it best to use a different kind of model that deals with binary outcomes?
Here is the data: CSV file