我觉得我错过了一些非常明显的东西,但经过一个小时的摆弄/谷歌搜索后,我无法让它发挥作用。代码:
#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)
optimx.check出错(par,optcfg $ ufn,optcfg $ ugr,optcfg $ uhess,lower,: 无法在初始参数下评估函数
我尝试了多种输入参数的方法。
答案 0 :(得分:8)
尝试
optimx(params, function(x) firstorder(x[1], x[2]))