def pene(geno, age, ind):
    """Define an age-dependent penetrance function"""
    if age < 16:
        return 0
    if ind.affected():
        return 1
    return (0.0, 0.001 * age, 0.001 * age)[sum(geno)]
