from pgmpy.inference import VariableElimination
from pgmpy.models import BayesianModel
import numpy as np
import pandas as pd
 values = pd.DataFrame(np.random.randint(low=0, high=2, size=(1000, 5)),
                      columns=['A', 'B', 'C', 'D', 'E'])
model = BayesianModel([('A', 'B'), ('C', 'B'), ('C', 'D'), ('B', 'E')])
model.fit(values)
inference = VariableElimination(model)

inference.query(variables=None)
Out[5]: 
{<Factor representing phi(D:2, C:2) at 0x7f2c3a07abe0>,
 <Factor representing phi(B:2, A:2, C:2) at 0x7f2c3a07ae48>,
 <Factor representing phi(A:2) at 0x7f2c3a07af98>,
 <Factor representing phi(C:2) at 0x7f2c3a07ab00>,
 <Factor representing phi(E:2, B:2) at 0x7f2c3a07af60>}

inference.query(variables=['A'])
Out[6]: {'A': <Factor representing phi(A:2) at 0x7f2c58718400>}