Я использую функцию 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):
Но, если я заставлю верхнюю и параметры нижнего предела, это хорошо:
b <- drm(Value ~ Dose, data = data, fct = LL.4(fixed = c(NA,100,0,NA)))
Это странно, потому что второй параметр должен быть ниже, а третий должен быть верхним, и я должен написать их в обратном порядке. Но результат в порядке:
Я новичок в этой библиотеке, но все примеры, которые я нашел, сделаны именно так, и я не смог найти причину, почему его инверсия кривой.