public void TTestPowerAnalysisConstructorTest2()
{
// Example from R's graphical manual
// http://rgm2.lab.nig.ac.jp/RGM2/func.php?rd_id=pwr:pwr.t2n.test
double effectSize = 0.6;
int n1 = 90;
int n2 = 60;
TwoSampleTTestPowerAnalysis target;
double actual, expected;
target = new TwoSampleTTestPowerAnalysis(TwoSampleHypothesis.FirstValueIsGreaterThanSecond)
{
Effect = effectSize,
Samples1 = n1,
Samples2 = n2,
};
target.ComputePower();
expected = 0.9737262;
actual = target.Power;
Assert.AreEqual(expected, actual, 1e-6);
target = new TwoSampleTTestPowerAnalysis(TwoSampleHypothesis.FirstValueIsSmallerThanSecond)
{
Effect = effectSize,
Samples1 = n1,
Samples2 = n2,
};
target.ComputePower();
expected = 0.0;
actual = target.Power;
Assert.AreEqual(expected, actual, 1e-6);
target = new TwoSampleTTestPowerAnalysis(TwoSampleHypothesis.ValuesAreDifferent)
{
Effect = effectSize,
Samples1 = n1,
Samples2 = n2,
};
target.ComputePower();
expected = 0.9470154;
actual = target.Power;
Assert.AreEqual(expected, actual, 1e-6);
}