## Portland cement data set is used.
data(pcd)
d<-0.05
aul(Y~X1+X2+X3+X4-1,d,data=pcd) # Model without the intercept is considered.
## To obtain the variation of MSE of Almost Unbiased Liu Estimator.
data(pcd)
d<-c(1:10/10)
plot(aul(Y~X1+X2+X3+X4-1,d,data=pcd),
main=c("Plot of MSE of Almost Unbiased Liu Estimator"),las=1)
Run the code above in your browser using DataLab