Recent reactor neutrino experiments have shown anomalous results in both the reactor flux and spectrum measurements. Compared the measurements, Reactor neutrino flux shows a 6% deficit while reactor netrino spectrum illustrates a bump-like structure at around 5 MeV region. In this presentation we will employ the methods of both theoretical model prediction and global neutrino data analysis to expore all the possible solutions to these reactor anomalies. We will discuss drawbacks of the model predictions within the standard model, and the possible new physics solution beyond the standard model.