Displaying 2 results from an estimated 2 matches for "dlbar".
Did you mean:
lbar
2007 Sep 09
2
What does it mean by "initial value not available"?
....0)))))))))))
fit.error=function(rec,lambda,lbar)
{sum((eval(apple)*1000-orange)^2/(orange^2))
}
fit.error.grr=function(rec,lambda, lbar)
{drec=sum(eval(D(apple,'rec'))*(eval(apple)*1000-orange)/(orange^2))
dlambda=sum(eval(D(apple,'lambda'))*(eval(apple)*1000-orange)/(orange^2))
dlbar=sum(eval(D(apple,'lbar'))*(eval(apple)*1000-orange)/(orange^2))
c(drec,dlambda,dlbar)
}
rr=c(4.33,4.22,4.27,4.43,4.43,4.44,4.45,4.65,4.77,4.77)
tot=rep(13319.17,10)
sh=rep(1553656,10)
sigmae=c(0.172239074,0.188209271,0.193703774,0.172659891,0.164427247,
0.24602361,0.173555309,0.186701165,...
2007 Sep 10
2
Are the error messages of ConstrOptim() consisten with each other?
...(eval(apple2)*1000-orange)^2/(orange^2))
}
fit.error.grr=function(rec1,lambda1,lbar1)
{rec=rec1
lambda=lambda1
lbar=lbar1
drec=sum(20000*eval(D(apple2,'rec'))*(eval(apple2)*10000-orange)/(orange^2))
dlambda=sum(20000*eval(D(apple2,'lambda'))*(eval(apple2)*10000-orange)/(orange^2))
dlbar=sum(20000*eval(D(apple2,'lbar'))*(eval(apple2)*10000-orange)/(orange^2))
c(drec,dlambda,dlbar)
}
ui=matrix(c(1,-1,0,0,0,0,0,0,1,-1,0,0,0,0,0,0,1,-1),6,3)
ci=c(0,-0.5,0,-2,0,-0.6)
constrOptim(c(0.5,0.3,0.5), f=fit.error, gr=fit.error.grr, ui=ui,ci=ci)
constrOptim(c(0.5,0.9,0.5), f=fit.err...