用pyevolve恢复优化

Question

我已经用Pyevolve做了一个优化，在看了结果之后，我想添加一些代以获得更好的收敛性。由于评估很长，我想知道我是否可以将优化恢复到上一代，并添加20多个代。一切都必须在数据库中设置，我希望这样他才能成为可能。

以下是我的GA属性(类似于第一个示例，但具有更复杂的评估函数):

    # Genome instance, 1D List of 6 elements
genome = G1DList.G1DList(6)
# Sets the range max and min of the 1D List
genome.setParams(rangemin=1, rangemax=15)
# The evaluator function (evaluation function)
genome.evaluator.set(eval_func)
# Genetic Algorithm Instance
ga=GSimpleGA.GSimpleGA(genome)
# Set the Roulette Wheel selector method, the number of generations and
# the termination criteria
ga.selector.set(Selectors.GRouletteWheel)
ga.setGenerations(50)
ga.setPopulationSize(10)
ga.terminationCriteria.set(GSimpleGA.ConvergenceCriteria)
# Sets the DB Adapter, the resetDB flag will make the Adapter recreate
# the database and erase all data every run, you should use this flag
# just in the first time, after the pyevolve.db was created, you can
# omit it.
sqlite_adapter = DBAdapters.DBSQLite(identify="F-Beam-Optimization", resetDB=True)
ga.setDBAdapter(sqlite_adapter)
# Do the evolution, with stats dump
# frequency of 5 generations
ga.evolve(freq_stats=2)

有人有这个想法吗?

Answer 1

嗨，在查看了Pyevolve的文档之后，似乎没有任何方法可以根据您在数据库中存储的内容恢复进化(奇怪的行为)。

如果你想实现这种类型的机制，你可以考虑偶尔pickle你的人口，并在Pyevolve中实现整个事情。

或者，您可以尝试DEAP，这是一个非常开放的框架，可以让您透明地查看和操作进化算法的各个方面。并且已经实现了一些检查点机制。

下面是你的代码在DEAP中的样子。

import random    
from deap import algorithms, base, creator, tools
# Create the needed types
creator.create("FitnessMax", base.Fitness, weights=(1.0,))
creator.create("Individual", list, fitness=creator.FitnessMax)
# Container for the evolutionary tools
toolbox = base.Toolbox()
toolbox.register("attr", random.random, 1, 15)
toolbox.register("individual", tools.initRepeat, creator.Individual, toolbox.attr, 6)
toolbox.register("population", tools.initRepeat, list, toolbox.individual)
# Operator registering
toolbox.register("evaluate", eval_func)
toolbox.register("mate", tools.cxTwoPoints)
toolbox.register("mutate", tools.mutGaussian, mu=0, sigma=1, indpb=0.05)
toolbox.register("select", tools.selTournament, tournsize=3)
population = toolbox.population(n=10)
stats = tools.Statistics(key=lambda ind: ind.fitness.values)
stats.register("Max", max)
checkpoint = tools.Checkpoint(population=population)
GEN, CXPB, MUTPB = 0, 0.5, 0.1
while stats.Max() < CONDITION:
    # Apply standard variation (crossover followed by mutation)
    offspring = algorithms.varSimple(toolbox, population, cxpb=CXPB, mutpb=MUTPB)
    # Evaluate the individuals
    fits = toolbox.map(toolbox.evaluate, offspring)
    for fit, ind in zip(fits, offspring):
        ind.fitness.values = fit
    # Select the fittest individuals
    offspring = [toolbox.clone(ind) for ind in toolbox.select(offspring, len(offspring)]
    # The "[:]" is important to not replace the label but what it contains
    population[:] = offspring
    stats.update(population)
    if GEN % 20 == 0:
        checkpoint.dump("my_checkpoint")
    GEN += 1

注意上面的代码没有经过测试。但它能满足你所有的要求。现在来看看如何加载检查点并重新启动进化。

checkpoint = tools.Checkpoint()
checkpoint.load("my_checkpoint.ems")
population = checkpoint["population"]
# Continue the evolution has in before

此外，DEAP有很好的文档，并且有超过25个不同的示例，可以帮助新用户快速入门，我还听说开发人员回答问题的速度非常快。