import flavio
from wilson import Wilson
with open('Scotogenic.yaml', 'r') as f:
    my_w = Wilson.load_wc(f)
print(flavio.np_prediction('BR(mu->egamma)', my_w))