A neural network (NN) aided model is proposed for the filtered reaction rate in moderate or intense low-oxygen dilution (MILD) combustion. The framework of the present model is based on the partially stirred reactor (PaSR) approach, and the fraction of the reactive structure appearing in the PaSR is predicted using different NN's, to consider both premixed and non-premixed conditions while allowing the use of imbalanced training data between premixed and non-premixed combustion direct numerical simulation (DNS) data. The key ingredient in the present model is the use of local combustion mode prediction performed by using another NN, which is developed in a previous study. The trained model was then assessed by using two unknown combustion DNS cases, which yields much higher dilution level (more intense MILD condition) and higher Karlovitz number than the DNS cases used as training data. The model performance assessment has been carried out by means of the Pearson's correlation coefficient and mean squared error. For both the present model and zeroth-order approximated reaction rate, the correlation coefficient with the target values shows relatively high values, suggesting that the trend of predicted field, by the present model and zeroth-order approximation, is well correlated with the actual reaction rate field. This suggests that the use of PaSR equation is promising if the fraction of the reactive structure is appropriately predicted, which is the objective in the present study. On the other hand, substantially lower mean squared error is observed for a range of filter sizes for the present model than that for the zeroth-order approximation. This suggests that the present filtered reaction rate model can account for the SGS contribution reasonably well.@en