Question: Error in drm() Fitting dose-response models using 4PL
0
gravatar for mschmidt
8 weeks ago by
mschmidt0
mschmidt0 wrote:

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.

error R fitting • 94 views
ADD COMMENTlink modified 8 weeks ago • written 8 weeks ago by mschmidt0
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1354 users visited in the last hour