代码如下:
require(DiceDesign)
require(DoE.wrapper)
set.seed(27662)
LHD=lhs.design(200,3,type="dmax",factor.names=list(Angle1=c(80,90),Angle2=c(90,100),Angle3=c(32,40)),
range=0.98,niter_max=2000)
LHDnew=maximinESE_LHS(LHD,T0=0.005*phiP(LHD,p=50),inner_it=200,J=50,it=4)
plot(LHDnew$design$Angle1,LHDnew$design$Angle2)
plot(LHDnew$critValues,type="l")
LHDopt=LHDnew$design
运行上述代码时,每次运行得到的LHDopt的结果都不一样,如何才能确保运行得出的LHDopt不变?