search for: parset2

...l minima. I'm searching for something faster and more intelligent than this example.] > 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....