## Portland cement data set is used.
data(pcd)
k<-0.1650
d<--0.1300
lte2(Y~X1+X2+X3+X4-1,k,d,data=pcd) # Model without the intercept is considered.
## To obtain the variation of MSE of Type (2) Liu Estimator.
data(pcd)
k<-c(0:10/10)
d<-c(-10:10/10)
msemat<-lte2(Y~X1+X2+X3+X4-1,k,d,data=pcd)
matplot(d,lte2(Y~X1+X2+X3+X4-1,k,d,data=pcd),type="l",ylab=c("MSE"),
main=c("Plot of MSE of Type (2) Liu Estimator"),las=1)
text(y=msemat[1,],x=d[1],labels=c(paste0("k=",k)))
## Use "press=TRUE" to obtain the variation of PRESS of Type (2) Liu Estimator.
Run the code above in your browser using DataCamp Workspace