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.
Try
optimx(params, function(x) firstorder(x[1], x[2]))
If you love us? You can donate to us via Paypal or buy me a coffee so we can maintain and grow! Thank you!
Donate Us With