In recent years, modeling neurons and neuronal collections with high accuracy have become central issues of neuroscience. The development of efficient algorithms for their simulation as well as the increase in computational power and parallelization need to keep up with the quantity and complexity of novel recordings and reconstructions reported by the experimental neuroscientists. The extraction of low-order equivalents that capture the essential aspects of the high-accuracy models is an essential part of the simulation process. The complexity of these models require the use of black-box data-oriented reduction approaches. We create a detailed model of the nervous system of a very known organism, C. Elegans, and show that it can be reduced using a modified data-driven model reduction method up to the order of 4 with very little loss in accuracy. The reduced model is able to predict the behaviour of the original for time ranges beyond the data used for the reduction.