hip<-read.table("hip.txt",col.names=c("y","age","sex","subj","time")) lmm.gibbs<-function (its=1000,p=5,beta.s=rep(0,p),beta.m=rep(0,p),beta.p=diag(rep(0.000001,p)),prec.s=1,prec.sh=0.001,prec.sc=0.001,precb.s=1,precb.sh=0.001,precb.sc=1) { y<-hip$y n<-length(y) subj<-unique(hip$subj) k<-length(subj) x<-cbind(rep(1,n),hip$age,hip$sex,ifelse(hip$time==2,1,0),ifelse(hip$time==3,1,0)) z<-NULL for (i in subj) z<-cbind(z,ifelse(hip$subj==i,1,0)) beta<-beta.s prec<-prec.s precb<-precb.s betasave<-NULL bsave<-NULL sigmasave<-NULL i<-0 while (i