xx = seq(-15,+15,length=1000) # the range of beta's to sample 

# the specific parameters:
betaHat = -5.0
lambda  = 1.0

postscript("../../WriteUp/Graphics/Chapter3/b_negative_ge_lambda.eps", onefile=FALSE, horizontal=FALSE)
plot( xx, ( xx - betaHat )^2 + lambda*abs(xx) )
dev.off()



betaHat = -1.0
lambda  = 5.0

postscript("../../WriteUp/Graphics/Chapter3/b_negative_le_lambda.eps", onefile=FALSE, horizontal=FALSE)
plot( xx, ( xx - betaHat )^2 + lambda*abs(xx) )
dev.off()



betaHat = 5.0
lambda  = 1.0

postscript("../../WriteUp/Graphics/Chapter3/b_positive_ge_lambda.eps", onefile=FALSE, horizontal=FALSE)
plot( xx, ( xx - betaHat )^2 + lambda*abs(xx) )
dev.off() 




betaHat = 1.0
lambda  = 5.0

postscript("../../WriteUp/Graphics/Chapter3/b_positive_le_lambda.eps", onefile=FALSE, horizontal=FALSE)
plot( xx, ( xx - betaHat )^2 + lambda*abs(xx) )
dev.off()