The SpHere INfrared Exoplanet (SHINE) project is a 500-star survey performed with SPHERE on the Very Large Telescope for the purpose ofdirectly detecting new substellar companions and understanding their formation and early evolution. Here we present an initial statistical analysisfor a subsample of 150 stars spanning spectral types from B to M that are representative of the full SHINE sample. Our goal is to constrain thefrequency of substellar companions with masses between 1 and 75MJupand semimajor axes between 5 and 300 au. For this purpose, we adoptdetection limits as a function of angular separation from the survey data for all stars converted into mass and projected orbital separation using theBEX-COND-hot evolutionary tracks and known distance to each system. Based on the results obtained for each star and on the 13 detections inthe sample, we use a Markov chain Monte Carlo tool to compare our observations to two different types of models. The first is a parametric modelbased on observational constraints, and the second type are numerical models that combine advanced core accretion and gravitational instabilityplanet population synthesis. Using the parametric model, we show that the frequencies of systems with at least one substellar companion are23.0+13.5−9.7%, 5.8+4.7−2.8%, and 12.6+12.9−7.1% for BA, FGK, and M stars, respectively. We also demonstrate that a planet-like formation pathway probablydominates the mass range from 1–75MJupfor companions around BA stars, while for M dwarfs, brown dwarf binaries dominate detections. Incontrast, a combination of binary star-like and planet-like formation is required to best fit the observations for FGK stars. Using our populationmodel and restricting our sample to FGK stars, we derive a frequency of 5.7+3.8−2.8%, consistent with predictions from the parametric model. Moregenerally, the frequency values that we derive are in excellent agreement with values obtained in previous studies.