Logo Questions Linux Laravel Mysql Ubuntu Git Menu
 

R- Optimx for exponential function with 2 parameters - cannot evaluate function at initial parameter values

I feel like I missed something very obvious but after an hour of fiddling/googling I cannot get this to work. Code:

#Generate data from exponential model

xdata<-seq_len(100)
ydata<-2*exp(-2*(xdata+rnorm(100)))

#Fit exponential model to data
firstorder<-function(C0,k){
 ynew<-C0*exp(-k*xdata)
 RMSE<-sum((ynew-ydata)^2,na.rm=TRUE)
 return(RMSE)
}

#Initial parameter values
params<-c(1,1)

#Optimize
optimx(params,firstorder)

Error in optimx.check(par, optcfg$ufn, optcfg$ugr, optcfg$uhess, lower, : Cannot evaluate function at initial parameters

I tried a variety of ways to input the parameters.

like image 653
Pinemangoes Avatar asked Mar 04 '14 15:03

Pinemangoes


1 Answers

Try

optimx(params, function(x) firstorder(x[1], x[2]))
like image 89
G. Grothendieck Avatar answered Oct 19 '22 05:10

G. Grothendieck