Displaying 1 result from an estimated 1 matches for "parval1".
Did you mean:
parval
2013 Feb 19
0
calcMin
...ple.]
> for(i in 1:9) {
+
fitpars=optim(par=parset1,fn=optimwrap,control=list(maxit=100),resid=res
id,caa=caa,na=na,vcode=vcode,maa=maa,ny=ny,nind=nind,qage=qage,selmod=se
lmod,oldagei=oldagei,vpaflag=vpaflag)
+ parset2=fitpars$par
+ parval2=fitpars$value
+ print(parval2)
+ print(parset2)
+ if((parval1-parval2)>1.0) {
+ parval1=parval2
+ parset1=parset2
+ }
+ if((parval1-parval2)<1.0 & parval2>0.1) {
+ parset1=par*pardown[i]
+ }
+ }
[1] 193.1500
[1] 43.84105 43.47576 42.91166 41.81009 42.57930
[1] 15.69411
[1] 13.35845 12.52650 12.28231 11.78171 11.29699
[1] 15.65625
[1] 11.91989 11....