Установочные кривые с пакетом DRC в R?

Question

Я пытаюсь согласовать кривые с пакетом DRC в R.

Пример:

 x_yrs<-c(2005, 2006, 2007, 2008, 2009, 2010, 2011, 2012, 2013, 2014, 
     2015, 2016, 2017)
 y<-c(1.89, 0.34, 0.47, 2.46, 2.13, 7.49, 47.24, 117.84, 202.8, 322.7, 
 540.72, 744.22, 1148.7)

 MaxPop<-110000
 Y_Adj<-y/MaxPop

 EV<-drm(y~ x_yrs,fct = LL.3(fixed = c(NA, NA, NA)))
 plot(EV, broken = TRUE, type = "all")

 EV<-drm(y~ x_yrs,fct = LL.5(fixed = c(NA, NA, NA, NA, NA)))
 plot(EV, broken = TRUE, type = "all")

 x_yrs_Adj<- x_yrs-2004
 EV<-drm(Y_Adj~ x_yrs_Adj,fct = LL.5(fixed = c(NA, NA, NA, NA, NA)))
 plot(EV, broken = TRUE, type = "all",xlim = c(0, 40), ylim = c(0, 1))

Я бы хотел, чтобы максимальное значение кривой было "1" или«MaxPop», т. Е. В качестве верхней асимптоты.

Как мне изменить модель drm для достижения этой цели?

Answer 1

" Я бы хотел, чтобы будущая численность населения достигла 110 000 человек. " Я не думаю, что будет возможно согласовать модель с этим ограничением на основе предоставленных вами данных.Ответ, который у вас есть на поддержку функции, даже близко не подходит к этой (потенциально?) Асимптотической области.Поэтому я думаю, что вам нужно переосмыслить свой подход.

Кроме того, в drc вы можете реализовать ограничения, указав значения для конкретных параметров через аргумент функции fixed.

EV <- drm(Y_Adj ~ x_yrs_Adj, fct = LL.5(fixed = c(NA, 0, 1, NA, NA)))

Вы можете узнать об отдельных параметрах, например, ?LL.5:

  LL.5(fixed = c(NA, NA, NA, NA, NA), names = c("b", "c", "d", "e", "f"), ...)

[...]

Логистика с пятью параметрамифункция задается выражением

      f(x) = c + \frac{d-c}{(1+\exp(b(\log(x)-\log(e))))^f}

Так что в этом случае мы устанавливаем c на ноль, а затем фиксируем d = 1.

Давайте покажем график

plot(EV, broken = TRUE, type = "all", xlim = c(0, 40000), ylim = c(0, 1))

Вы можете увидеть проблему здесь.Поскольку у вас нет поддержки значений x_yrs_Adj ближе к асимптотическому поведению функции, ваша подгонка (и результирующие оценочные параметры) будут плохими.