[R-es] Curva dosis-respuesta

[Olivier Nuñez]

Jorge,

la curva dosis-respuesta que sugiere tus datos tiene un patrón  
bastante atípico.
Típicamente la respuesta (o efecto) se estabiliza a partir de un  
cierto umbral (tipo sigmoide).
De ahí posiblemente las dificultades del modelo para ajustarse a tus  
datos.

Un saludo. Olivier

--  
____________________________________

Olivier G. Nuñez
Email: onunez en iberstat.es
Tel : +34 663 03 69 09
Web: http://www.iberstat.es

____________________________________




El 15/05/2012, a las 21:41, Jorge I Velez escribió:

> # datos
> d <- structure(list(y = c(509.75, 498.8629, 584.9069, 808.5421,  
> 679.3904,
> 597.3327, 561.4653, 496.4428, 474.4776), dose = c(2.27726728500976,
> 1.65500383670491, 2.93012117741345, 7.29822890763926,  
> 7.77257337407207,
> 8.42489485944128, 9.1495092419634, 12.1197337185771, 11.7126676116462
> )), .Names = c("y", "dose"), row.names = c(NA, -9L), class =  
> "data.frame")
>
> # curva
> # install.packages('drc')
> require(drc, quietly = TRUE, warn.conflicts = FALSE)
> with(d, plot(dose, y, las = 1))
> fit <- drm(y ~ dose, fct = LL.4(), control = drmc(maxIt = 100000,  
> relTol =
> 1e-4), data = d)
> # Error in optim(startVec, opfct, hessian = TRUE, method = optMethod,
> control = list(maxit = maxIt,  :
> #  non-finite finite-difference value [4]
> # Error in drmOpt(opfct, opdfct1, startVecSc, optMethod, constrained,
> warnVal,  :
> #  Convergence failed