## Not run: ------------------------------------
# dev.new()
# harray=seq(0.05,1.5,len=1000)
# ISEarray=ISE_fstar(harray,100)
# h_ISE=round(harray[which.min(ISEarray)],digits=4)
# dev.new()
# plot(harray,ISEarray,lwd=3,'l',xlab="h",ylab="ISE",main="ISE(h)",cex.main=2,cex.lab=1.7,
# cex.axis=1.7)
# legend(0.35,ISEarray[5],legend=c("n=100",paste("h_ISE=",h_ISE)),cex=1.8,bty="n")
## ---------------------------------------------
Run the code above in your browser using DataLab