In this study, machine learning modeling is utilized on a bioconvection flow problem. The two dimensional, time independent natural convection flow problem in the presence of oxytactic bacteria is considered inside a unit square cavity. The dimensionless governing equations in terms of stream function, temperature, oxygen concentration, bacteria density and vorticity are numerically solved by radial basis function collocation method. In variation of Rayleigh, bioconvection Rayleigh, Peclet, Lewis and buoyancy ratio parameter, average Nusselt, Sherwood and bacteria density along the heated wall as well as average bacteria density throughout the entire cavity are calculated, and organized into a matrix to form a dataset. This dataset is used for neural network modeling to model the outcomes as averages. The obtained mean squared error metric results show the good fit even by one layer neural network. A feature importance analysis also approves the importance of Rayleigh number on average Nusselt number, and Peclet number on all other averages.