Hi Wondering if anyone knows how to run mixed-effect Tobit model? Normally in cross-sectional data would use "proc qlim" but not sure if I can use this code in longitudinal data. My data consist of 4 fixed effects {x1, x2, x3, x4},  random effect intercept z0 and time varying random effect  z1 . Example of code using proc mixed if this helps. proc mixed data=data method=ml; class gender treatment; model y= x1 x2 x3 x4 /s; RANDOM intercept time /s; run; If possible I would like to avoid Bayesian as my knowledge is minimal in that field. Many thank, ... View more

Hi all I was wondering if someone can help me create the following macro. I am given dataset  call it "dat" with 100,000 observations with variables a b c d e, where b and c are categorical variables. 1) I want to split the data into groups each containing 5000 observations. 2) Run proc reg on each individually created dataset. 3) Create a table containing all of the parameter estimates. This is how far I have managed to get; %macro split(ndat=20); data %do i = 1 %to &ndat.; dat&i. %end; ; retain x; set dat nobs=nobs; if _n_ eq 1 then do; if mod(nobs,&ndat.) eq 0 then x=int(nobs/&ndat.); else x=int(nobs/&ndat.)+1; end; if _n_ le x then output dat1; %do i = 2 %to &ndat.; else if _n_ le (&i.*x) then output dat&i.; %end; run; %mend split; %split(ndat=20) ; /*creates 20 dat-sets I want to somehow implete the code below for each dataset*/ %macro test(data=); ods output ParameterEstimates=parm; /*this creates parameter estimates I want to combine into one table containing all 20 set of parameter estimates*/ proc reg data= &data; class b c ; model e= a b c d/ solution ss3; run; %mend test; ... View more