Differential D22163 Diff 60856 tests/stats/ttest/TTestTest.cpp

Changeset View

Standalone View

View Options

tests/stats/ttest/TTestTest.cpp

  • This file was added.
1​///***************************************************************************
2​// File : CorrelationTest.cpp
3​// Project : LabPlot
4​// Description : Tests for data correlation
5​// --------------------------------------------------------------------
6​// Copyright : (C) 2019 Devanshu Agarwal (agarwaldevanshu8@gmail.com)
7​// ***************************************************************************/
8​
9​///***************************************************************************
10​// * *
11​// * This program is free software; you can redistribute it and/or modify *
12​// * it under the terms of the GNU General Public License as published by *
13​// * the Free Software Foundation; either version 2 of the License, or *
14​// * (at your option) any later version. *
15​// * *
16​// * This program is distributed in the hope that it will be useful, *
17​// * but WITHOUT ANY WARRANTY; without even the implied warranty of *
18​// * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the *
19​// * GNU General Public License for more details. *
20​// * *
21​// * You should have received a copy of the GNU General Public License *
22​// * along with this program; if not, write to the Free Software *
23​// * Foundation, Inc., 51 Franklin Street, Fifth Floor, *
24​// * Boston, MA 02110-1301 USA *
25​// * *
26​// ***************************************************************************/
27​
28​#include "TTestTest.h"
29​#include "backend/hypothesisTest/HypothesisTest.h"
30​
31​#include "backend/core/AbstractColumn.h"
32​#include "backend/core/column/Column.h"
33​
34​void TTestTest::twoSampleIndependent() {
35​ QVector<double> col1Data = {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1};
36​ QVector<double> col2Data = {3, 1, 5, 4, 6, 4, 6, 2, 0, 5, 4, 5, 4, 3, 6, 6, 8, 5, 5, 4, 2, 5, 7, 5};
37​
38​ Column* col1 = new Column("col1", AbstractColumn::Numeric);
39​ Column* col2 = new Column("col2", AbstractColumn::Numeric);
40​
41​ col1->replaceValues(0, col1Data);
sgerlachUnsubmitted
Done

add comment: "first sample"

sgerlach: add comment: "first sample"
42​ col2->replaceValues(0, col2Data);
43​
44​ QVector<Column*> cols;
45​ cols << col1 << col2;
46​
47​ HypothesisTest tTest("Test");
48​ tTest.setColumns(cols);
49​
sgerlachUnsubmitted
Done

add comment "second sample"

sgerlach: add comment "second sample"
50​ HypothesisTest::Test test;
51​ test.type = HypothesisTest::Test::Type::TTest;
52​ test.subtype = HypothesisTest::Test::SubType::TwoSampleIndependent;
53​ test.tail = HypothesisTest::Test::Tail::Two;
54​
sgerlachUnsubmitted
Done

add space before

sgerlach: add space before
devanshuagarwalAuthorUnsubmitted
Done

Sorry, I didn't get. Add space before?

devanshuagarwal: Sorry, I didn't get. Add space before?
sgerlachUnsubmitted
Done

I meant to add an empty line to make it more readable.

sgerlach: I meant to add an empty line to make it more readable.
55​ bool categoricalVariable = true;
56​ bool equalVariance = true;
57​
58​ tTest.performTest(test, categoricalVariable, equalVariance);
59​ double tValue = tTest.statisticValue();
60​ double pValue = tTest.pValue();
61​
62​ qDebug() << "tValue is " << tValue;
63​ qDebug() << "pValue is: " << pValue;
64​
65​ FuzzyCompare(tValue, -1.713, 1.e-2);
66​ FuzzyCompare(pValue, 0.101, 1.e-2);
67​}
68​
69​
70​void TTestTest::twoSampleIndependent_data() {
71​}
72​
73​QTEST_MAIN(TTestTest)
sgerlachUnsubmitted
Done

Can you use the DEBUG() macro instead of qDebug()?
QDEBUG("tValue is " << tValue)

sgerlach: Can you use the DEBUG() macro instead of qDebug()? QDEBUG("tValue is " << tValue)
devanshuagarwalAuthorUnsubmitted
Done

QDEBUG and DEBUG macro are not working in my pc, may be due to some setting. I have to check for it. I will use the macro, once I will fix the setting.

devanshuagarwal: QDEBUG and DEBUG macro are not working in my pc, may be due to some setting. I have to check…
sgerlachUnsubmitted
Done

I thought that you already fixed it for you. If it's not working for you, please add a TODO-comment when using qDebug().

sgerlach: I thought that you already fixed it for you. If it's not working for you, please add a TODO…
sgerlachUnsubmitted
Done

using an absolute error of 0.01 depends heavily on tValue. Why is the error so big? Can you use relative errors instead?

sgerlach: using an absolute error of 0.01 depends heavily on tValue. Why is the error so big? Can you use…
devanshuagarwalAuthorUnsubmitted
Done

I tried using relative error. But, I am not getting what value of delta to use. If I use as small as 0.1 also then the test is not passing for values:

actual: 1.51815e-05
expected < 0.001

which is very more accurate than expected.

devanshuagarwal: I tried using relative error. But, I am not getting what value of delta to use. If I use as…
sgerlachUnsubmitted
Done

the problem here is that the expected value is not known more precise? An absolute error of 0.01 make not much sense here.
I would suggest to use a more precise expected value or test for actual < expected in this case.

sgerlach: the problem here is that the expected value is not known more precise? An absolute error of 0.
sgerlachUnsubmitted
Done

see above

sgerlach: see above
sgerlachUnsubmitted
Done

DEBUG()?

sgerlach: DEBUG()?
sgerlachUnsubmitted
Done

see above

sgerlach: see above
sgerlachUnsubmitted
Done

DEBUG()?

sgerlach: DEBUG()?
sgerlachUnsubmitted
Done

check comment

sgerlach: check comment