R回归函数的上限和下限

时间:2017-09-07 22:28:07

标签: r lower-bound upperbound

我试图在函数

中绑定以下变量
lower = c(Hyp.b = 0.01, Hyp.Di = .0000001),
upper = c(Hyp.b = 1.01, Hyp.Di = .1)

脚本如下:

Hyp.q.forward.fun = function( time, Hyp.qi, Hyp.b, Hyp.Di ){ # time in days
Hyp.q.theo = Hyp.qi*(1 + Hyp.b*Hyp.Di*time)^(-1/Hyp.b)
return(Hyp.q.theo)
}


residfun = function(x,x.days,y.prod){
 Hyp.qi = x[1]
 Hyp.b  = x[2]
 Hyp.Di = x[3]
 q.theo = (365.25/12)*Hyp.q.forward.fun(
   time=x.days, 
   Hyp.qi=Hyp.qi, 
   Hyp.b=Hyp.b, 
   Hyp.Di=Hyp.Di)
 #plot(x.days,y.prod); lines(x.days,q.theo)
 residual = sqrt(sum((q.theo-y.prod)^2))
 return(residual)
 }

我不确定我是否使用正确的方法来绑定这两个变量。任何帮助将不胜感激。

0 个答案:

没有答案