Background: I have the lmer model

lmer(log(log(Amplification)) ~ poly(Voltage, 3) + 
           (poly(Voltage, 3) | Serial_number), data = APD)

for a data frame APD. This contains lot of devices, separated by Serial_number. Each device's amplification was measured in the voltage range from 0 V to about 400 V. In total I have about 1500 devices (and nearly 80,000 observables).

Motivation: My goal is to get the corresponding voltage of a certain amplification for each device; Voltage(Amplification=150)?. With the help of extracting the coefficients of the polynomials I'm able to do an inverse regression and to calculate this voltage.

Problem: When I select the data range in different ways I calculate different voltages. This is more or less reasonable because the fits vary. How can I ensure which voltage is correct? When I take the whole data range the calculated (Amplification=150, Voltage(A=150))-point is far beside the measurement data. When I limit the data range then the point moves closer and closer to the measurement data. So obviously I can influence that strongly by selecting the data range.

Bottom line: My very goal is to find the very most reliable voltage for each device having an amplification of Amplification=150. A model which is often used in literature is


however I would have to invert the model in my case because I want to predict the Voltage given that Amplification=150.

The plots of the calculated data points:

Whole data range for fitting with lmer:

Single device complete data range

Limited data range with ignoring all measurement data below Amplification<100 and all higher than Amplification>200.

Single device

The origin curve of one single device:

Single subset/Serial_number

Minimal data set and script: https://files.fm/u/5yy22kkm

Thanks a lot in advance!

edit: I also posted this question on https://stackoverflow.com/questions/47119570/select-appropriate-data-range-for-mixed-model to gain more response.



Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Browse other questions tagged or ask your own question.