コード例 #1
0
def test_convex4():

    logger = getLogger('test_convex4')

    logger.info("Test  convex Poly(4)-Distribution")
    model = Models.Convex_pdf('C4', x, power=4, increasing=False, convex=True)

    result, f = model.fitTo(dataset, silent=True, refit=5)
    with wait(1), use_canvas('test_convex4'):
        model.draw(dataset)

    title = 'Positive decreasing convex polynomial'
    logger.info('%s \n%s' % (title, result.table(title=title, prefix='# ')))

    models.add(model)
コード例 #2
0
def test_convex4():

    logger.info("Test  convex Poly(4)-Distribution")
    model = Models.Convex_pdf('C4', x, power=4, increasing=False, convex=True)

    result, f = model.fitTo(dataset, silent=True)

    if 0 != result.status() or 3 != result.covQual():
        logger.warning('Fit is not perfect MIGRAD=%d QUAL=%d ' %
                       (result.status(), result.covQual()))
        print(result)
    else:
        for phi in model.phis:
            logger.info("\tConvex4:      phi= %-17s " % phi.ve())

    models.add(model)