Блог

Я пытаюсь вычислить кривую мощности (уравнение Фрейндлиха) в R, используя drm на основе, https://www.statforbiology.com/nonlinearregression/usefulequations . Вот код:

   model <- drm(numSpecies ~ Area, fct = DRC.powerCurve(),
                 data = speciesArea)
  

DRC.powerCurve() предоставляет оценку для первой производной (т. Е. Наклона), однако мне было интересно, возможно ли получить вторую производную? Согласно документации R (http://ftp.uni-bayreuth.de/math/statlib/R/CRAN/doc/packages/drc.pdf ) вы можете указать, какие параметры следует оценивать в операторе ‘fixed’

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

Однако, когда я пытаюсь сделать это в R с помощью DRC.PowerCurve устанавливая оценки параметров равными нулю, я получаю сообщение об ошибке «Неверный ‘фиксированный’ аргумент». Я также не уверен, какие начальные значения соответствуют второй производной.

Я не могу найти практически никакой информации о DRC.Кривая мощности, не говоря уже о получении второй производной. На самом деле, DRC.PowerCurve даже не указана в качестве потенциальной функции в документации R или в getMeanFunction() .