How can I identify the most predictive factors?
It basically lies within the fitted object. 

model.fit(X, y)
importances = model.feature_importances_


The importances are in the same order as the columns of features ($X$).