numsamples = 10;

dt = Current Data Table();
statistics = New Table( "statistics.jmp", New Column( "statistic" ) );

For( i = 1, i <= 2, i++,
	resultsTable = New Table( "Tables.jmp",
		New Column( "Sample1", numeric, formula( Random Binomial( 535, 0.0935 ) ) ),
		New Column( "Sample2", numeric, formula( Random Binomial( 588, 0.0935 ) ) ),
		New Column( "Diff", numeric, formula(:Sample1 / 535 - :Sample2 / 588)),
		New Column( "OR", numeric, formula(:Sample1 / (535 - :Sample1) /(:Sample2)*( 588 - :Sample2 )))
	);
	resultsTable << add rows( 1 );
	resultsTable << evalFormula;
	resultsTable << runFormulas;

	difference = Mean(:Diff);
	oddsratio = Mean(:OR);
	
	Wait( 1 );
	Close( resultsTable );
	Current Data Table( statistics );
	statistics << AddRow( {statistic = difference} );
	
);
	Wait( .1 );
For( i = 3, i <= numsamples, i++,
	resultsTable = New Table( "Tables.jmp",
		New Column( "Sample1", numeric, formula( Random Binomial( 535, 0.0935 ) ) ),
		New Column( "Sample2", numeric, formula( Random Binomial( 588, 0.0935 ) ) ),
		New Column( "Diff", numeric, formula(:Sample1 / 535 - :Sample2 / 588)),
		New Column( "OR", numeric, formula(:Sample1 / (535 - :Sample1) /(:Sample2)*( 588 - :Sample2 )))
	);
	resultsTable << add rows( 1 );
	resultsTable << evalFormula;
	resultsTable << runFormulas;
	difference = Col Mean(:Diff);
	oddsratio = Col Mean(:OR);
	Current Data Table( statistics );
	Close( resultsTable );
	Wait (0.1);
	statistics << AddRow( {statistic = difference} );
);