  We consider an analytic model of cosmic star formation which incorporatessupernova feedback, gas accretion and enriched outflows, reproducing thehistory of cosmic star formation, metallicity, supernovae type II rates and thefraction of baryons allocated to structures. We present a new statisticaltreatment of the available observational data on the star formation rate andmetallicity that accounts for the presence of possible systematics. We thenemploy a Bayesian Markov Chain Monte Carlo method to compare the predictions ofour model with observations and derive constraints on the 7 free parameters ofthe model. We find that the dust correction scheme one chooses to adopt for thestar formation data is critical in determining which scenario is favouredbetween a hierarchical star formation model, where star formation is prolongedby accretion, infall and merging, and a monolithic scenario, where starformation is rapid and efficient. We distinguish between these modes bydefining a characteristic minimum mass, M > 10^{11} M_solar, in our fiducialmodel, for early type galaxies where star formation occurs efficiently. Ourresults indicate that the hierarchical star formation model can achieve betteragreement with the data, but that this requires a high efficiency ofsupernova-driven outflows. In a monolithic model, our analysis points to theneed for a mechanism that drives metal-poor winds, perhaps in the form ofsupermassive black hole-induced outflows. Furthermore, the relative absence ofstar formation beyond z ~ 5 in the monolithic scenario requires an alternativemechanism to dwarf galaxies for reionizing the universe at z ~ 11, as requiredby observations of the microwave background. While the monolithic scenario isless favoured in terms of its quality-of-fit, it cannot yet be excluded.
We present a new statistical treatment of the available observational data on the star formation rate and metallicity that accounts for the presence of possible systematics. We then employ a Bayesian Markov Chain Monte Carlo method to compare the predictions of our model with observations and derive constraints on the 7 free parameters of the model.