public void TTestPowerAnalysisConstructorTest4()
{
// Example from http://www.ats.ucla.edu/stat/stata/dae/t_test_power2.htm,
// tested against G*Power results
double meanA = 0;
double meanB = 10;
double sdA = 15;
double sdB = 17;
double varA = sdA * sdA;
double varB = sdB * sdB;
{
var priori = TwoSampleTTestPowerAnalysis.GetSampleSize(10,
variance1: varA, variance2: varB, power: 0.8);
Assert.AreEqual(41, Math.Truncate(priori.Samples1));
Assert.AreEqual(41, Math.Truncate(priori.Samples2));
}
{
TwoSampleTTest test = new TwoSampleTTest(
meanA, varA, 30,
meanB, varB, 30);
Assert.AreEqual(0.661222, test.Analysis.Power, 1e-6);
}
{
TwoSampleTTest test = new TwoSampleTTest(
meanA, varA, 20,
meanB, varB, 40);
Assert.AreEqual(0.6102516, test.Analysis.Power, 1e-6);
}
{
var priori = TwoSampleTTestPowerAnalysis.GetSampleSize(10,
variance1: varA, variance2: varB, power: 0.8, alpha: 0.07);
Assert.AreEqual(37, Math.Truncate(priori.Samples1));
Assert.AreEqual(37, Math.Truncate(priori.Samples2));
}
}