Hi All,
Recently I'm doing my thesis with LHS, and I created the design by using the
folloing code.
"require(DiceDesign)
require(DoE.wrapper)
set.seed(27662)
LHD=lhs.design(200,7,type="dmax",factor.names=list(
Angle1=c(80,90),Angle2=c(90,100),
Angle3=c(32,40),Radial2=c(10,15),
Radial3=c(25,35),Length1=c(120,130),
Length2=c(0,10)),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‍"
As R language is new for me, Each time I run the code, the different LHDopt can
get. I don't know how can I avoid this situation.
Anybody could help me?
BR,
Dean
[[alternative HTML version deleted]]
______________________________________________
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
