In SAS, we can fit a two-way random effect model for pooled time series cross-sectional data using procedures PANEL or TSCSREG with the option /RANTWO.For an example, see http://rt.uits.iu.edu/visualization/analytics/docs/panel-docs/panel7.php

The two-way random effect model is specified as:y_it = B_0 + B_1*X1_it + u_i + v_t + e_it(i is group identifier, t is time identifier)

Can we fit such a model in other software packages such as HLM, SPSS, R, etc? If yes, how?
I