search for: rmultireg

...gorithm + # + Data=list(m=nalt,X=X[rowx1[i]:rowx2[i],],y=Y[rowy1[i]:rowy2[i],]) + bbar=t(theta)%*%z[i,] + prior=list(A=delta,betabar=bbar) + a=rmnlRwMetrop1(Data,prior,Mcmc=parameters,beta[i,]) + beta[i,]=a$betadraw + accept[i]=a$acceptr + } + mregout=rmultireg(beta,z,thetabar,A1,nu1,V) + delta=mregout$Sigma + theta=mregout$B + if (j%%keep==0){ + thetas[j/keep,]=as.vector(theta) + betastore[j/keep]=list(beta) + accepts[j/keep]=mean(accept) + ftime=proc.time()[3] + cat('Time taken by iteration: ',j, ' = ',round((ftime-itime)/60,2),'\n...