[R] DRC: Effective doses versus Predicted values
Iban
kalekume at gmail.com
Thu Aug 12 14:18:16 CEST 2010
Hi!
I want to use the DRC package in order to calculate the IC50 value of an
enzyme inhibition assay.
The problem is that the estimated ED50, is always out of the fitted curve.
In the example below, I had a ED50 value of 2.2896,
But when I predict the response level for this concentration I get a value
of 45.71 instead of the expected value of 50.
This is my data:
#Dose unit is concentration (mM); response unit is % of inhibition
dat <- data.frame(dose=c(0.480,0.957,2.870,4.780,6.700,9.570),
response=c(11.88171,28.62986,51.50847,65.51793,77.29431,82.83489))
model <- drm(response ~ dose, data=dat, fct= W2.3())
ED(model,50)
#get a ED50 valu of2.2896
predict(model, data.frame(dose=2.2896))
#get a response level of 45.714
plot(model)
points(2.2896,50, col="red")
#the ED50 value is out of the fitted curve
--
View this message in context: http://r.789695.n4.nabble.com/DRC-Effective-doses-versus-Predicted-values-tp2322613p2322613.html
Sent from the R help mailing list archive at Nabble.com.
More information about the R-help
mailing list