mdl1<-lme(y ~ x1,random=(~1|state),data=data,method="ML")
mdl2<-lme(y ~ x1 + x2,random=(~1|state),data=data,method="ML")
mdl3<-lme(y ~ x2,random=(~1|state),data=data,method="ML")
mdl4<-lme(y ~ x3,random=(~1|state),data=data,method="ML")
mdl5<-lme(y ~ x3 + x1,random=(~1|state),data=data,method="ML")
.
.
.
.
mdlx<-lme(y ~ xi,random=(~1|state),data=data,method="ML")