R help es - Jul 2009 - Problema con argumento en funcion nls

Hola,

Estamos escribiendo un wrapper para la función nls de la librería stats.
Tenemos un problema con uno de los argumentos weightsArgument que parece que
no lo coge la funcion nls.

Adjunto una versión del código que se puede ejecutar.

library(stats)
wrappernls <-
function(modelArgument,dataArgument,startArgument,weightsArgument)
{
F.mod <-try(nls(formula=modelArgument,data=dataArgument,start=startArgument,
weights=weightsArgument,
control=nls.control(maxiter = 500, tol = 1e-05, minFactor
1/102400)),silent=TRUE)
return(F.mod)
}

F<-c(0.4091867,0.4060938,0.4032078,0.4089090,0.4138126,0.4183426,0.4073004,0.4145457,0.4137699,
0.4161127,0.4228770,0.4231176,0.4295189,0.4290417,0.4348761,0.4517475,0.4899147,0.5463731,
0.6273890,0.7458752,0.8960531,1.0280455,1.1753147,1.3122100,1.4370375,1.5610782,1.6666418,
1.7771591,1.8853868,1.9628042,2.0407750,2.1179830,2.1647203,2.2147800,2.2848194,2.3442066,
2.3706858,2.4121310,2.4508073,2.4767710)
skipCycles=3;

y00<-max(F);
a0<- min(F)-y00;
x00<-which.min(abs(F-(a0/2+y00)));
derivativeAtc0= (F[x00+1]-F[x00-1])/2
b0=-4*x00*derivativeAtc0/a0
w<-c(rep(0,skipCycles), rep(1,length(F)-skipCycles))
x=1:length(F);

modelArg=F ~ y0+D*x+a/(1+(x/x0)^b)
dataArg=data.frame(F,x)
startArg=list(y0 = y00, a=a0, x0=x00, b=b0, D=0)
weightsArg=w

F.modWrapper <- wrappernls(modelArg,dataArg,startArg,weightsArg)
F.mod <-
nls(formula=modelArg,data=dataArg,start=startArg,weights=weightsArg,
control=nls.control(maxiter = 500, tol = 1e-05, minFactor = 1/102400))

*Y la salida es:*
> F.mod
Nonlinear regression model
model: F ~ y0 + D * x + a/(1 + (x/x0)^b)
data: dataArg
y0 a x0 b D
2.891357 -2.434011 25.604522 5.880144 -0.007322
weighted residual sum-of-squares: 0.02664

Number of iterations to convergence: 8
Achieved convergence tolerance: 3.517e-06
> F.modWrapper
[1] "Error in eval(expr, envir, enclos) : object
''weightsArgument'' not
found\n"
attr(,"class")
[1] "try-error"




Alguna idea de porqué pasa esto?
Gracias.

	[[alternative HTML version deleted]]

Hola, ¿qué tal?

Lee la ayuda de ?nls y verás que trae algunos comentarios sobre dónde
busca la función nls (algunos de) sus argumentos:

Variables in formula (and weights if not missing) are looked for first
in data, then the environment of formula and finally along the search
path. Functions in formula are searched for first in the environment
of formula and then along the search path.

En tu ejemplo, nls no busca su argumento weightsArgument "dentro" de
la función wrappernls sino fuera de ella, que es donde has definido la
fórmula.

Se me ocurren varias opciones:

1) Meter el peso dentro de los datos como columna adicional.

2) Recrear la fórmula dentro de la función para asignarle un
"environment" local.

3) Preguntarle a la fórmula cuál es su "environment" y acudir a él a
buscar los pesos.

Personalmente, trataría de implementar (2). Podrías hacerlo (y seguro
que hay formas mejores) definiendo, dentro de tu función

nueva.formula <- formula( deparse( formula.argumento ) )

Un saludo,

Carlos J. Gil Bellosta