// This returns data that follow a logit model with the
// supplied parameters. n is the number of observations.
// Note how the binary 0/1 variable is generated so that
// it really follows the logit model

logit_dgp(const n, const theta)
{
	decl k, x, y;
	
	k = sizer(theta); //sizer returns the number of rows of theta
	x = ones(n,1) ~ rann(n,k-1);

	y = (1 ./ (1 + exp(-x*theta)) .> ranu(n,1));
	
	return y ~ x;
	
}