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.
See Question&Answers more detail:
os 与恶龙缠斗过久,自身亦成为恶龙;凝视深渊过久,深渊将回以凝视…