Я хотел использовать функцию mselec из пакета drc для выбора наилучшей модели набора данных доза-реакция. Однако mselect не работает, если вы используете его внутри функции.
Работает следующий код:
library(drc)
ryegrass.m1 <- drm(rootl~conc, data = ryegrass, fct = LL.4())
mselect(ryegrass.m1,list(LL.5(), LN.4(), W1.4(), W2.4()))
Но не этот:
best.fit=function(data){
model1=drm(rootl~conc, data=data, fct=LL.4())
M1=drc::mselect(model1, list(LL.5(), LN.4(), W1.4(), W2.4()))
return(M1)
}
best.fit(ryegrass)
Я думаю, что проблема связана с окружением в R, но я не знаю, как это исправить. Кто-нибудь может мне помочь?