Обратная кривая подгонка с использованием Drm - PullRequest
0 голосов
/ 26 марта 2020

Я использую функцию drm библиотеки dr c для подбора кривой доза-эффект. Единица ответа, которую я использую - это процент. Это мой код:

data <- data.table(Dose = dose, Value = response)
b <- drm(Value ~ Dose, data = data, fct = LL.4(fixed = c(NA,NA,NA,NA)))

Поскольку я использую код из java, используя Rserve, я извлекаю coefs и печатаю данные с другой библиотекой.

double[] coefs = c.eval("coef(b)").asDoubles();

Проблема заключается в том, что он выполняет подбор кривой, что он правильный, но инвертированный (отрицательный HS): enter image description here

Но, если я заставлю верхнюю и параметры нижнего предела, это хорошо:

b <- drm(Value ~ Dose, data = data, fct = LL.4(fixed = c(NA,100,0,NA)))

Это странно, потому что второй параметр должен быть ниже, а третий должен быть верхним, и я должен написать их в обратном порядке. Но результат в порядке:

enter image description here

Я новичок в этой библиотеке, но все примеры, которые я нашел, сделаны именно так, и я не смог найти причину, почему его инверсия кривой.

...