%include data(baseball);

%missing(data=baseball);

*-- Transform variables;
data baseball;
	set baseball;
	if salary ^=. 
		then logsal = log(salary);
	years7 = min(years,7);
	trpc = (runsc + rbic + homerc) / years;
	trp = runs + rbi + homer;
	label logsal = 'log Salary'
		trpc='Total career runs/year'
		trp ='Total runs 1986'
		years7='Years, up to 7';

%missing(data=baseball, var=logsal years7 trpc batavgc,
	print=pmiss n nmiss);

title 'Original data';
proc reg data=baseball outest=stato;
	model logsal = years7 trpc batavgc ;

%miss(data=baseball, var=logsal years7 trpc batavgc,
	out=impute, danum=4);

title 'Imputed data';
proc reg data=impute outest=stats covout noprint;
	model logsal = years7 trpc batavgc;
	by dsnum;


%misscomb(data=stats, regdep=logsal);

data stats;
	set stato(in=ino) stats;
	if ino then dsnum=0;

proc print data=stats;
	where (_type_='PARMS');
	id dsnum;
	var INTERCEP YEARS7 TRPC BATAVGC _rmse_;
	format dsnum 5.0 intercep 6.3 years7 trpc batavgc _rmse_ 6.4;