## Portland cement data set is used.
data(pcd)
d<-0.05
r<-c(2.1930,1.1533,0.75850)
R<-c(1,0,0,0,0,1,0,0,0,0,1,0)
dpn<-c(0.0439,0.0029,0.0325)
delt<-c(0,0,0)
srliu(Y~X1+X2+X3+X4-1,r,R,dpn,delt,d,data=pcd)
# Model without the intercept is considered.
## To obtain the variation of MSE of Stochastic Restricted Liu Estimator.
data(pcd)
d<-c(0:100/100)
r<-c(2.1930,1.1533,0.75850)
R<-c(1,0,0,0,0,1,0,0,0,0,1,0)
dpn<-c(0.0439,0.0029,0.0325)
delt<-c(0,0,0)
plot(srliu(Y~X1+X2+X3+X4-1,r,R,dpn,delt,d,data=pcd),
main=c("Plot of MSE of Stochastic Restricted Liu Estimator"),las=1)
Run the code above in your browser using DataCamp Workspace