%include data(oralprac); title 'Two-Sample Multivariate Power Analysis'; title2 'Effects of delay on oral practice'; proc glm data=oral outstat=stats; class group; model listen speak read write = group / nouni ss3; manova h = group / short; *include macros(mpower); %mpower(data=stats, out=_out_, yvar=listen speak read write, alpha=.05);