public void SampleSizeTest1()
{
// Example from http://udel.edu/~mcdonald/statttest.html
// Computed using R's function power.t.test
double mean1 = 3.2;
double mean2 = 0;
double var1 = System.Math.Pow(4.3, 2);
double var2 = System.Math.Pow(4.3, 2);
double alpha = 0.05;
double power = 0.80;
TwoSampleTTest test = new TwoSampleTTest(
mean1: mean1, var1: var1, samples1: 10,
mean2: mean2, var2: var2, samples2: 10,
assumeEqualVariances: true, alternate: TwoSampleHypothesis.ValuesAreDifferent);
var target = (TwoSampleTTestPowerAnalysis)test.Analysis.Clone();
target.Power = power;
target.Size = alpha;
target.ComputeSamples(1);
double actual = target.Samples1;
double expected = 29.33682;
Assert.AreEqual(expected, actual, 1e-5);
}
}