Error in drm() Fitting dose-response models using 4PL
0
0
Entering edit mode
5.2 years ago
mschmidt ▴ 80

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 • 2.4k views
ADD COMMENT

Login before adding your answer.

Traffic: 1817 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6