[R] aproximate a titration kurve to the measure data.

[Bernardo Rangel Tura]

On Thu, 2009-10-01 at 10:18 -0700, awayguy wrote:
> yes, halo thank you.
> my measure data:
> v2 <- c(0, 2, 4, 6, 6.2, 6.4, 6.6, 6.8, 7, 7.2, 7.4, 7.6, 7.8, 8, 10, 12,
> 14)
> ph2 <- c(12.10, 11.94, 11.68, 11.11, 10.91, 10.74, 10.47, 9.71, 7.1, 4.24,
> 3.3, 3.08, 2.98, 2.86, 2.33, 2.11, 1.98)
> 
> with regards

Ok

temp<-data.frame(ph2,v2)
> drm(temp,, fct = LL.4())

A 'drc' model.

Call:
drm(formula = temp, fct = LL.4())

Coefficients:
b:(Intercept)  c:(Intercept)  d:(Intercept)  e:(Intercept)  
       39.420          2.425         11.487          7.002 

-- 
Bernardo Rangel Tura, M.D,MPH,Ph.D
National Institute of Cardiology
Brazil