#include <oxstd.h>

#import "../code/sim_autometrics"

main()
{
	decl cm = 1000;
	decl dgp = <2,7,8,9>, truesize = <1,3,3,5>;

	format(1000);
	
	for (decl i = 0; i < sizerc(dgp); ++i)
	{
		decl exp = new AutometricsExp("HP", dgp[i]);
		exp.Stepwise("Lasso", 0, -1, 3);	   // BIC
		exp.Run(cm);
		delete exp;
	}
	for (decl i = 0; i < sizerc(dgp); ++i)
	{
		decl exp = new AutometricsExp("HP", dgp[i]);
		exp.Stepwise("Lasso", 0, truesize[i], -1);	   // true size
		exp.Run(cm);
		delete exp;
	}
	for (decl i = 0; i < sizerc(dgp); ++i)
	{
		decl exp = new AutometricsExp("HP", dgp[i]);
		exp.Stepwise("Auto_size", 0, truesize[i], -1);	   // true size
		exp.Run(cm);
		delete exp;
	}
}