Method testHistogramGaussianML

tests/analysis/fit/FitTest.cpp:3831–3885

}

void FitTest::testHistogramGaussianML() {
	Spreadsheet spreadsheet(QStringLiteral("test"), false);
	AsciiFilter filter;

	const QString& fileName = QFINDTESTDATA(QLatin1String("data/Gaussian.dat"));

	auto properties = filter.properties();
	properties.headerEnabled = false;
	filter.setProperties(properties);
	filter.readDataFromFile(fileName, &spreadsheet, AbstractFileFilter::ImportMode::Replace);

	QCOMPARE(spreadsheet.rowCount(), 1000);
	QCOMPARE(spreadsheet.columnCount(), 1);

	Worksheet worksheet(QStringLiteral("test"), false);
	auto* plot = new CartesianPlot(QStringLiteral("plot"));
	worksheet.addChild(plot);

	auto* hist = new Histogram(QStringLiteral("Histogram"));
	plot->addChild(hist);
	hist->setDataColumn(spreadsheet.column(0));

	// Do the fit
	plot->addHistogramFit(hist, nsl_sf_stats_gaussian);

	auto fit = dynamic_cast<XYFitCurve*>(plot->child<XYFitCurve>(0));
	QVERIFY(fit != nullptr);

	QCOMPARE(fit->name(), i18n("Distribution Fit to '%1'", hist->name()));
	// get results
	const XYFitCurve::FitResult& fitResult = fit->fitResult();

	QCOMPARE(fitResult.available, true);
	QCOMPARE(fitResult.valid, true);

	WARN(std::setprecision(15) << fitResult.paramValues.at(0));
	QCOMPARE(fitResult.paramValues.at(0), 210.380459328);
	WARN(std::setprecision(15) << fitResult.paramValues.at(1));
	QCOMPARE(fitResult.paramValues.at(1), 0.999776858937);
	WARN(std::setprecision(15) << fitResult.errorValues.at(1));
	QCOMPARE(fitResult.errorValues.at(1), 0.0223507017166885);
	WARN(std::setprecision(15) << fitResult.paramValues.at(1) - fitResult.marginValues.at(1));
	QCOMPARE(fitResult.paramValues.at(1) - fitResult.marginValues.at(1), 0.955917043549699);
	WARN(std::setprecision(15) << fitResult.paramValues.at(1) + fitResult.marginValues.at(1));
	QCOMPARE(fitResult.paramValues.at(1) + fitResult.marginValues.at(1), 1.0436366743251);

	WARN(std::setprecision(15) << fitResult.paramValues.at(2));
	QCOMPARE(fitResult.paramValues.at(2), -0.0294045302042);
	WARN(std::setprecision(15) << fitResult.errorValues.at(2));
	QCOMPARE(fitResult.errorValues.at(2), 0.0316157202617);
	WARN(std::setprecision(15) << fitResult.paramValues.at(2) - fitResult.marginValues.at(2));
	QCOMPARE(fitResult.paramValues.at(2) - fitResult.marginValues.at(2), -0.0914455198610062);
	WARN(std::setprecision(15) << fitResult.paramValues.at(2) + fitResult.marginValues.at(2));
	QCOMPARE(fitResult.paramValues.at(2) + fitResult.marginValues.at(2), 0.0326364594526431);
}

void FitTest::testHistogramExponentialML() {
	Spreadsheet spreadsheet(QStringLiteral("test"), false);