Hi:
 
Is there a command in R (similar to hblm I guess) that would
allow me to fit a clustered version of the hierarchichal
bayes model, 
 
Y_ij = X_ij\beta + \delta_i + \epsilon_ij
 
the \epsilon_ij have known variances,
\delta_i ~ N(0,\tau^2)
etc..
 
Thanks!
 
 


		
---------------------------------


	[[alternative HTML version deleted]]


