[R-es] Curva dosis-respuesta

[Marcuzzi, Javier Rubén]

Estimado Jorge

No sabría la solución a su problema, leí su correo para aprender de usted, 
pensé en aumentar las interacciones para llegar a una convergencia válida, 
de la ayuda vi tres puntos, dos son palabras y una referencia. No se si le 
sirve, pero ...

niter:  the number of iterations that nlm used.

iterlim

http://cran.r-project.org/web/views/ChemPhys.html

Javier

-----Mensaje original----- 
From: Jorge I Velez
Sent: Tuesday, May 15, 2012 12:42 PM
To: R-help-es
Subject: [R-es] Curva dosis-respuesta

Buenos dias R-help-es,

Estoy interesado en estimar una curva dosis-respuesta para un conjunto de
datos y para ello, estoy utilizando la libreria "drm".  Hasta ahi todo bien.

Me gustaria automatizar algunas cosas y el primer paso para ello es la
estimacion del modelo.  Si la estimacion funciona, todo lo demas funciona;
de lo contrario, todo fallara.  Tengo algunas lineas que mitigan un poco el
problema, pero no cuando hay problemas de convergencia como el que se
ilustra a continuacion:

# 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('drm')
require(drm, 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

He cambiado un poco los parametros de control de la funcion drm, pero no ha
sido efectivo. Alguna otra sugerencia?.

Muchas gracias!
Jorge.-

[[alternative HTML version deleted]]

_______________________________________________
R-help-es mailing list
R-help-es en r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es