## Portland cement data set is used.
data(pcd)
k<-0.1650
d<--0.1300
aa<-c(0.958451,1.021155,0.857821,1.040296)
alte1(Y~X1+X2+X3+X4-1,k,d,aa,data=pcd) # Model without the intercept is considered.
## To obtain the variation of MSE of Type (1) Adjusted Liu Estimator.
data(pcd)
k<-c(0:10/10)
d<-c(-10:10/10)
aa<-c(0.958451,1.021155,0.857821,1.040296)
msemat<-alte1(Y~X1+X2+X3+X4-1,k,d,aa,data=pcd)
matplot(d,alte1(Y~X1+X2+X3+X4-1,k,d,aa,data=pcd),type="l",ylab=c("MSE"),
main=c("Plot of MSE of Type (1) Adjusted 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 (1) Adjusted Liu Estimator.
Run the code above in your browser using DataCamp Workspace