This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Error in drm() Fitting dose-response models using 4PL

I want to use a R script (the beginning below) to calculate ELISA results. My ELISA kit has 5 standards, and sample result is (the first standard is 'zero' but it need to be changed into non-0)

OD conc 1.16544997692108 0.000001 0.861099988222122 0.25 0.54805001616478 1 0.344400003552437 4 0.112999998033047 16

library("drc")
#####ELISA TEST
rm(list=ls())
library (drc)
###copy data for standards from excel (header – OD, conc) in windows use paste function
stdcrvdata <- read.table(pipe("pbpaste"), sep="\t", header=T) 
stdcrvdata$logconc <-log10(stdcrvdata$conc)# log10 from conc
plot(stdcrvdata$logconc, stdcrvdata$OD, main="log standard curve", xlab="x=log(conc)", ylab="y=OD")
fit<-drm(formula =   OD ~ logconc , data = stdcrvdata, fct = LL.4())

the last line of the script (fitting) results with error

Error in optim(startVec, opfct, hessian = TRUE, method = optMethod, control = list(maxit = maxIt, : initial value in 'vmmin' is not finite Error in drmOpt(opfct, opdfct1, startVecSc, optMethod, constrained, warnVal, : Convergence failed

how to fix that? The script works when used for 9-standard set.

r fitting error

0 answers

No answers yet.

Log in to answer this question.